diff --git a/src/teb/consphy.fh b/src/teb/consphy.fh
deleted file mode 100644
index 15350c19..00000000
--- a/src/teb/consphy.fh
+++ /dev/null
@@ -1,16 +0,0 @@
-! MC2 (Mesoscale Compressible Community) model.
-! Copyright (c) 2001-2020 Environment Canada (EC-RPN COMM Group).
-! Licensed under the MIT license.
-
-      REAL CPD, CPV, RGASD, RGASV, TRPL, TCDK, RAUW, EPS1, EPS2
-      REAL DELTA, CAPPA, TGL, CONSOL, GRAV, RAYT, STEFAN, PI
-      REAL OMEGA
-      REAL KNAMS, STLO, KARMAN, RIC, CHLC, CHLF
-      REAL T1S, T2S, AW, BW, AI, BI, SLP
-      LOGICAL INIT
-!
-      COMMON/CTESPHY/ INIT, CPD, CPV, RGASD, RGASV, TRPL, TCDK, RAUW,   
-     1                EPS1, EPS2, DELTA, CAPPA, TGL, CONSOL,            
-     1                GRAV, RAYT, STEFAN, PI, OMEGA,                    
-     1                KNAMS, STLO, KARMAN, RIC, CHLC, CHLF,             
-     1                T1S, T2S, AW, BW, AI, BI, SLP
diff --git a/src/teb/impnone.fh b/src/teb/impnone.fh
deleted file mode 100644
index 5128fd60..00000000
--- a/src/teb/impnone.fh
+++ /dev/null
@@ -1,5 +0,0 @@
-! MC2 (Mesoscale Compressible Community) model.
-! Copyright (c) 2001-2020 Environment Canada (EC-RPN COMM Group).
-! Licensed under the MIT license.
-
-      implicit none
diff --git a/src/teb/stabfunc2.fh b/src/teb/stabfunc2.fh
deleted file mode 100644
index 48a476d1..00000000
--- a/src/teb/stabfunc2.fh
+++ /dev/null
@@ -1,61 +0,0 @@
-C MC2 (Mesoscale Compressible Community) model.
-C Copyright (c) 2001-2020 Environment Canada (EC-RPN COMM Group).
-C Licensed under the MIT license.
-
-C   Internal function FMI
-C   Stability function for momentum in the unstable regime (ilmo<0)
-c   Reference: Delage Y. and Girard C. BLM 58 (19-31) Eq. 19
-c
-      REAL FUNCTION FMI(Z2,Z02,LZZ02,ILMO2,X,X0)
-#include "impnone.cdk"
-C
-      REAL, INTENT(IN ) :: Z2,Z02,LZZ02,ILMO2
-      REAL, INTENT(OUT) :: X,X0
-c
-      X =(1-CI*Z2 *BETA*ILMO2)**(0.16666666)
-      X0=(1-CI*Z02*BETA*ILMO2)**(0.16666666)
-      FMI=LZZ02+LOG((X0+1)**2*SQRT(X0**2-X0+1)*(X0**2+X0+1)**1.5
-     %               /((X+1)**2*SQRT(X**2-X+1)*(X**2+X+1)**1.5))
-     %              +RAC3*ATAN(RAC3*((X**2-1)*X0-(X0**2-1)*X)/
-     %              ((X0**2-1)*(X**2-1)+3*X*X0))
-c
-      RETURN
-      END FUNCTION FMI
-c
-C   Internal function FHI
-C   Stability function for heat and moisture in the unstable regime (ilmo<0)
-c   Reference: Delage Y. and Girard C. BLM 58 (19-31) Eq. 17
-c
-      REAL FUNCTION FHI(Z2,Z0T2,LZZ0T2,ILMO2,Y,Y0)
-#include "impnone.cdk"
-C
-      REAL, INTENT(IN ) :: Z2,Z0T2,LZZ0T2,ILMO2
-      REAL, INTENT(OUT) :: Y,Y0
-c
-      Y =(1-CI*Z2  *BETA*ILMO2)**(0.33333333)
-      Y0=(1-CI*Z0T2*BETA*ILMO2)**(0.33333333)
-      FHI=BETA*(LZZ0T2+1.5*LOG((Y0**2+Y0+1)/(Y**2+Y+1))+RAC3*
-     %        ATAN(RAC3*2*(Y-Y0)/((2*Y0+1)*(2*Y+1)+3)))
-c
-      RETURN
-      END FUNCTION FHI
-C
-C   Internal function psi
-C   Stability function for momentum in the stable regime (unsl>0)
-c   Reference :  Y. Delage, BLM, 82 (p23-48) (Eqs.33-37)
-c
-      REAL FUNCTION PSI(Z2,HI2,ILMO2)
-#include "impnone.cdk"
-C
-      REAL a,b,c,d
-      REAL, INTENT(IN ) :: ILMO2,Z2,HI2
-c
-      d = 4*AS*BETA*ILMO2
-      c = d*hi2 - hi2**2
-      b = d - 2*hi2
-      a = sqrt(1 + b*z2 - c*z2**2)
-      psi = 0.5 * (a-z2*hi2-log(1+b*z2*0.5+a)-
-     +            b/(2*sqrt(c))*asin((b-2*c*z2)/d))
-c
-      RETURN
-      END FUNCTION PSI
diff --git a/src/teb/surfcon.fh b/src/teb/surfcon.fh
deleted file mode 100644
index d813d3bf..00000000
--- a/src/teb/surfcon.fh
+++ /dev/null
@@ -1,6 +0,0 @@
-C MC2 (Mesoscale Compressible Community) model.
-C Copyright (c) 2001-2020 Environment Canada (EC-RPN COMM Group).
-C Licensed under the MIT license.
-
-      REAL AS,ASX,CI,BS,BETA,FACTN,HMIN,ANGMAX,RAC3
-      COMMON / SURFCON / AS,ASX,CI,BS,BETA,FACTN,HMIN,ANGMAX,RAC3