MITgcm |
|
|
|
|||
File indexing completed on 2022-06-17 05:09:17 UTC
view on githubraw file Latest commit cc60455f on 2022-06-16 14:55:16 UTC082e18c36c Jean*0001 c 0002 c 0003 c ================================================================== 0004 c HEADER exf_constants 0005 c ================================================================== 0006 c 0007 c o Header file for constants. 0008 c These include - numbers (e.g. 1, 2, 1/2, ...) 0009 c - physical constants (e.g. gravitational const.) 0010 c - empirical parameters 0011 c - control parameters (e.g. max. no of iteration) 0012 c 0013 c started: Patrick Heimbach heimbach@mit.edu 06-May-2000 0014 c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002 0015 c 0016 c ================================================================== 0017 c HEADER exf_constants 0018 c ================================================================== 0019 0020 c 1. numbers 0021 0022 c exf_half 0.5 0023 c exf_one 1.0 0024 c exf_two 2.0 0025 0026 _RL exf_half 0027 _RL exf_one 0028 _RL exf_two 0029 ec1f56ea0e Jean*0030 PARAMETER( 0031 & exf_half = 0.5 _d 0 , 0032 & exf_one = 1.0 _d 0 , 0033 & exf_two = 2.0 _d 0 082e18c36c Jean*0034 & ) 0035 ec1f56ea0e Jean*0036 c real exf_undef 0037 c PARAMETER( exf_undef = -9000. ) 082e18c36c Jean*0038 0039 c 2. physical constants 0040 75863b3332 Jean*0041 c stefanBoltzmann :: Stefan-Boltzmann constant [J*K^-4*m^-2*s^-1] 0042 c sigma = (2*pi^5*k^4)/(15*h^3*c^2) 0043 c karman :: von Karman constant 0044 _RL stefanBoltzmann 0045 _RL karman ec1f56ea0e Jean*0046 PARAMETER ( stefanBoltzmann = 5.670 _d -8 ) 0047 PARAMETER ( karman = 0.4 _d 0 ) 082e18c36c Jean*0048 20cf30c902 Patr*0049 c 3. empirical parameters 082e18c36c Jean*0050 0051 c To invert the relationship ustar = ustar(umagn) the following 0052 c parameterization is used: 75863b3332 Jean*0053 c 082e18c36c Jean*0054 c ustar**2 = umagn**2 * CDN(umagn) 75863b3332 Jean*0055 c 082e18c36c Jean*0056 c / cquadrag_1 * umagn**2 + cquadrag_2; 0 < u < 11 m/s 0057 c CDN(umagn) = 0058 c \ clindrag_1 * umagn + clindrag_2 ; u > 11 m/s 75863b3332 Jean*0059 c 082e18c36c Jean*0060 c clindrag_[n] - n = 1, 2 coefficients used to evaluate 0061 c LINEAR relationship of Large and Pond 1981 0062 c cquadrag_[n] - n = 1, 2 coefficients used to evaluate 0063 c quadratic relationship 0064 c u11 - u = 11 m/s wind speed 0065 c ustofu11 - ustar = 0.3818 m/s, corresponding to u = 11 m/s 0066 0067 _RL clindrag_1, clindrag_2 0068 _RL cquadrag_1, cquadrag_2 0069 _RL u11 0070 _RL ustofu11 0071 ec1f56ea0e Jean*0072 PARAMETER ( 0073 & ustofu11 = 0.381800 _d 0 , 0074 & u11 = 11. _d 0 , 0075 & clindrag_1 = 0.000065 _d 0 , 0076 & clindrag_2 = 0.000490 _d 0 , 082e18c36c Jean*0077 & cquadrag_1 = clindrag_1/u11/2 , 0078 & cquadrag_2 = clindrag_1*u11/2 + clindrag_2 0079 & ) 0080 75863b3332 Jean*0081 c 4. control parameters 0082 0083 c niter_bulk - Number of iterations to be performed for the 0084 c evaluation of the bulk surface fluxes. The ncom 0085 c model uses 2 hardwired interation steps (loop 0086 c unrolled). 0087 c ec1f56ea0e Jean*0088 INTEGER niter_bulk 0089 PARAMETER ( niter_bulk = 2 ) 75863b3332 Jean*0090 0091 C 5. other constants or parameters 0092 ec1f56ea0e Jean*0093 C COMMON /EXF_PARAM_R_2/ 75863b3332 Jean*0094 C cen2kel :: conversion of deg. Centigrade to Kelvin 0095 C gravity_mks :: gravitational acceleration [m/s^2] 0096 C atmrho :: mean atmospheric density [kg/m^3] 0097 C atmcp :: mean atmospheric specific heat [J/kg/K] 0098 C flamb :: latent heat of evaporation [J/kg] 0099 C flami :: latent heat of melting of pure ice [J/kg] 0100 C cvapor_[] :: Coeff to calculate Saturation Specific Humidity 0101 C see e.g. Gill (1982) p.41 Eq. (3.1.15) 0102 C humid_fac :: constant entering the evaluation of the virtual 0103 C temperature 0104 C gamma_blk :: adiabatic lapse rate 0105 C saltsat :: reduction of saturation vapor pressure over salt water ec1f56ea0e Jean*0106 C sstExtrapol :: extrapolation coeff from 1rst 2 levels up to surface 0107 C snow_emissivity :: longwave snow emissivity [-] (with pkg thsice/seaice) 75863b3332 Jean*0108 C-- to evaluate turbulent transfert coefficients: 0109 C cdrag_[n] :: n = 1,2,3 coefficients used to evaluate cc60455fbb Mart*0110 C drag coefficient, 0111 C For Large and Yeager (2009): extra coefficient n = 8 and 0112 C cdragMax :: maximum drag coefficient ... 0113 C umax :: ... at maximum wind 75863b3332 Jean*0114 C cstanton_[n] :: n = 1,2 coefficients used to evaluate 0115 C the Stanton number (stable/unstable cond.) 0116 C cdalton :: coefficient used to evaluate the Dalton number 0117 C zolmin :: minimum stability parameter 0118 C psim_fac :: coef used in turbulent fluxes calculation [-] 0119 C zref :: reference height 0120 C hu :: height of mean wind 0121 C ht :: height of mean temperature 0122 C hq :: height of mean rel. humidity 0123 C umin :: minimum absolute wind speed used to evaluate 0124 C drag coefficient [m/s] 0125 C exf_iceCd :: drag coefficient over sea-ice (fixed) 0126 C exf_iceCe :: transfert coeff. over sea-ice, for Evaporation (fixed) 0127 C exf_iceCh :: transfert coeff. over sea-ice, for Sens.Heating (fixed) 0128 C-- radiation: 0129 C exf_albedo :: Sea-water albedo 0130 C ocean_emissivity :: longwave ocean-surface emissivity [-] 0131 C ice_emissivity :: longwave seaice emissivity [-] (with pkg thsice/seaice) 0132 C snow_emissivity :: longwave snow emissivity [-] (with pkg thsice/seaice) 0133 0134 _RL cen2kel 0135 _RL gravity_mks 0136 _RL atmrho 0137 _RL atmcp 0138 _RL flamb, flami 0139 _RL cvapor_fac, cvapor_exp 0140 _RL cvapor_fac_ice, cvapor_exp_ice 0141 _RL humid_fac 0142 _RL gamma_blk 0143 _RL saltsat ec1f56ea0e Jean*0144 _RL sstExtrapol cc60455fbb Mart*0145 _RL cdrag_1, cdrag_2, cdrag_3, cdrag_8, cdragMax, umax 75863b3332 Jean*0146 _RL cstanton_1, cstanton_2 0147 _RL cdalton 0148 _RL zolmin 0149 _RL psim_fac 0150 _RL zref 0151 _RL hu 0152 _RL ht 0153 _RL hq 0154 _RL umin 0155 _RL exf_iceCd 0156 _RL exf_iceCe 0157 _RL exf_iceCh 0158 _RL exf_albedo 0159 _RL ocean_emissivity 0160 _RL ice_emissivity 0161 _RL snow_emissivity 082e18c36c Jean*0162 ec1f56ea0e Jean*0163 COMMON /EXF_PARAM_R_2/ 75863b3332 Jean*0164 & cen2kel, 0165 & gravity_mks, 20cf30c902 Patr*0166 & atmrho, 0167 & atmcp, 0168 & flamb, 0169 & flami, 75863b3332 Jean*0170 & cvapor_fac, cvapor_exp, 0171 & cvapor_fac_ice, cvapor_exp_ice, 20cf30c902 Patr*0172 & humid_fac, 0173 & gamma_blk, 0174 & saltsat, ec1f56ea0e Jean*0175 & sstExtrapol, cc60455fbb Mart*0176 & cdrag_1, cdrag_2, cdrag_3, cdrag_8, cdragMax, umax, 75863b3332 Jean*0177 & cstanton_1, cstanton_2, 0178 & cdalton, 0179 & zolmin, 20cf30c902 Patr*0180 & psim_fac, 75863b3332 Jean*0181 & zref, 20cf30c902 Patr*0182 & hu, 0183 & ht, 75863b3332 Jean*0184 & hq, 0185 & umin, 0186 & exf_iceCd, exf_iceCe, exf_iceCh, 0187 & exf_albedo, 0188 & ocean_emissivity, 0189 & ice_emissivity, 0190 & snow_emissivity
| [ 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 |
|