MITgcm |
** Warning **
Issuing rollback() due to DESTROY without explicit disconnect() of DBD::mysql::db handle dbname=MITgcm at /usr/local/share/lxr/lib/LXR/Common.pm line 1224.|
|
|
|||
File indexing completed on 2025-11-07 06:08:54 UTC
view on githubraw file Latest commit b7411f1a on 2025-11-06 19:05:26 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_snapshot_mdsio :: write snap-shot output using MDSIO 0141 C thSIce_mon_stdio :: write monitor to std-outp 0142 C thSIce_snapshot_mnc :: write snap-shot output using MNC 0143 C thSIce_mon_mnc :: write monitor to netcdf file 0144 C thSIce_pickup_read_mnc :: pickup read w/ MNC 0145 C thSIce_pickup_write_mnc :: pickup write w/ MNC 0146 C thSIce_pickup_write_mdsio :: pickup write w/ MDSIO 7269783f6f Jean*0147 COMMON / THSICE_PAR_L / 8a23e1b5d8 Jean*0148 & stepFwd_oceMxL, thSIce_skipThermo, ce354ad541 Jean*0149 & thSIce_calc_albNIR, b7411f1a84 Jean*0150 & thSIce_snapshot_mdsio, thSIce_mon_stdio, 0151 & thSIce_snapshot_mnc, thSIce_mon_mnc, 7269783f6f Jean*0152 & thSIce_pickup_read_mnc, df4e8f7bcf Ed H*0153 & thSIce_pickup_write_mdsio, 0154 & thSIce_pickup_write_mnc 87ea84cac6 Jean*0155 0156 LOGICAL stepFwd_oceMxL 8a23e1b5d8 Jean*0157 LOGICAL thSIce_skipThermo ce354ad541 Jean*0158 LOGICAL thSIce_calc_albNIR b7411f1a84 Jean*0159 LOGICAL thSIce_snapshot_mdsio, thSIce_mon_stdio 0160 LOGICAL thSIce_snapshot_mnc, thSIce_mon_mnc df4e8f7bcf Ed H*0161 LOGICAL thSIce_pickup_read_mnc 0162 LOGICAL thSIce_pickup_write_mdsio 0163 LOGICAL thSIce_pickup_write_mnc 87ea84cac6 Jean*0164 fc7306ba7d Jean*0165 C-- COMMON / THSICE_PAR_I / ice model (integer) parameters 49d6795014 Jean*0166 C startIceModel :: =1 : start ice model at nIter0 ; =0 : use pickup files 0167 C :: -1 : start from a small pickup (without Mix.Layer) 0168 C nitMaxTsf :: maximum Nb of iter to find Surface Temp (Trsf) 0169 C thSIceAdvScheme :: thSIce Advection scheme selector e40ff4922c Jean*0170 C thSIceBalanceAtmFW :: select balancing Fresh-Water flux from Atm+Land 0171 C :: =0 : none ; =1 : uniform ; =2 : scaled by Precip 7269783f6f Jean*0172 COMMON / THSICE_PAR_I / e40ff4922c Jean*0173 & startIceModel, nitMaxTsf, thSIceAdvScheme, thSIceBalanceAtmFW fc7306ba7d Jean*0174 0175 INTEGER startIceModel 0176 INTEGER nitMaxTsf 49d6795014 Jean*0177 INTEGER thSIceAdvScheme e40ff4922c Jean*0178 INTEGER thSIceBalanceAtmFW fc7306ba7d Jean*0179 0180 C-- COMMON / THSICE_PAR_R / ice model (real) parameters 2e83cb4d40 Jeff*0181 C thSIce_deltaT :: ice model time-step, seaice thicken/extend [s] 0182 C thSIce_dtTemp :: ice model time-step, solve4temp [s] 49d6795014 Jean*0183 C ocean_deltaT :: ocean mixed-layer time-step [s] 7c5790770d Davi*0184 C tauRelax_MxL :: Relaxation time scale for MixLayer T [s] 0185 C tauRelax_MxL_salt :: Relaxation time scale for MixLayer S [s] 87ea84cac6 Jean*0186 C hMxL_default :: default value for ocean MixLayer thickness [m] ba0b047096 Mart*0187 C sMxL_default :: default value for salinity in MixLayer [g/kg] 87ea84cac6 Jean*0188 C vMxL_default :: default value for ocean current velocity in MxL [m/s] 49d6795014 Jean*0189 C thSIce_diffK :: thickness (horizontal) diffusivity [m^2/s] 87ea84cac6 Jean*0190 C stressReduction :: reduction factor for wind-stress under sea-ice [0-1] 0191 C thSIce_diagFreq :: Frequency^-1 for diagnostic output [s] 0192 C thSIce_monFreq :: Frequency^-1 for monitor output [s] fc7306ba7d Jean*0193 COMMON / THSICE_PAR_R / 7c5790770d Davi*0194 & thSIce_deltaT, thSIce_dtTemp, ocean_deltaT, 0195 & tauRelax_MxL, tauRelax_MxL_salt, 87ea84cac6 Jean*0196 & hMxL_default, sMxL_default, vMxL_default, 49d6795014 Jean*0197 & thSIce_diffK, stressReduction, b7411f1a84 Jean*0198 & thSIce_diagFreq, thSIce_monFreq fc7306ba7d Jean*0199 7c5790770d Davi*0200 _RL thSIce_deltaT, thSIce_dtTemp, ocean_deltaT 0201 _RL tauRelax_MxL, tauRelax_MxL_salt 7269783f6f Jean*0202 _RL hMxL_default, sMxL_default, vMxL_default 49d6795014 Jean*0203 _RL thSIce_diffK, stressReduction b7411f1a84 Jean*0204 _RL thSIce_diagFreq, thSIce_monFreq fc7306ba7d Jean*0205 987916f561 Jean*0206 C-- COMMON / THSICE_PAR_C / ice model (character) parameters 0207 C thSIceFract_InitFile :: File name for initial ice fraction 0208 C thSIceThick_InitFile :: File name for initial ice thickness 0209 C thSIceSnowH_InitFile :: File name for initial snow thickness 0210 C thSIceSnowA_InitFile :: File name for initial snow Age 0211 C thSIceEnthp_InitFile :: File name for initial ice enthalpy 0212 C thSIceTsurf_InitFile :: File name for initial surf. temp 0213 COMMON / THSICE_PAR_C / 0214 & thSIceFract_InitFile, 0215 & thSIceThick_InitFile, 0216 & thSIceSnowH_InitFile, 0217 & thSIceSnowA_InitFile, 0218 & thSIceEnthp_InitFile, 7269783f6f Jean*0219 & thSIceTsurf_InitFile 987916f561 Jean*0220 CHARACTER*(MAX_LEN_FNAM) thSIceFract_InitFile 0221 CHARACTER*(MAX_LEN_FNAM) thSIceThick_InitFile 0222 CHARACTER*(MAX_LEN_FNAM) thSIceSnowH_InitFile 0223 CHARACTER*(MAX_LEN_FNAM) thSIceSnowA_InitFile 0224 CHARACTER*(MAX_LEN_FNAM) thSIceEnthp_InitFile 0225 CHARACTER*(MAX_LEN_FNAM) thSIceTsurf_InitFile 0226 fc7306ba7d Jean*0227 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| 0228 0229 #endif /* ALLOW_THSICE */ df4e8f7bcf Ed H*0230 0231 CEH3 ;;; Local Variables: *** 0232 CEH3 ;;; mode:fortran *** 0233 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 |
|