C $Header: /u/gcmpack/MITgcm/pkg/exf/EXF_CONSTANTS.h,v 1.5 2008/11/17 23:08:29 jmc Exp $
C $Name: $
c
c
c ==================================================================
c HEADER exf_constants
c ==================================================================
c
c o Header file for constants.
c These include - numbers (e.g. 1, 2, 1/2, ...)
c - physical constants (e.g. gravitational const.)
c - empirical parameters
c - control parameters (e.g. max. no of iteration)
c
c started: Patrick Heimbach heimbach@mit.edu 06-May-2000
c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
c
c ==================================================================
c HEADER exf_constants
c ==================================================================
c 1. numbers
c exf_half 0.5
c exf_one 1.0
c exf_two 2.0
_RL exf_half
_RL exf_one
_RL exf_two
PARAMETER(
& exf_half = 0.5 _d 0 ,
& exf_one = 1.0 _d 0 ,
& exf_two = 2.0 _d 0
& )
c real exf_undef
c PARAMETER( exf_undef = -9000. )
c 2. physical constants
c stefanBoltzmann :: Stefan-Boltzmann constant [J*K^-4*m^-2*s^-1]
c sigma = (2*pi^5*k^4)/(15*h^3*c^2)
c karman :: von Karman constant
_RL stefanBoltzmann
_RL karman
PARAMETER ( stefanBoltzmann = 5.670 _d -8 )
PARAMETER ( karman = 0.4 _d 0 )
c 3. empirical parameters
c To invert the relationship ustar = ustar(umagn) the following
c parameterization is used:
c
c ustar**2 = umagn**2 * CDN(umagn)
c
c / cquadrag_1 * umagn**2 + cquadrag_2; 0 < u < 11 m/s
c CDN(umagn) =
c \ clindrag_1 * umagn + clindrag_2 ; u > 11 m/s
c
c clindrag_[n] - n = 1, 2 coefficients used to evaluate
c LINEAR relationship of Large and Pond 1981
c cquadrag_[n] - n = 1, 2 coefficients used to evaluate
c quadratic relationship
c u11 - u = 11 m/s wind speed
c ustofu11 - ustar = 0.3818 m/s, corresponding to u = 11 m/s
_RL clindrag_1, clindrag_2
_RL cquadrag_1, cquadrag_2
_RL u11
_RL ustofu11
PARAMETER (
& ustofu11 = 0.381800 _d 0 ,
& u11 = 11. _d 0 ,
& clindrag_1 = 0.000065 _d 0 ,
& clindrag_2 = 0.000490 _d 0 ,
& cquadrag_1 = clindrag_1/u11/2 ,
& cquadrag_2 = clindrag_1*u11/2 + clindrag_2
& )
c 4. control parameters
c niter_bulk - Number of iterations to be performed for the
c evaluation of the bulk surface fluxes. The ncom
c model uses 2 hardwired interation steps (loop
c unrolled).
c
INTEGER niter_bulk
PARAMETER ( niter_bulk = 2 )
C 5. other constants or parameters
C COMMON /EXF_PARAM_R_2/
C cen2kel :: conversion of deg. Centigrade to Kelvin
C gravity_mks :: gravitational acceleration [m/s^2]
C atmrho :: mean atmospheric density [kg/m^3]
C atmcp :: mean atmospheric specific heat [J/kg/K]
C flamb :: latent heat of evaporation [J/kg]
C flami :: latent heat of melting of pure ice [J/kg]
C cvapor_[] :: Coeff to calculate Saturation Specific Humidity
C see e.g. Gill (1982) p.41 Eq. (3.1.15)
C humid_fac :: constant entering the evaluation of the virtual
C temperature
C gamma_blk :: adiabatic lapse rate
C saltsat :: reduction of saturation vapor pressure over salt water
C sstExtrapol :: extrapolation coeff from 1rst 2 levels up to surface
C snow_emissivity :: longwave snow emissivity [-] (with pkg thsice/seaice)
C-- to evaluate turbulent transfert coefficients:
C cdrag_[n] :: n = 1,2,3 coefficients used to evaluate
C drag coefficient
C cstanton_[n] :: n = 1,2 coefficients used to evaluate
C the Stanton number (stable/unstable cond.)
C cdalton :: coefficient used to evaluate the Dalton number
C zolmin :: minimum stability parameter
C psim_fac :: coef used in turbulent fluxes calculation [-]
C zref :: reference height
C hu :: height of mean wind
C ht :: height of mean temperature
C hq :: height of mean rel. humidity
C umin :: minimum absolute wind speed used to evaluate
C drag coefficient [m/s]
C exf_iceCd :: drag coefficient over sea-ice (fixed)
C exf_iceCe :: transfert coeff. over sea-ice, for Evaporation (fixed)
C exf_iceCh :: transfert coeff. over sea-ice, for Sens.Heating (fixed)
C-- radiation:
C exf_albedo :: Sea-water albedo
C ocean_emissivity :: longwave ocean-surface emissivity [-]
C ice_emissivity :: longwave seaice emissivity [-] (with pkg thsice/seaice)
C snow_emissivity :: longwave snow emissivity [-] (with pkg thsice/seaice)
_RL cen2kel
_RL gravity_mks
_RL atmrho
_RL atmcp
_RL flamb, flami
_RL cvapor_fac, cvapor_exp
_RL cvapor_fac_ice, cvapor_exp_ice
_RL humid_fac
_RL gamma_blk
_RL saltsat
_RL sstExtrapol
_RL cdrag_1, cdrag_2, cdrag_3
_RL cstanton_1, cstanton_2
_RL cdalton
_RL zolmin
_RL psim_fac
_RL zref
_RL hu
_RL ht
_RL hq
_RL umin
_RL exf_iceCd
_RL exf_iceCe
_RL exf_iceCh
_RL exf_albedo
_RL ocean_emissivity
_RL ice_emissivity
_RL snow_emissivity
COMMON /EXF_PARAM_R_2/
& cen2kel,
& gravity_mks,
& atmrho,
& atmcp,
& flamb,
& flami,
& cvapor_fac, cvapor_exp,
& cvapor_fac_ice, cvapor_exp_ice,
& humid_fac,
& gamma_blk,
& saltsat,
& sstExtrapol,
& cdrag_1, cdrag_2, cdrag_3,
& cstanton_1, cstanton_2,
& cdalton,
& zolmin,
& psim_fac,
& zref,
& hu,
& ht,
& hq,
& umin,
& exf_iceCd, exf_iceCe, exf_iceCh,
& exf_albedo,
& ocean_emissivity,
& ice_emissivity,
& snow_emissivity