Changeset 4830 for LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90
- Timestamp:
- Feb 22, 2024, 5:29:02 PM (7 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/lmdz_lscp_ini.F90
r4818 r4830 112 112 !$OMP THREADPRIVATE(coef_eva) 113 113 114 REAL, SAVE, PROTECTED :: coef_ eva_i! tuning coefficient ice precip sublimation115 !$OMP THREADPRIVATE(coef_ eva_i)114 REAL, SAVE, PROTECTED :: coef_sub ! tuning coefficient ice precip sublimation 115 !$OMP THREADPRIVATE(coef_sub) 116 116 117 117 REAL, SAVE, PROTECTED :: expo_eva=0.5 ! tuning coefficient liquid precip evaporation 118 118 !$OMP THREADPRIVATE(expo_eva) 119 119 120 REAL, SAVE, PROTECTED :: expo_ eva_i! tuning coefficient ice precip sublimation121 !$OMP THREADPRIVATE(expo_ eva_i)120 REAL, SAVE, PROTECTED :: expo_sub ! tuning coefficient ice precip sublimation 121 !$OMP THREADPRIVATE(expo_sub) 122 122 123 123 REAL, SAVE, PROTECTED :: cice_velo=1.645 ! factor in the ice fall velocity formulation … … 137 137 !$OMP THREADPRIVATE(ok_poprecip) 138 138 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 139 145 REAL, SAVE, PROTECTED :: rain_int_min=0.001 ! Minimum local rain intensity [mm/s] before the decrease in associated precipitation fraction 140 146 !$OMP THREADPRIVATE(rain_int_min) … … 155 161 !$OMP THREADPRIVATE(rho_rain) 156 162 157 REAL, SAVE, PROTECTED :: rho_snow =500.! A COMMENTER TODO [kg/m3]163 REAL, SAVE, PROTECTED :: rho_snow ! A COMMENTER TODO [kg/m3] 158 164 !$OMP THREADPRIVATE(rho_snow) 159 165 … … 164 170 !$OMP THREADPRIVATE(r_snow) 165 171 166 REAL, SAVE, PROTECTED :: tau_auto_snow_min=1 800. ! A COMMENTER TODO [s]172 REAL, SAVE, PROTECTED :: tau_auto_snow_min=100. ! A COMMENTER TODO [s] 167 173 !$OMP THREADPRIVATE(tau_auto_snow_min) 168 174 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] 170 176 !$OMP THREADPRIVATE(tau_auto_snow_max) 171 177 … … 182 188 !$OMP THREADPRIVATE(alpha_freez) 183 189 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 [-] 185 194 !$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) 186 213 !--End of the parameters for poprecip 187 214 … … 246 273 CALL getin_p('cld_lc_lsc',cld_lc_lsc) 247 274 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) 248 279 CALL getin_p('cld_tau_lsc',cld_tau_lsc) 249 280 CALL getin_p('cld_tau_con',cld_tau_con) … … 253 284 CALL getin_p('ffallv_lsc',ffallv_con) 254 285 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) 257 289 CALL getin_p('expo_eva',expo_eva) 258 expo_ eva_i=expo_eva259 CALL getin_p('expo_ eva_i',expo_eva_i)290 expo_sub=expo_eva 291 CALL getin_p('expo_sub',expo_sub) 260 292 CALL getin_p('iflag_autoconversion',iflag_autoconversion) 261 293 CALL getin_p('dist_liq',dist_liq) … … 267 299 CALL getin_p('gamma_col',gamma_col) 268 300 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) 269 313 270 314 … … 289 333 WRITE(lunout,*) 'lscp_ini, cld_lc_lsc', cld_lc_lsc 290 334 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 291 337 WRITE(lunout,*) 'lscp_ini, cld_tau_lsc', cld_tau_lsc 292 338 WRITE(lunout,*) 'lscp_ini, cld_tau_con', cld_tau_con … … 296 342 WRITE(lunout,*) 'lscp_ini, ffallv_con', ffallv_con 297 343 WRITE(lunout,*) 'lscp_ini, coef_eva', coef_eva 298 WRITE(lunout,*) 'lscp_ini, coef_ eva_i', coef_eva_i344 WRITE(lunout,*) 'lscp_ini, coef_sub', coef_sub 299 345 WRITE(lunout,*) 'lscp_ini, expo_eva', expo_eva 300 WRITE(lunout,*) 'lscp_ini, expo_ eva_i', expo_eva_i346 WRITE(lunout,*) 'lscp_ini, expo_sub', expo_sub 301 347 WRITE(lunout,*) 'lscp_ini, iflag_autoconversion', iflag_autoconversion 302 348 WRITE(lunout,*) 'lscp_ini, dist_liq', dist_liq … … 309 355 WRITE(lunout,*) 'lscp_ini, gamma_col:', gamma_col 310 356 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 311 363 312 364 … … 341 393 ENDIF 342 394 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 343 399 344 400 !AA Temporary initialisation
Note: See TracChangeset
for help on using the changeset viewer.