CTBTO FLEXPART WO4 (2015-10-15)
 All Classes Files Functions Variables
fpmetbinary_mod.F90
Go to the documentation of this file.
1 MODULE fpmetbinary_mod
2 
3  !*****************************************************************************
4  ! *
5  ! Contains data and routines for dumping and loading processed met *
6  ! fields. *
7  ! Authors Don Morton (Don.Morton@borealscicomp.com) *
8  ! Delia Arnold (deliona.arnold@gmail.com) *
9  ! *
10  ! 15 Sep 2015 *
11  ! *
12  ! Currently, the only data being dumped and loaded has data structures *
13  ! defined in com_mod.f90. In the future, perhaps it will be necessary *
14  ! to use data structures from other parts of the FLEXPART code system. *
15  ! *
16  ! *
17  ! Most of the data structures from com_mod.f90 that are dumped and *
18  ! loaded have a final dimension of size two, so that they may hold data *
19  ! from two met files. When we dump the contents into a .fp file, we *
20  ! need to specify which of the two to dump. Likewise, when we load *
21  ! from a .fp file, we need to specify which of the two possible indices *
22  ! to load into. *
23  ! *
24  ! Note that these routines need more robustness. For example, what *
25  ! what happens if the filename can't be read or written. Or, what *
26  ! happens if a read or write fails in any way. Right now, it's crash *
27  ! city. *
28  ! *
29  !*****************************************************************************
30 
31  USE com_mod
32  USE conv_mod
33 
34  IMPLICIT NONE
35 
36  ! Users may want to change these IO Unit values if they conflict with other parts
37  ! of code
38  INTEGER, PARAMETER :: IOUNIT_DUMP = 33, IOUNIT_LOAD = 34, &
39  IOUNIT_TEXTOUT = 35
40  PRIVATE iounit_dump, iounit_load, iounit_textout, fpio
41 
42 
43 CONTAINS
44 
45  !*****************************************************************************
46  ! *
47  ! Subroutines fpmetbinary_dump() and fpmetbinary_load() provide the *
48  ! public interface to *
49  ! this module functionality. I created the PRIVATE fpio() because I *
50  ! wanted all interactions with variables to be in one place. The read *
51  ! and write operations need to be done in exactly the same sequence, so *
52  ! I felt like keeping them in the same routine would at least allow for *
53  ! coders to more easily compare the two sequences than if they were *
54  ! separate. *
55  ! *
56  ! As mentioned above, the dumps and loads will, for most variables, *
57  ! need to refer to one of two index values for the last dimension of *
58  ! the array. *
59  ! *
60  !*****************************************************************************
61 
62 
63  SUBROUTINE fpmetbinary_dump(filename, cm_index)
64  CHARACTER(LEN=*), INTENT(IN) :: filename ! Full path for file
65  INTEGER, INTENT(IN) :: cm_index ! Index of last dimension in
66  ! most com_mod variables.
67  ! Should be 1 or 2
68 
69  INTEGER millisecs_start, millisecs_stop, count_rate, count_max
70 
71  CALL system_clock(millisecs_start, count_rate, count_max)
72  OPEN(iounit_dump, file=filename, action='WRITE', status='REPLACE', form="unformatted", access="stream")
73  CALL fpio(iounit_dump, 'DUMP', cm_index)
74  CLOSE(iounit_dump)
75  CALL system_clock(millisecs_stop, count_rate, count_max)
76 
77  !PRINT *, 'Dump walltime secs: ', (millisecs_stop-millisecs_start)/1000.0
78  END SUBROUTINE fpmetbinary_dump
79 
80  SUBROUTINE fpmetbinary_load(filename, cm_index)
81  CHARACTER(LEN=*), INTENT(IN) :: filename ! Full path for file
82  INTEGER, INTENT(IN) :: cm_index ! Index of last dimension in
83  ! most com_mod variables.
84  ! Should be 1 or 2
85 
86  INTEGER millisecs_start, millisecs_stop, count_rate, count_max
87 
88  CALL system_clock(millisecs_start, count_rate, count_max)
89  OPEN(iounit_load, file=filename, action='READ', status='OLD', form="unformatted", access="stream")
90  CALL fpio(iounit_load, 'LOAD', cm_index)
91  CLOSE(iounit_load)
92  CALL system_clock(millisecs_stop, count_rate, count_max)
93  !PRINT *, 'Load walltime secs: ', (millisecs_stop-millisecs_start)/1000.0
94  END SUBROUTINE fpmetbinary_load
95 
96  SUBROUTINE fpmetbinary_zero(cm_index)
97  INTEGER, INTENT(IN) :: cm_index ! Index of last dimension in
98  ! most com_mod variables.
99  ! Should be 1 or 2
100 
101 
102  ! Zeroes out, in local datastructures, the values dumped/loaded
103  ! This was written primarily as a testing mechanism.
104  ! The lines here correspond to READ and WRITE in the dump/load routines
105 
106  ! Scalar values
107  nx=0.0; ny=0.0; nxmin1=0.0; nymin1=0.0; nxfield=0.0
108  nuvz=0.0; nwz=0.0; nz=0.0; nmixz=0.0; nlev_ec=0.0
109  dx=0.0; dy=0.0; xlon0=0.0; ylat0=0.0; dxconst=0.0; dyconst=0.0
110 
111  ! Fixed fields, static in time
112  oro=0.0; excessoro=0.0; lsm=0.0; xlanduse=0.0; height=0.0
113 
114  ! 3d fields
115  uu(:,:,:,cm_index) = 0.0
116  vv(:,:,:,cm_index) = 0.0
117  uupol(:,:,:,cm_index) = 0.0
118  vvpol(:,:,:,cm_index) = 0.0
119  ww(:,:,:,cm_index) = 0.0
120  tt(:,:,:,cm_index) = 0.0
121  qv(:,:,:,cm_index) = 0.0
122  pv(:,:,:,cm_index) = 0.0
123  rho(:,:,:,cm_index) = 0.0
124  drhodz(:,:,:,cm_index) = 0.0
125  tth(:,:,:,cm_index) = 0.0
126  qvh(:,:,:,cm_index) = 0.0
127  pplev(:,:,:,cm_index) = 0.0
128  clouds(:,:,:,cm_index) = 0.0
129  cloudsh(:,:,cm_index) = 0.0
130 
131  ! 2d fields
132  ps(:,:,:,cm_index) = 0.0
133  sd(:,:,:,cm_index) = 0.0
134  msl(:,:,:,cm_index) = 0.0
135  tcc(:,:,:,cm_index) = 0.0
136  u10(:,:,:,cm_index) = 0.0
137  v10(:,:,:,cm_index) = 0.0
138  tt2(:,:,:,cm_index) = 0.0
139  td2(:,:,:,cm_index) = 0.0
140  lsprec(:,:,:,cm_index) = 0.0
141  convprec(:,:,:,cm_index) = 0.0
142  sshf(:,:,:,cm_index) = 0.0
143  ssr(:,:,:,cm_index) = 0.0
144  surfstr(:,:,:,cm_index) = 0.0
145  ustar(:,:,:,cm_index) = 0.0
146  wstar(:,:,:,cm_index) = 0.0
147  hmix(:,:,:,cm_index) = 0.0
148  tropopause(:,:,:,cm_index) = 0.0
149  oli(:,:,:,cm_index) = 0.0
150  diffk(:,:,:,cm_index) = 0.0
151  vdep(:,:,:,cm_index) = 0.0
152 
153  ! 1d fields
154  z0(:) = 0.0
155  akm(:) = 0.0
156  bkm(:) = 0.0
157  akz(:) = 0.0
158  bkz(:) = 0.0
159  aknew(:) = 0.0
160  bknew(:) = 0.0
161 
162  ! Nested, scalar values (for each nest)
163  nxn(:) = 0.0
164  nyn(:) = 0.0
165  dxn(:) = 0.0
166  dyn(:) = 0.0
167  xlon0n(:) = 0.0
168  ylat0n(:) = 0.0
169 
170  ! Nested fields, static in time
171  oron=0.0; excessoron=0.0; lsmn=0.0; xlandusen=0.0
172 
173  ! 3d nested fields
174  uun(:,:,:,cm_index,:) = 0.0
175  wwn(:,:,:,cm_index,:) = 0.0
176  ttn(:,:,:,cm_index,:) = 0.0
177  qvn(:,:,:,cm_index,:) = 0.0
178  pvn(:,:,:,cm_index,:) = 0.0
179  cloudsn(:,:,:,cm_index,:) = 0.0
180  cloudsnh(:,:,cm_index,:) = 0.0
181  rhon(:,:,:,cm_index,:) = 0.0
182  drhodzn(:,:,:,cm_index,:) = 0.0
183  tthn(:,:,:,cm_index,:) = 0.0
184  qvhn(:,:,:,cm_index,:) = 0.0
185 
186  ! 2d nested fields
187  psn(:,:,:,cm_index,:) = 0.0
188  sdn(:,:,:,cm_index,:) = 0.0
189  msln(:,:,:,cm_index,:) = 0.0
190  tccn(:,:,:,cm_index,:) = 0.0
191  u10n(:,:,:,cm_index,:) = 0.0
192  v10n(:,:,:,cm_index,:) = 0.0
193  tt2n(:,:,:,cm_index,:) = 0.0
194  td2n(:,:,:,cm_index,:) = 0.0
195  lsprecn(:,:,:,cm_index,:) = 0.0
196  convprecn(:,:,:,cm_index,:) = 0.0
197  sshfn(:,:,:,cm_index,:) = 0.0
198  ssrn(:,:,:,cm_index,:) = 0.0
199  surfstrn(:,:,:,cm_index,:) = 0.0
200  ustarn(:,:,:,cm_index,:) = 0.0
201  wstarn(:,:,:,cm_index,:) = 0.0
202  hmixn(:,:,:,cm_index,:) = 0.0
203  tropopausen(:,:,:,cm_index,:) = 0.0
204  olin(:,:,:,cm_index,:) = 0.0
205  diffkn(:,:,:,cm_index,:) = 0.0
206  vdepn(:,:,:,cm_index,:) = 0.0
207 
208  ! Auxiliary variables for nests
209  xresoln(:) = 0.0
210  yresoln(:) = 0.0
211  xln(:) = 0.0
212  yln(:) = 0.0
213  xrn(:) = 0.0
214  yrn(:) = 0.0
215 
216  ! Variables for polar stereographic projection
217  xglobal=.false.; sglobal=.false.; nglobal=.false.
218  switchnorthg=0.0; switchsouthg=0.0
219  southpolemap(:) = 0.0
220  northpolemap(:) = 0.0
221 
222  ! Variables declared in conv_mod (convection)
223  pconv(:) = 0.0
224  phconv(:) = 0.0
225  dpr(:) = 0.0
226  pconv_hpa(:) = 0.0
227  phconv_hpa(:) = 0.0
228  ft(:) = 0.0
229  fq(:) = 0.0
230  fmass(:,:) = 0.0
231  sub(:) = 0.0
232  fmassfrac(:,:) = 0.0
233  cbaseflux(:,:) = 0.0
234  cbasefluxn(:,:,:) = 0.0
235  tconv(:) = 0.0
236  qconv(:) = 0.0
237  qsconv(:) = 0.0
238  psconv=0.0; tt2conv=0.0; td2conv=0.0
239  nconvlev=0.0; nconvtop=0.0
240 
241  END SUBROUTINE fpmetbinary_zero
242 
243  SUBROUTINE fpio(iounit, op, cm_index)
244  IMPLICIT NONE
245  INTEGER, INTENT(IN) :: iounit
246  CHARACTER(LEN=4), INTENT(IN) :: op ! Operation - DUMP or LOAD
247  INTEGER, INTENT(IN) :: cm_index ! Index of last dimension in
248  ! most com_mod variables.
249  ! Should be 1 or 2
250  CHARACTER(LEN=128) :: errmesg
251 
252  if (op == 'DUMP') THEN
253 
254  ! Scalar values
255  WRITE(iounit) nx, ny, nxmin1, nymin1, nxfield
256  WRITE(iounit) nuvz, nwz, nz, nmixz, nlev_ec
257  WRITE(iounit) dx, dy, xlon0, ylat0, dxconst, dyconst
258 
259  ! Fixed fields, static in time
260  WRITE(iounit) oro, excessoro, lsm, xlanduse, height
261 
262  ! 3d fields
263  WRITE(iounit) uu(:,:,:,cm_index)
264  WRITE(iounit) vv(:,:,:,cm_index)
265  WRITE(iounit) uupol(:,:,:,cm_index)
266  WRITE(iounit) vvpol(:,:,:,cm_index)
267  WRITE(iounit) ww(:,:,:,cm_index)
268  WRITE(iounit) tt(:,:,:,cm_index)
269  WRITE(iounit) qv(:,:,:,cm_index)
270  WRITE(iounit) pv(:,:,:,cm_index)
271  WRITE(iounit) rho(:,:,:,cm_index)
272  WRITE(iounit) drhodz(:,:,:,cm_index)
273  WRITE(iounit) tth(:,:,:,cm_index)
274  WRITE(iounit) qvh(:,:,:,cm_index)
275  WRITE(iounit) pplev(:,:,:,cm_index)
276  WRITE(iounit) clouds(:,:,:,cm_index)
277  WRITE(iounit) cloudsh(:,:,cm_index)
278 
279  ! 2d fields
280  WRITE(iounit) ps(:,:,:,cm_index)
281  WRITE(iounit) sd(:,:,:,cm_index)
282  WRITE(iounit) msl(:,:,:,cm_index)
283  WRITE(iounit) tcc(:,:,:,cm_index)
284  WRITE(iounit) u10(:,:,:,cm_index)
285  WRITE(iounit) v10(:,:,:,cm_index)
286  WRITE(iounit) tt2(:,:,:,cm_index)
287  WRITE(iounit) td2(:,:,:,cm_index)
288  WRITE(iounit) lsprec(:,:,:,cm_index)
289  WRITE(iounit) convprec(:,:,:,cm_index)
290  WRITE(iounit) sshf(:,:,:,cm_index)
291  WRITE(iounit) ssr(:,:,:,cm_index)
292  WRITE(iounit) surfstr(:,:,:,cm_index)
293  WRITE(iounit) ustar(:,:,:,cm_index)
294  WRITE(iounit) wstar(:,:,:,cm_index)
295  WRITE(iounit) hmix(:,:,:,cm_index)
296  WRITE(iounit) tropopause(:,:,:,cm_index)
297  WRITE(iounit) oli(:,:,:,cm_index)
298  WRITE(iounit) diffk(:,:,:,cm_index)
299  WRITE(iounit) vdep(:,:,:,cm_index)
300 
301  ! 1d fields
302  WRITE(iounit) z0(:)
303  WRITE(iounit) akm(:)
304  WRITE(iounit) bkm(:)
305  WRITE(iounit) akz(:)
306  WRITE(iounit) bkz(:)
307  WRITE(iounit) aknew(:)
308  WRITE(iounit) bknew(:)
309 
310  ! Nested, scalar values (for each nest)
311  WRITE(iounit) nxn(:)
312  WRITE(iounit) nyn(:)
313  WRITE(iounit) dxn(:)
314  WRITE(iounit) dyn(:)
315  WRITE(iounit) xlon0n(:)
316  WRITE(iounit) ylat0n(:)
317 
318  ! Nested fields, static over time
319  WRITE(iounit) oron, excessoron, lsmn, xlandusen
320 
321  ! 3d nested fields
322  WRITE(iounit) uun(:,:,:,cm_index,:)
323  WRITE(iounit) vvn(:,:,:,cm_index,:)
324  WRITE(iounit) wwn(:,:,:,cm_index,:)
325  WRITE(iounit) ttn(:,:,:,cm_index,:)
326  WRITE(iounit) qvn(:,:,:,cm_index,:)
327  WRITE(iounit) pvn(:,:,:,cm_index,:)
328  WRITE(iounit) cloudsn(:,:,:,cm_index,:)
329  WRITE(iounit) cloudsnh(:,:,cm_index,:)
330  WRITE(iounit) rhon(:,:,:,cm_index,:)
331  WRITE(iounit) drhodzn(:,:,:,cm_index,:)
332  WRITE(iounit) tthn(:,:,:,cm_index,:)
333  WRITE(iounit) qvhn(:,:,:,cm_index,:)
334 
335  ! 2d nested fields
336  WRITE(iounit) psn(:,:,:,cm_index,:)
337  WRITE(iounit) sdn(:,:,:,cm_index,:)
338  WRITE(iounit) msln(:,:,:,cm_index,:)
339  WRITE(iounit) tccn(:,:,:,cm_index,:)
340  WRITE(iounit) u10n(:,:,:,cm_index,:)
341  WRITE(iounit) v10n(:,:,:,cm_index,:)
342  WRITE(iounit) tt2n(:,:,:,cm_index,:)
343  WRITE(iounit) td2n(:,:,:,cm_index,:)
344  WRITE(iounit) lsprecn(:,:,:,cm_index,:)
345  WRITE(iounit) convprecn(:,:,:,cm_index,:)
346  WRITE(iounit) sshfn(:,:,:,cm_index,:)
347  WRITE(iounit) ssrn(:,:,:,cm_index,:)
348  WRITE(iounit) surfstrn(:,:,:,cm_index,:)
349  WRITE(iounit) ustarn(:,:,:,cm_index,:)
350  WRITE(iounit) wstarn(:,:,:,cm_index,:)
351  WRITE(iounit) hmixn(:,:,:,cm_index,:)
352  WRITE(iounit) tropopausen(:,:,:,cm_index,:)
353  WRITE(iounit) olin(:,:,:,cm_index,:)
354  WRITE(iounit) diffkn(:,:,:,cm_index,:)
355  WRITE(iounit) vdepn(:,:,:,cm_index,:)
356 
357  ! Auxiliary variables for nests
358  WRITE(iounit) xresoln(:)
359  WRITE(iounit) yresoln(:)
360  WRITE(iounit) xln(:)
361  WRITE(iounit) yln(:)
362  WRITE(iounit) xrn(:)
363  WRITE(iounit) yrn(:)
364 
365  ! Variables for polar stereographic projection
366  WRITE(iounit) xglobal, sglobal, nglobal
367  WRITE(iounit) switchnorthg, switchsouthg
368  WRITE(iounit) southpolemap(:)
369  WRITE(iounit) northpolemap(:)
370 
371  ! Variables declared in conv_mod (convection)
372  WRITE(iounit) pconv(:)
373  WRITE(iounit) phconv(:)
374  WRITE(iounit) dpr(:)
375  WRITE(iounit) pconv_hpa(:)
376  WRITE(iounit) phconv_hpa(:)
377  WRITE(iounit) ft(:)
378  WRITE(iounit) fq(:)
379  WRITE(iounit) fmass(:,:)
380  WRITE(iounit) sub(:)
381  WRITE(iounit) fmassfrac(:,:)
382  WRITE(iounit) cbaseflux(:,:)
383  WRITE(iounit) cbasefluxn(:,:,:)
384  WRITE(iounit) tconv(:)
385  WRITE(iounit) qconv(:)
386  WRITE(iounit) qsconv(:)
387  WRITE(iounit) psconv, tt2conv, td2conv
388  WRITE(iounit) nconvlev, nconvtop
389 
390  ELSE IF (op == 'LOAD') THEN
391 
392  ! Scalar values
393  READ(iounit) nx, ny, nxmin1, nymin1, nxfield
394  READ(iounit) nuvz, nwz, nz, nmixz, nlev_ec
395  READ(iounit) dx, dy, xlon0, ylat0, dxconst, dyconst
396 
397  ! Fixed fields, static in time
398  READ(iounit) oro, excessoro, lsm, xlanduse, height
399 
400  ! 3d fields
401  READ(iounit) uu(:,:,:,cm_index)
402  READ(iounit) vv(:,:,:,cm_index)
403  READ(iounit) uupol(:,:,:,cm_index)
404  READ(iounit) vvpol(:,:,:,cm_index)
405  READ(iounit) ww(:,:,:,cm_index)
406  READ(iounit) tt(:,:,:,cm_index)
407  READ(iounit) qv(:,:,:,cm_index)
408  READ(iounit) pv(:,:,:,cm_index)
409  READ(iounit) rho(:,:,:,cm_index)
410  READ(iounit) drhodz(:,:,:,cm_index)
411  READ(iounit) tth(:,:,:,cm_index)
412  READ(iounit) qvh(:,:,:,cm_index)
413  READ(iounit) pplev(:,:,:,cm_index)
414  READ(iounit) clouds(:,:,:,cm_index)
415  READ(iounit) cloudsh(:,:,cm_index)
416 
417  ! 2d fields
418  READ(iounit) ps(:,:,:,cm_index)
419  READ(iounit) sd(:,:,:,cm_index)
420  READ(iounit) msl(:,:,:,cm_index)
421  READ(iounit) tcc(:,:,:,cm_index)
422  READ(iounit) u10(:,:,:,cm_index)
423  READ(iounit) v10(:,:,:,cm_index)
424  READ(iounit) tt2(:,:,:,cm_index)
425  READ(iounit) td2(:,:,:,cm_index)
426  READ(iounit) lsprec(:,:,:,cm_index)
427  READ(iounit) convprec(:,:,:,cm_index)
428  READ(iounit) sshf(:,:,:,cm_index)
429  READ(iounit) ssr(:,:,:,cm_index)
430  READ(iounit) surfstr(:,:,:,cm_index)
431  READ(iounit) ustar(:,:,:,cm_index)
432  READ(iounit) wstar(:,:,:,cm_index)
433  READ(iounit) hmix(:,:,:,cm_index)
434  READ(iounit) tropopause(:,:,:,cm_index)
435  READ(iounit) oli(:,:,:,cm_index)
436  READ(iounit) diffk(:,:,:,cm_index)
437  READ(iounit) vdep(:,:,:,cm_index)
438 
439  ! 1d fields
440  READ(iounit) z0(:)
441  READ(iounit) akm(:)
442  READ(iounit) bkm(:)
443  READ(iounit) akz(:)
444  READ(iounit) bkz(:)
445  READ(iounit) aknew(:)
446  READ(iounit) bknew(:)
447 
448 
449  ! Nested, scalar values (for each nest)
450  READ(iounit) nxn(:)
451  READ(iounit) nyn(:)
452  READ(iounit) dxn(:)
453  READ(iounit) dyn(:)
454  READ(iounit) xlon0n(:)
455  READ(iounit) ylat0n(:)
456 
457 
458  ! Nested fields, static over time
459  READ(iounit) oron, excessoron, lsmn, xlandusen
460 
461  ! 3d nested fields
462  READ(iounit) uun(:,:,:,cm_index,:)
463  READ(iounit) vvn(:,:,:,cm_index,:)
464  READ(iounit) wwn(:,:,:,cm_index,:)
465  READ(iounit) ttn(:,:,:,cm_index,:)
466  READ(iounit) qvn(:,:,:,cm_index,:)
467  READ(iounit) pvn(:,:,:,cm_index,:)
468  READ(iounit) cloudsn(:,:,:,cm_index,:)
469  READ(iounit) cloudsnh(:,:,cm_index,:)
470  READ(iounit) rhon(:,:,:,cm_index,:)
471  READ(iounit) drhodzn(:,:,:,cm_index,:)
472  READ(iounit) tthn(:,:,:,cm_index,:)
473  READ(iounit) qvhn(:,:,:,cm_index,:)
474 
475  ! 2d nested fields
476  READ(iounit) psn(:,:,:,cm_index,:)
477  READ(iounit) sdn(:,:,:,cm_index,:)
478  READ(iounit) msln(:,:,:,cm_index,:)
479  READ(iounit) tccn(:,:,:,cm_index,:)
480  READ(iounit) u10n(:,:,:,cm_index,:)
481  READ(iounit) v10n(:,:,:,cm_index,:)
482  READ(iounit) tt2n(:,:,:,cm_index,:)
483  READ(iounit) td2n(:,:,:,cm_index,:)
484  READ(iounit) lsprecn(:,:,:,cm_index,:)
485  READ(iounit) convprecn(:,:,:,cm_index,:)
486  READ(iounit) sshfn(:,:,:,cm_index,:)
487  READ(iounit) ssrn(:,:,:,cm_index,:)
488  READ(iounit) surfstrn(:,:,:,cm_index,:)
489  READ(iounit) ustarn(:,:,:,cm_index,:)
490  READ(iounit) wstarn(:,:,:,cm_index,:)
491  READ(iounit) hmixn(:,:,:,cm_index,:)
492  READ(iounit) tropopausen(:,:,:,cm_index,:)
493  READ(iounit) olin(:,:,:,cm_index,:)
494  READ(iounit) diffkn(:,:,:,cm_index,:)
495  READ(iounit) vdepn(:,:,:,cm_index,:)
496 
497  ! Auxiliary variables for nests
498  READ(iounit) xresoln(:)
499  READ(iounit) yresoln(:)
500  READ(iounit) xln(:)
501  READ(iounit) yln(:)
502  READ(iounit) xrn(:)
503  READ(iounit) yrn(:)
504 
505  ! Variables for polar stereographic projection
506  READ(iounit) xglobal, sglobal, nglobal
507  READ(iounit) switchnorthg, switchsouthg
508  READ(iounit) southpolemap(:)
509  READ(iounit) northpolemap(:)
510 
511  ! Variables declared in conv_mod (convection)
512  READ(iounit) pconv(:)
513  READ(iounit) phconv(:)
514  READ(iounit) dpr(:)
515  READ(iounit) pconv_hpa(:)
516  READ(iounit) phconv_hpa(:)
517  READ(iounit) ft(:)
518  READ(iounit) fq(:)
519  READ(iounit) fmass(:,:)
520  READ(iounit) sub(:)
521  READ(iounit) fmassfrac(:,:)
522  READ(iounit) cbaseflux(:,:)
523  READ(iounit) cbasefluxn(:,:,:)
524  READ(iounit) tconv(:)
525  READ(iounit) qconv(:)
526  READ(iounit) qsconv(:)
527  READ(iounit) psconv, tt2conv, td2conv
528  READ(iounit) nconvlev, nconvtop
529 
530  ELSE
531  stop 'fpio(): Illegal operation'
532 
533  ENDIF
534  END SUBROUTINE fpio
535 
536  SUBROUTINE fpmetbinary_filetext(filename, cm_index)
537 
538  ! This is a utility subroutine meant to be used for testing purposes.
539  ! It facilitates the text output of variables read in from the
540  ! specified .fp file. This routine will easily cause the program
541  ! to crash due memory allocation issues, particularly when you are
542  ! trying to text print 3d arrays in a single formatted statetment.
543 
544  CHARACTER(LEN=*), INTENT(IN) :: filename
545  INTEGER, INTENT(IN) :: cm_index ! Index of last dimension in
546  ! most com_mod variables.
547  ! Should be 1 or 2
548 
549  !OPEN(IOUNIT_TEXTOUT, file=filename, action='WRITE', status='REPLACE', &
550  ! form="formatted", access="stream")
551  OPEN(iounit_textout, file=filename, action='WRITE', &
552  form="formatted", access="APPEND")
553 
554  WRITE(iounit_textout, *) 'oro: ', oro
555  WRITE(iounit_textout, *) 'excessoro: ', excessoro
556  WRITE(iounit_textout, *) 'lsm: ', lsm
557  WRITE(iounit_textout, *) 'xlanduse: ', xlanduse
558  WRITE(iounit_textout, *) 'height: ', height
559 
560  WRITE(iounit_textout, *) 'uu: ', uu(:,:,:,cm_index)
561  WRITE(iounit_textout, *) 'vv: ', vv(:,:,:,cm_index)
562  WRITE(iounit_textout, *) 'uupol: ', uupol(:,:,:,cm_index)
563  WRITE(iounit_textout, *) 'vvpol: ', vvpol(:,:,:,cm_index)
564  WRITE(iounit_textout, *) 'ww: ', ww(:,:,:,cm_index)
565  WRITE(iounit_textout, *) 'tt: ', tt(:,:,:,cm_index)
566  WRITE(iounit_textout, *) 'qv: ', qv(:,:,:,cm_index)
567  WRITE(iounit_textout, *) 'pv: ', pv(:,:,:,cm_index)
568  WRITE(iounit_textout, *) 'rho: ', rho(:,:,:,cm_index)
569  WRITE(iounit_textout, *) 'drhodz: ', drhodz(:,:,:,cm_index)
570  WRITE(iounit_textout, *) 'tth: ', tth(:,:,:,cm_index)
571  WRITE(iounit_textout, *) 'qvh: ', qvh(:,:,:,cm_index)
572  WRITE(iounit_textout, *) 'pplev: ', pplev(:,:,:,cm_index)
573  WRITE(iounit_textout, *) 'clouds: ', clouds(:,:,:,cm_index)
574  WRITE(iounit_textout, *) 'cloudsh: ', cloudsh(:,:,cm_index)
575 
576 
577 
578 
579  CLOSE(iounit_textout)
580  END SUBROUTINE fpmetbinary_filetext
581 
582 
583 END MODULE fpmetbinary_mod
fpmetbinary_mod::fpmetbinary_zero
subroutine fpmetbinary_zero(cm_index)
Definition: fpmetbinary_mod.F90:96
fpmetbinary_mod::fpmetbinary_filetext
subroutine fpmetbinary_filetext(filename, cm_index)
Definition: fpmetbinary_mod.F90:536
fpmetbinary_mod::fpmetbinary_dump
subroutine fpmetbinary_dump(filename, cm_index)
Definition: fpmetbinary_mod.F90:63
fpmetbinary_mod::fpio
subroutine, private fpio(iounit, op, cm_index)
Definition: fpmetbinary_mod.F90:243
fpmetbinary_mod::fpmetbinary_load
subroutine fpmetbinary_load(filename, cm_index)
Definition: fpmetbinary_mod.F90:80
conv_mod
Definition: conv_mod.f90:11
com_mod
Definition: com_mod.f90:13
fpmetbinary_mod
Definition: fpmetbinary_mod.F90:1
Generated on Thu Oct 15 2015 06:20:58 for CTBTO FLEXPART WO4 (2015-10-15) by   doxygen 1.8.6