Changeset 8506
- Timestamp:
- 2017-09-07T12:28:42+02:00 (7 years ago)
- Location:
- branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3
- Files:
-
- 4 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r8505 r8506 176 176 INTEGER , PUBLIC :: jiceprt !: debug j-point 177 177 178 ! !!** ice-init namelist (namiceini) **179 ! -- iceist -- !180 LOGICAL , PUBLIC :: ln_limini ! initialization or not181 LOGICAL , PUBLIC :: ln_limini_file ! Ice initialization state from 2D netcdf file182 REAL(wp), PUBLIC :: rn_thres_sst ! threshold water temperature for initial sea ice183 REAL(wp), PUBLIC :: rn_hts_ini_n ! initial snow thickness in the north184 REAL(wp), PUBLIC :: rn_hts_ini_s ! initial snow thickness in the south185 REAL(wp), PUBLIC :: rn_hti_ini_n ! initial ice thickness in the north186 REAL(wp), PUBLIC :: rn_hti_ini_s ! initial ice thickness in the south187 REAL(wp), PUBLIC :: rn_ati_ini_n ! initial leads area in the north188 REAL(wp), PUBLIC :: rn_ati_ini_s ! initial leads area in the south189 REAL(wp), PUBLIC :: rn_smi_ini_n ! initial salinity190 REAL(wp), PUBLIC :: rn_smi_ini_s ! initial salinity191 REAL(wp), PUBLIC :: rn_tmi_ini_n ! initial temperature192 REAL(wp), PUBLIC :: rn_tmi_ini_s ! initial temperature193 194 178 ! !!** ice-dynamics namelist (namicedyn) ** 195 179 ! -- iceadv -- ! … … 252 236 REAL(wp), PUBLIC :: rn_simin !: minimum ice salinity [PSU] 253 237 254 ! !!** ice-mechanical redistribution namelist (namiceitdme)255 REAL(wp), PUBLIC :: rn_cs !: fraction of shearing energy contributing to ridging256 INTEGER , PUBLIC :: nn_partfun !: participation function: =0 Thorndike et al. (1975), =1 Lipscomb et al. (2007)257 REAL(wp), PUBLIC :: rn_gstar !: fractional area of young ice contributing to ridging258 REAL(wp), PUBLIC :: rn_astar !: equivalent of G* for an exponential participation function259 LOGICAL , PUBLIC :: ln_ridging !: ridging of ice or not260 REAL(wp), PUBLIC :: rn_hstar !: thickness that determines the maximal thickness of ridged ice261 REAL(wp), PUBLIC :: rn_por_rdg !: initial porosity of ridges (0.3 regular value)262 REAL(wp), PUBLIC :: rn_fsnowrdg !: fractional snow loss to the ocean during ridging263 REAL(wp), PUBLIC :: rn_fpondrdg !: fractional melt pond loss to the ocean during ridging264 LOGICAL , PUBLIC :: ln_rafting !: rafting of ice or not265 REAL(wp), PUBLIC :: rn_hraft !: threshold thickness (m) for rafting / ridging266 REAL(wp), PUBLIC :: rn_craft !: coefficient for smoothness of the hyperbolic tangent in rafting267 REAL(wp), PUBLIC :: rn_fsnowrft !: fractional snow loss to the ocean during ridging268 REAL(wp), PUBLIC :: rn_fpondrft !: fractional snow loss to the ocean during rafting269 270 238 ! MV MP 2016 271 239 ! !!** melt pond namelist (namicemp) -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceist.F90
r8498 r8506 46 46 INTEGER , PARAMETER :: jp_smi = 6 ! index of ice sali at T-point 47 47 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: si ! structure of input fields (file informations, fields read) 48 ! 49 ! ** namelist (namiceini) ** 50 LOGICAL :: ln_limini ! initialization or not 51 LOGICAL :: ln_limini_file ! Ice initialization state from 2D netcdf file 52 REAL(wp) :: rn_thres_sst ! threshold water temperature for initial sea ice 53 REAL(wp) :: rn_hts_ini_n ! initial snow thickness in the north 54 REAL(wp) :: rn_hts_ini_s ! initial snow thickness in the south 55 REAL(wp) :: rn_hti_ini_n ! initial ice thickness in the north 56 REAL(wp) :: rn_hti_ini_s ! initial ice thickness in the south 57 REAL(wp) :: rn_ati_ini_n ! initial leads area in the north 58 REAL(wp) :: rn_ati_ini_s ! initial leads area in the south 59 REAL(wp) :: rn_smi_ini_n ! initial salinity 60 REAL(wp) :: rn_smi_ini_s ! initial salinity 61 REAL(wp) :: rn_tmi_ini_n ! initial temperature 62 REAL(wp) :: rn_tmi_ini_s ! initial temperature 63 48 64 !!---------------------------------------------------------------------- 49 65 !! NEMO/ICE 4.0 , NEMO Consortium (2017) -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icerdgrft.F90
r8500 r8506 36 36 PUBLIC ice_rdgrft_alloc ! called by ice_init 37 37 38 !-----------------------------------------------------------------------39 38 ! Variables shared among ridging subroutines 40 !-----------------------------------------------------------------------41 39 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: asum ! sum of total ice and open water area 42 40 REAL(wp), ALLOCATABLE, DIMENSION(:,:) :: aksum ! ratio of area removed to area ridged … … 51 49 REAL(wp), PARAMETER :: krdgmin = 1.1_wp ! min ridge thickness multiplier 52 50 REAL(wp), PARAMETER :: kraft = 0.5_wp ! rafting multipliyer 53 !54 51 REAL(wp) :: zdrho ! 55 52 ! 53 ! ** namelist (namiceitdme) ** 54 REAL(wp) :: rn_cs !: fraction of shearing energy contributing to ridging 55 INTEGER :: nn_partfun !: participation function: =0 Thorndike et al. (1975), =1 Lipscomb et al. (2007) 56 REAL(wp) :: rn_gstar !: fractional area of young ice contributing to ridging 57 REAL(wp) :: rn_astar !: equivalent of G* for an exponential participation function 58 LOGICAL :: ln_ridging !: ridging of ice or not 59 REAL(wp) :: rn_hstar !: thickness that determines the maximal thickness of ridged ice 60 REAL(wp) :: rn_por_rdg !: initial porosity of ridges (0.3 regular value) 61 REAL(wp) :: rn_fsnowrdg !: fractional snow loss to the ocean during ridging 62 REAL(wp) :: rn_fpondrdg !: fractional melt pond loss to the ocean during ridging 63 LOGICAL :: ln_rafting !: rafting of ice or not 64 REAL(wp) :: rn_hraft !: threshold thickness (m) for rafting / ridging 65 REAL(wp) :: rn_craft !: coefficient for smoothness of the hyperbolic tangent in rafting 66 REAL(wp) :: rn_fsnowrft !: fractional snow loss to the ocean during ridging 67 REAL(wp) :: rn_fpondrft !: fractional snow loss to the ocean during rafting 56 68 ! 57 69 !!---------------------------------------------------------------------- -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icethd.F90
r8498 r8506 448 448 END DO 449 449 ! 450 ! Change thickness to volume 450 ! Change thickness to volume (replaces routine ice_var_eqv2glo) 451 451 v_i_1d(1:nidx) = ht_i_1d(1:nidx) * a_i_1d(1:nidx) 452 452 v_s_1d(1:nidx) = ht_s_1d(1:nidx) * a_i_1d(1:nidx)
Note: See TracChangeset
for help on using the changeset viewer.