107 |
! gnmax : borne superieure de gn (0.0233) |
! gnmax : borne superieure de gn (0.0233) |
108 |
! gninf : vrai si gn est en dessous de sa borne inferieure |
! gninf : vrai si gn est en dessous de sa borne inferieure |
109 |
! gnsup : vrai si gn est en dessus de sa borne superieure |
! gnsup : vrai si gn est en dessus de sa borne superieure |
|
! gm : drole d'objet bien utile |
|
110 |
! ri : nombre de Richardson |
! ri : nombre de Richardson |
111 |
! sn : coefficient de stabilite pour n |
! sn : coefficient de stabilite pour n |
112 |
! snq2 : premier terme du developement limite de sn en q2 |
! snq2 : premier terme du developement limite de sn en q2 |
118 |
REAL gnmax |
REAL gnmax |
119 |
LOGICAL gninf |
LOGICAL gninf |
120 |
LOGICAL gnsup |
LOGICAL gnsup |
|
REAL gm |
|
|
! REAL ri(klon, klev+1) |
|
121 |
REAL sn(klon, klev+1) |
REAL sn(klon, klev+1) |
122 |
REAL snq2(klon, klev+1) |
REAL snq2(klon, klev+1) |
123 |
REAL sm(klon, klev+1) |
REAL sm(klon, klev+1) |
323 |
|
|
324 |
gn=-long(igrid, ilev)**2 / q2(igrid, ilev) & |
gn=-long(igrid, ilev)**2 / q2(igrid, ilev) & |
325 |
* n2(igrid, ilev) |
* n2(igrid, ilev) |
|
gm=long(igrid, ilev)**2 / q2(igrid, ilev) & |
|
|
* m2(igrid, ilev) |
|
|
|
|
326 |
gninf=.false. |
gninf=.false. |
327 |
gnsup=.false. |
gnsup=.false. |
328 |
long(igrid, ilev)=long(igrid, ilev) |
long(igrid, ilev)=long(igrid, ilev) |