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 |
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) |
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) |