/[lmdze]/trunk/libf/phylmd/conema3.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/conema3.f90

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 12 by guez, Mon Jul 21 16:05:07 2008 UTC
# Line 62  c Line 62  c
62        INTEGER ntra,ntrac !number of tracers; if no tracer transport        INTEGER ntra,ntrac !number of tracers; if no tracer transport
63                           ! is needed, set ntra = 1 (or 0)                           ! is needed, set ntra = 1 (or 0)
64        PARAMETER (ntrac=nqmx-2)        PARAMETER (ntrac=nqmx-2)
65        REAL dtime        REAL, intent(in):: dtime
66  c  c
67        REAL d_t2(klon,klev), d_q2(klon,klev) ! sbl        REAL d_t2(klon,klev), d_q2(klon,klev) ! sbl
68        REAL d_u2(klon,klev), d_v2(klon,klev) ! sbl        REAL d_u2(klon,klev), d_v2(klon,klev) ! sbl
# Line 70  c Line 70  c
70        REAL em_d_u2(klev), em_d_v2(klev)     ! sbl          REAL em_d_u2(klev), em_d_v2(klev)     ! sbl  
71  c  c
72        REAL, intent(in):: paprs(klon,klev+1)        REAL, intent(in):: paprs(klon,klev+1)
73        real pplay(klon,klev)        real, intent(in):: pplay(klon,klev)
74        REAL t(klon,klev), q(klon,klev), d_t(klon,klev), d_q(klon,klev)        REAL t(klon,klev), q(klon,klev), d_t(klon,klev), d_q(klon,klev)
75        REAL u(klon,klev), v(klon,klev), tra(klon,klev,ntra)        REAL u(klon,klev), v(klon,klev), tra(klon,klev,ntra)
76        REAL d_u(klon,klev), d_v(klon,klev), d_tra(klon,klev,ntra)        REAL d_u(klon,klev), d_v(klon,klev), d_tra(klon,klev,ntra)
# Line 249  c       CALL cv_thermo(iflagcon) Line 249  c       CALL cv_thermo(iflagcon)
249    
250  c -- appel schema de nuages:  c -- appel schema de nuages:
251    
 c       CALL CLOUDS_SUB_LS(klev,em_q,em_qs,em_t  
 c    i          ,em_p,em_ph,dtime,em_qcondc  
 c    o          ,em_cldf,em_cldq,em_pradj,em_ftadj,em_fradj,em_ifc)  
   
252          do k = 1, klev          do k = 1, klev
253           cldf(i,k)  = em_cldf(k)  ! cloud fraction (0-1)           cldf(i,k)  = em_cldf(k)  ! cloud fraction (0-1)
254           cldq(i,k)  = em_cldq(k)  ! in-cloud water content (kg/kg)           cldq(i,k)  = em_cldq(k)  ! in-cloud water content (kg/kg)
# Line 356  c  precip. Line 352  c  precip.
352              if (ktop(i)-kbas(i).gt.0.and.              if (ktop(i)-kbas(i).gt.0.and.
353       s         l.ge.kbas(i).and.l.le.ktop(i)) then       s         l.ge.kbas(i).and.l.le.ktop(i)) then
354                 qcond_incld(i,l)=rain(i)*8.e4                 qcond_incld(i,l)=rain(i)*8.e4
 c    s         *(pplay(i,l      )-paprs(i,ktop(i)+1))  
355       s         /(pplay(i,kbas(i))-pplay(i,ktop(i)))       s         /(pplay(i,kbas(i))-pplay(i,ktop(i)))
356  c    s         **2  c    s         **2
357              else              else
Legend:
Removed from v.3  
changed lines
  Added in v.12
  ViewVC Help
Powered by ViewVC 1.1.21