- Timestamp:
- Nov 17, 2023, 11:24:09 AM (20 months ago)
- Location:
- trunk/LMDZ.MARS
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/changelog.txt
r3126 r3127 4340 4340 * THe choice of the isotherm for adsorption can be now done by setting the integer choice_ads in the callphys.def choice_ads = 1 adsorption rate is computed with the H2O thermal speed; choice_ads = 2 adsorption rate is computed based on exeperimental resutls, choice_ads =3 no adsorption 4341 4341 4342 == 17/11/2023 == CS 4343 * Bug fixes for improvedcloud_mod.F : 4344 - Initialization of dMice when low ccn number 4345 - Initialization of zdq at first sub time step (spenttime=0) 4346 - Move 'coef' to fix the formula for zimicro -
trunk/LMDZ.MARS/libf/phymars/improvedclouds_mod.F
r3092 r3127 516 516 endif 517 517 518 ELSE 519 ! Initialization of dMice when it's not computed 520 dMice=0 ! no condensation/sublimation to account for 518 521 ENDIF !of if Nccn>1 519 522 … … 570 573 IF (cloud_adapt_ts) then 571 574 c Estimation of how much is actually condensing/subliming 572 zdq=(dMice/spenttime)*(ptimestep-spenttime) 575 IF (spenttime.ne.0) then 576 zdq=(dMice/spenttime)*(ptimestep-spenttime) 577 ELSE 578 !Initialization for spenttime=0 579 zdq=zpotcond(ig,l)*((ptimestep-spenttime)/ 580 & ptimestep) 581 ENDIF 573 582 c Threshold with powerlaw (sanity check) 574 583 zdq=min(abs(zdq), … … 689 698 alpha=1.81846504e+11 690 699 beta=1.54550140e+00 691 zimicro=ceiling( min(max(coef*alpha*abs(potcond)**beta,5.),7000.))700 zimicro=ceiling(coef*min(max(alpha*abs(potcond)**beta,5.),7000.)) 692 701 693 702 END SUBROUTINE adapt_imicro
Note: See TracChangeset
for help on using the changeset viewer.