20 |
REAL, intent(in):: pqen(klon, klev) ! humidite (environnement) |
REAL, intent(in):: pqen(klon, klev) ! humidite (environnement) |
21 |
REAL, intent(in):: pqsen(klon, klev) ! humidite saturante (environnement) |
REAL, intent(in):: pqsen(klon, klev) ! humidite saturante (environnement) |
22 |
REAL, intent(in):: pgeo(klon, klev) ! geopotentiel (g * metre) |
REAL, intent(in):: pgeo(klon, klev) ! geopotentiel (g * metre) |
23 |
REAL paph(klon, klev+1) ! pression aux demi-niveaux |
REAL, intent(in):: paph(klon, klev+1) ! pression aux demi-niveaux |
24 |
REAL pgeoh(klon, klev) ! geopotentiel aux demi-niveaux |
REAL pgeoh(klon, klev) ! geopotentiel aux demi-niveaux |
25 |
REAL ptenh(klon, klev) ! temperature aux demi-niveaux |
REAL ptenh(klon, klev) ! temperature aux demi-niveaux |
26 |
REAL pqenh(klon, klev) ! humidite aux demi-niveaux |
REAL pqenh(klon, klev) ! humidite aux demi-niveaux |
72 |
ENDDO |
ENDDO |
73 |
|
|
74 |
icall=0 |
icall=0 |
75 |
CALL flxadjtq(paph(1, k), ptenh(1, k), pqsenh(1, k), llflag, icall) |
CALL flxadjtq(paph(:, k), ptenh(1, k), pqsenh(1, k), llflag, icall) |
76 |
|
|
77 |
DO i = 1, klon |
DO i = 1, klon |
78 |
pqenh(i, k)=MIN(pqen(i, k-1), pqsen(i, k-1)) & |
pqenh(i, k)=MIN(pqen(i, k-1), pqsen(i, k-1)) & |