- Timestamp:
- Jul 24, 2024, 4:23:34 PM (4 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_qsat.F90
r5116 r5117 24 24 REAL Tbef(klon),DT(klon) 25 25 REAL tdelta,qsatbef,zcor,qlbef,zdelta,zcvm5,dqsat,num,denom,dqsat_dT 26 logicalZsat26 LOGICAL Zsat 27 27 REAL RLvCp 28 28 … … 49 49 50 50 DO ig=1,klon 51 if(active(ig)) THEN51 IF (active(ig)) THEN 52 52 Tbef(ig)=ztemp(ig) 53 53 zdelta=MAX(0.,SIGN(1.,RTT-Tbef(ig))) … … 65 65 ! On ne condense pas mais on dit que le qsat est le qta 66 66 DO ig=1,klon 67 if(active(ig)) THEN68 if (0.<abs(DT(ig)).and.abs(DT(ig))<=DDT0) THEN67 IF (active(ig)) THEN 68 IF (0.<abs(DT(ig)).AND.abs(DT(ig))<=DDT0) THEN 69 69 zqsat(ig)=zqta(ig) 70 70 endif … … 75 75 afaire(:)=abs(DT(:))>DDT0 76 76 do ig=1,klon 77 if(afaire(ig)) THEN77 IF (afaire(ig)) THEN 78 78 Tbef(ig)=Tbef(ig)+DT(ig) 79 79 zdelta=MAX(0.,SIGN(1.,RTT-Tbef(ig)))
Note: See TracChangeset
for help on using the changeset viewer.