Changeset 1307 for trunk/LMDZ.MARS/libf/phymars/updaterad.F90
- Timestamp:
- Jul 9, 2014, 7:26:01 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/LMDZ.MARS/libf/phymars/updaterad.F90
r1266 r1307 23 23 real, parameter :: r3icemax = 125.e-12 ! ie ricemax = 500 microns 24 24 real, parameter :: ricemax = 500.e-6 25 26 real, parameter :: qice_threshold = 1.e-15 ! 1.e-1027 28 25 29 26 … … 79 76 !! Indeed, one can obtain ricemin for small but non negligible qice values, and therefore hugely opaque clouds. 80 77 81 if (qice .ge. qice_threshold) then 82 83 rhocloud = (qice*rho_ice + qccn_true*rho_dust) / (qice + qccn_true) 84 rhocloud = min(max(rhocloud,rho_ice),rho_dust) 85 86 rice = (qice + qccn_true) * 0.75 / pi / rhocloud / nccn_true 78 79 rhocloud = (qice*rho_ice + qccn_true*rho_dust) / (qice + qccn_true) 80 rhocloud = min(max(rhocloud,rho_ice),rho_dust) 81 82 rice = (qice + qccn_true) * 0.75 / pi / rhocloud / nccn_true 87 83 88 if (rice .le. r3icemin) then 89 rice = ricemin 90 else if (rice .ge. r3icemax) then 91 rice = ricemax 92 else 93 rice = rice**(1./3.) ! here rice is always positive 94 endif 95 96 else 97 84 if (rice .le. r3icemin) then 98 85 rice = ricemin 99 rhocloud = rho_dust 100 101 endif 86 else if (rice .ge. r3icemax) then 87 rice = ricemax 88 else 89 rice = rice**(1./3.) ! here rice is always positive 90 endif 91 102 92 103 93 end subroutine updaterice_micro
Note: See TracChangeset
for help on using the changeset viewer.