|
||||
File indexing completed on 2021-04-08 05:12:33 UTC
view on githubraw file Latest commit ba0b0470 on 2021-04-08 01:06:32 UTCfc7306ba7d Jean*0001 #ifdef ALLOW_THSICE 0002 C *==========================================================* 0003 C | THSICE_PARAMS.h 0004 C | o Header file for Therm_SeaIce package parameters: 0005 C | - basic parameter ( I/O frequency, etc ...) 0006 C | - physical constants (used in therm_SeaIce pkg) 0007 C *==========================================================* 0008 0009 C---------------------------------------------------------------------------- 0010 C.. Common blocks for almost everything that the sea ice model passes around. 0011 C---------------------------------------------------------------------------- 0012 0013 C-- COMMON / THSICE_PHYSPAR_R / physical (real) parameter 0014 C.. densities c2c29a499f Jean*0015 C rhos :: density of snow [kg/m^3] 0016 C rhoi :: density of ice [kg/m^3] 0017 C rhosw :: density of seawater [kg/m^3] 0018 C rhofw :: density of fresh water [kg/m^3] e6b6bab319 Davi*0019 C floodFac :: flooding factor = (rhosw-rhoi)/rhos [dimensionless] fc7306ba7d Jean*0020 C.. specific heats c2c29a499f Jean*0021 C cpIce :: specific heat of fresh ice [J/kg/K] 0022 C cpWater :: specific heat of water [J/kg/K] fc7306ba7d Jean*0023 C .. thermal conductivity. QQ check units c2c29a499f Jean*0024 C kIce :: thermal conductivity of pure ice [W/m/K] 0025 C kSnow :: thermal conductivity of snow [W/m/K] fc7306ba7d Jean*0026 C .. heat transfer coefficient c2c29a499f Jean*0027 C bMeltCoef :: base-melting heat transfer coefficient 0028 C (between ice & water) [no unit] fc7306ba7d Jean*0029 C .. latent heat c2c29a499f Jean*0030 C Lfresh :: latent heat of melting of pure ice [J/kg] fcd60511e1 Jean*0031 C .. Enthalpy c2c29a499f Jean*0032 C qsnow :: snow enthalpy [J/kg] fc7306ba7d Jean*0033 C .. Albedo c854f591d1 Jean*0034 C albColdSnow :: albedo of cold (=dry) new snow (Tsfc < tempSnowAlb) fcd60511e1 Jean*0035 C albWarmSnow :: albedo of warm (=wet) new snow (Tsfc = 0) c2c29a499f Jean*0036 C tempSnowAlb :: temperature transition from ColdSnow to WarmSnow Alb. [oC] fcd60511e1 Jean*0037 C albOldSnow :: albedo of old snow (snowAge > 35.d) 0038 C albIceMax :: max albedo of bare ice (thick ice) 0039 C albIceMin :: minimum ice albedo (very thin ice) 0040 C hAlbIce :: ice thickness for albedo transition: thin/thick ice albedo 0041 C hAlbSnow :: snow thickness for albedo transition: snow/ice albedo 0042 C hNewSnowAge :: new snow thickness that refresh the snow-age (by 1/e) c2c29a499f Jean*0043 C snowAgTime :: snow aging time scale [s] fcd60511e1 Jean*0044 C .. Solar parameters c2c29a499f Jean*0045 C i0swFrac :: fraction of penetrating solar rad 0046 C ksolar :: bulk solar abs coeff of sea ice [m^-1] 0047 C dhSnowLin :: half slope of linear distribution of snow thickness within 0048 C the grid-cell (from hSnow-dhSnow to hSnow+dhSnow, if full 0049 C ice & snow cover) [m] ; (only used for SW radiation). fc7306ba7d Jean*0050 C .. Salinity c2c29a499f Jean*0051 C saltIce :: salinity of ice [g/kg] 0052 C S_winton :: Winton salinity of ice [g/kg] fc7306ba7d Jean*0053 C .. melting c2c29a499f Jean*0054 C Tf0kel :: Freezing temp of fresh water in Kelvin = 273.15 ba0b047096 Mart*0055 C mu_Tf :: linear dependence of melting temperature on salinity [oC/(g/kg)] c2c29a499f Jean*0056 C Tf(sea-water) = -mu_Tf * S 0057 C Tmlt1 :: Winton melting temperature: Tmlt1 = -mu_Tf * S_winton 0058 C Terrmax :: Temperature convergence criteria [oC] fc7306ba7d Jean*0059 C .. Min/Max c2c29a499f Jean*0060 C hIceMin :: Minimum ice thickness [m] 0061 C hiMax :: Maximum ice thickness [m] 0062 C hsMax :: Maximum snow thickness [m] fc7306ba7d Jean*0063 C .. for fractional ice c2c29a499f Jean*0064 C iceMaskMax :: maximum Ice fraction (=1 for no fractional ice) 0065 C iceMaskMin :: mimimum Ice fraction (=1 for no fractional ice) 0066 C fracEnFreez :: fraction of energy going to lateral freezing (vs height increase) 0067 C fracEnMelt :: fraction of energy going to lateral melting (vs height decrease) 0068 C (=0 for no fract. ice) 0069 C hThinIce :: ice height above which fracEnMelt/Freez are applied [m] 0070 C (=hIceMin for no fractional ice) 0071 C hThickIce :: ice height below which fracEnMelt/Freez are applied [m] 0072 C (=large for no fractional ice) 0073 C hNewIceMax :: new ice maximum thickness [m] fc7306ba7d Jean*0074 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| 0075 7269783f6f Jean*0076 COMMON / THSICE_PHYSPAR_R / e6b6bab319 Davi*0077 & rhos, rhoi, rhosw, rhofw, floodFac, c2c29a499f Jean*0078 & cpIce, cpWater, 0079 & kIce, kSnow, 0080 & bMeltCoef, Lfresh, qsnow, c854f591d1 Jean*0081 & albColdSnow, albWarmSnow, tempSnowAlb, 0082 & albOldSnow, hNewSnowAge, snowAgTime, fcd60511e1 Jean*0083 & albIceMax, albIceMin, hAlbIce, hAlbSnow, c2c29a499f Jean*0084 & i0swFrac, ksolar, dhSnowLin, 0085 & saltIce, S_winton, mu_Tf, 0086 & Tf0kel, Tmlt1, Terrmax, 0087 & hIceMin, hiMax, hsMax, 0088 & iceMaskMax, iceMaskMin, 0089 & fracEnMelt, fracEnFreez, 0090 & hThinIce, hThickIce, hNewIceMax fc7306ba7d Jean*0091 0092 _RL rhos 0093 _RL rhoi 0094 _RL rhosw 0095 _RL rhofw e6b6bab319 Davi*0096 _RL floodFac c2c29a499f Jean*0097 _RL cpIce 0098 _RL cpWater 0099 _RL kIce 0100 _RL kSnow 0101 _RL bMeltCoef fc7306ba7d Jean*0102 _RL Lfresh 0103 _RL qsnow fcd60511e1 Jean*0104 _RL albColdSnow 0105 _RL albWarmSnow c854f591d1 Jean*0106 _RL tempSnowAlb fcd60511e1 Jean*0107 _RL albOldSnow c854f591d1 Jean*0108 _RL hNewSnowAge 0109 _RL snowAgTime fcd60511e1 Jean*0110 _RL albIceMax 0111 _RL albIceMin 0112 _RL hAlbIce 0113 _RL hAlbSnow c2c29a499f Jean*0114 _RL i0swFrac fc7306ba7d Jean*0115 _RL ksolar c2c29a499f Jean*0116 _RL dhSnowLin 0117 _RL saltIce fc7306ba7d Jean*0118 _RL S_winton 0119 _RL mu_Tf 0120 _RL Tf0kel 0121 _RL Tmlt1 0122 _RL Terrmax c2c29a499f Jean*0123 _RL hIceMin fc7306ba7d Jean*0124 _RL hiMax 0125 _RL hsMax c2c29a499f Jean*0126 _RL iceMaskMax 0127 _RL iceMaskMin 0128 _RL fracEnMelt 0129 _RL fracEnFreez 0130 _RL hThinIce 0131 _RL hThickIce 0132 _RL hNewIceMax fc7306ba7d Jean*0133 0134 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| 0135 87ea84cac6 Jean*0136 C-- COMMON / THSICE_PAR_L / ice model (logical) parameters df4e8f7bcf Ed H*0137 C stepFwd_oceMxL :: step forward mixed-layer T & S (slab-ocean) 8a23e1b5d8 Jean*0138 C thSIce_skipThermo :: by-pass seaice thermodynamics ce354ad541 Jean*0139 C thSIce_calc_albNIR :: calculate Near Infra-Red Albedo df4e8f7bcf Ed H*0140 C thSIce_tave_mdsio :: write TimeAverage output using MDSIO 0141 C thSIce_snapshot_mdsio :: write snap-shot output using MDSIO 0142 C thSIce_mon_stdio :: write monitor to std-outp 0143 C thSIce_tave_mnc :: write TimeAverage output using MNC 0144 C thSIce_snapshot_mnc :: write snap-shot output using MNC 0145 C thSIce_mon_mnc :: write monitor to netcdf file 0146 C thSIce_pickup_read_mnc :: pickup read w/ MNC 0147 C thSIce_pickup_write_mnc :: pickup write w/ MNC 0148 C thSIce_pickup_write_mdsio :: pickup write w/ MDSIO 7269783f6f Jean*0149 COMMON / THSICE_PAR_L / 8a23e1b5d8 Jean*0150 & stepFwd_oceMxL, thSIce_skipThermo, ce354ad541 Jean*0151 & thSIce_calc_albNIR, df4e8f7bcf Ed H*0152 & thSIce_tave_mdsio, thSIce_snapshot_mdsio, thSIce_mon_stdio, 0153 & thSIce_tave_mnc, thSIce_snapshot_mnc, thSIce_mon_mnc, 7269783f6f Jean*0154 & thSIce_pickup_read_mnc, df4e8f7bcf Ed H*0155 & thSIce_pickup_write_mdsio, 0156 & thSIce_pickup_write_mnc 87ea84cac6 Jean*0157 0158 LOGICAL stepFwd_oceMxL 8a23e1b5d8 Jean*0159 LOGICAL thSIce_skipThermo ce354ad541 Jean*0160 LOGICAL thSIce_calc_albNIR df4e8f7bcf Ed H*0161 LOGICAL thSIce_tave_mdsio, thSIce_snapshot_mdsio, thSIce_mon_stdio 0162 LOGICAL thSIce_tave_mnc, thSIce_snapshot_mnc, thSIce_mon_mnc 0163 LOGICAL thSIce_pickup_read_mnc 0164 LOGICAL thSIce_pickup_write_mdsio 0165 LOGICAL thSIce_pickup_write_mnc 87ea84cac6 Jean*0166 fc7306ba7d Jean*0167 C-- COMMON / THSICE_PAR_I / ice model (integer) parameters 49d6795014 Jean*0168 C startIceModel :: =1 : start ice model at nIter0 ; =0 : use pickup files 0169 C :: -1 : start from a small pickup (without Mix.Layer) 0170 C nitMaxTsf :: maximum Nb of iter to find Surface Temp (Trsf) 0171 C thSIceAdvScheme :: thSIce Advection scheme selector e40ff4922c Jean*0172 C thSIceBalanceAtmFW :: select balancing Fresh-Water flux from Atm+Land 0173 C :: =0 : none ; =1 : uniform ; =2 : scaled by Precip 7269783f6f Jean*0174 COMMON / THSICE_PAR_I / e40ff4922c Jean*0175 & startIceModel, nitMaxTsf, thSIceAdvScheme, thSIceBalanceAtmFW fc7306ba7d Jean*0176 0177 INTEGER startIceModel 0178 INTEGER nitMaxTsf 49d6795014 Jean*0179 INTEGER thSIceAdvScheme e40ff4922c Jean*0180 INTEGER thSIceBalanceAtmFW fc7306ba7d Jean*0181 0182 C-- COMMON / THSICE_PAR_R / ice model (real) parameters 2e83cb4d40 Jeff*0183 C thSIce_deltaT :: ice model time-step, seaice thicken/extend [s] 0184 C thSIce_dtTemp :: ice model time-step, solve4temp [s] 49d6795014 Jean*0185 C ocean_deltaT :: ocean mixed-layer time-step [s] 7c5790770d Davi*0186 C tauRelax_MxL :: Relaxation time scale for MixLayer T [s] 0187 C tauRelax_MxL_salt :: Relaxation time scale for MixLayer S [s] 87ea84cac6 Jean*0188 C hMxL_default :: default value for ocean MixLayer thickness [m] ba0b047096 Mart*0189 C sMxL_default :: default value for salinity in MixLayer [g/kg] 87ea84cac6 Jean*0190 C vMxL_default :: default value for ocean current velocity in MxL [m/s] 49d6795014 Jean*0191 C thSIce_diffK :: thickness (horizontal) diffusivity [m^2/s] 87ea84cac6 Jean*0192 C stressReduction :: reduction factor for wind-stress under sea-ice [0-1] 0193 C thSIce_taveFreq :: Frequency^-1 for time-Aver. output [s] 0194 C thSIce_diagFreq :: Frequency^-1 for diagnostic output [s] 0195 C thSIce_monFreq :: Frequency^-1 for monitor output [s] fc7306ba7d Jean*0196 COMMON / THSICE_PAR_R / 7c5790770d Davi*0197 & thSIce_deltaT, thSIce_dtTemp, ocean_deltaT, 0198 & tauRelax_MxL, tauRelax_MxL_salt, 87ea84cac6 Jean*0199 & hMxL_default, sMxL_default, vMxL_default, 49d6795014 Jean*0200 & thSIce_diffK, stressReduction, fc7306ba7d Jean*0201 & thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq 0202 7c5790770d Davi*0203 _RL thSIce_deltaT, thSIce_dtTemp, ocean_deltaT 0204 _RL tauRelax_MxL, tauRelax_MxL_salt 7269783f6f Jean*0205 _RL hMxL_default, sMxL_default, vMxL_default 49d6795014 Jean*0206 _RL thSIce_diffK, stressReduction fc7306ba7d Jean*0207 _RL thSIce_taveFreq, thSIce_diagFreq, thSIce_monFreq 0208 987916f561 Jean*0209 C-- COMMON / THSICE_PAR_C / ice model (character) parameters 0210 C thSIceFract_InitFile :: File name for initial ice fraction 0211 C thSIceThick_InitFile :: File name for initial ice thickness 0212 C thSIceSnowH_InitFile :: File name for initial snow thickness 0213 C thSIceSnowA_InitFile :: File name for initial snow Age 0214 C thSIceEnthp_InitFile :: File name for initial ice enthalpy 0215 C thSIceTsurf_InitFile :: File name for initial surf. temp 0216 COMMON / THSICE_PAR_C / 0217 & thSIceFract_InitFile, 0218 & thSIceThick_InitFile, 0219 & thSIceSnowH_InitFile, 0220 & thSIceSnowA_InitFile, 0221 & thSIceEnthp_InitFile, 7269783f6f Jean*0222 & thSIceTsurf_InitFile 987916f561 Jean*0223 CHARACTER*(MAX_LEN_FNAM) thSIceFract_InitFile 0224 CHARACTER*(MAX_LEN_FNAM) thSIceThick_InitFile 0225 CHARACTER*(MAX_LEN_FNAM) thSIceSnowH_InitFile 0226 CHARACTER*(MAX_LEN_FNAM) thSIceSnowA_InitFile 0227 CHARACTER*(MAX_LEN_FNAM) thSIceEnthp_InitFile 0228 CHARACTER*(MAX_LEN_FNAM) thSIceTsurf_InitFile 0229 fc7306ba7d Jean*0230 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| 0231 0232 #endif /* ALLOW_THSICE */ df4e8f7bcf Ed H*0233 0234 CEH3 ;;; Local Variables: *** 0235 CEH3 ;;; mode:fortran *** 0236 CEH3 ;;; End: ***
[ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
This page was automatically generated from https://github.com/MITgcm/MITgcm by the 2.2.1-MITgcm-0.1 LXR engine. The LXR team |