/[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/libf/phylmd/CV3_routines/cv3_undilute2.f90 revision 47 by guez, Fri Jul 1 15:00:48 2011 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
7              use cvparam3              use cv3_param_m
8              use cvthermo              use cvthermo
9        implicit none        implicit none
10    
# Line 28  Line 28 
28    
29    
30  ! inputs:  ! inputs:
31        integer 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.47  
changed lines
  Added in v.178
  ViewVC Help
Powered by ViewVC 1.1.21