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

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

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

revision 1.2 by dimitri, Sat Dec 28 10:11:11 2002 UTC revision 1.3 by dimitri, Tue Feb 18 05:33:54 2003 UTC
# Line 12  c                      - empirical param Line 12  c                      - empirical param
12  c                      - control parameters (e.g. max. no of iteration)  c                      - control parameters (e.g. max. no of iteration)
13  c  c
14  c     started: Patrick Heimbach heimbach@mit.edu  06-May-2000  c     started: Patrick Heimbach heimbach@mit.edu  06-May-2000
15  c     changed: Dimitris Menemenlis menemenlis@jpl.nasa.gov 20-Dec-2002  c     mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
16  c  c
17  c     ==================================================================  c     ==================================================================
18  c     HEADER exf_constants  c     HEADER exf_constants
# Line 29  c     exf_two    2.0 Line 29  c     exf_two    2.0
29        _RL exf_two        _RL exf_two
30    
31        parameter(        parameter(
32       &              exf_half =  0.5,       &              exf_half =  0.5 _d 0 ,
33       &              exf_one  =  1.0,       &              exf_one  =  1.0 _d 0 ,
34       &              exf_two  =  2.0       &              exf_two  =  2.0 _d 0
35       &         )       &         )
36    
   
37        real       exf_undef        real       exf_undef
38        parameter( exf_undef = -9000. )        parameter( exf_undef = -9000. )
39    
# Line 43  c     2. physical constants Line 42  c     2. physical constants
42  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
43  c     is identical to "gravity" used in MITgcmUV  c     is identical to "gravity" used in MITgcmUV
44  c     needs to be marmonized through common constants.h file  c     needs to be marmonized through common constants.h file
45        _RL     gravity_mks        _RL         gravity_mks
46        parameter (gravity_mks = 9.81D0)        parameter ( gravity_mks = 9.81 _d 0 )
47  #endif  #endif
48    
49  c     3. empirical parameters  c     3. empirical parameters
50    
51        _RL         climtempfreeze        _RL         climtempfreeze
52        parameter ( climtempfreeze = -1.9 )        parameter ( climtempfreeze = -1.9 _d 0 )
53    
54  #ifdef ALLOW_BULKFORMULAE  #ifdef ALLOW_BULKFORMULAE
55    
56  c     atmrho       - mean atmospheric density [kg/(m*3)]  c     atmrho       - mean atmospheric density [kg/(m*3)]
57  c     atmcp        - mean atmospheric specific heat [ ? ]  c     atmcp        - mean atmospheric specific heat [J/kg/deg K]
58  c     flamb        - latent heat of evaporation [ ? ]  c     flamb        - latent heat of evaporation [J/kg]
59  c     cdrag_[n]    - n = 1,2,3 coefficients used to evaluate  c     cdrag_[n]    - n = 1,2,3 coefficients used to evaluate
60  c                    drag coefficient  c                    drag coefficient
61  c     cstanton_[n] - n = 1,2   coefficients used to evaluate  c     cstanton_[n] - n = 1,2   coefficients used to evaluate
# Line 98  c     hq           - height of mean rel. Line 97  c     hq           - height of mean rel.
97        _RL ht        _RL ht
98        _RL hq        _RL hq
99    
100        parameter ( cdrag_1     =       0.0027000,        parameter ( cdrag_1     =       0.0027000 _d 0 ,
101       &            cdrag_2     =       0.0001420,       &            cdrag_2     =       0.0001420 _d 0 ,
102       &            cdrag_3     =       0.0000764,       &            cdrag_3     =       0.0000764 _d 0 ,
103       &            cstanton_1  =       0.0327000,       &            cstanton_1  =       0.0327000 _d 0 ,
104       &            cstanton_2  =       0.0180000,       &            cstanton_2  =       0.0180000 _d 0 ,
105       &            cdalton     =       0.0346000,       &            cdalton     =       0.0346000 _d 0 ,
106       &            atmrho      =       1.200,       &            atmrho      =       1.200     _d 0 ,
107       &            atmcp       =    1005.000,       &            atmcp       =    1005.000     _d 0 ,
108       &            flamb       = 2500000.000,       &            flamb       = 2500000.000     _d 0 ,
109       &            umin        =       0.500,       &            umin        =       0.500     _d 0 ,
110       &            zolmin      =    -100.000,       &            zolmin      =    -100.000     _d 0 ,
111       &            zref        =      10.000,       &            zref        =      10.000     _d 0 ,
112       &            karman      =       0.400,       &            karman      =       0.400     _d 0 ,
113       &            cvapor_fac  =  640380.000,       &            cvapor_fac  =  640380.000     _d 0 ,
114       &            cvapor_exp  =    5107.400,       &            cvapor_exp  =    5107.400     _d 0 ,
115       &            humid_fac   =       0.606,       &            humid_fac   =       0.606     _d 0 ,
116       &            gamma_blk   =       0.010,       &            gamma_blk   =       0.010     _d 0 ,
117       &            saltsat     =       0.980,       &            saltsat     =       0.980     _d 0 ,
118       &            psim_fac    =       5.000,       &            psim_fac    =       5.000     _d 0 ,
119       &            cen2kel     =     273.150,       &            cen2kel     =     273.150     _d 0 ,
120       &            hu          =      10.000,       &            hu          =      10.000     _d 0 ,
121       &            ht          =       2.000,       &            ht          =       2.000     _d 0 ,
122       &            hq          =       2.000       &            hq          =       2.000     _d 0
123       &          )       &          )
124    
125    
# Line 148  c       ustofu11     - ustar = 0.3818 m/ Line 147  c       ustofu11     - ustar = 0.3818 m/
147        _RL ustofu11        _RL ustofu11
148    
149        parameter (        parameter (
150       &            ustofu11    =         0.381800 ,       &            ustofu11    =         0.381800 _d 0 ,
151       &            u11         =        11.       ,       &            u11         =        11.       _d 0 ,
152       &            clindrag_1  =         0.000065 ,       &            clindrag_1  =         0.000065 _d 0 ,
153       &            clindrag_2  =         0.000490 ,       &            clindrag_2  =         0.000490 _d 0 ,
154       &            cquadrag_1  = clindrag_1/u11/2 ,       &            cquadrag_1  = clindrag_1/u11/2 ,
155       &            cquadrag_2  = clindrag_1*u11/2 + clindrag_2       &            cquadrag_2  = clindrag_1*u11/2 + clindrag_2
156       &          )       &          )
# Line 163  c       ustofu11     - ustar = 0.3818 m/ Line 162  c       ustofu11     - ustar = 0.3818 m/
162        parameter ( czol = hu*karman*gravity_mks )        parameter ( czol = hu*karman*gravity_mks )
163  #endif  #endif
164    
 #endif  
   
165  c     4. control parameters  c     4. control parameters
166    
   
 #ifdef ALLOW_BULKFORMULAE  
   
167  c     niter_bulk   - Number of iterations to be performed for the  c     niter_bulk   - Number of iterations to be performed for the
168  c                    evaluation of the bulk surface fluxes. The ncom  c                    evaluation of the bulk surface fluxes. The ncom
169  c                    model uses 2 hardwired interation steps (loop  c                    model uses 2 hardwired interation steps (loop
# Line 178  c Line 172  c
172        integer     niter_bulk        integer     niter_bulk
173        parameter ( niter_bulk = 2 )        parameter ( niter_bulk = 2 )
174    
175  #endif  #endif /* ALLOW_BULKFORMULAE */
   
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22