Changeset 8422 for branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3
- Timestamp:
- 2017-08-08T15:58:05+02:00 (7 years ago)
- Location:
- branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3
- Files:
-
- 8 added
- 8 deleted
- 6 edited
-
ice.F90 (modified) (6 diffs)
-
iceerr1.F90 (modified) (2 diffs)
-
iceerr2.F90 (modified) (2 diffs)
-
iceistate.F90 (modified) (1 diff)
-
iceitd.F90 (added)
-
icestp.F90 (modified) (4 diffs)
-
icethd.F90 (added)
-
icethd_da.F90 (added)
-
icethd_dh.F90 (added)
-
icethd_dif.F90 (added)
-
icethd_ent.F90 (added)
-
icethd_lac.F90 (added)
-
icethd_sal.F90 (added)
-
limitd_th.F90 (deleted)
-
limthd.F90 (deleted)
-
limthd_da.F90 (deleted)
-
limthd_dh.F90 (deleted)
-
limthd_dif.F90 (deleted)
-
limthd_ent.F90 (deleted)
-
limthd_lac.F90 (deleted)
-
limthd_sal.F90 (deleted)
-
limvar.F90 (modified) (2 diffs)
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/ice.F90
r8413 r8422 77 77 !! - ice_adv : ice advection | 78 78 !! - ice_rdgrft : ice ridging and rafting | 79 !! - lim_thd : ice halo-thermodynamics | 80 !! - lim_itd_th : thermodynamic changes in ice thickness distribution | 81 !! and creation of new ice | 79 !! - ice_thd : ice halo-thermodynamics and creation of new ice | 80 !! - ice_itd : thermodynamic changes in ice thickness distribution | 82 81 !! | 83 82 !! See the associated routines for more information | … … 150 149 !! bv_i ! - ! relative brine volume | ??? | 151 150 !!===================================================================== 152 153 LOGICAL, PUBLIC :: con_i = .false. ! switch for conservation test154 151 155 152 !!-------------------------------------------------------------------------- … … 221 218 222 219 ! !!** ice-thermodynamics namelist (namicethd) ** 223 ! -- limthd_dif -- !220 ! -- icethd_dif -- ! 224 221 REAL(wp), PUBLIC :: rn_kappa_i !: coef. for the extinction of radiation Grenfell et al. (2006) [1/m] 225 222 INTEGER , PUBLIC :: nn_ice_thcon !: thermal conductivity: =0 Untersteiner (1964) ; =1 Pringle et al (2007) … … 227 224 INTEGER , PUBLIC :: nn_monocat !: virtual ITD mono-category parameterizations (1-4) or not (0) 228 225 REAL(wp), PUBLIC :: rn_cdsn !: thermal conductivity of the snow [W/m/K] 229 ! -- limthd_dh -- !226 ! -- icethd_dh -- ! 230 227 LOGICAL , PUBLIC :: ln_limdH !: activate ice thickness change from growing/melting (T) or not (F) 231 228 REAL(wp), PUBLIC :: rn_betas !: coef. for partitioning of snowfall between leads and sea ice 232 ! -- limthd_da -- !229 ! -- icethd_da -- ! 233 230 LOGICAL , PUBLIC :: ln_limdA !: activate lateral melting param. (T) or not (F) 234 231 REAL(wp), PUBLIC :: rn_beta !: coef. beta for lateral melting param. 235 232 REAL(wp), PUBLIC :: rn_dmin !: minimum floe diameter for lateral melting param. 236 ! -- limthd_lac -- !233 ! -- icethd_lac -- ! 237 234 LOGICAL , PUBLIC :: ln_limdO !: activate ice growth in open-water (T) or not (F) 238 235 REAL(wp), PUBLIC :: rn_hnewice !: thickness for new ice formation (m) … … 241 238 REAL(wp), PUBLIC :: rn_vfrazb !: threshold drift speed for collection of bottom frazil ice 242 239 REAL(wp), PUBLIC :: rn_Cfrazb !: squeezing coefficient for collection of bottom frazil ice 243 ! -- limitd_th-- !240 ! -- iceitd -- ! 244 241 REAL(wp), PUBLIC :: rn_himin !: minimum ice thickness 245 242 … … 366 363 367 364 ! heat flux associated with ice-ocean mass exchange 368 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hfx_thd !: ice-ocean heat flux from thermo processes ( limthd_dh) [W.m-2]365 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hfx_thd !: ice-ocean heat flux from thermo processes (icethd_dh) [W.m-2] 369 366 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hfx_dyn !: ice-ocean heat flux from ridging [W.m-2] 370 367 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: hfx_res !: residual heat flux due to correction of ice thickness [W.m-2] -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceerr1.F90
r8420 r8422 17 17 USE ice 18 18 USE ice1D ! LIM thermodynamic sea-ice variables 19 USE limitd_th19 USE iceitd 20 20 USE limvar 21 21 USE icecons ! conservation tests … … 104 104 ! Rebin categories with thickness out of bounds 105 105 !---------------------------------------------------- 106 IF ( jpl > 1 ) CALL lim_itd_th_reb106 IF ( jpl > 1 ) CALL ice_itd_reb 107 107 108 108 !----------------- -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceerr2.F90
r8420 r8422 16 16 USE phycst ! physical constants 17 17 USE ice 18 USE ice1D ! LIM thermodynamic sea-ice variables19 USE limitd_th18 USE ice1D ! LIM thermodynamic sea-ice variables 19 USE iceitd 20 20 USE limvar 21 21 USE icecons ! conservation tests … … 118 118 ! Rebin categories with thickness out of bounds 119 119 !---------------------------------------------------- 120 IF ( jpl > 1 ) CALL lim_itd_th_reb120 IF ( jpl > 1 ) CALL ice_itd_reb 121 121 122 122 !----------------- -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/iceistate.F90
r8413 r8422 307 307 ! This case below should not be used if (ht_s/ht_i) is ok in namelist 308 308 ! In case snow load is in excess that would lead to transformation from snow to ice 309 ! Then, transfer the snow excess into the ice (different from limthd_dh)309 ! Then, transfer the snow excess into the ice (different from icethd_dh) 310 310 zdh = MAX( 0._wp, ( rhosn * ht_s(ji,jj,jl) + ( rhoic - rau0 ) * ht_i(ji,jj,jl) ) * r1_rau0 ) 311 311 ! recompute ht_i, ht_s avoiding out of bounds values -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/icestp.F90
r8420 r8422 34 34 USE icerhg ! Ice rheology 35 35 USE iceadv ! Ice advection 36 USE limthd ! Ice thermodynamics36 USE icethd ! Ice thermodynamics 37 37 USE icerdgrft ! Ice ridging/rafting 38 38 USE iceupdate ! sea surface boundary condition … … 48 48 ! END MV MP 2016 49 49 USE iceistate ! LIM initial state 50 USE limthd_sal ! LIM ice thermodynamics: salinity50 USE icethd_sal ! LIM ice thermodynamics: salinity 51 51 ! 52 52 USE c1d ! 1D vertical configuration … … 190 190 !----------------------------! 191 191 ! --- zap this if no ice thermo --- ! 192 IF( ln_limthd ) CALL lim_thd( kt ) ! -- Ice thermodynamics192 IF( ln_limthd ) CALL ice_thd( kt ) ! -- Ice thermodynamics 193 193 194 194 ! MV MP 2016 … … 271 271 CALL ice_itd_init ! ice thickness distribution initialization 272 272 ! 273 CALL lim_thd_init ! set ice thermodynics parameters274 ! 275 CALL lim_thd_sal_init ! set ice salinity parameters273 CALL ice_thd_init ! set ice thermodynics parameters 274 ! 275 CALL ice_thd_sal_init ! set ice salinity parameters 276 276 277 277 ! MV MP 2016 -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/LIM_SRC_3/limvar.F90
r8420 r8422 431 431 SUBROUTINE lim_var_salprof1d 432 432 !!------------------------------------------------------------------- 433 !! *** ROUTINE lim_ thd_salprof1d ***433 !! *** ROUTINE lim_var_salprof1d *** 434 434 !! 435 435 !! ** Purpose : 1d computation of the sea ice salinity profile … … 760 760 zht_s(ji,jl) = zht_i(ji,jl) * ( zhts(ji) / zhti(ji) ) 761 761 ! In case snow load is in excess that would lead to transformation from snow to ice 762 ! Then, transfer the snow excess into the ice (different from limthd_dh)762 ! Then, transfer the snow excess into the ice (different from icethd_dh) 763 763 zdh = MAX( 0._wp, ( rhosn * zht_s(ji,jl) + ( rhoic - rau0 ) * zht_i(ji,jl) ) * r1_rau0 ) 764 764 ! recompute ht_i, ht_s avoiding out of bounds values
Note: See TracChangeset
for help on using the changeset viewer.
