Ignore:
Timestamp:
Mar 16, 2018, 7:28:53 PM (7 years ago)
Author:
musat
Message:

Seuil maximum nombre de noyaux de condensation cdnc_max
passe a physiq.def (utile pour l'experience de type aquaplanete)

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/IPSLCM6.0.15/libf/phylmd/newmicro.F90

    r3248 r3282  
    166166  REAL,SAVE :: cdnc_min_m3
    167167  !$OMP THREADPRIVATE(cdnc_min,cdnc_min_m3)
     168  REAL,SAVE :: cdnc_max=-1.
     169  REAL,SAVE :: cdnc_max_m3
     170  !$OMP THREADPRIVATE(cdnc_max,cdnc_max_m3)
    168171 
    169172  ! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     
    183186    cdnc_min_m3=cdnc_min*1E6
    184187    IF (cdnc_min_m3<0.) cdnc_min_m3=20.E6 ! astuce pour retrocompatibilite
     188    write(lunout,*)'cdnc_min=', cdnc_min_m3/1.E6
     189    call getin_p('cdnc_max',cdnc_max)
     190    cdnc_max_m3=cdnc_max*1E6
     191    IF (cdnc_max_m3<0.) cdnc_max_m3=1000.E6 ! astuce pour retrocompatibilite
     192    write(lunout,*)'cdnc_max=', cdnc_max_m3/1.E6
    185193  ENDIF
    186194
     
    251259        cdnc(i, k) = 10.**(bl95_b0+bl95_b1*log(max(mass_solu_aero(i,k), &
    252260          1.E-4))/log(10.))*1.E6 !-m-3
    253         cdnc(i, k) = min(1000.E6, max(cdnc_min_m3,cdnc(i,k)))
     261        cdnc(i, k) = min(cdnc_max_m3, max(cdnc_min_m3,cdnc(i,k)))
    254262
    255263        ! --pre-industrial case
    256264        cdnc_pi(i, k) = 10.**(bl95_b0+bl95_b1*log(max(mass_solu_aero_pi(i,k), &
    257265          1.E-4))/log(10.))*1.E6 !-m-3
    258         cdnc_pi(i, k) = min(1000.E6, max(cdnc_min_m3,cdnc_pi(i,k)))
     266        cdnc_pi(i, k) = min(cdnc_max_m3, max(cdnc_min_m3,cdnc_pi(i,k)))
    259267
    260268        ! --present-day case
Note: See TracChangeset for help on using the changeset viewer.