/[lmdze]/trunk/Sources/phylmd/CV30_routines/cv30_undilute1.f
ViewVC logotype

Diff of /trunk/Sources/phylmd/CV30_routines/cv30_undilute1.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

trunk/phylmd/CV3_routines/cv3_undilute1.f revision 103 by guez, Fri Aug 29 13:00:05 2014 UTC trunk/Sources/phylmd/CV30_routines/cv3_undilute1.f revision 184 by guez, Wed Mar 16 14:50:46 2016 UTC
# Line 33  Line 33 
33  ! local variables:  ! local variables:
34        integer i, k        integer i, k
35        integer icb1(len), icbs(len), icbsmax2 ! convect3        integer icb1(len), icbs(len), icbsmax2 ! convect3
36        real tg, qg, alv, s, ahg, tc, denom, es, rg        real tg, qg, alv, s, ahg, tc, denom, es
37        real ah0(len), cpp(len)        real ah0(len), cpp(len)
38        real tnk(len), qnk(len), gznk(len), ticb(len), gzicb(len)        real tnk(len), qnk(len), gznk(len), ticb(len), gzicb(len)
39        real qsicb(len) ! convect3        real qsicb(len) ! convect3
# Line 176  Line 176 
176           clw(i,icbs(i))=qnk(i)-qg           clw(i,icbs(i))=qnk(i)-qg
177           clw(i,icbs(i))=max(0.0,clw(i,icbs(i)))           clw(i,icbs(i))=max(0.0,clw(i,icbs(i)))
178    
          rg=qg/(1.-qnk(i))  
 ! ori         tvp(i,icb(i))=tp(i,icb(i))*(1.+rg*epsi)  
179  ! convect3: (qg utilise au lieu du vrai mixing ratio rg)  ! convect3: (qg utilise au lieu du vrai mixing ratio rg)
180           tvp(i,icbs(i))=tp(i,icbs(i))*(1.+qg/eps-qnk(i)) !whole thing           tvp(i,icbs(i))=tp(i,icbs(i))*(1.+qg/eps-qnk(i)) !whole thing
181    
# Line 275  Line 273 
273           clw(i,icb(i)+1)=qnk(i)-qg           clw(i,icb(i)+1)=qnk(i)-qg
274           clw(i,icb(i)+1)=max(0.0,clw(i,icb(i)+1))           clw(i,icb(i)+1)=max(0.0,clw(i,icb(i)+1))
275    
          rg=qg/(1.-qnk(i))  
 ! ori         tvp(i,icb(i))=tp(i,icb(i))*(1.+rg*epsi)  
276  ! convect3: (qg utilise au lieu du vrai mixing ratio rg)  ! convect3: (qg utilise au lieu du vrai mixing ratio rg)
277           tvp(i,icb(i)+1)=tp(i,icb(i)+1)*(1.+qg/eps-qnk(i)) !whole thing           tvp(i,icb(i)+1)=tp(i,icb(i)+1)*(1.+qg/eps-qnk(i)) !whole thing
278    
Legend:
Removed from v.103  
changed lines
  Added in v.184
  ViewVC Help
Powered by ViewVC 1.1.21