pkg/ebm/EBM.h

2c03616886 Jean*0001 CBOP
                0002 C     !ROUTINE: EBM.h
                0003 C     !INTERFACE:
                0004 C     #include EBM.h
                0005
                0006 C     !DESCRIPTION:
                0007 C     Header file defining EBM pkg parameters and variables
                0008 CEOP
                0009
                0010 C     ZonalMeanSST :: zonal mean sea surface temperature
                0011 C     CountX  :: number of ocean points in each latitude band
                0012 C     TmlS    :: Average mid-latitude temperature in the Southern Hemisphere
                0013 C     TmlN    ::     same in Northern Hemisphere
                0014 C     TtS     :: Average tropical temperature in the Southern Hemisphere
                0015 C     TtN     ::     same in Northern Hemisphere
                0016 C     winPert :: weather  patterns added to the background wind.
                0017 C     Run     :: Runoff into the ocean [unit = m/s unless scale_runoff <> 1]
8831cd24b3 Jean*0018       COMMON /EBM_FLD_RL/
d8206d87ee Patr*0019      &               ZonalMeanSST, CountX,
                0020      &               TmlN, TmlS, TtN, TtS
4e8c9bded2 Patr*0021       _RL ZonalMeanSST(1-OLy:sNy+OLy,nSy)
                0022       _RL CountX(1-OLy:sNy+OLy,nSy)
d8206d87ee Patr*0023       _RL TmlS, TmlN, TtS, TtN
                0024
8831cd24b3 Jean*0025       COMMON /EBM_FLD_RS/
d8206d87ee Patr*0026      &               Run,
                0027      &               winPert
                0028       _RS Run      (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0029       _RS winPert  (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
                0030
2c03616886 Jean*0031 C     RunoffFile  :: Runoff input file name
d8206d87ee Patr*0032       COMMON /EBM_PARAM_C/
                0033      &               RunoffFile
                0034       CHARACTER*(MAX_LEN_FNAM) RunoffFile
                0035
2c03616886 Jean*0036 C     latBnd           :: Latitude boundaries used in EBM
                0037 C     tauThetaZonRelax :: time-scale [s] for relaxation towards Zon.Aver. SST
d8206d87ee Patr*0038       COMMON /EBM_PARAM_RL/
2c03616886 Jean*0039      &               latBnd,
d8206d87ee Patr*0040      &               tauThetaZonRelax,
                0041      &               scale_runoff
2c03616886 Jean*0042       _RL latBnd(3)
d8206d87ee Patr*0043       _RL tauThetaZonRelax
                0044       _RL scale_runoff
                0045
2c03616886 Jean*0046 C--   Physical Constants
d8206d87ee Patr*0047       _RL t_mlt, lv, cp, rho_air
8831cd24b3 Jean*0048       PARAMETER ( t_mlt = 273.15 _d 0 )
                0049       PARAMETER ( lv = 2.5 _d 6 )
                0050       PARAMETER ( cp = 1004. _d 0 )
                0051       PARAMETER ( rho_air = 1.27 _d 0 )
2c03616886 Jean*0052
                0053 C--   Constant parameters for EBM computation
                0054 C-    Shortwave and albedo parameters
d8206d87ee Patr*0055       _RL Q0, Q2, A0, A2, A4
8831cd24b3 Jean*0056       PARAMETER ( Q0 = 1365. _d 0 , Q2 = -0.482 _d 0 )
                0057       PARAMETER ( A0 = 0.322 _d 0 , A2 = 0.231 _d 0 , A4 = 0.086 _d 0  )
2c03616886 Jean*0058 C-    Longwave parameters
d8206d87ee Patr*0059       _RL LW0, LW1
8831cd24b3 Jean*0060       PARAMETER ( LW0 = 195. _d 0 , LW1 = 2.78 _d 0 )
2c03616886 Jean*0061 C-    Heat transport parameters
d8206d87ee Patr*0062       _RL H1, H3, H5
95514d96b2 Jean*0063       PARAMETER ( H1 = 3.866 _d 0 , H3 = -2.851 _d 0, H5 = -1.016 _d 0 )
2c03616886 Jean*0064 C-    Freshwater flux parameters
d8206d87ee Patr*0065       _RL F1, F2, F3, F4, F5
8831cd24b3 Jean*0066       PARAMETER ( F1 = 2.092 _d 0 , F2 = 5.796 _d 0 , F3 = 8.472 _d 0 )
                0067       PARAMETER ( F4 = 7.728 _d 0 , F5 = 2.362 _d 0  )
2c03616886 Jean*0068 C-    Parameters used to calculate the transport efficiency
d8206d87ee Patr*0069       _RL trans_eff, Hw, Nw, Tw, At, dz, htil, tau
8831cd24b3 Jean*0070       PARAMETER ( trans_eff = 2. _d 0 , Hw = 8. _d 3, Nw = 0.013 _d 0 )
                0071       PARAMETER ( Tw = 300. _d 0 , At = 3. _d 0 , dz = 450. _d 0 )
                0072       PARAMETER ( htil = 0.8 _d 0 , tau = 5. _d -3)
2c03616886 Jean*0073 C-    Climate change parameter (W/(m^2 y) )
d8206d87ee Patr*0074       _RL DLW
8831cd24b3 Jean*0075       PARAMETER (DLW = 0.06 _d 0 )
d8206d87ee Patr*0076
2c03616886 Jean*0077 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|