Changeset 3112 for LMDZ6/trunk/libf/phylmd/newmicro.F90
- Timestamp:
- Dec 5, 2017, 12:53:00 PM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/newmicro.F90
r2596 r3112 10 10 USE dimphy 11 11 USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, & 12 reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra 12 reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra, icc3dcon, icc3dstra 13 13 USE phys_state_var_mod, ONLY: rnebcon, clwcon 14 14 USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14) … … 586 586 lcc3dstra(i, k) = lcc3dstra(i, k) - lcc3dcon(i, k) ! eau liquide stratiforme 587 587 lcc3dstra(i, k) = max(lcc3dstra(i,k), 0.0) 588 ! 589 icc3dcon(i, k) = rnebcon(i, k)*(1-phase3d(i, k))*clwcon(i, k) ! glace 590 ! convective 591 icc3dstra(i, k) = pclc(i, k)*pqlwp(i, k)*(1-phase3d(i, k)) 592 icc3dstra(i, k) = icc3dstra(i, k) - icc3dcon(i, k) ! glace stratiforme 593 icc3dstra(i, k) = max( icc3dstra(i, k), 0.0) 594 ! 588 595 ! Compute cloud droplet radius as above in meter 589 596 radius = 1.1*((pqlwp(i,k)*pplay(i,k)/(rd*t(i,k)))/(4./3*rpi*1000.* & … … 626 633 IF (lcc3dcon(i,k)<=0.0) lcc3dcon(i, k) = 0.0 627 634 IF (lcc3dstra(i,k)<=0.0) lcc3dstra(i, k) = 0.0 635 IF (icc3dcon(i,k)<=0.0) icc3dcon(i, k) = 0.0 636 IF (icc3dstra(i,k)<=0.0) icc3dstra(i, k) = 0.0 628 637 END DO 629 638 IF (reffclwtop(i)<=0.0) reffclwtop(i) = 0.0
Note: See TracChangeset
for help on using the changeset viewer.