3 |
C |
C |
4 |
C CountX - number of ocean points in each latitude band |
C CountX - number of ocean points in each latitude band |
5 |
C ZonalMeanSST - zonal mean sea surface temperature |
C ZonalMeanSST - zonal mean sea surface temperature |
6 |
C TmlS - Average mid-latitude temperature in the |
C TmlS - Average mid-latitude temperature in the |
7 |
c Southern Hemisphere |
c Southern Hemisphere |
8 |
C TmlN - Northern |
C TmlN - Northern |
9 |
C TtS - Average tropical temperature in the Southern Hemisphere |
C TtS - Average tropical temperature in the Southern Hemisphere |
10 |
C TtN - Northern |
C TtN - Northern |
11 |
C winPert - weather patterns added to the background wind. |
C winPert - weather patterns added to the background wind. |
12 |
COMMON /EBM_FLD_RL/ |
COMMON /EBM_FLD_RL/ |
13 |
& ZonalMeanSST, CountX, |
& ZonalMeanSST, CountX, |
14 |
& TmlN, TmlS, TtN, TtS |
& TmlN, TmlS, TtN, TtS |
15 |
_RL ZonalMeanSST(1-OLy:sNy+OLy,nSy) |
_RL ZonalMeanSST(1-OLy:sNy+OLy,nSy) |
16 |
_RL CountX(1-OLy:sNy+OLy,nSy) |
_RL CountX(1-OLy:sNy+OLy,nSy) |
17 |
_RL TmlS, TmlN, TtS, TtN |
_RL TmlS, TmlN, TtS, TtN |
18 |
|
|
19 |
COMMON /EBM_FLD_RS/ |
COMMON /EBM_FLD_RS/ |
20 |
& Run, |
& Run, |
21 |
& winPert |
& winPert |
22 |
_RS Run (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
_RS Run (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
36 |
|
|
37 |
C Constant parameters |
C Constant parameters |
38 |
_RL t_mlt, lv, cp, rho_air |
_RL t_mlt, lv, cp, rho_air |
39 |
PARAMETER (t_mlt = 273.15, lv = 2.5e6, cp = 1004., rho_air = 1.27) |
PARAMETER ( t_mlt = 273.15 _d 0 ) |
40 |
|
PARAMETER ( lv = 2.5 _d 6 ) |
41 |
|
PARAMETER ( cp = 1004. _d 0 ) |
42 |
|
PARAMETER ( rho_air = 1.27 _d 0 ) |
43 |
C sin(lat) and Legendre polynomials |
C sin(lat) and Legendre polynomials |
44 |
cph We will make these three (i,j) arrays to |
cph We will make these three (i,j) arrays to |
45 |
cph avoid AD recomputations |
cph avoid AD recomputations |
49 |
C Shortwave and albedo parameters |
C Shortwave and albedo parameters |
50 |
_RL SW(1-OLy:sNy+OLy,nSy) |
_RL SW(1-OLy:sNy+OLy,nSy) |
51 |
_RL Q0, Q2, A0, A2, A4 |
_RL Q0, Q2, A0, A2, A4 |
52 |
PARAMETER ( Q0 = 1365., Q2 = -0.482) |
PARAMETER ( Q0 = 1365. _d 0 , Q2 = -0.482 _d 0 ) |
53 |
PARAMETER ( A0 = 0.322, A2 = 0.231, A4 = 0.086 ) |
PARAMETER ( A0 = 0.322 _d 0 , A2 = 0.231 _d 0 , A4 = 0.086 _d 0 ) |
54 |
C Longwave parameters |
C Longwave parameters |
55 |
_RL LW(1-OLy:sNy+OLy,nSy) |
_RL LW(1-OLy:sNy+OLy,nSy) |
56 |
_RL LW0, LW1 |
_RL LW0, LW1 |
57 |
PARAMETER ( LW0 = 195., LW1 = 2.78 ) |
PARAMETER ( LW0 = 195. _d 0 , LW1 = 2.78 _d 0 ) |
58 |
C Heat transport parameters |
C Heat transport parameters |
59 |
_RL Hd(1-OLy:sNy+OLy,nSy), Hd35(2) |
_RL Hd(1-OLy:sNy+OLy,nSy), Hd35(2) |
60 |
_RL H1, H3, H5 |
_RL H1, H3, H5 |
61 |
PARAMETER ( H1 = 3.866, H3 = -2.851, H5 = -1.016 ) |
PARAMETER ( H1 = 3.866 _d 0 , H3 = -2.851 _d 0 , H5 = -1.016 _d 0 ) |
62 |
C Freshwater flux parameters |
C Freshwater flux parameters |
63 |
_RL Fw(1-OLy:sNy+OLy,nSy), Fw35(2) |
_RL Fw(1-OLy:sNy+OLy,nSy), Fw35(2) |
64 |
_RL F1, F2, F3, F4, F5 |
_RL F1, F2, F3, F4, F5 |
65 |
PARAMETER ( F1 = 2.092, F2 = 5.796, F3 = 8.472, |
PARAMETER ( F1 = 2.092 _d 0 , F2 = 5.796 _d 0 , F3 = 8.472 _d 0 ) |
66 |
& F4 = 7.728, F5 = 2.362 ) |
PARAMETER ( F4 = 7.728 _d 0 , F5 = 2.362 _d 0 ) |
67 |
C Temperature parameterization |
C Temperature parameterization |
68 |
_RL T(1-OLy:sNy+OLy,nSy) |
_RL T(1-OLy:sNy+OLy,nSy) |
69 |
_RL T_var(4), T0(2), T2(2), T35(2), DTDy35(2) |
_RL T_var(4), T0(2), T2(2), T35(2), DTDy35(2) |
71 |
_RL Cl, Cf, Cs, C |
_RL Cl, Cf, Cs, C |
72 |
_RL gamma, kappa, De |
_RL gamma, kappa, De |
73 |
_RL trans_eff, Hw, Nw, Tw, At, dz, htil, tau |
_RL trans_eff, Hw, Nw, Tw, At, dz, htil, tau |
74 |
PARAMETER (trans_eff = 2.0, Hw = 8.e3, Nw = 0.013, |
PARAMETER ( trans_eff = 2. _d 0 , Hw = 8. _d 3, Nw = 0.013 _d 0 ) |
75 |
& Tw = 300., At = 3., dz = 450., htil = 0.8, tau = 5.e-3) |
PARAMETER ( Tw = 300. _d 0 , At = 3. _d 0 , dz = 450. _d 0 ) |
76 |
|
PARAMETER ( htil = 0.8 _d 0 , tau = 5. _d -3) |
77 |
C Climate change parameter (W/(m^2 y) ) |
C Climate change parameter (W/(m^2 y) ) |
78 |
_RL DLW |
_RL DLW |
79 |
PARAMETER (DLW = 0.06) |
PARAMETER (DLW = 0.06 _d 0 ) |
80 |
C Latitude boundaries used |
C Latitude boundaries used |
81 |
_RL lat(3) |
_RL lat(3) |
82 |
DATA lat / 0.0, 35.0, 85.0 / |
DATA lat / 0. _d 0 , 35. _d 0 , 85. _d 0 / |
83 |
|
|