Changeset 992


Ignore:
Timestamp:
Aug 7, 2008, 6:43:31 PM (16 years ago)
Author:
Laurent Fairhead
Message:

Corrections il y avait incohérence, puisque l'on homogénéise les tendances
totales et pas celles dues aux seules descentes précipitantes. JYG
LF

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ4/trunk/libf/phylmd/cv3_routines.F

    r991 r992  
    30063006       csum(il)=0.0
    30073007       dsum(il)=0.0
    3008       enddo
    3009 
     3008        esum(il)=0.0
     3009        fsum(il)=0.0
     3010        gsum(il)=0.0
     3011        hsum(il)=0.0
     3012      enddo
     3013c
     3014      do i=1,nl
     3015       do il=1,ncum
     3016          th_wake(il,i)=t_wake(il,i)*(1000.0/p(il,i))**rdcp
     3017       enddo
     3018      enddo
     3019c
    30103020      do i=1,nl
    30113021       do il=1,ncum
    30123022        if (i.le.(icb(il)-1) .and. iflag(il) .le. 1) then
    3013       asum(il)=asum(il)+ft(il,i)*(ph(il,i)-ph(il,i+1))
    3014       bsum(il)=bsum(il)+fr(il,i)*(lv(il,i)+(cl-cpd)*(t(il,i)-t(il,1)))
     3023cjyg  Saturated part : use T profile
     3024      asum(il)=asum(il)+(ft(il,i)-ftd(il,i))*(ph(il,i)-ph(il,i+1))
     3025      bsum(il)=bsum(il)+(fr(il,i)-fqd(il,i))
     3026     :              *(lv(il,i)+(cl-cpd)*(t(il,i)-t(il,1)))
    30153027     :                  *(ph(il,i)-ph(il,i+1))
    30163028      csum(il)=csum(il)+(lv(il,i)+(cl-cpd)*(t(il,i)-t(il,1)))
    30173029     :                      *(ph(il,i)-ph(il,i+1))
    30183030      dsum(il)=dsum(il)+t(il,i)*(ph(il,i)-ph(il,i+1))/th(il,i)
     3031cjyg  Unsaturated part : use T_wake profile
     3032      esum(il)=esum(il)+ftd(il,i)*(ph(il,i)-ph(il,i+1))
     3033      fsum(il)=fsum(il)+fqd(il,i)
     3034     :              *(lv(il,i)+(cl-cpd)*(t_wake(il,i)-t_wake(il,1)))
     3035     :                  *(ph(il,i)-ph(il,i+1))
     3036      gsum(il)=gsum(il)+(lv(il,i)+(cl-cpd)*(t_wake(il,i)-t_wake(il,1)))
     3037     :                      *(ph(il,i)-ph(il,i+1))
     3038      hsum(il)=hsum(il)+t_wake(il,i)
     3039     ;                      *(ph(il,i)-ph(il,i+1))/th_wake(il,i)
    30193040        endif
    30203041       enddo
     
    30253046       do il=1,ncum
    30263047        if (i.le.(icb(il)-1) .and. iflag(il) .le. 1) then
    3027          ft(il,i)=asum(il)*t(il,i)/(th(il,i)*dsum(il))
    3028          fr(il,i)=bsum(il)/csum(il)
     3048         ftd(il,i)=esum(il)*t_wake(il,i)/(th_wake(il,i)*hsum(il))
     3049         fqd(il,i)=fsum(il)/gsum(il)
     3050         ft(il,i)=ftd(il,i)+asum(il)*t(il,i)/(th(il,i)*dsum(il))
     3051         fr(il,i)=fqd(il,i)+bsum(il)/csum(il)
    30293052        endif
    30303053       enddo
    30313054      enddo
    3032 
    30333055c
    30343056c   ***           reset counter and return           ***
Note: See TracChangeset for help on using the changeset viewer.