1 |
|
|
module blowing_snow_ini_mod |
2 |
|
|
|
3 |
|
|
implicit none |
4 |
|
|
|
5 |
|
|
save |
6 |
|
|
integer :: prt_level=0,lunout |
7 |
|
|
real RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG |
8 |
|
|
real coef_eva_bs, expo_eva_bs, fallv_bs, zeta_bs |
9 |
|
|
real prt_bs, pbst_bs, qbst_bs |
10 |
|
|
integer :: niter_bs, iflag_saltation_bs |
11 |
|
|
|
12 |
|
|
!$OMP THREADPRIVATE(prt_level, lunout) |
13 |
|
|
!$OMP THREADPRIVATE(RCPD, RLSTT, RLVTT, RLMLT, RVTMP2, RTT,RD,RG) |
14 |
|
|
!$OMP THREADPRIVATE(coef_eva_bs, expo_eva_bs, fallv_bs, zeta_bs) |
15 |
|
|
!$OMP THREADPRIVATE(pbst_bs, prt_bs, qbst_bs) |
16 |
|
|
!$OMP THREADPRIVATE(niter_bs, iflag_saltation_bs) |
17 |
|
|
|
18 |
|
|
contains |
19 |
|
|
|
20 |
|
1 |
subroutine blowing_snow_ini(prt_level_in,lunout_in, & |
21 |
|
|
RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in,& |
22 |
|
|
RVTMP2_in, RTT_in,RD_in,RG_in) |
23 |
|
|
|
24 |
|
|
USE ioipsl_getin_p_mod, ONLY : getin_p |
25 |
|
|
|
26 |
|
|
integer, intent(in) :: prt_level_in,lunout_in |
27 |
|
|
real, intent(in) :: RCPD_in, RLSTT_in, RLVTT_in, RLMLT_in |
28 |
|
|
real, intent(in) :: RVTMP2_in, RTT_in, RD_in, RG_in |
29 |
|
|
|
30 |
|
|
|
31 |
|
1 |
prt_level=prt_level_in |
32 |
|
1 |
RG=RG_in |
33 |
|
1 |
RD=RD_in |
34 |
|
1 |
RCPD=RCPD_in |
35 |
|
1 |
RLVTT=RLVTT_in |
36 |
|
1 |
RLSTT=RLSTT_in |
37 |
|
1 |
RLMLT=RLMLT_in |
38 |
|
1 |
RTT=RTT_in |
39 |
|
|
RG=RG_in |
40 |
|
1 |
RVTMP2=RVTMP2_in |
41 |
|
1 |
lunout=lunout_in |
42 |
|
|
|
43 |
|
|
|
44 |
|
1 |
qbst_bs= 0.001 |
45 |
|
1 |
CALL getin_p('qbst_bs',qbst_bs) |
46 |
|
|
|
47 |
|
1 |
pbst_bs= 0.01 |
48 |
|
1 |
CALL getin_p('pbst_bs',pbst_bs) |
49 |
|
|
|
50 |
|
1 |
prt_bs= 0.0003 |
51 |
|
1 |
CALL getin_p('prt_bs',prt_bs) |
52 |
|
|
|
53 |
|
|
|
54 |
|
1 |
zeta_bs= 3. |
55 |
|
1 |
CALL getin_p('zeta_bs',zeta_bs) |
56 |
|
|
|
57 |
|
1 |
fallv_bs = 0.1 |
58 |
|
1 |
CALL getin_p('fallv_bs',fallv_bs) |
59 |
|
|
|
60 |
|
1 |
coef_eva_bs = 2.e-5 |
61 |
|
1 |
CALL getin_p('coef_eva_bs',coef_eva_bs) |
62 |
|
|
|
63 |
|
1 |
expo_eva_bs = 0.5 |
64 |
|
1 |
CALL getin_p('expo_eva_bs',expo_eva_bs) |
65 |
|
|
|
66 |
|
1 |
niter_bs = 5 |
67 |
|
1 |
CALL getin_p('niter_bs',niter_bs) |
68 |
|
|
|
69 |
|
1 |
iflag_saltation_bs=0 |
70 |
|
1 |
CALL getin_p('iflag_saltation_bs',iflag_saltation_bs) |
71 |
|
|
|
72 |
|
|
|
73 |
|
1 |
end subroutine blowing_snow_ini |
74 |
|
|
|
75 |
|
|
end module blowing_snow_ini_mod |