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

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

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

trunk/phylmd/CV3_routines/cv3_undilute2.f revision 97 by guez, Fri Apr 25 14:58:31 2014 UTC trunk/Sources/phylmd/CV3_routines/cv3_undilute2.f revision 178 by guez, Fri Mar 11 18:47:26 2016 UTC
# Line 1  Line 1 
1    
2        SUBROUTINE cv3_undilute2(nloc,ncum,nd,icb,icbs,nk &        SUBROUTINE cv3_undilute2(nloc,ncum,nd,icb,icbs,nk &
3                               ,tnk,qnk,gznk,t,q,qs,gz &                               ,tnk,qnk,gznk,t,qs,gz &
4                               ,p,h,tv,lv,pbase,buoybase,plcl &                               ,p,h,tv,lv,pbase,buoybase,plcl &
5                               ,inb,tp,tvp,clw,hp,ep,sigp,buoy)                               ,inb,tp,tvp,clw,hp,ep,sigp,buoy)
6        use conema3_m        use conema3_m
# Line 30  Line 30 
30  ! inputs:  ! inputs:
31        integer, intent(in):: ncum, nd, nloc        integer, intent(in):: ncum, nd, nloc
32        integer icb(nloc), icbs(nloc), nk(nloc)        integer icb(nloc), icbs(nloc), nk(nloc)
33        real t(nloc,nd), q(nloc,nd), qs(nloc,nd), gz(nloc,nd)        real t(nloc,nd), qs(nloc,nd), gz(nloc,nd)
34        real p(nloc,nd)        real p(nloc,nd)
35        real tnk(nloc), qnk(nloc), gznk(nloc)        real tnk(nloc), qnk(nloc), gznk(nloc)
36        real lv(nloc,nd), tv(nloc,nd), h(nloc,nd)        real lv(nloc,nd), tv(nloc,nd), h(nloc,nd)
# Line 44  Line 44 
44    
45  ! local variables:  ! local variables:
46        integer i, k        integer i, k
47        real tg,qg,ahg,alv,s,tc,es,denom,rg,tca,elacrit        real tg,qg,ahg,alv,s,tc,es,denom
48        real by, defrac, pden        real pden
49        real ah0(nloc), cape(nloc), capem(nloc), byp(nloc)        real ah0(nloc)
       logical lcape(nloc)  
50    
51  !=====================================================================  !=====================================================================
52  ! --- SOME INITIALIZATIONS  ! --- SOME INITIALIZATIONS
# Line 142  Line 141 
141    
142                 clw(i,k)=qnk(i)-qg                 clw(i,k)=qnk(i)-qg
143                 clw(i,k)=max(0.0,clw(i,k))                 clw(i,k)=max(0.0,clw(i,k))
                rg=qg/(1.-qnk(i))  
 ! ori               tvp(i,k)=tp(i,k)*(1.+rg*epsi)  
144  ! convect3: (qg utilise au lieu du vrai mixing ratio rg):  ! convect3: (qg utilise au lieu du vrai mixing ratio rg):
145                 tvp(i,k)=tp(i,k)*(1.+qg/eps-qnk(i)) ! whole thing                 tvp(i,k)=tp(i,k)*(1.+qg/eps-qnk(i)) ! whole thing
146              endif              endif
Legend:
Removed from v.97  
changed lines
  Added in v.178
  ViewVC Help
Powered by ViewVC 1.1.21