Changeset 4835 for LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_ini.F90
- Timestamp:
- Feb 29, 2024, 7:42:12 PM (3 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/trunk/libf/phylmd/lmdz_blowing_snow_ini.F90
r4821 r4835 3 3 implicit none 4 4 5 real, save, protected :: RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG 6 real, save, protected :: coef_eva_bs, expo_eva_bs, fallv_bs, zeta_bs 7 real, save, protected :: prt_bs, pbst_bs, qbst_bs 8 9 integer, save, protected :: iflag_saltation_bs 5 real, save, protected :: RCPD, RV, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG, RPI 6 real, save, protected :: coef_sub_bs, fallv_bs, zeta_bs, c_esalt_bs 7 real, save, protected :: prt_bs, pbst_bs, qbst_bs, r_bs 8 integer, save, protected :: iflag_saltation_bs, iflag_sedim_bs, iflag_sublim_bs 10 9 11 !$OMP THREADPRIVATE(RCPD, R LSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG)12 !$OMP THREADPRIVATE(coef_ eva_bs, expo_eva_bs, fallv_bs, zeta_bs)10 !$OMP THREADPRIVATE(RCPD, RV, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG, RPI) 11 !$OMP THREADPRIVATE(coef_sub_bs, fallv_bs, r_bs, zeta_bs, c_esalt_bs) 13 12 !$OMP THREADPRIVATE(pbst_bs, prt_bs, qbst_bs) 14 !$OMP THREADPRIVATE(iflag_s altation_bs)13 !$OMP THREADPRIVATE(iflag_sedim_bs, iflag_sublim_bs) 15 14 16 real, save, protected :: tbsmelt=278.15 ! parameter to calculate melting fraction of BS sedimentation 17 real, save, protected :: taumeltbs0=1800.0 ! Melting time scale of blowing snow at 273.15K 18 real, save, protected :: qbsmin=1.E-10 ! Minimum blowing snow specific content 15 real, save, protected :: tbsmelt=278.15 ! parameter to calculate melting fraction of BS sedimentation 16 real, save, protected :: taumeltbs0=600.0 ! Melting time scale of blowing snow at 273.15K 17 real, save, protected :: qbmin=1.E-10 ! Minimum blowing snow specific content 18 !$OMP THREADPRIVATE(tbsmelt, taumeltbs0, qbmin) 19 19 20 !$OMP THREADPRIVATE(tbsmelt, taumeltbs0, qbsmin) 20 real, save, protected :: tau_dens0_bs=864000. ! 10 days by default, in s 21 real, save, protected :: tau_densmin_bs= 21600. ! 1/4 days according to in situ obs by C. Amory during blowing snow + 22 ! Marshall et al. 1999 (snow densification during rain) 23 real, save, protected :: tau_eqsalt_bs= 10. ! Nemoto and Nishimura 2004 show that steady-state saltation is achieved within a time tau_eqsalt of about 10s 24 real, save, protected :: rhofresh_bs = 300.0 ! fresh snow density kg/m3 25 real, save, protected :: rhohard_bs = 450.0 ! hard snow density kg/m3 26 real, save, protected :: rhoice_bs = 920.0 ! ice density kg/m3 27 real, save, protected :: rhobs=900.0 ! blowing snow density (kg/m3) following Bintanja et al. 2001 part I 28 !$OMP THREADPRIVATE(rhoice_bs, rhofresh_bs, rhohard_bs, tau_dens0_bs, tau_densmin_bs, tau_eqsalt_bs, rhobs) 21 29 22 30 … … 24 32 25 33 subroutine blowing_snow_ini(RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in,& 26 RVTMP2_in, RTT_in,RD_in,RG_in )34 RVTMP2_in, RTT_in,RD_in,RG_in, RV_in, RPI_in) 27 35 28 36 USE ioipsl_getin_p_mod, ONLY : getin_p 29 37 30 real, intent(in) :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in 31 real, intent(in) :: RVTMP2_in, RTT_in, RD_in, RG_in 38 real, intent(in) :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in, RPI_in 39 real, intent(in) :: RVTMP2_in, RTT_in, RD_in, RG_in, RV_in 32 40 33 41 34 42 RG=RG_in 35 43 RD=RD_in 44 RV=RV_in 36 45 RCPD=RCPD_in 37 46 RLVTT=RLVTT_in … … 41 50 RG=RG_in 42 51 RVTMP2=RVTMP2_in 52 RPI=RPI_in 43 53 54 c_esalt_bs= 3.25 55 CALL getin_p('c_esalt_bs',c_esalt_bs) 44 56 45 57 qbst_bs= 0.001 … … 58 70 CALL getin_p('fallv_bs',fallv_bs) 59 71 60 coef_ eva_bs = 2.e-561 CALL getin_p('coef_ eva_bs',coef_eva_bs)72 coef_sub_bs = 0.1 73 CALL getin_p('coef_sub_bs',coef_sub_bs) 62 74 63 expo_eva_bs = 0.564 CALL getin_p(' expo_eva_bs',expo_eva_bs)75 iflag_sublim_bs=1 76 CALL getin_p('iflag_sublim_bs',iflag_sublim_bs) 65 77 66 iflag_s altation_bs=067 CALL getin_p('iflag_s altation_bs',iflag_saltation_bs)78 iflag_sedim_bs=1 79 CALL getin_p('iflag_sedim_bs',iflag_sedim_bs) 68 80 81 r_bs=150.0e-6 82 CALL getin_p('r_bs',r_bs) 69 83 70 84 end subroutine blowing_snow_ini
Note: See TracChangeset
for help on using the changeset viewer.