| 1 | ! |
|---|
| 2 | ! $Id: surface_data.f90 5662 2025-05-20 14:24:41Z rkazeroni $ |
|---|
| 3 | ! |
|---|
| 4 | MODULE surface_data |
|---|
| 5 | |
|---|
| 6 | IMPLICIT NONE |
|---|
| 7 | |
|---|
| 8 | REAL, PARAMETER :: calice=1.0/(5.1444e+06*0.15) |
|---|
| 9 | REAL, PARAMETER :: calsno=1./(2.3867e+06*.15) |
|---|
| 10 | |
|---|
| 11 | LOGICAL, SAVE :: ok_veget ! true for use of vegetation model ORCHIDEE |
|---|
| 12 | !$OMP THREADPRIVATE(ok_veget) |
|---|
| 13 | |
|---|
| 14 | CHARACTER(len=10), SAVE :: type_veget ! orchidee/y/bucket/n/betaclim |
|---|
| 15 | !$OMP THREADPRIVATE(type_veget) |
|---|
| 16 | |
|---|
| 17 | CHARACTER(len=6), SAVE :: type_ocean ! force/slab/couple |
|---|
| 18 | !$OMP THREADPRIVATE(type_ocean) |
|---|
| 19 | |
|---|
| 20 | !GG |
|---|
| 21 | INTEGER, SAVE :: iflag_seaice |
|---|
| 22 | !$OMP THREADPRIVATE(iflag_seaice) |
|---|
| 23 | INTEGER, SAVE :: iflag_seaice_alb |
|---|
| 24 | !$OMP THREADPRIVATE(iflag_seaice_alb) |
|---|
| 25 | INTEGER, SAVE :: iflag_leads |
|---|
| 26 | !$OMP THREADPRIVATE(iflag_leads) |
|---|
| 27 | |
|---|
| 28 | !For sea-ice |
|---|
| 29 | REAL,SAVE :: sice_cond |
|---|
| 30 | !$OMP THREADPRIVATE(sice_cond) |
|---|
| 31 | REAL,SAVE :: sisno_cond |
|---|
| 32 | !$OMP THREADPRIVATE(sisno_cond) |
|---|
| 33 | REAL,SAVE :: sisno_den |
|---|
| 34 | !$OMP THREADPRIVATE(sisno_den) |
|---|
| 35 | REAL,SAVE :: sisno_min |
|---|
| 36 | !$OMP THREADPRIVATE(sisno_min) |
|---|
| 37 | REAL,SAVE :: sithick_min |
|---|
| 38 | !$OMP THREADPRIVATE(sithick_min) |
|---|
| 39 | REAL,SAVE :: sisno_wfact |
|---|
| 40 | !$OMP THREADPRIVATE(sisno_wfact) |
|---|
| 41 | REAL,SAVE :: amax_n |
|---|
| 42 | !$OMP THREADPRIVATE(amax_n) |
|---|
| 43 | REAL,SAVE :: amax_s |
|---|
| 44 | !$OMP THREADPRIVATE(amax_s) |
|---|
| 45 | REAL,SAVE :: rn_alb_sdry |
|---|
| 46 | !$OMP THREADPRIVATE(rn_alb_sdry) |
|---|
| 47 | REAL,SAVE :: rn_alb_smlt |
|---|
| 48 | !$OMP THREADPRIVATE(rn_alb_smlt) |
|---|
| 49 | REAL,SAVE :: rn_alb_idry |
|---|
| 50 | !$OMP THREADPRIVATE(rn_alb_idry) |
|---|
| 51 | REAL,SAVE :: rn_alb_imlt |
|---|
| 52 | !$OMP THREADPRIVATE(rn_alb_imlt) |
|---|
| 53 | REAL,SAVE :: si_pen_frac |
|---|
| 54 | !$OMP THREADPRIVATE(si_pen_frac) |
|---|
| 55 | REAL,SAVE :: si_pen_ext |
|---|
| 56 | !$OMP THREADPRIVATE(si_pen_ext) |
|---|
| 57 | REAL,SAVE :: fseaN |
|---|
| 58 | !$OMP THREADPRIVATE(fseaN) |
|---|
| 59 | REAL,SAVE :: fseaS |
|---|
| 60 | !$OMP THREADPRIVATE(fseaS) |
|---|
| 61 | !GG |
|---|
| 62 | |
|---|
| 63 | ! if type_ocean=couple : version_ocean=opa8 ou nemo |
|---|
| 64 | ! if type_ocean=slab : version_ocean=sicOBS or sicINT or sicNO |
|---|
| 65 | CHARACTER(len=6), SAVE :: version_ocean |
|---|
| 66 | !$OMP THREADPRIVATE(version_ocean) |
|---|
| 67 | |
|---|
| 68 | ! Pas de temps couplage atm/oce (en secondes) |
|---|
| 69 | REAL, SAVE :: t_coupl |
|---|
| 70 | !$OMP THREADPRIVATE(t_coupl) |
|---|
| 71 | |
|---|
| 72 | ! FOR INLANDSIS: |
|---|
| 73 | !=============== |
|---|
| 74 | |
|---|
| 75 | ! 1 for coupling with INLANDSIS |
|---|
| 76 | INTEGER, SAVE :: landice_opt ! 1 for coupling with INLANDSIS |
|---|
| 77 | !$OMP THREADPRIVATE(landice_opt) |
|---|
| 78 | |
|---|
| 79 | ! temperature calculation options within the soil and at the surface |
|---|
| 80 | INTEGER, SAVE :: iflag_tsurf_inlandsis,iflag_temp_inlandsis |
|---|
| 81 | !$OMP THREADPRIVATE(iflag_tsurf_inlandsis,iflag_temp_inlandsis) |
|---|
| 82 | |
|---|
| 83 | ! flags for albedo and roughness calc. |
|---|
| 84 | INTEGER, SAVE :: iflag_albcalc,iflag_z0m_snow |
|---|
| 85 | !$OMP THREADPRIVATE(iflag_albcalc,iflag_z0m_snow) |
|---|
| 86 | |
|---|
| 87 | ! with or without snow module/ blowing snow, ascii outfile |
|---|
| 88 | LOGICAL, SAVE :: SnoMod,BloMod,ok_outfor |
|---|
| 89 | !$OMP THREADPRIVATE(SnoMod,BloMod,ok_outfor) |
|---|
| 90 | |
|---|
| 91 | ! activate slush, korlyakov snow density, RN z0h calc. |
|---|
| 92 | LOGICAL, SAVE :: is_ok_slush,is_ok_density_kotlyakov,is_ok_z0h_rn |
|---|
| 93 | !$OMP THREADPRIVATE(is_ok_slush,is_ok_density_kotlyakov,is_ok_z0h_rn) |
|---|
| 94 | |
|---|
| 95 | ! activate detection snow/ice layers and option XF discrtet/option runoff AC |
|---|
| 96 | LOGICAL, SAVE :: ok_zsn_ii,discret_xf,opt_runoff_ac |
|---|
| 97 | !$OMP THREADPRIVATE(ok_zsn_ii,discret_xf, opt_runoff_ac) |
|---|
| 98 | |
|---|
| 99 | ! value of z0m snow when prescribed and albedo correction term |
|---|
| 100 | REAL, SAVE :: prescribed_z0m_snow,correc_alb |
|---|
| 101 | !$OMP THREADPRIVATE(prescribed_z0m_snow, correc_alb) |
|---|
| 102 | |
|---|
| 103 | ! value of sphericity [0-99] and snow grain size [e-4m] for polar buffer snow |
|---|
| 104 | ! layer |
|---|
| 105 | REAL, SAVE :: buf_sph_pol,buf_siz_pol |
|---|
| 106 | !$OMP THREADPRIVATE(buf_sph_pol,buf_siz_pol) |
|---|
| 107 | |
|---|
| 108 | |
|---|
| 109 | |
|---|
| 110 | END MODULE surface_data |
|---|