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

Contents of /MITgcm/pkg/exf/exf_fields.h

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


Revision 1.6 - (show annotations) (download)
Fri May 23 18:37:31 2003 UTC (21 years ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint50h_post, checkpoint51e_post, checkpoint51b_post, checkpoint51c_post, checkpoint50i_post, checkpoint51, checkpoint51b_pre, checkpoint50g_post, checkpoint50f_post, checkpoint51d_post, checkpoint51a_post
Changes since 1.5: +14 -16 lines
File MIME type: text/plain
o exf totally restructured
  - replaced all exf_set_ routines by one generic exf_set_gen.F
    and made sure it's still differentiable
  - split exf_getffields.F into new exf_bulkformulae.F
  - exf_inscal_/exf_outscal_ suffixes consistent with field names
  - updated adjoint/makefile and relevant verif. accordingly

1 c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_fields.h,v 1.1.6.6 2003/05/23 16:18:57 heimbach Exp $
2 c
3 c
4 c ==================================================================
5 c HEADER exf_fields
6 c ==================================================================
7 c
8 c o Header file for the surface flux data.
9 c
10 c started: Ralf.Giering@FastOpt.de 25-Mai-2000
11 c changed: field swap in adj. mode; heimbach@mit.edu 10-Jan-2002
12 c included runoff D. Stammer, Nov. 25, 2001
13 c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
14 c
15 c ==================================================================
16 c HEADER exf_fields
17 c ==================================================================
18 c
19 c
20 c Field definitions, units, and sign conventions:
21 c ===============================================
22 c
23 c ustress :: Zonal surface wind stress in N/m^2
24 c > 0 for increase in uVel, which is west to
25 c east for cartesian and spherical polar grids
26 c Typical range: -0.5 < ustress < 0.5
27 c Southwest C-grid U point
28 c Input field
29 c
30 c vstress :: Meridional surface wind stress in N/m^2
31 c > 0 for increase in vVel, which is south to
32 c north for cartesian and spherical polar grids
33 c Typical range: -0.5 < vstress < 0.5
34 c Southwest C-grid V point
35 c Input field
36 c
37 c hflux :: Net upward surface heat flux excluding shortwave in W/m^2
38 c hflux = latent + sensible + lwflux
39 c > 0 for decrease in theta (ocean cooling)
40 c Typical range: -250 < hflux < 600
41 c Southwest C-grid tracer point
42 c Input field
43 c
44 c sflux :: Net upward freshwater flux in m/s
45 c sflux = evap - precip - runoff
46 c > 0 for increase in salt (ocean salinity)
47 c Typical range: -1e-7 < sflux < 1e-7
48 c Southwest C-grid tracer point
49 c Input field
50 c
51 c swflux :: Net upward shortwave radiation in W/m^2
52 c swflux = - ( swdown - ice and snow absorption - reflected )
53 c > 0 for decrease in theta (ocean cooling)
54 c Typical range: -350 < swflux < 0
55 c Southwest C-grid tracer point
56 c Input field
57 c
58 c uwind :: Surface (10-m) zonal wind velocity in m/s
59 c > 0 for increase in uVel, which is west to
60 c east for cartesian and spherical polar grids
61 c Typical range: -10 < uwind < 10
62 c Southwest C-grid U point
63 c Input or input/output field
64 c
65 c vwind :: Surface (10-m) meridional wind velocity in m/s
66 c > 0 for increase in vVel, which is south to
67 c north for cartesian and spherical polar grids
68 c Typical range: -10 < vwind < 10
69 c Southwest C-grid V point
70 c Input or input/output field
71 c
72 c atemp :: Surface (2-m) air temperature in deg K
73 c Typical range: 200 < atemp < 300
74 c Southwest C-grid tracer point
75 c Input or input/output field
76 c
77 c aqh :: Surface (2m) specific humidity in kg/kg
78 c Typical range: 0 < aqh < 0.02
79 c Southwest C-grid tracer point
80 c Input or input/output field
81 c
82 c lwflux :: Net upward longwave radiation in W/m^2
83 c lwflux = - ( lwdown - ice and snow absorption - emitted )
84 c > 0 for decrease in theta (ocean cooling)
85 c Typical range: -20 < lwflux < 170
86 c Southwest C-grid tracer point
87 c Input field
88 c
89 c evap :: Evaporation in m/s
90 c > 0 for increase in salt (ocean salinity)
91 c Typical range: 0 < evap < 2.5e-7
92 c Southwest C-grid tracer point
93 c Input, input/output, or output field
94 c
95 c precip :: Precipitation in m/s
96 c > 0 for decrease in salt (ocean salinity)
97 c Typical range: 0 < precip < 5e-7
98 c Southwest C-grid tracer point
99 c Input or input/output field
100 c
101 c runoff :: River and glacier runoff in m/s
102 c > 0 for decrease in salt (ocean salinity)
103 c Typical range: 0 < runoff < ????
104 c Southwest C-grid tracer point
105 c Input or input/output field
106 c !!! WATCH OUT: Default exf_inscal_runoff !!!
107 c !!! in exf_readparms.F is not 1.0 !!!
108 c
109 c swdown :: Downward shortwave radiation in W/m^2
110 c > 0 for increase in theta (ocean warming)
111 c Typical range: 0 < swdown < 450
112 c Southwest C-grid tracer point
113 c Input/output field
114 c
115 c lwdown :: Downward longwave radiation in W/m^2
116 c > 0 for increase in theta (ocean warming)
117 c Typical range: 50 < lwdown < 450
118 c Southwest C-grid tracer point
119 c Input/output field
120 c
121 c apressure :: Atmospheric pressure field in N/m^2
122 c > 0 for ????
123 c Typical range: ???? < apressure < ????
124 c Southwest C-grid tracer point
125 c Input field
126 c
127 c
128 c NOTES:
129 c ======
130 c
131 c Input and output units and sign conventions can be customized
132 c using variables exf_inscal_* and exf_outscal_*, which are set
133 c by exf_readparms.F
134 c
135 c Output fields fu, fv, Qnet, Qsw, and EmPmR are
136 c defined in FFIELDS.h
137 c
138 c #ifndef SHORTWAVE_HEATING, hflux includes shortwave,
139 c that is, hflux = latent + sensible + lwflux +swflux
140 c
141 c If (EXFwindOnBgrid .EQ. .TRUE.), uwind and vwind are
142 c defined on northeast B-grid U and V points, respectively.
143 c
144 c Arrays *0 and *1 below are used for temporal interpolation.
145 c
146
147 common /exf_stress_r/ ustress, vstress
148 _RL ustress (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
149 _RL vstress (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
150 common /exfl_ustress_r/ ustress0, ustress1
151 _RL ustress0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
152 _RL ustress1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
153 common /exfl_vstress_r/ vstress0, vstress1
154 _RL vstress0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
155 _RL vstress1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
156
157 #ifdef ALLOW_ATM_WIND
158 common /exf_atm_wind_r/ uwind, vwind
159 _RL uwind (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
160 _RL vwind (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
161 common /exfl_uwind_r/ uwind0, uwind1
162 _RL uwind0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
163 _RL uwind1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
164 common /exfl_vwind_r/ vwind0, vwind1
165 _RL vwind0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
166 _RL vwind1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
167 #endif
168
169 common /exf_hsflux_r/ hflux, sflux
170 _RL hflux (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
171 _RL sflux (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
172 common /exfl_hflux_r/ hflux0, hflux1
173 _RL hflux0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
174 _RL hflux1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
175 common /exfl_sflux_r/ sflux0, sflux1
176 _RL sflux0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
177 _RL sflux1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
178
179 #if defined(ALLOW_ATM_TEMP) || defined(EXF_READ_EVAP)
180 common /exf_evap/ evap
181 _RL evap (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
182 #endif
183
184 #ifdef ALLOW_ATM_TEMP
185 common /exf_atm_temp_r/ atemp, aqh, lwflux, precip
186 _RL atemp (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
187 _RL aqh (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
188 _RL lwflux (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
189 _RL precip (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
190 common /exfl_atemp_r/ atemp0, atemp1
191 _RL atemp0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
192 _RL atemp1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
193 common /exfl_aqh_r/ aqh0, aqh1
194 _RL aqh0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
195 _RL aqh1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
196 common /exfl_lwflux_r/ lwflux0, lwflux1
197 _RL lwflux0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
198 _RL lwflux1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
199 common /exfl_precip_r/ precip0, precip1
200 _RL precip0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
201 _RL precip1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
202 #endif
203
204 #if defined(ALLOW_ATM_TEMP) || defined(SHORTWAVE_HEATING)
205 common /exf_swflux_r/ swflux
206 _RL swflux (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
207 common /exfl_swflux_r/ swflux0, swflux1
208 _RL swflux0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
209 _RL swflux1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
210 #endif
211
212 #if defined(ALLOW_ATM_TEMP) || defined(EXF_READ_EVAP)
213 common /exfl_evap_r/ evap0, evap1
214 _RL evap0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
215 _RL evap1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
216 #endif
217
218 #if defined (ALLOW_RUNOFF) || defined (ALLOW_SEAICE)
219 common /exfl_runoff_r/ runoff
220 _RL runoff (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
221 #endif
222
223 #ifdef ALLOW_DOWNWARD_RADIATION
224 common /exf_rad_down_r/
225 & swdown, lwdown, swdown0, swdown1, lwdown0, lwdown1
226 _RL swdown (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
227 _RL lwdown (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
228 _RL swdown0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
229 _RL swdown1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
230 _RL lwdown0 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
231 _RL lwdown1 (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
232 #endif
233
234 #ifdef ATMOSPHERIC_LOADING
235 common /exf_apressure_r/ apressure, apressure0, apressure1
236 _RL apressure (1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
237 _RL apressure0(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
238 _RL apressure1(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
239 #endif

  ViewVC Help
Powered by ViewVC 1.1.22