| 1 |
C $Header: /u/gcmpack/MITgcm/pkg/ebm/EBM.h,v 1.3 2004/05/21 21:45:35 heimbach Exp $ |
| 2 |
C $Name: $ |
| 3 |
C |
| 4 |
C CountX - number of ocean points in each latitude band |
| 5 |
C ZonalMeanSST - zonal mean sea surface temperature |
| 6 |
C TmlS - Average mid-latitude temperature in the |
| 7 |
c Southern Hemisphere |
| 8 |
C TmlN - Northern |
| 9 |
C TtS - Average tropical temperature in the Southern Hemisphere |
| 10 |
C TtN - Northern |
| 11 |
C winPert - weather patterns added to the background wind. |
| 12 |
COMMON /EBM_FLD_RL/ |
| 13 |
& ZonalMeanSST, CountX, |
| 14 |
& TmlN, TmlS, TtN, TtS |
| 15 |
_RL ZonalMeanSST(1-OLy:sNy+OLy,nSy) |
| 16 |
_RL CountX(1-OLy:sNy+OLy,nSy) |
| 17 |
_RL TmlS, TmlN, TtS, TtN |
| 18 |
|
| 19 |
COMMON /EBM_FLD_RS/ |
| 20 |
& Run, |
| 21 |
& winPert |
| 22 |
_RS Run (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
| 23 |
_RS winPert (1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy) |
| 24 |
|
| 25 |
COMMON /EBM_PARAM_C/ |
| 26 |
& RunoffFile |
| 27 |
CHARACTER*(MAX_LEN_FNAM) RunoffFile |
| 28 |
|
| 29 |
COMMON /EBM_PARAM_RL/ |
| 30 |
& tauThetaZonRelax, |
| 31 |
& lambdaThetaZonRelax, |
| 32 |
& scale_runoff |
| 33 |
_RL tauThetaZonRelax |
| 34 |
_RL lambdaThetaZonRelax |
| 35 |
_RL scale_runoff |
| 36 |
|
| 37 |
C Constant parameters |
| 38 |
_RL t_mlt, lv, cp, rho_air |
| 39 |
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 |
| 44 |
cph We will make these three (i,j) arrays to |
| 45 |
cph avoid AD recomputations |
| 46 |
_RL S(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSy) |
| 47 |
_RL P2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSy) |
| 48 |
_RL P4(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSy) |
| 49 |
C Shortwave and albedo parameters |
| 50 |
_RL SW(1-OLy:sNy+OLy,nSy) |
| 51 |
_RL Q0, Q2, A0, A2, A4 |
| 52 |
PARAMETER ( Q0 = 1365. _d 0 , Q2 = -0.482 _d 0 ) |
| 53 |
PARAMETER ( A0 = 0.322 _d 0 , A2 = 0.231 _d 0 , A4 = 0.086 _d 0 ) |
| 54 |
C Longwave parameters |
| 55 |
_RL LW(1-OLy:sNy+OLy,nSy) |
| 56 |
_RL LW0, LW1 |
| 57 |
PARAMETER ( LW0 = 195. _d 0 , LW1 = 2.78 _d 0 ) |
| 58 |
C Heat transport parameters |
| 59 |
_RL Hd(1-OLy:sNy+OLy,nSy), Hd35(2) |
| 60 |
_RL H1, H3, H5 |
| 61 |
PARAMETER ( H1 = 3.866 _d 0 , H3 = -2.851 _d 0 , H5 = -1.016 _d 0 ) |
| 62 |
C Freshwater flux parameters |
| 63 |
_RL Fw(1-OLy:sNy+OLy,nSy), Fw35(2) |
| 64 |
_RL F1, F2, F3, F4, F5 |
| 65 |
PARAMETER ( F1 = 2.092 _d 0 , F2 = 5.796 _d 0 , F3 = 8.472 _d 0 ) |
| 66 |
PARAMETER ( F4 = 7.728 _d 0 , F5 = 2.362 _d 0 ) |
| 67 |
C Temperature parameterization |
| 68 |
_RL T(1-OLy:sNy+OLy,nSy) |
| 69 |
_RL T_var(4), T0(2), T2(2), T35(2), DTDy35(2) |
| 70 |
C Parameters used to calculate the transport efficiency |
| 71 |
_RL Cl, Cf, Cs, C |
| 72 |
_RL gamma, kappa, De |
| 73 |
_RL trans_eff, Hw, Nw, Tw, At, dz, htil, tau |
| 74 |
PARAMETER ( trans_eff = 2. _d 0 , Hw = 8. _d 3, Nw = 0.013 _d 0 ) |
| 75 |
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) ) |
| 78 |
_RL DLW |
| 79 |
PARAMETER (DLW = 0.06 _d 0 ) |
| 80 |
C Latitude boundaries used |
| 81 |
_RL lat(3) |
| 82 |
DATA lat / 0. _d 0 , 35. _d 0 , 85. _d 0 / |
| 83 |
|
Parent Directory
| 
Revision Log
| 
Revision Graph