/[MITgcm]/MITgcm/pkg/exf/exf_constants.h
ViewVC logotype

Annotation of /MITgcm/pkg/exf/exf_constants.h

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

Revision 1.3 - (hide annotations) (download)
Tue Feb 18 05:33:54 2003 UTC (21 years, 4 months ago) by dimitri
Branch: MAIN
CVS Tags: checkpoint48f_post, checkpoint51k_post, checkpoint51l_post, checkpoint51j_post, checkpoint52l_pre, checkpoint48i_post, checkpoint52e_pre, hrcube4, checkpoint52j_post, checkpoint51o_pre, checkpoint50e_post, checkpoint52e_post, checkpoint50c_post, checkpoint51n_pre, checkpoint52d_pre, checkpoint51f_pre, checkpoint48h_post, checkpoint50c_pre, branchpoint-genmake2, checkpoint52j_pre, checkpoint50h_post, branch-netcdf, checkpoint50d_pre, checkpoint51r_post, checkpoint52b_pre, checkpoint51i_post, checkpoint50b_pre, checkpoint51e_post, checkpoint51b_post, checkpoint51l_pre, checkpoint52m_post, checkpoint51c_post, hrcube5, checkpoint49, checkpoint50i_post, checkpoint51o_post, checkpoint48g_post, checkpoint51q_post, checkpoint52l_post, checkpoint52k_post, checkpoint51, checkpoint50, checkpoint52, checkpoint50d_post, checkpoint52d_post, checkpoint51b_pre, checkpoint52a_post, checkpoint52b_post, checkpoint52f_post, checkpoint52c_post, checkpoint51h_pre, checkpoint50g_post, checkpoint51g_post, ecco_c52_e35, checkpoint51f_post, checkpoint50b_post, checkpoint50f_post, checkpoint50f_pre, checkpoint52a_pre, checkpoint51d_post, checkpoint50a_post, checkpoint51m_post, checkpoint51t_post, checkpoint52i_post, checkpoint51a_post, checkpoint50e_pre, checkpoint51p_post, checkpoint51n_post, checkpoint51i_pre, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, hrcube_1, hrcube_2, hrcube_3, checkpoint51s_post
Branch point for: netcdf-sm0, branch-genmake2, branch-nonh, tg2-branch, checkpoint51n_branch
Changes since 1.2: +38 -45 lines
File MIME type: text/plain 
Merging from release1_p12:
o Modifications for using pkg/exf with pkg/seaice
  - improved description of the various forcing configurations
  - added basic radiation bulk formulae to pkg/exf
  - units/sign fix for evap computation in exf_getffields.F
  - updated verification/global_with_exf/results/output.txt
o Added pkg/sbo for computing IERS Special Bureau for the Oceans
  (SBO) core products, including oceanic mass, center-of-mass,
  angular, and bottom pressure (see pkg/sbo/README.sbo).
o Lower bound for viscosity/diffusivity in pkg/kpp/kpp_routines.F
  to avoid negative values in shallow regions.
  - updated verification/natl_box/results/output.txt
  - updated verification/lab_sea/results/output.txt
o MPI gather, scatter: eesupp/src/gather_2d.F and scatter_2d.F
o Added useSingleCpuIO option (see PARAMS.h).
o Updated useSingleCpuIO option in mdsio_writefield.F to
  work with multi-field files, e.g., for single-file pickup.
o pkg/seaice:
  - bug fix in growth.F: QNET for no shortwave case
  - added HeffFile for specifying initial sea-ice thickness
  - changed SEAICE_EXTERNAL_FLUXES wind stress implementation
o Added missing /* */ to CPP comments in pkg/seaice, pkg/exf,
  kpp_transport_t.F, forward_step.F, and the_main_loop.F
o pkg/seaice:
  - adjoint-friendly modifications
  - added a SEAICE_WRITE_PICKUP at end of the_model_main.F
1 dimitri 1.3 c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_constants.h,v 1.1.4.3 2003/01/30 05:55:53 dimitri Exp $
2 heimbach 1.1 c
3     c
4     c ==================================================================
5     c HEADER exf_constants
6     c ==================================================================
7     c
8     c o Header file for constants.
9     c These include - numbers (e.g. 1, 2, 1/2, ...)
10     c - physical constants (e.g. gravitational const.)
11     c - empirical parameters
12     c - control parameters (e.g. max. no of iteration)
13     c
14     c started: Patrick Heimbach heimbach@mit.edu 06-May-2000
15 dimitri 1.3 c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
16 heimbach 1.1 c
17     c ==================================================================
18     c HEADER exf_constants
19     c ==================================================================
20    
21     c 1. numbers
22    
23     c exf_half 0.5
24     c exf_one 1.0
25     c exf_two 2.0
26    
27     _RL exf_half
28     _RL exf_one
29     _RL exf_two
30    
31     parameter(
32 dimitri 1.3 & exf_half = 0.5 _d 0 ,
33     & exf_one = 1.0 _d 0 ,
34     & exf_two = 2.0 _d 0
35 heimbach 1.1 & )
36    
37     real exf_undef
38     parameter( exf_undef = -9000. )
39    
40     c 2. physical constants
41    
42     #ifdef ALLOW_ATM_TEMP
43     c is identical to "gravity" used in MITgcmUV
44     c needs to be marmonized through common constants.h file
45 dimitri 1.3 _RL gravity_mks
46     parameter ( gravity_mks = 9.81 _d 0 )
47 heimbach 1.1 #endif
48    
49     c 3. empirical parameters
50    
51     _RL climtempfreeze
52 dimitri 1.3 parameter ( climtempfreeze = -1.9 _d 0 )
53 heimbach 1.1
54     #ifdef ALLOW_BULKFORMULAE
55    
56     c atmrho - mean atmospheric density [kg/(m*3)]
57 dimitri 1.3 c atmcp - mean atmospheric specific heat [J/kg/deg K]
58     c flamb - latent heat of evaporation [J/kg]
59 heimbach 1.1 c cdrag_[n] - n = 1,2,3 coefficients used to evaluate
60     c drag coefficient
61     c cstanton_[n] - n = 1,2 coefficients used to evaluate
62     c the Stanton number (stable/unstable cond.)
63     c dalton - coefficient used to evaluate the Dalton number
64     c umin - minimum absolute wind speed used to evaluate
65     c drag coefficient [m/s]
66     c zolmin - minimum stability parameter
67     c zref - reference height
68     c
69     c karman - von Karman constant
70     c cvapor - see e.g. Gill (1982) p.41 Eq. (3.1.15)
71     c humid_fac - constant entering the evaluation of the virtual
72     c temperature
73     c gamma_blk - adiabatic lapse rate
74     c saltsat - reduction of saturation vapor pressure over salt water
75     c psim_fac -
76     c cen2kel - conversion of deg. Centigrade to Kelvin
77     c hu - height of mean wind
78     c ht - height of mean temperature
79     c hq - height of mean rel. humidity
80    
81     _RL atmrho, atmcp
82     _RL flamb
83     _RL cdrag_1, cdrag_2, cdrag_3
84     _RL cstanton_1, cstanton_2
85     _RL cdalton
86     _RL umin
87     _RL zolmin
88     _RL zref
89     _RL karman
90     _RL cvapor_fac, cvapor_exp
91     _RL humid_fac
92     _RL gamma_blk
93     _RL saltsat
94     _RL psim_fac
95     _RL cen2kel
96     _RL hu
97     _RL ht
98     _RL hq
99    
100 dimitri 1.3 parameter ( cdrag_1 = 0.0027000 _d 0 ,
101     & cdrag_2 = 0.0001420 _d 0 ,
102     & cdrag_3 = 0.0000764 _d 0 ,
103     & cstanton_1 = 0.0327000 _d 0 ,
104     & cstanton_2 = 0.0180000 _d 0 ,
105     & cdalton = 0.0346000 _d 0 ,
106     & atmrho = 1.200 _d 0 ,
107     & atmcp = 1005.000 _d 0 ,
108     & flamb = 2500000.000 _d 0 ,
109     & umin = 0.500 _d 0 ,
110     & zolmin = -100.000 _d 0 ,
111     & zref = 10.000 _d 0 ,
112     & karman = 0.400 _d 0 ,
113     & cvapor_fac = 640380.000 _d 0 ,
114     & cvapor_exp = 5107.400 _d 0 ,
115     & humid_fac = 0.606 _d 0 ,
116     & gamma_blk = 0.010 _d 0 ,
117     & saltsat = 0.980 _d 0 ,
118     & psim_fac = 5.000 _d 0 ,
119     & cen2kel = 273.150 _d 0 ,
120     & hu = 10.000 _d 0 ,
121     & ht = 2.000 _d 0 ,
122     & hq = 2.000 _d 0
123 heimbach 1.1 & )
124    
125    
126     #ifndef ALLOW_ATM_WIND
127     #ifdef ALLOW_ATM_TEMP
128     c To invert the relationship ustar = ustar(umagn) the following
129     c parameterization is used:
130     c
131     c ustar**2 = umagn**2 * CDN(umagn)
132     c
133     c / cquadrag_1 * umagn**2 + cquadrag_2; 0 < u < 11 m/s
134     c CDN(umagn) =
135     c \ clindrag_1 * umagn + clindrag_2 ; u > 11 m/s
136     c
137     c clindrag_[n] - n = 1, 2 coefficients used to evaluate
138     c LINEAR relationship of Large and Pond 1981
139     c cquadrag_[n] - n = 1, 2 coefficients used to evaluate
140     c quadratic relationship
141     c u11 - u = 11 m/s wind speed
142     c ustofu11 - ustar = 0.3818 m/s, corresponding to u = 11 m/s
143    
144     _RL clindrag_1, clindrag_2
145     _RL cquadrag_1, cquadrag_2
146     _RL u11
147     _RL ustofu11
148    
149     parameter (
150 dimitri 1.3 & ustofu11 = 0.381800 _d 0 ,
151     & u11 = 11. _d 0 ,
152     & clindrag_1 = 0.000065 _d 0 ,
153     & clindrag_2 = 0.000490 _d 0 ,
154 heimbach 1.1 & cquadrag_1 = clindrag_1/u11/2 ,
155     & cquadrag_2 = clindrag_1*u11/2 + clindrag_2
156     & )
157     #endif
158     #endif
159    
160     #ifdef ALLOW_ATM_TEMP
161     _RL czol
162     parameter ( czol = hu*karman*gravity_mks )
163     #endif
164    
165     c 4. control parameters
166    
167     c niter_bulk - Number of iterations to be performed for the
168     c evaluation of the bulk surface fluxes. The ncom
169     c model uses 2 hardwired interation steps (loop
170     c unrolled).
171     c
172     integer niter_bulk
173     parameter ( niter_bulk = 2 )
174    
175 dimitri 1.3 #endif /* ALLOW_BULKFORMULAE */
  ViewVC Help
Powered by ViewVC 1.1.22