/[MITgcm]/MITgcm/pkg/ebm/EBM.h
ViewVC logotype

Diff of /MITgcm/pkg/ebm/EBM.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.3 by heimbach, Fri May 21 21:45:35 2004 UTC revision 1.4 by jmc, Thu Aug 6 13:52:52 2009 UTC
# Line 3  C $Name$ Line 3  C $Name$
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)
# Line 36  C     winPert - weather  patterns added Line 36  C     winPert - weather  patterns added
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
# Line 46  cph avoid AD recomputations Line 49  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)
# Line 68  C     Parameters used to calculate the t Line 71  C     Parameters used to calculate the t
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    
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22