[4380] | 1 | module lscp_ini_mod |
---|
| 2 | |
---|
| 3 | implicit none |
---|
| 4 | |
---|
| 5 | ! PARAMETERS for lscp: |
---|
| 6 | !-------------------- |
---|
| 7 | |
---|
| 8 | REAL, SAVE :: seuil_neb=0.001 ! cloud fraction threshold: a cloud really exists when exceeded |
---|
| 9 | !$OMP THREADPRIVATE(seuil_neb) |
---|
| 10 | |
---|
| 11 | INTEGER, SAVE :: ninter=5 ! number of iterations to calculate autoconversion to precipitation |
---|
| 12 | !$OMP THREADPRIVATE(ninter) |
---|
| 13 | |
---|
| 14 | INTEGER,SAVE :: iflag_evap_prec=1 ! precipitation evaporation flag. 0: nothing, 1: "old way", |
---|
| 15 | ! 2: Max cloud fraction above to calculate the max of reevaporation |
---|
| 16 | ! 4: LTP'method i.e. evaporation in the clear-sky fraction of the mesh only |
---|
| 17 | !$OMP THREADPRIVATE(iflag_evap_prec) |
---|
| 18 | |
---|
| 19 | REAL t_coup ! temperature threshold which determines the phase |
---|
| 20 | PARAMETER (t_coup=234.0) ! for which the saturation vapor pressure is calculated |
---|
| 21 | |
---|
| 22 | REAL DDT0 ! iteration parameter |
---|
| 23 | PARAMETER (DDT0=.01) |
---|
| 24 | |
---|
| 25 | REAL ztfondue ! parameter to calculate melting fraction of precipitation |
---|
| 26 | PARAMETER (ztfondue=278.15) |
---|
| 27 | |
---|
| 28 | REAL, SAVE :: rain_int_min=0.001 ! Minimum local rain intensity [mm/s] before the decrease in associated precipitation fraction |
---|
| 29 | !$OMP THREADPRIVATE(rain_int_min) |
---|
| 30 | |
---|
| 31 | REAL, SAVE :: a_tr_sca(4) ! Variables for tracers temporary: alpha parameter for scavenging, 4 possible scavenging processes |
---|
| 32 | !$OMP THREADPRIVATE(a_tr_sca) |
---|
| 33 | |
---|
| 34 | INTEGER, SAVE :: iflag_mpc_bl=0 ! flag to activate boundary layer mixed phase cloud param |
---|
| 35 | !$OMP THREADPRIVATE(iflag_mpc_bl) |
---|
| 36 | |
---|
[4412] | 37 | LOGICAL, SAVE :: ok_radocond_snow=.false. ! take into account the mass of ice precip in the cloud ice content seen by radiation |
---|
| 38 | !$OMP THREADPRIVATE(ok_radocond_snow) |
---|
[4380] | 39 | |
---|
[4420] | 40 | LOGICAL, SAVE :: ok_debug_autoconversion=.true. ! removes a bug in the autoconversion process |
---|
| 41 | !$OMP THREADPRIVATE(ok_debug_autoconversion) |
---|
[4380] | 42 | |
---|
[4420] | 43 | |
---|
[4380] | 44 | CONTAINS |
---|
| 45 | |
---|
| 46 | SUBROUTINE lscp_ini(dtime,ok_ice_sursat) |
---|
| 47 | |
---|
| 48 | |
---|
| 49 | USE ioipsl_getin_p_mod, ONLY : getin_p |
---|
| 50 | USE print_control_mod, ONLY: prt_level, lunout |
---|
| 51 | USE ice_sursat_mod, ONLY: ice_sursat_init |
---|
| 52 | |
---|
| 53 | REAL, INTENT(IN) :: dtime |
---|
| 54 | LOGICAL, INTENT(IN) :: ok_ice_sursat |
---|
| 55 | |
---|
| 56 | CALL getin_p('ninter',ninter) |
---|
| 57 | CALL getin_p('iflag_evap_prec',iflag_evap_prec) |
---|
| 58 | CALL getin_p('seuil_neb',seuil_neb) |
---|
| 59 | CALL getin_p('rain_int_min',rain_int_min) |
---|
| 60 | CALL getin_p('iflag_mpc_bl',iflag_mpc_bl) |
---|
[4420] | 61 | CALL getin_p('ok_radocond_snow',ok_radocond_snow) |
---|
| 62 | CALL getin_p('ok_debug_autoconversion',ok_debug_autoconversion) |
---|
[4380] | 63 | WRITE(lunout,*) 'lscp, ninter:', ninter |
---|
| 64 | WRITE(lunout,*) 'lscp, iflag_evap_prec:', iflag_evap_prec |
---|
| 65 | WRITE(lunout,*) 'lscp, seuil_neb:', seuil_neb |
---|
| 66 | WRITE(lunout,*) 'lscp, rain_int_min:', rain_int_min |
---|
| 67 | WRITE(lunout,*) 'lscp, iflag_mpc_bl:', iflag_mpc_bl |
---|
[4412] | 68 | WRITE(lunout,*) 'lscp, ok_radocond_snow:', ok_radocond_snow |
---|
[4420] | 69 | WRITE(lunout,*) 'lscp, ok_debug_autoconversion:', ok_debug_autoconversion |
---|
| 70 | |
---|
[4380] | 71 | ! check for precipitation sub-time steps |
---|
| 72 | IF (ABS(dtime/REAL(ninter)-360.0).GT.0.001) THEN |
---|
| 73 | WRITE(lunout,*) 'lscp: it is not expected, see Z.X.Li', dtime |
---|
| 74 | WRITE(lunout,*) 'I would prefer a 6 min sub-timestep' |
---|
| 75 | ENDIF |
---|
| 76 | |
---|
| 77 | |
---|
| 78 | !AA Temporary initialisation |
---|
| 79 | a_tr_sca(1) = -0.5 |
---|
| 80 | a_tr_sca(2) = -0.5 |
---|
| 81 | a_tr_sca(3) = -0.5 |
---|
| 82 | a_tr_sca(4) = -0.5 |
---|
| 83 | |
---|
| 84 | IF (ok_ice_sursat) CALL ice_sursat_init() |
---|
| 85 | |
---|
| 86 | |
---|
| 87 | |
---|
| 88 | end subroutine lscp_ini |
---|
| 89 | |
---|
| 90 | end module lscp_ini_mod |
---|