Changeset 8422
- Timestamp:
- 2017-08-08T15:58:05+02:00 (7 years ago)
- Location:
- branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO
- Files:
-
- 8 added
- 8 deleted
- 10 edited
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 -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/BDY/bdyice.F90
r8411 r8422 119 119 ! ----------------- 120 120 ! In case a) snow load would be in excess or b) ice is coming into a warmer environment that would lead to 121 ! very large transformation from snow to ice (see limthd_dh.F90)122 123 ! Then, a) transfer the snow excess into the ice (different from limthd_dh)121 ! very large transformation from snow to ice (see icethd_dh.F90) 122 123 ! Then, a) transfer the snow excess into the ice (different from icethd_dh) 124 124 zdh = MAX( 0._wp, ( rhosn * ht_s(ji,jj,jl) + ( rhoic - rau0 ) * ht_i(ji,jj,jl) ) * r1_rau0 ) 125 125 ! Or, b) transfer all the snow into ice (if incoming ice is likely to melt as it comes into a warmer environment) -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/LBC/lib_mpp.F90
r8226 r8422 145 145 146 146 ! variables used in case of sea-ice 147 INTEGER, PUBLIC :: ncomm_ice !: communicator made by the processors with sea-ice (public so that it can be freed in limthd)147 INTEGER, PUBLIC :: ncomm_ice !: communicator made by the processors with sea-ice (public so that it can be freed in icethd) 148 148 INTEGER :: ngrp_iworld ! group ID for the world processors (for rheology) 149 149 INTEGER :: ngrp_ice ! group ID for the ice processors (for rheology) -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk.F90
r8321 r8422 41 41 #if defined key_lim3 42 42 USE ice , ONLY : u_ice, v_ice, jpl, a_i_b, at_i_b 43 USE limthd_dh ! for CALL lim_thd_snwblow43 USE icethd_dh ! for CALL ice_thd_snwblow 44 44 #endif 45 45 USE sbcblk_algo_ncar ! => turb_ncar : NCAR - CORE (Large & Yeager, 2009) … … 789 789 ! --- evaporation minus precipitation --- ! 790 790 zsnw(:,:) = 0._wp 791 CALL lim_thd_snwblow( (1.-at_i_b(:,:)), zsnw ) ! snow distribution over ice after wind blowing791 CALL ice_thd_snwblow( (1.-at_i_b(:,:)), zsnw ) ! snow distribution over ice after wind blowing 792 792 emp_oce(:,:) = ( 1._wp - at_i_b(:,:) ) * zevap(:,:) - ( tprecip(:,:) - sprecip(:,:) ) - sprecip(:,:) * (1._wp - zsnw ) 793 793 emp_ice(:,:) = SUM( a_i_b(:,:,:) * evap_ice(:,:,:), dim=3 ) - sprecip(:,:) * zsnw -
branches/2017/dev_r8183_ICEMODEL/NEMOGCM/NEMO/OPA_SRC/SBC/sbccpl.F90
r8321 r8422 40 40 #endif 41 41 #if defined key_lim3 42 USE limthd_dh ! for CALL lim_thd_snwblow42 USE icethd_dh ! for CALL ice_thd_snwblow 43 43 #endif 44 44 ! … … 1621 1621 #if defined key_lim3 1622 1622 ! zsnw = snow fraction over ice after wind blowing (=picefr if no blowing) 1623 zsnw(:,:) = 0._wp ; CALL lim_thd_snwblow( ziceld, zsnw )1623 zsnw(:,:) = 0._wp ; CALL ice_thd_snwblow( ziceld, zsnw ) 1624 1624 1625 1625 ! --- evaporation minus precipitation corrected (because of wind blowing on snow) --- !
Note: See TracChangeset
for help on using the changeset viewer.