MODULE ice !!====================================================================== !! *** MODULE ice *** !! LIM-3 Sea Ice physics: diagnostics variables of ice defined in memory !!===================================================================== !! History : 3.0 ! 2008-03 (M. Vancoppenolle) original code LIM-3 !! 4.0 ! 2011-02 (G. Madec) dynamical allocation !!---------------------------------------------------------------------- #if defined key_lim3 !!---------------------------------------------------------------------- !! 'key_lim3' LIM-3 sea-ice model !!---------------------------------------------------------------------- USE in_out_manager ! I/O manager USE lib_mpp ! MPP library IMPLICIT NONE PRIVATE PUBLIC ice_alloc ! Called in sbc_lim_init !!====================================================================== !! LIM3 by the use of sweat, agile fingers and sometimes brain juice, !! was developed in Louvain-la-Neuve by : !! * Martin Vancoppenolle (UCL-ASTR, Belgium) !! * Sylvain Bouillon (UCL-ASTR, Belgium) !! * Miguel Angel Morales Maqueda (NOC-L, UK) !! !! Based on extremely valuable earlier work by !! * Thierry Fichefet !! * Hugues Goosse !! !! The following persons also contributed to the code in various ways !! * Gurvan Madec, Claude Talandier, Christian Ethe (LOCEAN, France) !! * Xavier Fettweis (UCL-ASTR), Ralph Timmermann (AWI, Germany) !! * Bill Lipscomb (LANL), Cecilia Bitz (UWa) !! and Elisabeth Hunke (LANL), USA. !! !! For more info, the interested user is kindly invited to consult the following references !! For model description and validation : !! * Vancoppenolle et al., Ocean Modelling, 2008a. !! * Vancoppenolle et al., Ocean Modelling, 2008b. !! For a specific description of EVP : !! * Bouillon et al., Ocean Modelling 2009. !! !! Or the reference manual, that should be available by 2011 !!====================================================================== !! | !! I C E S T A T E V A R I A B L E S | !! | !! Introduction : | !! -------------- | !! Every ice-covered grid cell is characterized by a series of state | !! variables. To account for unresolved spatial variability in ice | !! thickness, the ice cover in divided in ice thickness categories. | !! | !! Sea ice state variables depend on the ice thickness category | !! | !! Those variables are divided into two groups | !! * Extensive (or global) variables. | !! These are the variables that are transported by all means | !! * Intensive (or equivalent) variables. | !! These are the variables that are either physically more | !! meaningful and/or used in ice thermodynamics | !! | !! Routines in limvar.F90 perform conversions | !! - lim_var_glo2eqv : from global to equivalent variables | !! - lim_var_eqv2glo : from equivalent to global variables | !! | !! For various purposes, the sea ice state variables have sometimes | !! to be aggregated over all ice thickness categories. This operation | !! is done in : | !! - lim_var_agg | !! | !! in icestp.F90, the routines that compute the changes in the ice | !! state variables are called | !! - lim_dyn : ice dynamics | !! - lim_trp : ice transport | !! - lim_itd_me : mechanical redistribution (ridging and rafting) | !! - lim_thd : ice halo-thermodynamics | !! - lim_itd_th : thermodynamic changes in ice thickness distribution | !! and creation of new ice | !! | !! See the associated routines for more information | !! | !! List of ice state variables : | !! ----------------------------- | !! | !!-------------|-------------|---------------------------------|-------| !! name in | name in | meaning | units | !! 2D routines | 1D routines | | | !!-------------|-------------|---------------------------------|-------| !! | !! ******************************************************************* | !! *** Dynamical variables (prognostic) *** | !! ******************************************************************* | !! | !! u_ice | - | Comp. U of the ice velocity | m/s | !! v_ice | - | Comp. V of the ice velocity | m/s | !! | !! ******************************************************************* | !! *** Category dependent state variables (prognostic) *** | !! ******************************************************************* | !! | !! ** Global variables | !!-------------|-------------|---------------------------------|-------| !! a_i | a_i_1d | Ice concentration | | !! v_i | - | Ice volume per unit area | m | !! v_s | - | Snow volume per unit area | m | !! smv_i | - | Sea ice salt content | ppt.m | !! oa_i ! - ! Sea ice areal age content | s | !! e_i ! - ! Ice enthalpy | J/m2 | !! - ! q_i_1d ! Ice enthalpy per unit vol. | J/m3 | !! e_s ! - ! Snow enthalpy | J/m2 | !! - ! q_s_1d ! Snow enthalpy per unit vol. | J/m3 | !! | !!-------------|-------------|---------------------------------|-------| !! | !! ** Equivalent variables | !!-------------|-------------|---------------------------------|-------| !! | !! ht_i | ht_i_1d | Ice thickness | m | !! ht_s ! ht_s_1d | Snow depth | m | !! sm_i ! sm_i_1d | Sea ice bulk salinity ! ppt | !! s_i ! s_i_1d | Sea ice salinity profile ! ppt | !! o_i ! - | Sea ice Age ! s | !! t_i ! t_i_1d | Sea ice temperature ! K | !! t_s ! t_s_1d | Snow temperature ! K | !! t_su ! t_su_1d | Sea ice surface temperature ! K | !! | !! notes: the ice model only sees a bulk (i.e., vertically averaged) | !! salinity, except in thermodynamic computations, for which | !! the salinity profile is computed as a function of bulk | !! salinity | !! | !! the sea ice surface temperature is not associated to any | !! heat content. Therefore, it is not a state variable and | !! does not have to be advected. Nevertheless, it has to be | !! computed to determine whether the ice is melting or not | !! | !! ******************************************************************* | !! *** Category-summed state variables (diagnostic) *** | !! ******************************************************************* | !! at_i | at_i_1d | Total ice concentration | | !! vt_i | - | Total ice vol. per unit area | m | !! vt_s | - | Total snow vol. per unit ar. | m | !! smt_i | - | Mean sea ice salinity | ppt | !! tm_i | - | Mean sea ice temperature | K | !! et_i ! - ! Total ice enthalpy | J/m2 | !! et_s ! - ! Total snow enthalpy | J/m2 | !! bv_i ! - ! relative brine volume | ??? | !!===================================================================== LOGICAL, PUBLIC :: con_i = .false. ! switch for conservation test !!-------------------------------------------------------------------------- !! * Share Module variables !!-------------------------------------------------------------------------- ! !!** ice-generic parameters namelist (namicerun) ** INTEGER , PUBLIC :: jpl !: number of ice categories INTEGER , PUBLIC :: nlay_i !: number of ice layers INTEGER , PUBLIC :: nlay_s !: number of snow layers REAL(wp) , PUBLIC :: rn_amax_n !: maximum ice concentration Northern hemisphere REAL(wp) , PUBLIC :: rn_amax_s !: maximum ice concentration Southern hemisphere CHARACTER(len=32) , PUBLIC :: cn_icerst_in !: suffix of ice restart name (input) CHARACTER(len=32) , PUBLIC :: cn_icerst_out !: suffix of ice restart name (output) CHARACTER(len=256), PUBLIC :: cn_icerst_indir !: ice restart input directory CHARACTER(len=256), PUBLIC :: cn_icerst_outdir!: ice restart output directory LOGICAL , PUBLIC :: ln_limthd !: flag for ice thermo (T) or not (F) LOGICAL , PUBLIC :: ln_limdyn !: flag for ice dynamics (T) or not (F) INTEGER , PUBLIC :: nn_limdyn !: flag for ice dynamics REAL(wp) , PUBLIC :: rn_uice !: prescribed u-vel (case nn_limdyn=0) REAL(wp) , PUBLIC :: rn_vice !: prescribed v-vel (case nn_limdyn=0) ! !!** ice-diagnostics namelist (namicediag) ** LOGICAL , PUBLIC :: ln_limdiachk !: flag for ice diag (T) or not (F) LOGICAL , PUBLIC :: ln_limdiahsb !: flag for ice diag (T) or not (F) LOGICAL , PUBLIC :: ln_limctl !: flag for sea-ice points output (T) or not (F) INTEGER , PUBLIC :: iiceprt !: debug i-point INTEGER , PUBLIC :: jiceprt !: debug j-point ! !!** ice-init namelist (namiceini) ** ! -- limistate -- ! LOGICAL , PUBLIC :: ln_limini ! initialization or not LOGICAL , PUBLIC :: ln_limini_file ! Ice initialization state from 2D netcdf file REAL(wp), PUBLIC :: rn_thres_sst ! threshold water temperature for initial sea ice REAL(wp), PUBLIC :: rn_hts_ini_n ! initial snow thickness in the north REAL(wp), PUBLIC :: rn_hts_ini_s ! initial snow thickness in the south REAL(wp), PUBLIC :: rn_hti_ini_n ! initial ice thickness in the north REAL(wp), PUBLIC :: rn_hti_ini_s ! initial ice thickness in the south REAL(wp), PUBLIC :: rn_ati_ini_n ! initial leads area in the north REAL(wp), PUBLIC :: rn_ati_ini_s ! initial leads area in the south REAL(wp), PUBLIC :: rn_smi_ini_n ! initial salinity REAL(wp), PUBLIC :: rn_smi_ini_s ! initial salinity REAL(wp), PUBLIC :: rn_tmi_ini_n ! initial temperature REAL(wp), PUBLIC :: rn_tmi_ini_s ! initial temperature ! !!** ice-thickness distribution namelist (namiceitd) ** INTEGER , PUBLIC :: nn_catbnd !: categories distribution following: tanh function (1), or h^(-alpha) function (2) REAL(wp), PUBLIC :: rn_himean !: mean thickness of the domain (used to compute the distribution, nn_itdshp = 2 only) ! !!** ice-dynamics namelist (namicedyn) ** ! -- limtrp & limadv -- ! INTEGER , PUBLIC :: nn_limadv !: choose the advection scheme (-1=Prather ; 0=Ultimate-Macho) INTEGER , PUBLIC :: nn_limadv_ord !: choose the order of the advection scheme (if Ultimate-Macho) ! -- limitd_me -- ! INTEGER , PUBLIC :: nn_icestr !: ice strength parameterization (0=Hibler79 1=Rothrock75) REAL(wp), PUBLIC :: rn_pe_rdg !: ridging work divided by pot. energy change in ridging, nn_icestr = 1 REAL(wp), PUBLIC :: rn_pstar !: determines ice strength, Hibler JPO79 REAL(wp), PUBLIC :: rn_crhg !: determines changes in ice strength LOGICAL , PUBLIC :: ln_icestr_bvf !: use brine volume to diminish ice strength ! -- limdyn & limrhg -- ! REAL(wp), PUBLIC :: rn_cio !: drag coefficient for oceanic stress REAL(wp), PUBLIC :: rn_creepl !: creep limit : has to be under 1.0e-9 REAL(wp), PUBLIC :: rn_ecc !: eccentricity of the elliptical yield curve INTEGER , PUBLIC :: nn_nevp !: number of iterations for subcycling REAL(wp), PUBLIC :: rn_relast !: ratio => telast/rdt_ice (1/3 or 1/9 depending on nb of subcycling nevp) LOGICAL , PUBLIC :: ln_landfast !: landfast ice parameterization (T or F) REAL(wp), PUBLIC :: rn_gamma !: fraction of ocean depth that ice must reach to initiate landfast ice REAL(wp), PUBLIC :: rn_icebfr !: maximum bottom stress per unit area of contact (landfast ice) REAL(wp), PUBLIC :: rn_lfrelax !: relaxation time scale (s-1) to reach static friction (landfast ice) ! !!** ice-diffusion namelist (namicehdf) ** INTEGER , PUBLIC :: nn_ahi0 !: sea-ice hor. eddy diffusivity coeff. (3 ways of calculation) REAL(wp), PUBLIC :: rn_ahi0_ref !: sea-ice hor. eddy diffusivity coeff. (m2/s) ! !!** ice-thermodynamics namelist (namicethd) ** ! -- limthd_dif -- ! REAL(wp), PUBLIC :: rn_kappa_i !: coef. for the extinction of radiation Grenfell et al. (2006) [1/m] REAL(wp), PUBLIC :: nn_conv_dif !: maximal number of iterations for heat diffusion REAL(wp), PUBLIC :: rn_terr_dif !: maximal tolerated error (C) for heat diffusion INTEGER , PUBLIC :: nn_ice_thcon !: thermal conductivity: =0 Untersteiner (1964) ; =1 Pringle et al (2007) LOGICAL , PUBLIC :: ln_it_qnsice !: iterate surface flux with changing surface temperature or not (F) INTEGER , PUBLIC :: nn_monocat !: virtual ITD mono-category parameterizations (1) or not (0) REAL(wp), PUBLIC :: rn_cdsn !: thermal conductivity of the snow [W/m/K] ! -- limthd_dh -- ! LOGICAL , PUBLIC :: ln_limdH !: activate ice thickness change from growing/melting (T) or not (F) REAL(wp), PUBLIC :: rn_betas !: coef. for partitioning of snowfall between leads and sea ice ! -- limthd_da -- ! LOGICAL , PUBLIC :: ln_limdA !: activate lateral melting param. (T) or not (F) REAL(wp), PUBLIC :: rn_beta !: coef. beta for lateral melting param. REAL(wp), PUBLIC :: rn_dmin !: minimum floe diameter for lateral melting param. ! -- limthd_lac -- ! LOGICAL , PUBLIC :: ln_limdO !: activate ice growth in open-water (T) or not (F) REAL(wp), PUBLIC :: rn_hnewice !: thickness for new ice formation (m) LOGICAL , PUBLIC :: ln_frazil !: use of frazil ice collection as function of wind (T) or not (F) REAL(wp), PUBLIC :: rn_maxfrazb !: maximum portion of frazil ice collecting at the ice bottom REAL(wp), PUBLIC :: rn_vfrazb !: threshold drift speed for collection of bottom frazil ice REAL(wp), PUBLIC :: rn_Cfrazb !: squeezing coefficient for collection of bottom frazil ice ! -- limitd_th -- ! REAL(wp), PUBLIC :: rn_himin !: minimum ice thickness ! !!** ice-salinity namelist (namicesal) ** LOGICAL , PUBLIC :: ln_limdS !: activate gravity drainage and flushing (T) or not (F) INTEGER , PUBLIC :: nn_icesal !: salinity configuration used in the model ! ! 1 - constant salinity in both space and time ! ! 2 - prognostic salinity (s(z,t)) ! ! 3 - salinity profile, constant in time REAL(wp), PUBLIC :: rn_icesal !: bulk salinity (ppt) in case of constant salinity REAL(wp), PUBLIC :: rn_sal_gd !: restoring salinity for gravity drainage [PSU] REAL(wp), PUBLIC :: rn_time_gd !: restoring time constant for gravity drainage (= 20 days) [s] REAL(wp), PUBLIC :: rn_sal_fl !: restoring salinity for flushing [PSU] REAL(wp), PUBLIC :: rn_time_fl !: restoring time constant for gravity drainage (= 10 days) [s] REAL(wp), PUBLIC :: rn_simax !: maximum ice salinity [PSU] REAL(wp), PUBLIC :: rn_simin !: minimum ice salinity [PSU] ! !!** ice-mechanical redistribution namelist (namiceitdme) REAL(wp), PUBLIC :: rn_cs !: fraction of shearing energy contributing to ridging INTEGER , PUBLIC :: nn_partfun !: participation function: =0 Thorndike et al. (1975), =1 Lipscomb et al. (2007) REAL(wp), PUBLIC :: rn_gstar !: fractional area of young ice contributing to ridging REAL(wp), PUBLIC :: rn_astar !: equivalent of G* for an exponential participation function LOGICAL , PUBLIC :: ln_ridging !: ridging of ice or not REAL(wp), PUBLIC :: rn_hstar !: thickness that determines the maximal thickness of ridged ice REAL(wp), PUBLIC :: rn_por_rdg !: initial porosity of ridges (0.3 regular value) REAL(wp), PUBLIC :: rn_fsnowrdg !: fractional snow loss to the ocean during ridging REAL(wp), PUBLIC :: rn_fpondrdg !: fractional melt pond loss to the ocean during ridging LOGICAL , PUBLIC :: ln_rafting !: rafting of ice or not REAL(wp), PUBLIC :: rn_hraft !: threshold thickness (m) for rafting / ridging REAL(wp), PUBLIC :: rn_craft !: coefficient for smoothness of the hyperbolic tangent in rafting REAL(wp), PUBLIC :: rn_fsnowrft !: fractional snow loss to the ocean during ridging REAL(wp), PUBLIC :: rn_fpondrft !: fractional snow loss to the ocean during rafting ! MV MP 2016 ! !!** melt pond namelist (namicemp) LOGICAL , PUBLIC :: ln_pnd !: activate ponds or not LOGICAL , PUBLIC :: ln_pnd_rad !: ponds radiatively active or not LOGICAL , PUBLIC :: ln_pnd_fw !: ponds active wrt meltwater or not INTEGER , PUBLIC :: nn_pnd_scheme !: type of melt pond scheme: =0 prescribed, =1 empirical, =2 topographic REAL(wp), PUBLIC :: rn_apnd !: prescribed pond fraction (0