Ignore:
Timestamp:
Feb 22, 2024, 5:29:02 PM (7 months ago)
Author:
evignon
Message:

changements suite à l'atelier nuage d'aujourd'hui.

File:
1 edited

Legend:

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

    r4818 r4830  
    112112  !$OMP THREADPRIVATE(coef_eva)
    113113
    114   REAL, SAVE, PROTECTED :: coef_eva_i                       ! tuning coefficient ice precip sublimation
    115   !$OMP THREADPRIVATE(coef_eva_i)
     114  REAL, SAVE, PROTECTED :: coef_sub                        ! tuning coefficient ice precip sublimation
     115  !$OMP THREADPRIVATE(coef_sub)
    116116
    117117  REAL, SAVE, PROTECTED :: expo_eva=0.5                     ! tuning coefficient liquid precip evaporation
    118118  !$OMP THREADPRIVATE(expo_eva)
    119119
    120   REAL, SAVE, PROTECTED :: expo_eva_i                       ! tuning coefficient ice precip sublimation
    121   !$OMP THREADPRIVATE(expo_eva_i)
     120  REAL, SAVE, PROTECTED :: expo_sub                       ! tuning coefficient ice precip sublimation
     121  !$OMP THREADPRIVATE(expo_sub)
    122122
    123123  REAL, SAVE, PROTECTED :: cice_velo=1.645                  ! factor in the ice fall velocity formulation
     
    137137  !$OMP THREADPRIVATE(ok_poprecip)
    138138
     139  REAL, SAVE, PROTECTED :: cld_lc_lsc_snow               ! snow autoconversion coefficient, stratiform rain
     140  !$OMP THREADPRIVATE(cld_lc_lsc_snow)
     141
     142  REAL, SAVE, PROTECTED :: cld_lc_con_snow                ! snow autoconversion coefficient, convective rain
     143  !$OMP THREADPRIVATE(cld_lc_con_snow)
     144
    139145  REAL, SAVE, PROTECTED :: rain_int_min=0.001               ! Minimum local rain intensity [mm/s] before the decrease in associated precipitation fraction
    140146  !$OMP THREADPRIVATE(rain_int_min)
     
    155161  !$OMP THREADPRIVATE(rho_rain)
    156162
    157   REAL, SAVE, PROTECTED :: rho_snow=500.                   ! A COMMENTER TODO [kg/m3]
     163  REAL, SAVE, PROTECTED :: rho_snow                        ! A COMMENTER TODO [kg/m3]
    158164  !$OMP THREADPRIVATE(rho_snow)
    159165
     
    164170  !$OMP THREADPRIVATE(r_snow)
    165171
    166   REAL, SAVE, PROTECTED :: tau_auto_snow_min=1800.          ! A COMMENTER TODO [s]
     172  REAL, SAVE, PROTECTED :: tau_auto_snow_min=100.          ! A COMMENTER TODO [s]
    167173  !$OMP THREADPRIVATE(tau_auto_snow_min)
    168174
    169   REAL, SAVE, PROTECTED :: tau_auto_snow_max=7200.          ! A COMMENTER TODO [s]
     175  REAL, SAVE, PROTECTED :: tau_auto_snow_max=1000.          ! A COMMENTER TODO [s]
    170176  !$OMP THREADPRIVATE(tau_auto_snow_max)
    171177
     
    182188  !$OMP THREADPRIVATE(alpha_freez)
    183189
    184   REAL, SAVE, PROTECTED :: gamma_freez=0.1               ! A COMMENTER TODO [m-3.s-1]
     190  REAL, SAVE, PROTECTED :: beta_freez=0.1                 ! A COMMENTER TODO [m-3.s-1]
     191  !$OMP THREADPRIVATE(beta_freez)
     192
     193  REAL, SAVE, PROTECTED :: gamma_freez=1.                 ! A COMMENTER TODO [-]
    185194  !$OMP THREADPRIVATE(gamma_freez)
     195
     196  REAL, SAVE, PROTECTED :: rain_fallspeed=4.              ! A COMMENTER TODO [m/s]
     197  !$OMP THREADPRIVATE(rain_fallspeed)
     198
     199  REAL, SAVE, PROTECTED :: rain_fallspeed_clr              ! A COMMENTER TODO [m/s]
     200  !$OMP THREADPRIVATE(rain_fallspeed_clr)
     201
     202  REAL, SAVE, PROTECTED :: rain_fallspeed_cld             ! A COMMENTER TODO [m/s]
     203  !$OMP THREADPRIVATE(rain_fallspeed_cld)
     204
     205  REAL, SAVE, PROTECTED :: snow_fallspeed=1.             ! A COMMENTER TODO [m/s]
     206  !$OMP THREADPRIVATE(snow_fallspeed)
     207
     208  REAL, SAVE, PROTECTED :: snow_fallspeed_clr             ! A COMMENTER TODO [m/s]
     209  !$OMP THREADPRIVATE(snow_fallspeed_clr)
     210
     211  REAL, SAVE, PROTECTED :: snow_fallspeed_cld             ! A COMMENTER TODO [m/s]
     212  !$OMP THREADPRIVATE(snow_fallspeed_cld)
    186213  !--End of the parameters for poprecip
    187214
     
    246273    CALL getin_p('cld_lc_lsc',cld_lc_lsc)
    247274    CALL getin_p('cld_lc_con',cld_lc_con)
     275    cld_lc_lsc_snow=cld_lc_lsc
     276    cld_lc_con_snow=cld_lc_con
     277    CALL getin_p('cld_lc_lsc_snow',cld_lc_lsc_snow)
     278    CALL getin_p('cld_lc_con_snow',cld_lc_con_snow)
    248279    CALL getin_p('cld_tau_lsc',cld_tau_lsc)
    249280    CALL getin_p('cld_tau_con',cld_tau_con)
     
    253284    CALL getin_p('ffallv_lsc',ffallv_con)
    254285    CALL getin_p('coef_eva',coef_eva)
    255     coef_eva_i=coef_eva
    256     CALL getin_p('coef_eva_i',coef_eva_i)
     286    coef_sub=coef_eva
     287    CALL getin_p('coef_eva_i',coef_sub)
     288    CALL getin_p('coef_sub',coef_sub)
    257289    CALL getin_p('expo_eva',expo_eva)
    258     expo_eva_i=expo_eva
    259     CALL getin_p('expo_eva_i',expo_eva_i)
     290    expo_sub=expo_eva
     291    CALL getin_p('expo_sub',expo_sub)
    260292    CALL getin_p('iflag_autoconversion',iflag_autoconversion)
    261293    CALL getin_p('dist_liq',dist_liq)
     
    267299    CALL getin_p('gamma_col',gamma_col)
    268300    CALL getin_p('gamma_rim',gamma_rim)
     301    CALL getin_p('gamma_freez',gamma_freez)
     302    CALL getin_p('r_snow',r_snow)
     303    CALL getin_p('rain_fallspeed',rain_fallspeed)
     304    rain_fallspeed_clr=rain_fallspeed
     305    rain_fallspeed_cld=rain_fallspeed
     306    CALL getin_p('rain_fallspeed_clr',rain_fallspeed_clr)
     307    CALL getin_p('rain_fallspeed_cld',rain_fallspeed_cld)
     308    CALL getin_p('snow_fallspeed',snow_fallspeed)
     309    snow_fallspeed_clr=snow_fallspeed
     310    snow_fallspeed_cld=snow_fallspeed
     311    CALL getin_p('snow_fallspeed_clr',snow_fallspeed_clr)
     312    CALL getin_p('snow_fallspeed_cld',snow_fallspeed_cld)
    269313
    270314
     
    289333    WRITE(lunout,*) 'lscp_ini, cld_lc_lsc', cld_lc_lsc
    290334    WRITE(lunout,*) 'lscp_ini, cld_lc_con', cld_lc_con
     335    WRITE(lunout,*) 'lscp_ini, cld_lc_lsc_snow', cld_lc_lsc_snow
     336    WRITE(lunout,*) 'lscp_ini, cld_lc_con_snow', cld_lc_con_snow
    291337    WRITE(lunout,*) 'lscp_ini, cld_tau_lsc', cld_tau_lsc
    292338    WRITE(lunout,*) 'lscp_ini, cld_tau_con', cld_tau_con
     
    296342    WRITE(lunout,*) 'lscp_ini, ffallv_con', ffallv_con
    297343    WRITE(lunout,*) 'lscp_ini, coef_eva', coef_eva
    298     WRITE(lunout,*) 'lscp_ini, coef_eva_i', coef_eva_i
     344    WRITE(lunout,*) 'lscp_ini, coef_sub', coef_sub
    299345    WRITE(lunout,*) 'lscp_ini, expo_eva', expo_eva
    300     WRITE(lunout,*) 'lscp_ini, expo_eva_i', expo_eva_i
     346    WRITE(lunout,*) 'lscp_ini, expo_sub', expo_sub
    301347    WRITE(lunout,*) 'lscp_ini, iflag_autoconversion', iflag_autoconversion
    302348    WRITE(lunout,*) 'lscp_ini, dist_liq', dist_liq
     
    309355    WRITE(lunout,*) 'lscp_ini, gamma_col:', gamma_col
    310356    WRITE(lunout,*) 'lscp_ini, gamma_rim:', gamma_rim
     357    WRITE(lunout,*) 'lscp_ini, gamma_freez:', gamma_freez
     358    WRITE(lunout,*) 'lscp_ini, r_snow:', r_snow
     359    WRITE(lunout,*) 'lscp_ini, rain_fallspeed_clr:', rain_fallspeed_clr
     360    WRITE(lunout,*) 'lscp_ini, rain_fallspeed_cld:', rain_fallspeed_cld
     361    WRITE(lunout,*) 'lscp_ini, snow_fallspeed_clr:', snow_fallspeed_clr
     362    WRITE(lunout,*) 'lscp_ini, snow_fallspeed_cld:', snow_fallspeed_cld
    311363
    312364
     
    341393    ENDIF
    342394
     395    !--Initialisations of constants depending on other constants
     396    !--rho_snow formula from r_snow (Brandes et al. 2007 - JAMC)
     397    rho_snow = 1.e3 * 0.178 * (r_snow * 2. * 1000.)**(-0.922)
     398
    343399
    344400    !AA Temporary initialisation
Note: See TracChangeset for help on using the changeset viewer.