LMDZ
qsat.F90
Go to the documentation of this file.
1 subroutine qsat (dq, q, e, p, t, r)
2 
3 #include "tsmbkind.h"
4 
5 implicit none
6 
7 real_b :: dq, q, e, p, t, r, tf, a, b
8 
9 tf=273.16_jprb
10 if (t.gt.tf) then
11  a = 17.269_jprb
12  b = 35.86_jprb
13 else
14  a = 21.874_jprb
15  b = 7.66_jprb
16 end if
17 e = 6.11_jprb * exp(a*(t-tf)/(t-b))
18 q = 0.622_jprb * r*e/(p-r*e)
19 dq=a*p*q*(tf-b)
20 dq=dq/((p-r*e)*(t-b)*(t-b))
21 
22 !print 9001, P,T,R,A,B,E,Q
23 9001 format (1x,'QSAT ',5f10.3,2e12.5)
24 
25 !----------------------------------------
26 return
27 end subroutine qsat
r
!$Id mode_top_bound COMMON comconstr r
Definition: comconst.h:7
dq
do llm!au dessus on relaxe vers profil init!on fait l hypothese que dans ce il n y a plus d eau liq au dessus!donc la relaxation en thetal et qt devient relaxation en tempe et qv l dq1 relax dq(l, 1)
qsat
subroutine qsat(dq, q, e, p, t, r)
Definition: qsat.F90:2
x
!$Header!c c INCLUDE fxyprim h c c c Fonctions in line c c REAL fyprim REAL rj c c il faut la calculer avant d appeler ces fonctions c c c Fonctions a changer selon x(x) et y(y) choisis.c-----------------------------------------------------------------c c.....ici