Changeset 3117 for LMDZ6/trunk/libf/phylmd/newmicro.F90
- Timestamp:
- Dec 8, 2017, 2:04:01 PM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/newmicro.F90
r3113 r3117 10 10 USE dimphy 11 11 USE phys_local_var_mod, ONLY: scdnc, cldncl, reffclwtop, lcc, reffclws, & 12 reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra , icc3dcon, icc3dstra, zfice12 reffclwc, cldnvi, lcc3d, lcc3dcon, lcc3dstra 13 13 USE phys_state_var_mod, ONLY: rnebcon, clwcon 14 14 USE icefrac_lsc_mod ! computes ice fraction (JBM 3/14) … … 152 152 REAL zrho(klon, klev) !--rho pour la couche 153 153 REAL dh(klon, klev) !--dz pour la couche 154 REAL zfice(klon, klev) 154 155 REAL rad_chaud(klon, klev) !--rayon pour les nuages chauds 155 156 REAL rad_chaud_pi(klon, klev) !--rayon pour les nuages chauds pre-industriels … … 585 586 lcc3dstra(i, k) = lcc3dstra(i, k) - lcc3dcon(i, k) ! eau liquide stratiforme 586 587 lcc3dstra(i, k) = max(lcc3dstra(i,k), 0.0) 587 !588 icc3dcon(i, k) = rnebcon(i, k)*(1-phase3d(i, k))*clwcon(i, k) ! glace589 ! convective590 icc3dstra(i, k) = pclc(i, k)*pqlwp(i, k)*(1-phase3d(i, k))591 icc3dstra(i, k) = icc3dstra(i, k) - icc3dcon(i, k) ! glace stratiforme592 icc3dstra(i, k) = max( icc3dstra(i, k), 0.0)593 !594 588 ! Compute cloud droplet radius as above in meter 595 589 radius = 1.1*((pqlwp(i,k)*pplay(i,k)/(rd*t(i,k)))/(4./3*rpi*1000.* & … … 632 626 IF (lcc3dcon(i,k)<=0.0) lcc3dcon(i, k) = 0.0 633 627 IF (lcc3dstra(i,k)<=0.0) lcc3dstra(i, k) = 0.0 634 IF (icc3dcon(i,k)<=0.0) icc3dcon(i, k) = 0.0635 IF (icc3dstra(i,k)<=0.0) icc3dstra(i, k) = 0.0636 628 END DO 637 629 IF (reffclwtop(i)<=0.0) reffclwtop(i) = 0.0
Note: See TracChangeset
for help on using the changeset viewer.