--- trunk/Sources/phylmd/vdif_kcay.f 2015/06/16 15:23:29 145 +++ trunk/Sources/phylmd/vdif_kcay.f 2016/03/22 16:31:39 188 @@ -107,7 +107,6 @@ ! gnmax : borne superieure de gn (0.0233) ! gninf : vrai si gn est en dessous de sa borne inferieure ! gnsup : vrai si gn est en dessus de sa borne superieure - ! gm : drole d'objet bien utile ! ri : nombre de Richardson ! sn : coefficient de stabilite pour n ! snq2 : premier terme du developement limite de sn en q2 @@ -119,8 +118,6 @@ REAL gnmax LOGICAL gninf LOGICAL gnsup - REAL gm - ! REAL ri(klon, klev+1) REAL sn(klon, klev+1) REAL snq2(klon, klev+1) REAL sm(klon, klev+1) @@ -150,10 +147,7 @@ REAL termq3m2 ! q2min : borne inferieure de q2 - ! q2max : borne superieure de q2 - REAL q2min - REAL q2max ! knmin : borne inferieure de kn ! kmmin : borne inferieure de km @@ -177,7 +171,6 @@ PARAMETER (knmin=1.E-5) PARAMETER (kmmin=1.E-5) PARAMETER (q2min=1.e-5) - PARAMETER (q2max=1.E+2) PARAMETER (nlay=klev) PARAMETER (nlev=klev+1) @@ -326,9 +319,6 @@ gn=-long(igrid, ilev)**2 / q2(igrid, ilev) & * n2(igrid, ilev) - gm=long(igrid, ilev)**2 / q2(igrid, ilev) & - * m2(igrid, ilev) - gninf=.false. gnsup=.false. long(igrid, ilev)=long(igrid, ilev)