/[MITgcm]/MITgcm/pkg/bulk_force/BULKF_PARAMS.h
ViewVC logotype

Contents of /MITgcm/pkg/bulk_force/BULKF_PARAMS.h

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

Revision 1.2 - (show annotations) (download)
Sun Jan 22 15:51:34 2006 UTC (18 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint58a_post, checkpoint58b_post
Changes since 1.1: +75 -17 lines
File MIME type: text/plain 
- add option and parameters to use AIM surface-flux formulae.
- Change loading part: S/R BULKF_FIELDS_LOAD only take care of bulkf_fields ;
  others forcing fields (file-name defined in PARM05, parameter file "data")
  are loaded from S/R EXTERNAL_FIELDS_LOAD, whether or not pkg bulk-force is used.
- initialise all bulkf_fields in bulkf_init.F ; do in-lining of exf_bulkcdn.F ;
- use the right EXCH call for uwind,vwind (to work on CS-grid)
- re-arrange header files (move parameters from BULKF.h to BULKF_PARAMS.h)
  and parameters (note: calcWindStress replaces .NOT.readwindstress).
1 C $Header: /u/gcmpack/MITgcm/pkg/bulk_force/BULKF_PARAMS.h,v 1.1 2003/11/23 01:36:55 jmc Exp $
2 C $Name: $
3
4 #ifdef ALLOW_BULK_FORCE
5 C *==========================================================*
6 C | BULK_PARAMS.h
7 C | o Header file for BULK_FORC package parameters:
8 C | - basic parameter ( I/O frequency, etc ...)
9 C | - physical constants (not used in therm_seaice pkg)
10 C *==========================================================*
11
12 C-- COMMON / BULK_PHYS_CONST / physical constants
13 C.. densities
14 C rhoa :: density of air (kg/m^3)
15 C rhofw :: density of fresh water (kg/m^3)
16 C.. specific heats
17 C cpair :: specific heat of air (J/kg/K)
18 C cpwv :: specific heat of water vapour (J/kg/KC)
19 C .. latent heat
20 C Lvap :: latent heat of vaporizn at 0 C (J/kg)
21 C Lfresh :: latent heat of melting of pure ice (J/kg)
22 C .. melting
23 C Tf0kel :: Freezing temp of fresh ice in Kelvin = 273.15
24 C .. wind drag
25 C drag_[n] :: n = 1,2,3 coefficients used to evaluate drag coefficient
26 C .. constants
27 C stefan :: Stefan-Boltzmann constant (W/m^2 K^4)
28 C xkar :: Von Karman constant QQ units?
29 C Rvap :: gas constant for H2O vapor (J/kg/K)
30 C.. Miscellaneous
31 C p0 :: surface pressure (mb)
32 C.. Combinations used for efficiency
33 C Qcoef :: another constant for latent heat flux
34 C .. for bulk formula
35 C humid_fac :: const. for the evaluation of the virtual temp.
36 C saltsat :: reduction of sat. vapor pressure over salt water
37 C gamma_blk :: adiabatic lapse rate
38 C .. for Vince bulk formula QQ check units
39 C Lvap_ice :: latent heat from sublimation
40 C Rgas :: gas constant for dry air
41 C Sha ::
42 C .. emissivities
43 C atm_emissivity ::
44 C ocean_emissivity ::
45 C snow_emissivity ::
46 C ice_emissivity ::
47 C .. for BULKF_FORMULA_AIM
48 C FWIND0 :: ratio of near-sfc wind to lowest-level wind
49 C CHS :: heat exchange coefficient over sea
50 C VGUST :: wind speed for sub-grid-scale gusts
51 C DTHETA :: Potential temp. gradient for stability correction
52 C dTstab :: potential temp. increment for stability function derivative
53 C FSTAB :: Amplitude of stability correction (fraction)
54 C .. Albedo
55 C ocean_albedo :: ocean surface albedo [0-1]
56
57 _RL rhoa
58 _RL rhosw
59 _RL rhofw
60 _RL cpair
61 _RL cpwv
62 _RL Lvap
63 _RL Lfresh
64 _RL Tf0kel
65 _RL cdrag_1, cdrag_2, cdrag_3
66 _RL stefan
67 _RL xkar
68 _RL Rvap
69 _RL p0
70 _RL Qcoef
71 _RL humid_fac
72 _RL saltsat
73 _RL gamma_blk
74 _RL Lvap_ice
75 _RL Rgas
76 _RL Sha
77 _RL atm_emissivity
78 _RL ocean_emissivity
79 _RL snow_emissivity
80 _RL ice_emissivity
81 #ifdef ALLOW_FORMULA_AIM
82 _RL FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB
83 #endif
84 _RL ocean_albedo
85
86 COMMON / BULK_PHYS_CONST /
87 & rhoa, rhosw, rhofw,
88 & cpair, cpwv,
89 & Lvap, Lfresh,
90 & Tf0kel,
91 & cdrag_1, cdrag_2, cdrag_3,
92 & stefan, xkar, Rvap,
93 & p0, Qcoef,
94 & humid_fac, saltsat, gamma_blk,
95 & Lvap_ice, Rgas, Sha,
96 & atm_emissivity, ocean_emissivity,
97 & snow_emissivity, ice_emissivity,
98 #ifdef ALLOW_FORMULA_AIM
99 & FWIND0, CHS, VGUST, DTHETA, dTstab, FSTAB,
100 #endif
101 & ocean_albedo
102
103 C-- COMMON / BULK_PAR_R / real parameter
104
105 _RL blk_taveFreq
106
107 COMMON / BULK_PAR_R /
108 & blk_taveFreq
109
110 C-- COMMON / BULK_PAR_L / Logical parameters
111 C .. for BULKF_FORMULA_AIM
112 C calcWindStress :: True to calculate Wind-Stress from surface wind
113 C useFluxFormula_AIM :: set to T when using AIM flux formula rather
114 C than the default formula (LANL)
115
116 LOGICAL calcWindStress
117 LOGICAL useFluxFormula_AIM
118
119 COMMON / BULK_PAR_L /
120 & calcWindStress,
121 & useFluxFormula_AIM
122
123 C-- COMMON / BULK_PAR_C / Character string parameters
124 CHARACTER*(MAX_LEN_FNAM) AirTempFile
125 CHARACTER*(MAX_LEN_FNAM) AirHumidityFile
126 CHARACTER*(MAX_LEN_FNAM) RainFile
127 CHARACTER*(MAX_LEN_FNAM) SolarFile
128 CHARACTER*(MAX_LEN_FNAM) LongwaveFile
129 CHARACTER*(MAX_LEN_FNAM) UWindFile
130 CHARACTER*(MAX_LEN_FNAM) VWindFile
131 CHARACTER*(MAX_LEN_FNAM) WSpeedFile
132 CHARACTER*(MAX_LEN_FNAM) RunoffFile
133 CHARACTER*(MAX_LEN_FNAM) QnetFile
134 CHARACTER*(MAX_LEN_FNAM) EmPFile
135 CHARACTER*(MAX_LEN_FNAM) CloudFile
136 CHARACTER*(MAX_LEN_FNAM) SnowFile
137 CHARACTER*(MAX_LEN_FNAM) airPotTempFile
138
139 COMMON / BULK_PAR_C /
140 & AirTempFile, AirHumidityFile, RainFile,
141 & SolarFile, LongwaveFile, UWindFile, VWindFile,
142 & WSpeedFile, RunoffFile,
143 & QnetFile, EmPFile, CloudFile, SnowFile,
144 & airPotTempFile
145
146 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
147
148 #endif /*ALLOW_BULK_FORCE*/
  ViewVC Help
Powered by ViewVC 1.1.22