source: LMDZ6/trunk/libf/phylmd/blowing_snow_ini_mod.F90 @ 4608

Last change on this file since 4608 was 4529, checked in by evignon, 19 months ago

nouvelle formulation du flux d'emission de neige soufflee

File size: 1.9 KB
Line 
1module blowing_snow_ini_mod
2
3implicit none
4
5save
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      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         prt_level=prt_level_in
32         RG=RG_in
33         RD=RD_in
34         RCPD=RCPD_in
35         RLVTT=RLVTT_in
36         RLSTT=RLSTT_in
37         RLMLT=RLMLT_in
38         RTT=RTT_in
39         RG=RG_in
40         RVTMP2=RVTMP2_in
41         lunout=lunout_in
42
43
44         qbst_bs= 0.001
45         CALL getin_p('qbst_bs',qbst_bs)
46
47         pbst_bs= 0.01
48         CALL getin_p('pbst_bs',pbst_bs)
49
50         prt_bs= 0.0003
51         CALL getin_p('prt_bs',prt_bs)
52
53
54         zeta_bs= 3.
55         CALL getin_p('zeta_bs',zeta_bs)
56
57         fallv_bs = 0.1
58         CALL getin_p('fallv_bs',fallv_bs)
59
60         coef_eva_bs =  2.e-5
61         CALL getin_p('coef_eva_bs',coef_eva_bs)
62
63         expo_eva_bs = 0.5
64         CALL getin_p('expo_eva_bs',expo_eva_bs)
65
66         niter_bs = 5
67         CALL getin_p('niter_bs',niter_bs)
68
69         iflag_saltation_bs=0
70         CALL getin_p('iflag_saltation_bs',iflag_saltation_bs)
71
72
73      end subroutine blowing_snow_ini
74
75end module blowing_snow_ini_mod
Note: See TracBrowser for help on using the repository browser.