source: LMDZ6/branches/blowing_snow/libf/phylmd/blowing_snow_ini_mod.F90 @ 5450

Last change on this file since 5450 was 4485, checked in by evignon, 21 months ago

premier commit pour l'ajout de la neige soufflee sur la nouvelle branche

File size: 1.8 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
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)
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
70      end subroutine blowing_snow_ini
71
72end module blowing_snow_ini_mod
Note: See TracBrowser for help on using the repository browser.