Ignore:
Timestamp:
Jul 9, 2014, 7:26:01 PM (10 years ago)
Author:
tnavarro
Message:

Changed two microphysics thresholds for timesteps used in Mesoscale and LES

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/phymars/updaterad.F90

    r1266 r1307  
    2323real, parameter :: r3icemax = 125.e-12 ! ie ricemax  = 500 microns
    2424real, parameter :: ricemax  = 500.e-6
    25 
    26 real, parameter :: qice_threshold  = 1.e-15 ! 1.e-10
    27 
    2825
    2926
     
    7976!! Indeed, one can obtain ricemin for small but non negligible qice values, and therefore hugely opaque clouds.
    8077 
    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
     79rhocloud = (qice*rho_ice + qccn_true*rho_dust) / (qice + qccn_true)
     80rhocloud = min(max(rhocloud,rho_ice),rho_dust)
     81
     82rice = (qice + qccn_true) * 0.75 / pi / rhocloud / nccn_true
    8783 
    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 
     84if (rice .le. r3icemin) then
    9885  rice = ricemin
    99   rhocloud = rho_dust
    100 
    101 endif
     86else if (rice .ge. r3icemax) then
     87  rice = ricemax
     88else
     89  rice = rice**(1./3.) ! here rice is always positive
     90endif
     91
    10292
    10393end subroutine updaterice_micro
Note: See TracChangeset for help on using the changeset viewer.