/[lmdze]/trunk/phylmd/CV30_routines/cv30_yield.f
ViewVC logotype

Diff of /trunk/phylmd/CV30_routines/cv30_yield.f

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

trunk/phylmd/CV3_routines/cv3_yield.f revision 105 by guez, Thu Sep 4 10:40:24 2014 UTC trunk/Sources/phylmd/CV30_routines/cv30_yield.f revision 186 by guez, Mon Mar 21 15:36:26 2016 UTC
# Line 1  Line 1 
1  module cv3_yield_m  module cv30_yield_m
2    
3    implicit none    implicit none
4    
5  contains  contains
6    
7    SUBROUTINE cv3_yield(nloc,ncum,nd,na  &    SUBROUTINE cv30_yield(nloc,ncum,nd,na  &
8         ,icb,inb,delt &         ,icb,inb,delt &
9         ,t,rr,u,v,gz,p,ph,h,hp,lv,cpn,th &         ,t,rr,u,v,gz,p,ph,h,hp,lv,cpn,th &
10         ,ep,clw,m,tp,mp,rp,up,vp &         ,ep,clw,m,tp,mp,rp,up,vp &
# Line 14  contains Line 14  contains
14         ,iflag,precip,VPrecip,ft,fr,fu,fv &         ,iflag,precip,VPrecip,ft,fr,fu,fv &
15         ,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd)         ,upwd,dnwd,dnwd0,ma,mike,tls,tps,qcondc,wd)
16      use conema3_m      use conema3_m
17      use cv3_param_m      use cv30_param_m
18      use cvthermo      use cvthermo
19      use cvflag      use cvflag
20    
# Line 54  contains Line 54  contains
54      real rat, awat, delti      real rat, awat, delti
55      real ax, bx, cx, dx      real ax, bx, cx, dx
56      real cpinv, rdcp, dpinv      real cpinv, rdcp, dpinv
57      real lvcp(nloc,na), mke(nloc,na)      real lvcp(nloc,na)
58      real am(nloc), work(nloc), ad(nloc), amp1(nloc)      real am(nloc), work(nloc), ad(nloc), amp1(nloc)
59  !!!      real up1(nloc), dn1(nloc)  !!!      real up1(nloc), dn1(nloc)
60      real up1(nloc,nd,nd), dn1(nloc,nd,nd)      real up1(nloc,nd,nd), dn1(nloc,nd,nd)
# Line 230  contains Line 230  contains
230         enddo         enddo
231         if(num1.le.0) cycle         if(num1.le.0) cycle
232    
233         call zilch(amp1,ncum)         amp1(:ncum) = 0.
234         call zilch(ad,ncum)         ad(:ncum) = 0.
235    
236         do  k=i+1,nl+1         do  k=i+1,nl+1
237            do  il=1,ncum            do  il=1,ncum
# Line 627  contains Line 627  contains
627      !ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc      !ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
628    
629      do i=1,nd      do i=1,nd
        do il=1,ncum  
           mke(il,i)=upwd(il,i)+dnwd(il,i)  
        enddo  
     enddo  
   
     do i=1,nd  
630         DO  il=1,ncum         DO  il=1,ncum
631            rdcp=(rrd*(1.-rr(il,i))-rr(il,i)*rrv) &            rdcp=(rrd*(1.-rr(il,i))-rr(il,i)*rrv) &
632                 /(cpd*(1.-rr(il,i))+rr(il,i)*cpv)                 /(cpd*(1.-rr(il,i))+rr(il,i)*cpv)
# Line 703  contains Line 697  contains
697         enddo                                               ! cld         enddo                                               ! cld
698      enddo                                                ! cld      enddo                                                ! cld
699    
700    end SUBROUTINE cv3_yield    end SUBROUTINE cv30_yield
701    
702  end module cv3_yield_m  end module cv30_yield_m
Legend:
Removed from v.105  
changed lines
  Added in v.186
  ViewVC Help
Powered by ViewVC 1.1.21