C $Header: /u/gcmpack/MITgcm/pkg/ebm/EBM.h,v 1.6 2012/02/16 02:30:34 jmc Exp $ C $Name: $ CBOP C !ROUTINE: EBM.h C !INTERFACE: C #include EBM.h C !DESCRIPTION: C Header file defining EBM pkg parameters and variables CEOP C ZonalMeanSST :: zonal mean sea surface temperature C CountX :: number of ocean points in each latitude band C TmlS :: Average mid-latitude temperature in the Southern Hemisphere C TmlN :: same in Northern Hemisphere C TtS :: Average tropical temperature in the Southern Hemisphere C TtN :: same in Northern Hemisphere C winPert :: weather patterns added to the background wind. C Run :: Runoff into the ocean [unit = m/s unless scale_runoff <> 1] COMMON /EBM_FLD_RL/ & ZonalMeanSST, CountX, & TmlN, TmlS, TtN, TtS _RL ZonalMeanSST(1-OLy:sNy+OLy,nSy) _RL CountX(1-OLy:sNy+OLy,nSy) _RL TmlS, TmlN, TtS, TtN COMMON /EBM_FLD_RS/ & Run, & winPert _RS Run (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) _RS winPert (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) C RunoffFile :: Runoff input file name COMMON /EBM_PARAM_C/ & RunoffFile CHARACTER*(MAX_LEN_FNAM) RunoffFile C latBnd :: Latitude boundaries used in EBM C tauThetaZonRelax :: time-scale [s] for relaxation towards Zon.Aver. SST COMMON /EBM_PARAM_RL/ & latBnd, & tauThetaZonRelax, & scale_runoff _RL latBnd(3) _RL tauThetaZonRelax _RL scale_runoff C-- Physical Constants _RL t_mlt, lv, cp, rho_air PARAMETER ( t_mlt = 273.15 _d 0 ) PARAMETER ( lv = 2.5 _d 6 ) PARAMETER ( cp = 1004. _d 0 ) PARAMETER ( rho_air = 1.27 _d 0 ) C-- Constant parameters for EBM computation C- Shortwave and albedo parameters _RL Q0, Q2, A0, A2, A4 PARAMETER ( Q0 = 1365. _d 0 , Q2 = -0.482 _d 0 ) PARAMETER ( A0 = 0.322 _d 0 , A2 = 0.231 _d 0 , A4 = 0.086 _d 0 ) C- Longwave parameters _RL LW0, LW1 PARAMETER ( LW0 = 195. _d 0 , LW1 = 2.78 _d 0 ) C- Heat transport parameters _RL H1, H3, H5 PARAMETER ( H1 = 3.866 _d 0 , H3 = -2.851 _d 0, H5 = -1.016 _d 0 ) C- Freshwater flux parameters _RL F1, F2, F3, F4, F5 PARAMETER ( F1 = 2.092 _d 0 , F2 = 5.796 _d 0 , F3 = 8.472 _d 0 ) PARAMETER ( F4 = 7.728 _d 0 , F5 = 2.362 _d 0 ) C- Parameters used to calculate the transport efficiency _RL trans_eff, Hw, Nw, Tw, At, dz, htil, tau PARAMETER ( trans_eff = 2. _d 0 , Hw = 8. _d 3, Nw = 0.013 _d 0 ) PARAMETER ( Tw = 300. _d 0 , At = 3. _d 0 , dz = 450. _d 0 ) PARAMETER ( htil = 0.8 _d 0 , tau = 5. _d -3) C- Climate change parameter (W/(m^2 y) ) _RL DLW PARAMETER (DLW = 0.06 _d 0 ) C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|