/[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 10 by guez, Fri Apr 18 14:45:53 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
# Line 608  C          The interval on which CAPE is Line 608  C          The interval on which CAPE is
608           CAPE=CAPE+RD*BUOY(I-1)*DELTAP/P(I-1)           CAPE=CAPE+RD*BUOY(I-1)*DELTAP/P(I-1)
609           DCAPE=RD*BUOY(I-1)*DELTAP/P(I-1)           DCAPE=RD*BUOY(I-1)*DELTAP/P(I-1)
610           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(*,*) '############################################'  
611    
 c fin test sb  
612           CAPE=AMAX1(0.0,CAPE)           CAPE=AMAX1(0.0,CAPE)
613  C  C
614           SIGOLD=SIG(I)           SIGOLD=SIG(I)
# Line 639  C Line 631  C
631           AMU=0.5*(SIG(I)+SIGOLD)*W           AMU=0.5*(SIG(I)+SIGOLD)*W
632           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)
633    
 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 ---------  
   
634           W0(I)=W           W0(I)=W
635     98 CONTINUE     98 CONTINUE
636        W0(ICB)=0.5*W0(ICB+1)        W0(ICB)=0.5*W0(ICB+1)
Legend:
Removed from v.3  
changed lines
  Added in v.10
  ViewVC Help
Powered by ViewVC 1.1.21