1 |
C $Header: /u/gcmpack/MITgcm/pkg/ebm/EBM.h,v 1.5 2011/08/29 19:38:44 jmc Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
CBOP |
5 |
C !ROUTINE: EBM.h |
6 |
C !INTERFACE: |
7 |
C #include EBM.h |
8 |
|
9 |
C !DESCRIPTION: |
10 |
C Header file defining EBM pkg parameters and variables |
11 |
CEOP |
12 |
|
13 |
C ZonalMeanSST :: zonal mean sea surface temperature |
14 |
C CountX :: number of ocean points in each latitude band |
15 |
C TmlS :: Average mid-latitude temperature in the Southern Hemisphere |
16 |
C TmlN :: same in Northern Hemisphere |
17 |
C TtS :: Average tropical temperature in the Southern Hemisphere |
18 |
C TtN :: same in Northern Hemisphere |
19 |
C winPert :: weather patterns added to the background wind. |
20 |
C Run :: Runoff into the ocean [unit = m/s unless scale_runoff <> 1] |
21 |
COMMON /EBM_FLD_RL/ |
22 |
& ZonalMeanSST, CountX, |
23 |
& TmlN, TmlS, TtN, TtS |
24 |
_RL ZonalMeanSST(1-OLy:sNy+OLy,nSy) |
25 |
_RL CountX(1-OLy:sNy+OLy,nSy) |
26 |
_RL TmlS, TmlN, TtS, TtN |
27 |
|
28 |
COMMON /EBM_FLD_RS/ |
29 |
& Run, |
30 |
& winPert |
31 |
_RS Run (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
32 |
_RS winPert (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
33 |
|
34 |
C RunoffFile :: Runoff input file name |
35 |
COMMON /EBM_PARAM_C/ |
36 |
& RunoffFile |
37 |
CHARACTER*(MAX_LEN_FNAM) RunoffFile |
38 |
|
39 |
C latBnd :: Latitude boundaries used in EBM |
40 |
C tauThetaZonRelax :: time-scale [s] for relaxation towards Zon.Aver. SST |
41 |
COMMON /EBM_PARAM_RL/ |
42 |
& latBnd, |
43 |
& tauThetaZonRelax, |
44 |
& scale_runoff |
45 |
_RL latBnd(3) |
46 |
_RL tauThetaZonRelax |
47 |
_RL scale_runoff |
48 |
|
49 |
C-- Physical Constants |
50 |
_RL t_mlt, lv, cp, rho_air |
51 |
PARAMETER ( t_mlt = 273.15 _d 0 ) |
52 |
PARAMETER ( lv = 2.5 _d 6 ) |
53 |
PARAMETER ( cp = 1004. _d 0 ) |
54 |
PARAMETER ( rho_air = 1.27 _d 0 ) |
55 |
|
56 |
C-- Constant parameters for EBM computation |
57 |
C- Shortwave and albedo parameters |
58 |
_RL Q0, Q2, A0, A2, A4 |
59 |
PARAMETER ( Q0 = 1365. _d 0 , Q2 = -0.482 _d 0 ) |
60 |
PARAMETER ( A0 = 0.322 _d 0 , A2 = 0.231 _d 0 , A4 = 0.086 _d 0 ) |
61 |
C- Longwave parameters |
62 |
_RL LW0, LW1 |
63 |
PARAMETER ( LW0 = 195. _d 0 , LW1 = 2.78 _d 0 ) |
64 |
C- Heat transport parameters |
65 |
_RL H1, H3, H5 |
66 |
PARAMETER ( H1 = 3.866 _d 0 , H3 = -2.851 _d 0, H5 = -1.016 _d 0 ) |
67 |
C- Freshwater flux parameters |
68 |
_RL F1, F2, F3, F4, F5 |
69 |
PARAMETER ( F1 = 2.092 _d 0 , F2 = 5.796 _d 0 , F3 = 8.472 _d 0 ) |
70 |
PARAMETER ( F4 = 7.728 _d 0 , F5 = 2.362 _d 0 ) |
71 |
C- Parameters used to calculate the transport efficiency |
72 |
_RL trans_eff, Hw, Nw, Tw, At, dz, htil, tau |
73 |
PARAMETER ( trans_eff = 2. _d 0 , Hw = 8. _d 3, Nw = 0.013 _d 0 ) |
74 |
PARAMETER ( Tw = 300. _d 0 , At = 3. _d 0 , dz = 450. _d 0 ) |
75 |
PARAMETER ( htil = 0.8 _d 0 , tau = 5. _d -3) |
76 |
C- Climate change parameter (W/(m^2 y) ) |
77 |
_RL DLW |
78 |
PARAMETER (DLW = 0.06 _d 0 ) |
79 |
|
80 |
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |