Ignore:
Timestamp:
Sep 27, 2019, 11:17:19 PM (5 years ago)
Author:
fhourdin
Message:

Sortie de nouveaux parametres pour controler la distribution des nuages.

File:
1 edited

Legend:

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

    r3495 r3570  
    885885      REAL,SAVE :: sigma1s_factor=1.1
    886886      REAL,SAVE :: sigma1s_power=0.6
     887      REAL,SAVE :: sigma2s_factor=0.09
     888      REAL,SAVE :: sigma2s_power=0.5
    887889      REAL,SAVE :: cloudth_ratqsmin=-1.
    888       !$OMP THREADPRIVATE(sigma1s_factor,sigma1s_power,cloudth_ratqsmin)
     890      !$OMP THREADPRIVATE(sigma1s_factor,sigma1s_power,sigma2s_factor,sigma2s_power,cloudth_ratqsmin)
    889891      INTEGER, SAVE :: iflag_cloudth_vert_noratqs=0
    890892      !$OMP THREADPRIVATE(iflag_cloudth_vert_noratqs)
     
    951953        CALL getin_p('cloudth_sigma1s_power',sigma1s_power)
    952954        WRITE(*,*) 'cloudth_sigma1s_power = ', sigma1s_power
     955        ! Factor used in the calculation of sigma2s
     956        CALL getin_p('cloudth_sigma2s_factor',sigma2s_factor)
     957        WRITE(*,*) 'cloudth_sigma2s_factor = ', sigma2s_factor
     958        ! Power used in the calculation of sigma2s
     959        CALL getin_p('cloudth_sigma2s_power',sigma2s_power)
     960        WRITE(*,*) 'cloudth_sigma2s_power = ', sigma2s_power
    953961        ! Minimum value for the environmental air subgrid water distrib
    954962        CALL getin_p('cloudth_ratqsmin',cloudth_ratqsmin)
     
    10191027      ENDIF
    10201028      sigma1s = sigma1s_fraca + sigma1s_ratqs
    1021       sigma2s=(0.09*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**0.5))+0.002*zqta(ind1,ind2)
     1029      sigma2s=(sigma2s_factor*(((sth-senv)**2)**0.5)/((fraca(ind1,ind2)+0.02)**sigma2s_power))+0.002*zqta(ind1,ind2)
    10221030!      tests
    10231031!      sigma1s=(0.92**0.5)*(fraca(ind1,ind2)**0.5)/(1-fraca(ind1,ind2))*((sth-senv)**2)**0.5+ratqs(ind1,ind2)*po(ind1)
Note: See TracChangeset for help on using the changeset viewer.