/[lmdze]/trunk/libf/phylmd/convect3.f
ViewVC logotype

Diff of /trunk/libf/phylmd/convect3.f

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 14  cccc     *     DTVPDT1,DTVPDQ1,DPLCLDT,D Line 14  cccc     *     DTVPDT1,DTVPDQ1,DPLCLDT,D
14  C  C
15  C    ***  THE PARAMETER NA SHOULD IN GENERAL EQUAL ND   ***  C    ***  THE PARAMETER NA SHOULD IN GENERAL EQUAL ND   ***
16  C  C
17  c#################################################################  
18  cFleur       Introduction des traceurs dans convect3 le 6 juin 200  cFleur       Introduction des traceurs dans convect3 le 6 juin 200
19  c#################################################################  
20        use dimens_m        use dimens_m
21        use dimphy        use dimphy
22        use YOMCST        use YOMCST
23    
24          real, intent(in):: dtime, DELT
25        PARAMETER (NA=60)        PARAMETER (NA=60)
26    
27        integer NTRAC        integer NTRAC
# Line 608  C          The interval on which CAPE is Line 610  C          The interval on which CAPE is
610           CAPE=CAPE+RD*BUOY(I-1)*DELTAP/P(I-1)           CAPE=CAPE+RD*BUOY(I-1)*DELTAP/P(I-1)
611           DCAPE=RD*BUOY(I-1)*DELTAP/P(I-1)           DCAPE=RD*BUOY(I-1)*DELTAP/P(I-1)
612           DLNP=DELTAP/P(I-1)           DLNP=DELTAP/P(I-1)
 Cjyg2  
 c sb3d         print *,'buoy,dlnp,dcape,cape',buoy(i-1),dlnp,dcape,cape  
 c test sb:  
 c@       write(*,*) '############################################'  
 c@         write(*,*) 'cape,rrd,buoy,deltap,p,pbase,ph:'  
 c@     :     ,cape,rd,buoy(i-1),deltap,p(i-1),pbase,ph(i)  
 c@       write(*,*) '############################################'  
613    
 c fin test sb  
614           CAPE=AMAX1(0.0,CAPE)           CAPE=AMAX1(0.0,CAPE)
615  C  C
616           SIGOLD=SIG(I)           SIGOLD=SIG(I)
# Line 639  C Line 633  C
633           AMU=0.5*(SIG(I)+SIGOLD)*W           AMU=0.5*(SIG(I)+SIGOLD)*W
634           M(I)=AMU*0.007*P(I)*(PH(I)-PH(I+1))/TV(I)           M(I)=AMU*0.007*P(I)*(PH(I)-PH(I+1))/TV(I)
635    
 c --------- test sb:  
 c       write(*,*) '############################################'  
 c       write(*,*) 'k,amu,buoy(k-1),deltap,w,beta,fac,cape,w0(k)'  
 c       write(*,*) i,amu,buoy(i-1),deltap  
 c     :           ,w,beta,fac,cape,w0(i)  
 c       write(*,*) '############################################'  
 c ---------  
   
636           W0(I)=W           W0(I)=W
637     98 CONTINUE     98 CONTINUE
638        W0(ICB)=0.5*W0(ICB+1)        W0(ICB)=0.5*W0(ICB+1)
Legend:
Removed from v.3  
changed lines
  Added in v.12
  ViewVC Help
Powered by ViewVC 1.1.21