--- trunk/libf/phylmd/convect3.f 2008/03/31 13:58:05 9 +++ trunk/libf/phylmd/convect3.f 2008/04/18 14:45:53 10 @@ -14,9 +14,9 @@ C C *** THE PARAMETER NA SHOULD IN GENERAL EQUAL ND *** C -c################################################################# + cFleur Introduction des traceurs dans convect3 le 6 juin 200 -c################################################################# + use dimens_m use dimphy use YOMCST @@ -608,15 +608,7 @@ CAPE=CAPE+RD*BUOY(I-1)*DELTAP/P(I-1) DCAPE=RD*BUOY(I-1)*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(*,*) '############################################' -c fin test sb CAPE=AMAX1(0.0,CAPE) C SIGOLD=SIG(I) @@ -639,14 +631,6 @@ AMU=0.5*(SIG(I)+SIGOLD)*W M(I)=AMU*0.007*P(I)*(PH(I)-PH(I+1))/TV(I) -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 --------- - W0(I)=W 98 CONTINUE W0(ICB)=0.5*W0(ICB+1)