- Timestamp:
- Jul 23, 2024, 3:29:36 PM (8 weeks ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/lmdz_thermcell_qsat.F90
r5082 r5103 2 2 CONTAINS 3 3 4 subroutinethermcell_qsat(klon,active,pplev,ztemp,zqta,zqsat)4 SUBROUTINE thermcell_qsat(klon,active,pplev,ztemp,zqta,zqsat) 5 5 implicit none 6 6 … … 37 37 38 38 RLvCp = RLVTT/RCPD 39 tout_converge=. false.40 afaire(:)=. false.39 tout_converge=.FALSE. 40 afaire(:)=.FALSE. 41 41 DT(:)=0. 42 42 … … 48 48 !==================================================================== 49 49 50 doig=1,klon50 DO ig=1,klon 51 51 if (active(ig)) then 52 52 Tbef(ig)=ztemp(ig) … … 60 60 zqsat(ig)=qsatbef 61 61 endif 62 enddo 62 END DO 63 63 64 64 ! Traitement du cas ou il y a condensation mais faible 65 65 ! On ne condense pas mais on dit que le qsat est le qta 66 doig=1,klon66 DO ig=1,klon 67 67 if (active(ig)) then 68 68 if (0.<abs(DT(ig)).and.abs(DT(ig))<=DDT0) then … … 70 70 endif 71 71 endif 72 enddo 72 END DO 73 73 74 doiter=1,1074 DO iter=1,10 75 75 afaire(:)=abs(DT(:))>DDT0 76 76 do ig=1,klon … … 93 93 endif 94 94 enddo 95 enddo 95 END DO 96 96 97 97 return
Note: See TracChangeset
for help on using the changeset viewer.