1 |
c $Header: /u/gcmpack/MITgcm/pkg/exf/exf_param.h,v 1.22 2005/02/21 05:32:55 dimitri Exp $ |
2 |
c |
3 |
c |
4 |
c ================================================================== |
5 |
c HEADER exf_param |
6 |
c ================================================================== |
7 |
c |
8 |
c o Header file for the surface flux data. Used by the external |
9 |
c forcing package. |
10 |
c |
11 |
c started: Christian Eckert eckert@mit.edu 30-Jun-1999 |
12 |
c |
13 |
c changed: Christian Eckert eckert@mit.edu 14-Jan-2000 |
14 |
c - Restructured the original version in order to have a |
15 |
c better interface to the MITgcmUV. |
16 |
c |
17 |
c Christian Eckert eckert@mit.edu 12-Feb-2000 |
18 |
c - Changed some variables names (package prefix: exf_) |
19 |
c |
20 |
c Patrick Heimbach, heimbach@mit.edu 04-May-2000 |
21 |
c - included exf_iprec, exf_yftype to enable easy |
22 |
c switch between 32bit/64 bit data format |
23 |
c |
24 |
c Patrick Heimbach, heimbach@mit.edu 01-May-2001 |
25 |
c - added obcs parameters |
26 |
c |
27 |
c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002 |
28 |
c |
29 |
c ================================================================== |
30 |
c HEADER exf_param |
31 |
c ================================================================== |
32 |
|
33 |
c year in seconds |
34 |
_RL year2sec |
35 |
|
36 |
c Calendar data. |
37 |
_RL repeatPeriod |
38 |
|
39 |
c Sea-water albedo |
40 |
_RL exf_albedo |
41 |
|
42 |
c Maximum absolute windstress, used to reset unreastically high |
43 |
c data values |
44 |
_RL windstressmax |
45 |
|
46 |
integer hfluxstartdate1 |
47 |
integer hfluxstartdate2 |
48 |
_RL hfluxstartdate |
49 |
_RL hfluxperiod |
50 |
_RL hfluxconst |
51 |
_RL hflux_exfremo_intercept |
52 |
_RL hflux_exfremo_slope |
53 |
character*1 hfluxmask |
54 |
parameter( hfluxmask = 's' ) |
55 |
|
56 |
integer atempstartdate1 |
57 |
integer atempstartdate2 |
58 |
_RL atempstartdate |
59 |
_RL atempperiod |
60 |
_RL atempconst |
61 |
_RL atemp_exfremo_intercept |
62 |
_RL atemp_exfremo_slope |
63 |
character*1 atempmask |
64 |
parameter( atempmask = 's' ) |
65 |
|
66 |
integer aqhstartdate1 |
67 |
integer aqhstartdate2 |
68 |
_RL aqhstartdate |
69 |
_RL aqhperiod |
70 |
_RL aqhconst |
71 |
_RL aqh_exfremo_intercept |
72 |
_RL aqh_exfremo_slope |
73 |
character*1 aqhmask |
74 |
parameter( aqhmask = 's' ) |
75 |
|
76 |
integer sfluxstartdate1 |
77 |
integer sfluxstartdate2 |
78 |
_RL sfluxstartdate |
79 |
_RL sfluxperiod |
80 |
_RL sfluxconst |
81 |
_RL sflux_exfremo_intercept |
82 |
_RL sflux_exfremo_slope |
83 |
character*1 sfluxmask |
84 |
parameter( sfluxmask = 's' ) |
85 |
|
86 |
integer evapstartdate1 |
87 |
integer evapstartdate2 |
88 |
_RL evapstartdate |
89 |
_RL evapperiod |
90 |
_RL evapconst |
91 |
_RL evap_exfremo_intercept |
92 |
_RL evap_exfremo_slope |
93 |
character*1 evapmask |
94 |
parameter( evapmask = 's' ) |
95 |
|
96 |
integer precipstartdate1 |
97 |
integer precipstartdate2 |
98 |
_RL precipstartdate |
99 |
_RL precipperiod |
100 |
_RL precipconst |
101 |
_RL precip_exfremo_intercept |
102 |
_RL precip_exfremo_slope |
103 |
character*1 precipmask |
104 |
parameter( precipmask = 's' ) |
105 |
|
106 |
integer runoffstartdate1 |
107 |
integer runoffstartdate2 |
108 |
_RL runoffstartdate |
109 |
_RL runoffperiod |
110 |
_RL runoffconst |
111 |
_RL runoff_exfremo_intercept |
112 |
_RL runoff_exfremo_slope |
113 |
character*1 runoffmask |
114 |
parameter( runoffmask = 's' ) |
115 |
|
116 |
integer ustressstartdate1 |
117 |
integer ustressstartdate2 |
118 |
_RL ustressstartdate |
119 |
_RL ustressperiod |
120 |
_RL ustressconst |
121 |
_RL ustress_exfremo_intercept |
122 |
_RL ustress_exfremo_slope |
123 |
character*1 ustressmask |
124 |
parameter( ustressmask = 'u' ) |
125 |
|
126 |
integer vstressstartdate1 |
127 |
integer vstressstartdate2 |
128 |
_RL vstressstartdate |
129 |
_RL vstressperiod |
130 |
_RL vstressconst |
131 |
_RL vstress_exfremo_intercept |
132 |
_RL vstress_exfremo_slope |
133 |
character*1 vstressmask |
134 |
parameter( vstressmask = 'v' ) |
135 |
|
136 |
integer uwindstartdate1 |
137 |
integer uwindstartdate2 |
138 |
_RL uwindstartdate |
139 |
_RL uwindperiod |
140 |
_RL uwindconst |
141 |
_RL uwind_exfremo_intercept |
142 |
_RL uwind_exfremo_slope |
143 |
character*1 uwindmask |
144 |
parameter( uwindmask = 's' ) |
145 |
|
146 |
integer vwindstartdate1 |
147 |
integer vwindstartdate2 |
148 |
_RL vwindstartdate |
149 |
_RL vwindperiod |
150 |
_RL vwindconst |
151 |
_RL vwind_exfremo_intercept |
152 |
_RL vwind_exfremo_slope |
153 |
character*1 vwindmask |
154 |
parameter( vwindmask = 's' ) |
155 |
|
156 |
integer swfluxstartdate1 |
157 |
integer swfluxstartdate2 |
158 |
_RL swfluxstartdate |
159 |
_RL swfluxperiod |
160 |
_RL swfluxconst |
161 |
_RL swflux_exfremo_intercept |
162 |
_RL swflux_exfremo_slope |
163 |
character*1 swfluxmask |
164 |
parameter( swfluxmask = 's' ) |
165 |
|
166 |
integer lwfluxstartdate1 |
167 |
integer lwfluxstartdate2 |
168 |
_RL lwfluxstartdate |
169 |
_RL lwfluxperiod |
170 |
_RL lwfluxconst |
171 |
_RL lwflux_exfremo_intercept |
172 |
_RL lwflux_exfremo_slope |
173 |
character*1 lwfluxmask |
174 |
parameter( lwfluxmask = 's' ) |
175 |
|
176 |
integer swdownstartdate1 |
177 |
integer swdownstartdate2 |
178 |
_RL swdownstartdate |
179 |
_RL swdownperiod |
180 |
_RL swdownconst |
181 |
_RL swdown_exfremo_intercept |
182 |
_RL swdown_exfremo_slope |
183 |
character*1 swdownmask |
184 |
parameter( swdownmask = 's' ) |
185 |
|
186 |
integer lwdownstartdate1 |
187 |
integer lwdownstartdate2 |
188 |
_RL lwdownstartdate |
189 |
_RL lwdownperiod |
190 |
_RL lwdownconst |
191 |
_RL lwdown_exfremo_intercept |
192 |
_RL lwdown_exfremo_slope |
193 |
character*1 lwdownmask |
194 |
parameter( lwdownmask = 's' ) |
195 |
|
196 |
integer apressurestartdate1 |
197 |
integer apressurestartdate2 |
198 |
_RL apressurestartdate |
199 |
_RL apressureperiod |
200 |
_RL apressureconst |
201 |
_RL apressure_exfremo_intercept |
202 |
_RL apressure_exfremo_slope |
203 |
character*1 apressuremask |
204 |
parameter( apressuremask = 's' ) |
205 |
|
206 |
integer obcsNstartdate1 |
207 |
integer obcsNstartdate2 |
208 |
_RL obcsNstartdate |
209 |
_RL obcsNperiod |
210 |
|
211 |
integer obcsSstartdate1 |
212 |
integer obcsSstartdate2 |
213 |
_RL obcsSstartdate |
214 |
_RL obcsSperiod |
215 |
|
216 |
integer obcsEstartdate1 |
217 |
integer obcsEstartdate2 |
218 |
_RL obcsEstartdate |
219 |
_RL obcsEperiod |
220 |
|
221 |
integer obcsWstartdate1 |
222 |
integer obcsWstartdate2 |
223 |
_RL obcsWstartdate |
224 |
_RL obcsWperiod |
225 |
|
226 |
c File names. |
227 |
character*(128) hfluxfile |
228 |
character*(128) atempfile |
229 |
character*(128) aqhfile |
230 |
character*(128) evapfile |
231 |
character*(128) precipfile |
232 |
character*(128) sfluxfile |
233 |
character*(128) runofffile |
234 |
character*(128) ustressfile |
235 |
character*(128) vstressfile |
236 |
character*(128) uwindfile |
237 |
character*(128) vwindfile |
238 |
character*(128) swfluxfile |
239 |
character*(128) lwfluxfile |
240 |
character*(128) swdownfile |
241 |
character*(128) lwdownfile |
242 |
character*(128) apressurefile |
243 |
|
244 |
C useExfYearlyFields :: when set, automatically add extension |
245 |
C _YEAR to input file names |
246 |
C twoDigitYear :: when set, use 2-digit year extension YR |
247 |
C instead of _YEAR for useExfYearlyFields |
248 |
logical useExfYearlyFields, twoDigitYear |
249 |
logical useExfCheckRange |
250 |
|
251 |
common /exf_param_l/ |
252 |
& useExfYearlyFields, twoDigitYear, |
253 |
& useExfCheckRange |
254 |
common /exf_param_i/ |
255 |
& hfluxstartdate1, hfluxstartdate2, |
256 |
& atempstartdate1, atempstartdate2, |
257 |
& aqhstartdate1, aqhstartdate2, |
258 |
& sfluxstartdate1, sfluxstartdate2, |
259 |
& evapstartdate1, evapstartdate2, |
260 |
& runoffstartdate1, runoffstartdate2, |
261 |
& precipstartdate1, precipstartdate2, |
262 |
& ustressstartdate1, ustressstartdate2, |
263 |
& vstressstartdate1, vstressstartdate2, |
264 |
& uwindstartdate1, uwindstartdate2, |
265 |
& vwindstartdate1, vwindstartdate2, |
266 |
& swfluxstartdate1, swfluxstartdate2, |
267 |
& lwfluxstartdate1, lwfluxstartdate2, |
268 |
& swdownstartdate1, swdownstartdate2, |
269 |
& lwdownstartdate1, lwdownstartdate2, |
270 |
& obcsNstartdate1, obcsNstartdate2, |
271 |
& obcsSstartdate1, obcsSstartdate2, |
272 |
& obcsEstartdate1, obcsEstartdate2, |
273 |
& obcsWstartdate1, obcsWstartdate2, |
274 |
& apressurestartdate1,apressurestartdate2 |
275 |
|
276 |
common /exf_param_r/ |
277 |
& year2sec, windstressmax, |
278 |
& repeatPeriod, exf_albedo, |
279 |
& hfluxperiod, hfluxstartdate, |
280 |
& atempperiod, atempstartdate, |
281 |
& aqhperiod, aqhstartdate, |
282 |
& sfluxperiod, sfluxstartdate, |
283 |
& evapperiod, evapstartdate, |
284 |
& precipperiod, precipstartdate, |
285 |
& runoffperiod, runoffstartdate, |
286 |
& ustressperiod, ustressstartdate, |
287 |
& vstressperiod, vstressstartdate, |
288 |
& uwindperiod, uwindstartdate, |
289 |
& vwindperiod, vwindstartdate, |
290 |
& swfluxperiod, swfluxstartdate, |
291 |
& lwfluxperiod, lwfluxstartdate, |
292 |
& swdownperiod, swdownstartdate, |
293 |
& lwdownperiod, lwdownstartdate, |
294 |
& obcsNperiod, obcsNstartdate, |
295 |
& obcsSperiod, obcsSstartdate, |
296 |
& obcsEperiod, obcsEstartdate, |
297 |
& obcsWperiod, obcsWstartdate, |
298 |
& apressureperiod, apressurestartdate, |
299 |
& hfluxconst, |
300 |
& atempconst, |
301 |
& aqhconst, |
302 |
& sfluxconst, |
303 |
& evapconst, |
304 |
& precipconst, |
305 |
& runoffconst, |
306 |
& ustressconst, |
307 |
& vstressconst, |
308 |
& uwindconst, |
309 |
& vwindconst, |
310 |
& swfluxconst, |
311 |
& lwfluxconst, |
312 |
& swdownconst, |
313 |
& lwdownconst, |
314 |
& apressureconst |
315 |
|
316 |
common /exf_param_trend_removal/ |
317 |
& hflux_exfremo_intercept, |
318 |
& atemp_exfremo_intercept, |
319 |
& aqh_exfremo_intercept, |
320 |
& sflux_exfremo_intercept, |
321 |
& evap_exfremo_intercept, |
322 |
& precip_exfremo_intercept, |
323 |
& runoff_exfremo_intercept, |
324 |
& ustress_exfremo_intercept, |
325 |
& vstress_exfremo_intercept, |
326 |
& uwind_exfremo_intercept, |
327 |
& vwind_exfremo_intercept, |
328 |
& swflux_exfremo_intercept, |
329 |
& lwflux_exfremo_intercept, |
330 |
& swdown_exfremo_intercept, |
331 |
& lwdown_exfremo_intercept, |
332 |
& apressure_exfremo_intercept, |
333 |
& hflux_exfremo_slope, |
334 |
& atemp_exfremo_slope, |
335 |
& aqh_exfremo_slope, |
336 |
& sflux_exfremo_slope, |
337 |
& evap_exfremo_slope, |
338 |
& precip_exfremo_slope, |
339 |
& runoff_exfremo_slope, |
340 |
& ustress_exfremo_slope, |
341 |
& vstress_exfremo_slope, |
342 |
& uwind_exfremo_slope, |
343 |
& vwind_exfremo_slope, |
344 |
& swflux_exfremo_slope, |
345 |
& lwflux_exfremo_slope, |
346 |
& swdown_exfremo_slope, |
347 |
& lwdown_exfremo_slope, |
348 |
& apressure_exfremo_slope |
349 |
|
350 |
common /exf_param_c/ |
351 |
& hfluxfile, |
352 |
& atempfile, |
353 |
& aqhfile, |
354 |
& sfluxfile, |
355 |
& evapfile, |
356 |
& precipfile, |
357 |
& runofffile, |
358 |
& ustressfile, |
359 |
& vstressfile, |
360 |
& uwindfile, |
361 |
& vwindfile, |
362 |
& swfluxfile, |
363 |
& lwfluxfile, |
364 |
& swdownfile, |
365 |
& lwdownfile, |
366 |
& apressurefile |
367 |
|
368 |
c file precision and field type |
369 |
|
370 |
common /exf_param_type/ |
371 |
& exf_iprec, |
372 |
& exf_yftype |
373 |
|
374 |
integer exf_iprec |
375 |
character*(2) exf_yftype |
376 |
|
377 |
c exf_inscal_* input scaling factors |
378 |
c exf_offset_atemp input air temperature offset |
379 |
c (for conversion from C to K, if needed) |
380 |
c exf_outscale_* output scaling factors |
381 |
|
382 |
_RL exf_inscal_hflux |
383 |
_RL exf_inscal_sflux |
384 |
_RL exf_inscal_ustress |
385 |
_RL exf_inscal_vstress |
386 |
_RL exf_inscal_uwind |
387 |
_RL exf_inscal_vwind |
388 |
_RL exf_inscal_swflux |
389 |
_RL exf_inscal_lwflux |
390 |
_RL exf_inscal_precip |
391 |
_RL exf_inscal_sst |
392 |
_RL exf_inscal_sss |
393 |
_RL exf_inscal_atemp |
394 |
_RL exf_offset_atemp |
395 |
_RL exf_inscal_aqh |
396 |
_RL exf_inscal_evap |
397 |
_RL exf_inscal_apressure |
398 |
_RL exf_inscal_runoff |
399 |
_RL exf_inscal_swdown |
400 |
_RL exf_inscal_lwdown |
401 |
|
402 |
_RL exf_outscal_hflux |
403 |
_RL exf_outscal_sflux |
404 |
_RL exf_outscal_ustress |
405 |
_RL exf_outscal_vstress |
406 |
_RL exf_outscal_swflux |
407 |
_RL exf_outscal_sst |
408 |
_RL exf_outscal_sss |
409 |
_RL exf_outscal_apressure |
410 |
|
411 |
common /exf_param_scal/ |
412 |
& exf_inscal_hflux |
413 |
& , exf_inscal_sflux |
414 |
& , exf_inscal_ustress |
415 |
& , exf_inscal_vstress |
416 |
& , exf_inscal_uwind |
417 |
& , exf_inscal_vwind |
418 |
& , exf_inscal_swflux |
419 |
& , exf_inscal_lwflux |
420 |
& , exf_inscal_precip |
421 |
& , exf_inscal_sst |
422 |
& , exf_inscal_sss |
423 |
& , exf_inscal_atemp |
424 |
& , exf_offset_atemp |
425 |
& , exf_inscal_aqh |
426 |
& , exf_inscal_evap |
427 |
& , exf_inscal_apressure |
428 |
& , exf_inscal_runoff |
429 |
& , exf_inscal_swdown |
430 |
& , exf_inscal_lwdown |
431 |
& , exf_outscal_hflux |
432 |
& , exf_outscal_sflux |
433 |
& , exf_outscal_ustress |
434 |
& , exf_outscal_vstress |
435 |
& , exf_outscal_swflux |
436 |
& , exf_outscal_sst |
437 |
& , exf_outscal_sss |
438 |
& , exf_outscal_apressure |
439 |
|
440 |
#ifndef USE_EXF_INTERPOLATION |
441 |
c-- set dummy dimension 1 |
442 |
integer MAX_LAT_INC |
443 |
parameter(MAX_LAT_INC = 1) |
444 |
#else |
445 |
c for lat interpolation, arraysize currently set to 2176 max data values |
446 |
integer MAX_LAT_INC |
447 |
parameter(MAX_LAT_INC = 2176) |
448 |
_RL ustress_lon0, ustress_lon_inc |
449 |
_RL ustress_lat0, ustress_lat_inc(MAX_LAT_INC) |
450 |
INTEGER ustress_nlon, ustress_nlat |
451 |
_RL vstress_lon0, vstress_lon_inc |
452 |
_RL vstress_lat0, vstress_lat_inc(MAX_LAT_INC) |
453 |
INTEGER vstress_nlon, vstress_nlat |
454 |
_RL hflux_lon0, hflux_lon_inc |
455 |
_RL hflux_lat0, hflux_lat_inc(MAX_LAT_INC) |
456 |
INTEGER hflux_nlon, hflux_nlat |
457 |
_RL sflux_lon0, sflux_lon_inc |
458 |
_RL sflux_lat0, sflux_lat_inc(MAX_LAT_INC) |
459 |
INTEGER sflux_nlon, sflux_nlat |
460 |
_RL swflux_lon0, swflux_lon_inc |
461 |
_RL swflux_lat0, swflux_lat_inc(MAX_LAT_INC) |
462 |
INTEGER swflux_nlon, swflux_nlat |
463 |
_RL runoff_lon0, runoff_lon_inc |
464 |
_RL runoff_lat0, runoff_lat_inc(MAX_LAT_INC) |
465 |
INTEGER runoff_nlon, runoff_nlat |
466 |
_RL atemp_lon0, atemp_lon_inc |
467 |
_RL atemp_lat0, atemp_lat_inc(MAX_LAT_INC) |
468 |
INTEGER atemp_nlon, atemp_nlat |
469 |
_RL aqh_lon0, aqh_lon_inc |
470 |
_RL aqh_lat0, aqh_lat_inc(MAX_LAT_INC) |
471 |
INTEGER aqh_nlon, aqh_nlat |
472 |
_RL evap_lon0, evap_lon_inc |
473 |
_RL evap_lat0, evap_lat_inc(MAX_LAT_INC) |
474 |
INTEGER evap_nlon, evap_nlat |
475 |
_RL precip_lon0, precip_lon_inc |
476 |
_RL precip_lat0, precip_lat_inc(MAX_LAT_INC) |
477 |
INTEGER precip_nlon, precip_nlat |
478 |
_RL uwind_lon0, uwind_lon_inc |
479 |
_RL uwind_lat0, uwind_lat_inc(MAX_LAT_INC) |
480 |
INTEGER uwind_nlon, uwind_nlat |
481 |
_RL vwind_lon0, vwind_lon_inc |
482 |
_RL vwind_lat0, vwind_lat_inc(MAX_LAT_INC) |
483 |
INTEGER vwind_nlon, vwind_nlat |
484 |
_RL lwflux_lon0, lwflux_lon_inc |
485 |
_RL lwflux_lat0, lwflux_lat_inc(MAX_LAT_INC) |
486 |
INTEGER lwflux_nlon, lwflux_nlat |
487 |
_RL swdown_lon0, swdown_lon_inc |
488 |
_RL swdown_lat0, swdown_lat_inc(MAX_LAT_INC) |
489 |
INTEGER swdown_nlon, swdown_nlat |
490 |
_RL lwdown_lon0, lwdown_lon_inc |
491 |
_RL lwdown_lat0, lwdown_lat_inc(MAX_LAT_INC) |
492 |
INTEGER lwdown_nlon, lwdown_nlat |
493 |
_RL apressure_lon0,apressure_lon_inc |
494 |
_RL apressure_lat0,apressure_lat_inc(MAX_LAT_INC) |
495 |
INTEGER apressure_nlon,apressure_nlat |
496 |
|
497 |
common /exf_interpolation/ |
498 |
& ustress_lon0, ustress_lon_inc, |
499 |
& ustress_lat0, ustress_lat_inc, |
500 |
& ustress_nlon, ustress_nlat, |
501 |
& vstress_lon0, vstress_lon_inc, |
502 |
& vstress_lat0, vstress_lat_inc, |
503 |
& vstress_nlon, vstress_nlat, |
504 |
& hflux_lon0, hflux_lon_inc, |
505 |
& hflux_lat0, hflux_lat_inc, |
506 |
& hflux_nlon, hflux_nlat, |
507 |
& sflux_lon0, sflux_lon_inc, |
508 |
& sflux_lat0, sflux_lat_inc, |
509 |
& sflux_nlon, sflux_nlat, |
510 |
& swflux_lon0, swflux_lon_inc, |
511 |
& swflux_lat0, swflux_lat_inc, |
512 |
& swflux_nlon, swflux_nlat, |
513 |
& runoff_lon0, runoff_lon_inc, |
514 |
& runoff_lat0, runoff_lat_inc, |
515 |
& runoff_nlon, runoff_nlat, |
516 |
& atemp_lon0, atemp_lon_inc, |
517 |
& atemp_lat0, atemp_lat_inc, |
518 |
& atemp_nlon, atemp_nlat, |
519 |
& aqh_lon0, aqh_lon_inc, |
520 |
& aqh_lat0, aqh_lat_inc, |
521 |
& aqh_nlon, aqh_nlat, |
522 |
& evap_lon0, evap_lon_inc, |
523 |
& evap_lat0, evap_lat_inc, |
524 |
& evap_nlon, evap_nlat, |
525 |
& precip_lon0, precip_lon_inc, |
526 |
& precip_lat0, precip_lat_inc, |
527 |
& precip_nlon, precip_nlat, |
528 |
& uwind_lon0, uwind_lon_inc, |
529 |
& uwind_lat0, uwind_lat_inc, |
530 |
& uwind_nlon, uwind_nlat, |
531 |
& vwind_lon0, vwind_lon_inc, |
532 |
& vwind_lat0, vwind_lat_inc, |
533 |
& vwind_nlon, vwind_nlat, |
534 |
& lwflux_lon0, lwflux_lon_inc, |
535 |
& lwflux_lat0, lwflux_lat_inc, |
536 |
& lwflux_nlon, lwflux_nlat, |
537 |
& swdown_lon0, swdown_lon_inc, |
538 |
& swdown_lat0, swdown_lat_inc, |
539 |
& swdown_nlon, swdown_nlat, |
540 |
& lwdown_lon0, lwdown_lon_inc, |
541 |
& lwdown_lat0, lwdown_lat_inc, |
542 |
& lwdown_nlon, lwdown_nlat, |
543 |
& apressure_lon0,apressure_lon_inc, |
544 |
& apressure_lat0,apressure_lat_inc, |
545 |
& apressure_nlon,apressure_nlat |
546 |
#endif |