Ignore:
Timestamp:
Jul 19, 2017, 12:19:26 PM (7 years ago)
Author:
fhourdin
Message:

Un lievre leve, sans doute.
Un tout petit seuil sur les largeurs de distributions sous maille
du schéma de nuage, qui produit un gros effet sur le contraste
strato-cumulus / cumulus ...

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/cloudth_mod.F90

    r2959 r2960  
    876876      REAL,SAVE :: sigma1s_factor=1.1
    877877      REAL,SAVE :: sigma1s_power=0.6
    878       !$OMP THREADPRIVATE(sigma1s_factor,sigma1s_power)
     878      REAL,SAVE :: cloudth_ratqsmin=-1.
     879      !$OMP THREADPRIVATE(sigma1s_factor,sigma1s_power,cloudth_ratqsmin)
    879880      INTEGER, SAVE :: iflag_cloudth_vert_noratqs=0
    880881      !$OMP THREADPRIVATE(iflag_cloudth_vert_noratqs)
     
    929930        CALL getin_p('cloudth_sigma1s_power',sigma1s_power)
    930931        WRITE(*,*) 'cloudth_sigma1s_power = ', sigma1s_power
     932        ! Minimum value for the environmental air subgrid water distrib
     933        CALL getin_p('cloudth_ratqsmin',cloudth_ratqsmin)
     934        WRITE(*,*) 'cloudth_ratqsmin = ', cloudth_ratqsmin
    931935        ! Remove the dependency to ratqs from the variance of the vertical PDF
    932936        CALL getin_p('iflag_cloudth_vert_noratqs',iflag_cloudth_vert_noratqs)
     
    988992     &                (1-fraca(ind1,ind2))*((sth-senv)**2)**0.5
    989993!     sigma1s_fraca = (1.1**0.5)*(fraca(ind1,ind2)**0.6)/(1-fraca(ind1,ind2))*((sth-senv)**2)**0.5
    990       sigma1s_ratqs = ratqs(ind1,ind2)*po(ind1)
     994      IF (cloudth_ratqsmin>0.) THEN
     995         sigma1s_ratqs = cloudth_ratqsmin*po(ind1)
     996      ELSE
     997         sigma1s_ratqs = ratqs(ind1,ind2)*po(ind1)
     998      ENDIF
    991999      sigma1s = sigma1s_fraca + sigma1s_ratqs
    9921000      sigma2s=(0.09*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**0.5))+0.002*zqta(ind1,ind2)
Note: See TracChangeset for help on using the changeset viewer.