ROMS
Loading...
Searching...
No Matches
def_norm.F
Go to the documentation of this file.
1#include "cppdefs.h"
2 MODULE def_norm_mod
3#ifdef FOUR_DVAR
4!
5!git $Id$
6!================================================== Hernan G. Arango ===
7! Copyright (c) 2002-2025 The ROMS Group !
8! Licensed under a MIT/X style license !
9! See License_ROMS.md !
10!=======================================================================
11! !
12! This routine creates error covariance normalization factors file !
13! using either the standard NetCDF library or the Parallel-IO (PIO) !
14! library. It defines its dimensions, attributes, and variables. !
15! !
16! The NetCDF files are used for variational data assimilation. Four !
17! different files can be generates for initial, model, opem boundary, !
18! and surface forcing error covariances. !
19! !
20!=======================================================================
21!
22 USE mod_param
23 USE mod_parallel
24# ifdef BIOLOGY
25 USE mod_biology
26# endif
27 USE mod_fourdvar
28 USE mod_iounits
29 USE mod_ncparam
30 USE mod_scalars
31# if defined SEDIMENT || defined BBL_MODEL
32 USE mod_sediment
33# endif
34!
35 USE def_dim_mod, ONLY : def_dim
36 USE def_info_mod, ONLY : def_info
37 USE def_var_mod, ONLY : def_var
38 USE strings_mod, ONLY : founderror
39 USE wrt_info_mod, ONLY : wrt_info
40!
41 implicit none
42!
43 PUBLIC :: def_norm
44 PRIVATE :: def_norm_nf90
45# if defined PIO_LIB && defined DISTRIBUTE
46 PRIVATE :: def_norm_pio
47# endif
48!
49 CONTAINS
50!
51!***********************************************************************
52 SUBROUTINE def_norm (ng, model, ifile)
53!***********************************************************************
54!
55! Imported variable declarations.
56!
57 integer, intent(in) :: ng, model, ifile
58!
59! Local variable declarations.
60!
61 character (len=*), parameter :: myfile = &
62 & __FILE__
63!
64!-----------------------------------------------------------------------
65! Create a new history file according to IO type.
66!-----------------------------------------------------------------------
67!
68 SELECT CASE (nrm(ifile,ng)%IOtype)
69 CASE (io_nf90)
70 CALL def_norm_nf90 (ng, model, ifile)
71
72# if defined PIO_LIB && defined DISTRIBUTE
73 CASE (io_pio)
74 CALL def_norm_pio (ng, model, ifile)
75# endif
76 CASE DEFAULT
77 IF (master) WRITE (stdout,10) nrm(ifile,ng)%IOtype
78 exit_flag=3
79 END SELECT
80 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
81!
82 10 FORMAT (' DEF_NORM - Illegal output file type, io_type = ',i0, &
83 & /,12x,'Check KeyWord ''OUT_LIB'' in ''roms.in''.')
84!
85 RETURN
86 END SUBROUTINE def_norm
87!
88!***********************************************************************
89 SUBROUTINE def_norm_nf90 (ng, model, ifile)
90!***********************************************************************
91!
92 USE mod_netcdf
93!
94! Imported variable declarations.
95!
96 integer, intent(in) :: ng, model, ifile
97!
98! Local variable declarations.
99!
100 logical :: ldefine, got_var(nv)
101!
102 integer, parameter :: natt = 25
103
104 integer :: i, j, nvd3, nvd4
105 integer :: recdim, status, varid
106# ifdef ADJUST_BOUNDARY
107 integer :: iorjdim, ifield
108# endif
109 integer :: dimids(ndimid)
110 integer :: t2dgrd(3), u2dgrd(3), v2dgrd(3)
111# ifdef ADJUST_BOUNDARY
112 integer :: t2dobc(3)
113# endif
114# ifdef SOLVE3D
115 integer :: itrc
116
117 integer :: t3dgrd(4), u3dgrd(4), v3dgrd(4)
118# ifdef ADJUST_BOUNDARY
119 integer :: t3dobc(4)
120# endif
121# endif
122!
123 real(r8) :: aval(6)
124!
125 character (len=60 ) :: text
126 character (len=256) :: ncname
127 character (len=MaxLen) :: vinfo(natt)
128
129 character (len=*), parameter :: myfile = &
130 & __FILE__//", def_norm_nf90"
131!
132 sourcefile=myfile
133!
134!=======================================================================
135! Create a new background covariace normalization file.
136!=======================================================================
137!
138 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
139 ncname=nrm(ifile,ng)%name
140!
141 define : IF (ldefnrm(ifile,ng)) THEN
142 CALL netcdf_create (ng, itlm, trim(ncname), nrm(ifile,ng)%ncid)
143 IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
144 IF (master) WRITE (stdout,10) trim(ncname)
145 RETURN
146 END IF
147!
148!-----------------------------------------------------------------------
149! Define the dimensions of staggered fields.
150!-----------------------------------------------------------------------
151!
152 dimids=0
153!
154 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xi_rho', &
155 & iobounds(ng)%xi_rho, dimids( 1))
156 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
157
158 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xi_u', &
159 & iobounds(ng)%xi_u, dimids( 2))
160 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
161
162 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xi_v', &
163 & iobounds(ng)%xi_v, dimids( 3))
164 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
165
166 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xi_psi', &
167 & iobounds(ng)%xi_psi, dimids( 4))
168 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
169
170 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'eta_rho', &
171 & iobounds(ng)%eta_rho, dimids( 5))
172 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
173
174 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'eta_u', &
175 & iobounds(ng)%eta_u, dimids( 6))
176 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
177
178 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'eta_v', &
179 & iobounds(ng)%eta_v, dimids( 7))
180 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
181
182 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'eta_psi', &
183 & iobounds(ng)%eta_psi, dimids( 8))
184 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
185
186# ifdef ADJUST_BOUNDARY
187 IF (ifile.eq.3) THEN
188 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'IorJ', &
189 & iobounds(ng)%IorJ, iorjdim)
190 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
191 END IF
192# endif
193
194# if defined WRITE_WATER && defined MASKING
195 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xy_rho', &
196 & iobounds(ng)%xy_rho, dimids(17))
197 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
198
199 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xy_u', &
200 & iobounds(ng)%xy_u, dimids(18))
201 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
202
203 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xy_v', &
204 & iobounds(ng)%xy_v, dimids(19))
205 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
206# endif
207
208# ifdef SOLVE3D
209# if defined WRITE_WATER && defined MASKING
210 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xyz_rho', &
211 & iobounds(ng)%xy_rho*n(ng), dimids(20))
212 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
213
214 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xyz_u', &
215 & iobounds(ng)%xy_u*n(ng), dimids(21))
216 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
217
218 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xyz_v', &
219 & iobounds(ng)%xy_v*n(ng), dimids(22))
220 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
221
222 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xyz_w', &
223 & iobounds(ng)%xy_rho*(n(ng)+1), dimids(23))
224 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
225# endif
226
227 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 's_rho', &
228 & n(ng), dimids( 9))
229 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
230
231 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 's_w', &
232 & n(ng)+1, dimids(10))
233 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
234
235 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'tracer', &
236 & nt(ng), dimids(11))
237 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
238
239# ifdef SEDIMENT
240 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'NST', &
241 & nst, dimids(32))
242 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
243
244 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'Nbed', &
245 & nbed, dimids(16))
246 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
247
248# if defined WRITE_WATER && defined MASKING
249 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'xybed', &
250 & iobounds(ng)%xy_rho*nbed, dimids(24))
251 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
252# endif
253# endif
254
255# ifdef ECOSIM
256 status=def_dim(ng, inlm, nrm(ifile,ng)%ncid, ncname, 'Nbands', &
257 & nbands, dimids(33))
258 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
259
260 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'Nphy', &
261 & nphy, dimids(25))
262 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
263
264 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'Nbac', &
265 & nbac, dimids(26))
266 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
267
268 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'Ndom', &
269 & ndom, dimids(27))
270 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
271
272 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'Nfec', &
273 & nfec, dimids(28))
274 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
275# endif
276# endif
277
278 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'boundary',&
279 & 4, dimids(14))
280 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
281
282# ifdef FOUR_DVAR
283 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, 'Nstate', &
284 & nstatevar(ng), dimids(29))
285 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
286# endif
287
288 status=def_dim(ng, itlm, nrm(ifile,ng)%ncid, ncname, &
289 & trim(adjustl(vname(5,idtime))), &
290 & nf90_unlimited, dimids(12))
291 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
292
293 recdim=dimids(12)
294!
295! Set number of dimensions for output variables.
296!
297# if defined WRITE_WATER && defined MASKING
298 nvd3=2
299 nvd4=2
300# else
301 nvd3=3
302 nvd4=4
303# endif
304!
305! Define dimension vectors for staggered tracer type variables.
306!
307# if defined WRITE_WATER && defined MASKING
308 t2dgrd(1)=dimids(17)
309 t2dgrd(2)=dimids(12)
310# ifdef SOLVE3D
311 t3dgrd(1)=dimids(20)
312 t3dgrd(2)=dimids(12)
313# endif
314# else
315 t2dgrd(1)=dimids( 1)
316 t2dgrd(2)=dimids( 5)
317 t2dgrd(3)=dimids(12)
318# ifdef SOLVE3D
319 t3dgrd(1)=dimids( 1)
320 t3dgrd(2)=dimids( 5)
321 t3dgrd(3)=dimids( 9)
322 t3dgrd(4)=dimids(12)
323# endif
324# endif
325# ifdef ADJUST_BOUNDARY
326 IF (ifile.eq.3) THEN
327 t2dobc(1)=iorjdim
328 t2dobc(2)=dimids(14)
329 t2dobc(3)=dimids(12)
330# ifdef SOLVE3D
331 t3dobc(1)=iorjdim
332 t3dobc(2)=dimids( 9)
333 t3dobc(3)=dimids(14)
334 t3dobc(4)=dimids(12)
335# endif
336 END IF
337# endif
338!
339! Define dimension vectors for staggered u-momentum type variables.
340!
341# if defined WRITE_WATER && defined MASKING
342 u2dgrd(1)=dimids(18)
343 u2dgrd(2)=dimids(12)
344# ifdef SOLVE3D
345 u3dgrd(1)=dimids(21)
346 u3dgrd(2)=dimids(12)
347# endif
348# else
349 u2dgrd(1)=dimids( 2)
350 u2dgrd(2)=dimids( 6)
351 u2dgrd(3)=dimids(12)
352# ifdef SOLVE3D
353 u3dgrd(1)=dimids( 2)
354 u3dgrd(2)=dimids( 6)
355 u3dgrd(3)=dimids( 9)
356 u3dgrd(4)=dimids(12)
357# endif
358# endif
359!
360! Define dimension vectors for staggered v-momentum type variables.
361!
362# if defined WRITE_WATER && defined MASKING
363 v2dgrd(1)=dimids(19)
364 v2dgrd(2)=dimids(12)
365# ifdef SOLVE3D
366 v3dgrd(1)=dimids(22)
367 v3dgrd(2)=dimids(12)
368# endif
369# else
370 v2dgrd(1)=dimids( 3)
371 v2dgrd(2)=dimids( 7)
372 v2dgrd(3)=dimids(12)
373# ifdef SOLVE3D
374 v3dgrd(1)=dimids( 3)
375 v3dgrd(2)=dimids( 7)
376 v3dgrd(3)=dimids( 9)
377 v3dgrd(4)=dimids(12)
378# endif
379# endif
380!
381! Initialize unlimited time record dimension.
382!
383 nrm(ifile,ng)%Rindex=0
384!
385! Initialize local information variable arrays.
386!
387 DO i=1,natt
388 DO j=1,len(vinfo(1))
389 vinfo(i)(j:j)=' '
390 END DO
391 END DO
392 DO i=1,6
393 aval(i)=0.0_r8
394 END DO
395!
396!-----------------------------------------------------------------------
397! Define time-recordless information variables.
398!-----------------------------------------------------------------------
399!
400 CALL def_info (ng, itlm, nrm(ifile,ng)%ncid, ncname, dimids)
401 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
402!
403!-----------------------------------------------------------------------
404! Define background covariance normalization variables.
405!-----------------------------------------------------------------------
406!
407! Define model time.
408!
409 vinfo( 1)=vname(1,idtime)
410 vinfo( 2)=vname(2,idtime)
411 WRITE (vinfo( 3),'(a,a)') 'seconds since ', trim(rclock%string)
412 vinfo( 4)=trim(rclock%calendar)
413 vinfo(14)=vname(4,idtime)
414 vinfo(21)=vname(6,idtime)
415 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
416 & nrm(ifile,ng)%Vid(idtime), &
417 & nf_tout, 1, (/recdim/), aval, vinfo, ncname, &
418 & setparaccess = .true.)
419 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
420!
421!-----------------------------------------------------------------------
422! Initial conditions or model error covariance normalization factors.
423!-----------------------------------------------------------------------
424!
425 IF ((ifile.eq.1).or.(ifile.eq.2)) THEN
426 IF (ifile.eq.1) THEN
427 text='initial conditions error covariance normalization'
428 ELSE IF (ifile.eq.2) THEN
429 text='model error covariance normalization'
430 END IF
431!
432! Define free-surface normalization factor.
433!
434 vinfo( 1)=vname(1,idfsur)
435 WRITE (vinfo( 2),20) trim(vname(2,idfsur)), trim(text)
436 vinfo( 3)='nondimensional'
437 vinfo(14)=vname(4,idfsur)
438 vinfo(16)=vname(1,idtime)
439# if defined WRITE_WATER && defined MASKING
440 vinfo(20)='mask_rho'
441# endif
442 vinfo(21)=vname(6,idfsur)
443 vinfo(22)='coordinates'
444 aval(5)=real(iinfo(1,idfsur,ng),r8)
445 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
446 & nrm(ifile,ng)%Vid(idfsur), &
447 & nf_fout, nvd3, t2dgrd, aval, vinfo, ncname)
448 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
449!
450! Define 2D U-momentum normalization factor.
451!
452 vinfo( 1)=vname(1,idubar)
453 WRITE (vinfo( 2),20) trim(vname(2,idubar)), trim(text)
454 vinfo( 3)='nondimensional'
455 vinfo(14)=vname(4,idubar)
456 vinfo(16)=vname(1,idtime)
457# if defined WRITE_WATER && defined MASKING
458 vinfo(20)='mask_u'
459# endif
460 vinfo(21)=vname(6,idubar)
461 vinfo(22)='coordinates'
462 aval(5)=real(iinfo(1,idubar,ng),r8)
463 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
464 & nrm(ifile,ng)%Vid(idubar), &
465 & nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
466 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
467!
468! Define 2D V-momentum normalization factor.
469!
470 vinfo( 1)=vname(1,idvbar)
471 WRITE (vinfo( 2),20) trim(vname(2,idvbar)), trim(text)
472 vinfo( 3)='nondimensional'
473 vinfo(14)=vname(4,idvbar)
474 vinfo(16)=vname(1,idtime)
475# if defined WRITE_WATER && defined MASKING
476 vinfo(20)='mask_v'
477# endif
478 vinfo(21)=vname(6,idvbar)
479 vinfo(22)='coordinates'
480 aval(5)=real(iinfo(1,idvbar,ng),r8)
481 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
482 & nrm(ifile,ng)%Vid(idvbar), &
483 & nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
484 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
485
486# ifdef SOLVE3D
487!
488! Define 3D U-momentum normalization factor.
489!
490 vinfo( 1)=vname(1,iduvel)
491 WRITE (vinfo( 2),20) trim(vname(2,iduvel)), trim(text)
492 vinfo( 3)='nondimensional'
493 vinfo(14)=vname(4,iduvel)
494 vinfo(16)=vname(1,idtime)
495# if defined WRITE_WATER && defined MASKING
496 vinfo(20)='mask_u'
497# endif
498 vinfo(21)=vname(6,iduvel)
499 vinfo(22)='coordinates'
500 aval(5)=real(iinfo(1,iduvel,ng),r8)
501 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
502 & nrm(ifile,ng)%Vid(iduvel), &
503 & nf_fout, nvd4, u3dgrd, aval, vinfo, ncname)
504 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
505!
506! Define 3D V-momentum normalization factor.
507!
508 vinfo( 1)=vname(1,idvvel)
509 WRITE (vinfo( 2),20) trim(vname(2,idvvel)), trim(text)
510 vinfo( 3)='nondimensional'
511 vinfo(14)=vname(4,idvvel)
512 vinfo(16)=vname(1,idtime)
513# if defined WRITE_WATER && defined MASKING
514 vinfo(20)='mask_v'
515# endif
516 vinfo(21)=vname(6,idvvel)
517 vinfo(22)='coordinates'
518 aval(5)=real(iinfo(1,idvvel,ng),r8)
519 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
520 & nrm(ifile,ng)%Vid(idvvel), &
521 & nf_fout, nvd4, v3dgrd, aval, vinfo, ncname)
522 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
523!
524! Define tracer type normalization factors.
525!
526 DO itrc=1,nt(ng)
527 vinfo( 1)=vname(1,idtvar(itrc))
528 WRITE (vinfo( 2),20) trim(vname(2,idtvar(itrc))), trim(text)
529 vinfo( 3)='nondimensional'
530 vinfo(14)=vname(4,idtvar(itrc))
531 vinfo(16)=vname(1,idtime)
532# ifdef SEDIMENT
533 DO i=1,nst
534 IF (itrc.eq.idsed(i)) THEN
535 WRITE (vinfo(19),30) 1000.0_r8*sd50(i,ng)
536 END IF
537 END DO
538# endif
539# if defined WRITE_WATER && defined MASKING
540 vinfo(20)='mask_rho'
541# endif
542 vinfo(21)=vname(6,idtvar(itrc))
543 vinfo(22)='coordinates'
544 aval(5)=real(r3dvar,r8)
545 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
546 & nrm(ifile,ng)%Vid(idtvar(itrc)), &
547 & nf_fout, nvd4, t3dgrd, aval, vinfo, ncname)
548 END DO
549 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
550# endif
551
552# ifdef ADJUST_BOUNDARY
553!
554!-----------------------------------------------------------------------
555! Boundary conditions error covariance normalization factors.
556!-----------------------------------------------------------------------
557!
558 ELSE IF (ifile.eq.3) THEN
559 text='error covariance normalization'
560!
561! Define free-surface open boundaries.
562!
563 IF (any(lobc(:,isfsur,ng))) THEN
564 ifield=idsbry(isfsur)
565 vinfo( 1)=vname(1,ifield)
566 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
567 vinfo( 3)=vname(3,ifield)
568 vinfo(14)=vname(4,ifield)
569 vinfo(16)=vname(1,idtime)
570 vinfo(21)=vname(6,ifield)
571 aval(5)=real(iinfo(1,ifield,ng),r8)
572 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
573 & nrm(ifile,ng)%Vid(ifield), &
574 & nf_fout, 3, t2dobc, aval, vinfo, ncname, &
575 & setfillval = .false.)
576 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
577 END IF
578!
579! Define 2D U-momentum component open boundaries.
580!
581 IF (any(lobc(:,isubar,ng))) THEN
582 ifield=idsbry(isubar)
583 vinfo( 1)=vname(1,ifield)
584 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
585 vinfo( 3)=vname(3,ifield)
586 vinfo(14)=vname(4,ifield)
587 vinfo(16)=vname(1,idtime)
588 vinfo(21)=vname(6,ifield)
589 aval(5)=real(iinfo(1,ifield,ng),r8)
590 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
591 & nrm(ifile,ng)%Vid(ifield), &
592 & nf_fout, 3, t2dobc, aval, vinfo, ncname, &
593 & setfillval = .false.)
594 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
595 END IF
596!
597! Define 2D V-momentum component open boundaries.
598!
599 IF (any(lobc(:,isvbar,ng))) THEN
600 ifield=idsbry(isvbar)
601 vinfo( 1)=vname(1,ifield)
602 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
603 vinfo( 3)=vname(3,ifield)
604 vinfo(14)=vname(4,ifield)
605 vinfo(16)=vname(1,idtime)
606 vinfo(21)=vname(6,ifield)
607 aval(5)=real(iinfo(1,ifield,ng),r8)
608 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
609 & nrm(ifile,ng)%Vid(ifield), &
610 & nf_fout, 3, t2dobc, aval, vinfo, ncname, &
611 & setfillval = .false.)
612 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
613 END IF
614
615# ifdef SOLVE3D
616!
617! Define 3D U-momentum component open boundaries.
618!
619 IF (any(lobc(:,isuvel,ng))) THEN
620 ifield=idsbry(isuvel)
621 vinfo( 1)=vname(1,ifield)
622 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
623 vinfo( 3)=vname(3,ifield)
624 vinfo(14)=vname(4,ifield)
625 vinfo(16)=vname(1,idtime)
626 vinfo(21)=vname(6,ifield)
627 aval(5)=real(iinfo(1,ifield,ng),r8)
628 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
629 & nrm(ifile,ng)%Vid(ifield), &
630 & nf_fout, 4, t3dobc, aval, vinfo, ncname, &
631 & setfillval = .false.)
632 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
633 END IF
634!
635! Define 3D V-momentum component open boundaries.
636!
637 IF (any(lobc(:,isvvel,ng))) THEN
638 ifield=idsbry(isvvel)
639 vinfo( 1)=vname(1,ifield)
640 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
641 vinfo( 3)=vname(3,ifield)
642 vinfo(14)=vname(4,ifield)
643 vinfo(16)=vname(1,idtime)
644 vinfo(21)=vname(6,ifield)
645 aval(5)=real(iinfo(1,ifield,ng),r8)
646 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
647 & nrm(ifile,ng)%Vid(ifield), &
648 & nf_fout, 4, t3dobc, aval, vinfo, ncname, &
649 & setfillval = .false.)
650 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
651 END IF
652!
653! Define tracer type variables open boundaries.
654!
655 DO itrc=1,nt(ng)
656 IF (any(lobc(:,istvar(itrc),ng))) THEN
657 ifield=idsbry(istvar(itrc))
658 vinfo( 1)=vname(1,ifield)
659 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
660 vinfo( 3)=vname(3,ifield)
661 vinfo(14)=vname(4,ifield)
662 vinfo(16)=vname(1,idtime)
663# ifdef SEDIMENT
664 DO i=1,nst
665 IF (itrc.eq.idsed(i)) THEN
666 WRITE (vinfo(19),30) 1000.0_r8*sd50(i,ng)
667 END IF
668 END DO
669# endif
670 vinfo(21)=vname(6,ifield)
671 aval(5)=real(iinfo(1,ifield,ng),r8)
672 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
673 & nrm(ifile,ng)%Vid(ifield), &
674 & nf_fout, 4, t3dobc, aval, vinfo, ncname, &
675 & setfillval = .false.)
676 IF (founderror(exit_flag, noerror, &
677 & __line__, myfile)) RETURN
678 END IF
679 END DO
680# endif
681# endif
682
683# if defined ADJUST_WSTRESS || defined ADJUST_STFLUX
684!
685!-----------------------------------------------------------------------
686! Surface forcing error covariance normalization factors.
687!-----------------------------------------------------------------------
688!
689 ELSE IF (ifile.eq.4) THEN
690 text='error covariance normalization'
691
692# ifdef ADJUST_WSTRESS
693!
694! Define surface U-momentum stress normalization factors.
695!
696 vinfo( 1)=vname(1,idusms)
697 WRITE (vinfo( 2),20) trim(vname(2,idusms)), trim(text)
698 vinfo( 3)='nondimensional'
699 vinfo(14)=vname(4,idusms)
700 vinfo(16)=vname(1,idtime)
701# if defined WRITE_WATER && defined MASKING
702 vinfo(20)='mask_u'
703# endif
704 vinfo(21)=vname(6,idusms)
705 vinfo(22)='coordinates'
706 aval(5)=real(u2dvar,r8)
707 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
708 & nrm(ifile,ng)%Vid(idusms), &
709 & nf_fout, nvd3, u2dgrd, aval, vinfo, ncname)
710 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
711!
712! Define surface V-momentum stress normalization factors.
713!
714 vinfo( 1)=vname(1,idvsms)
715 WRITE (vinfo( 2),20) trim(vname(2,idvsms)), trim(text)
716 vinfo( 2)=vname(2,idvsms)
717 vinfo( 3)='nondimensional'
718 vinfo(14)=vname(4,idvsms)
719 vinfo(16)=vname(1,idtime)
720# if defined WRITE_WATER && defined MASKING
721 vinfo(20)='mask_v'
722# endif
723 vinfo(21)=vname(6,idvsms)
724 vinfo(22)='coordinates'
725 aval(5)=real(v2dvar,r8)
726 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
727 & nrm(ifile,ng)%Vid(idvsms), &
728 & nf_fout, nvd3, v2dgrd, aval, vinfo, ncname)
729 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
730# endif
731
732# if defined ADJUST_STFLUX && defined SOLVE3D
733!
734! Define surface tracer fluxes.
735!
736 DO itrc=1,nt(ng)
737 IF (lstflux(itrc,ng)) THEN
738 vinfo( 1)=vname(1,idtsur(itrc))
739 WRITE (vinfo( 2),20) trim(vname(2,idtsur(itrc))), &
740 & trim(text)
741 vinfo( 3)='nondimensional'
742 vinfo(14)=vname(4,idtsur(itrc))
743 vinfo(16)=vname(1,idtime)
744# if defined WRITE_WATER && defined MASKING
745 vinfo(20)='mask_rho'
746# endif
747 vinfo(21)=vname(6,idtsur(itrc))
748 vinfo(22)='coordinates'
749 aval(5)=real(r2dvar,r8)
750 status=def_var(ng, itlm, nrm(ifile,ng)%ncid, &
751 & nrm(ifile,ng)%Vid(idtsur(itrc)), nf_fout, &
752 & nvd3, t2dgrd, aval, vinfo, ncname)
753 IF (founderror(exit_flag, noerror, &
754 & __line__, myfile)) RETURN
755 END IF
756 END DO
757# endif
758# endif
759 END IF
760!
761!-----------------------------------------------------------------------
762! Leave definition mode.
763!-----------------------------------------------------------------------
764!
765 CALL netcdf_enddef (ng, model, ncname, nrm(ifile,ng)%ncid)
766 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
767!
768!-----------------------------------------------------------------------
769! Write out time-recordless, information variables. Deactive file
770! creation switch.
771!-----------------------------------------------------------------------
772!
773 CALL wrt_info (ng, model, nrm(ifile,ng)%ncid, ncname)
774 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
775 ldefnrm(ifile,ng)=.false.
776 END IF define
777!
778!=======================================================================
779! Open an existing normalization file, check its contents, and
780! prepare for appending data.
781!=======================================================================
782!
783 query : IF (.not.ldefnrm(ifile,ng)) THEN
784 ncname=nrm(ifile,ng)%name
785!
786! Open normalization file for read/write.
787!
788 CALL netcdf_open (ng, model, ncname, 1, nrm(ifile,ng)%ncid)
789 IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
790 WRITE (stdout,40) trim(ncname)
791 RETURN
792 END IF
793!
794! Inquire about the dimensions and check for consistency.
795!
796 CALL netcdf_check_dim (ng, model, ncname, &
797 & ncid = nrm(ifile,ng)%ncid)
798 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
799!
800! Inquire about the variables.
801!
802 CALL netcdf_inq_var (ng, model, ncname, &
803 & ncid = nrm(ifile,ng)%ncid)
804 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
805!
806! Initialize logical switches.
807!
808 DO i=1,nv
809 got_var(i)=.false.
810 END DO
811!
812! Scan variable list from input NetCDF and activate switches for
813! normalization variables. Get variable IDs.
814!
815 DO i=1,n_var
816 IF (trim(var_name(i)).eq.trim(vname(1,idtime))) THEN
817 got_var(idtime)=.true.
818 nrm(ifile,ng)%Vid(idtime)=var_id(i)
819 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idfsur))) THEN
820 got_var(idfsur)=.true.
821 nrm(ifile,ng)%Vid(idfsur)=var_id(i)
822 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubar))) THEN
823 got_var(idubar)=.true.
824 nrm(ifile,ng)%Vid(idubar)=var_id(i)
825 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbar))) THEN
826 got_var(idvbar)=.true.
827 nrm(ifile,ng)%Vid(idvbar)=var_id(i)
828# ifdef ADJUST_BOUNDARY
829 ELSE IF (trim(var_name(i)).eq. &
830 & trim(vname(1,idsbry(isfsur)))) THEN
831 got_var(idsbry(isfsur))=.true.
832 nrm(ifile,ng)%Vid(idsbry(isfsur))=var_id(i)
833 ELSE IF (trim(var_name(i)).eq. &
834 & trim(vname(1,idsbry(isubar)))) THEN
835 got_var(idsbry(isubar))=.true.
836 nrm(ifile,ng)%Vid(idsbry(isubar))=var_id(i)
837 ELSE IF (trim(var_name(i)).eq. &
838 & trim(vname(1,idsbry(isvbar)))) THEN
839 got_var(idsbry(isvbar))=.true.
840 nrm(ifile,ng)%Vid(idsbry(isvbar))=var_id(i)
841# endif
842# ifdef SOLVE3D
843 ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduvel))) THEN
844 got_var(iduvel)=.true.
845 nrm(ifile,ng)%Vid(iduvel)=var_id(i)
846 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvvel))) THEN
847 got_var(idvvel)=.true.
848 nrm(ifile,ng)%Vid(idvvel)=var_id(i)
849# ifdef ADJUST_BOUNDARY
850 ELSE IF (trim(var_name(i)).eq. &
851 & trim(vname(1,idsbry(isuvel)))) THEN
852 got_var(idsbry(isuvel))=.true.
853 nrm(ifile,ng)%Vid(idsbry(isuvel))=var_id(i)
854 ELSE IF (trim(var_name(i)).eq. &
855 & trim(vname(1,idsbry(isvvel)))) THEN
856 got_var(idsbry(isvvel))=.true.
857 nrm(ifile,ng)%Vid(idsbry(isvvel))=var_id(i)
858# endif
859# endif
860# ifdef ADJUST_WSTRESS
861 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idusms))) THEN
862 got_var(idusms)=.true.
863 nrm(ifile,ng)%Vid(idusms)=var_id(i)
864 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvsms))) THEN
865 got_var(idvsms)=.true.
866 nrm(ifile,ng)%Vid(idvsms)=var_id(i)
867# endif
868 END IF
869# ifdef SOLVE3D
870 DO itrc=1,nt(ng)
871 IF (trim(var_name(i)).eq.trim(vname(1,idtvar(itrc)))) THEN
872 got_var(idtvar(itrc))=.true.
873 nrm(ifile,ng)%Vid(idtvar(itrc))=var_id(i)
874# ifdef ADJUST_BOUNDARY
875 ELSE IF (trim(var_name(i)).eq. &
876 & trim(vname(1,idsbry(istvar(itrc))))) THEN
877 got_var(idsbry(istvar(itrc)))=.true.
878 nrm(ifile,ng)%Vid(idsbry(istvar(itrc)))=var_id(i)
879# endif
880# ifdef ADJUST_STFLUX
881 ELSE IF (trim(var_name(i)).eq. &
882 & trim(vname(1,idtsur(itrc)))) THEN
883 got_var(idtsur(itrc))=.true.
884 nrm(ifile,ng)%Vid(idtsur(itrc))=var_id(i)
885# endif
886 END IF
887 END DO
888# endif
889 END DO
890!
891! Check if initialization variables are available in input NetCDF
892! file.
893!
894 IF (.not.got_var(idtime)) THEN
895 IF (master) WRITE (stdout,50) trim(vname(1,idtime)), &
896 & trim(ncname)
897 exit_flag=3
898 RETURN
899 END IF
900 IF (.not.got_var(idfsur).and.(ifile.le.2)) THEN
901 IF (master) WRITE (stdout,50) trim(vname(1,idfsur)), &
902 & trim(ncname)
903 exit_flag=3
904 RETURN
905 END IF
906 IF (.not.got_var(idubar).and.(ifile.le.2)) THEN
907 IF (master) WRITE (stdout,50) trim(vname(1,idubar)), &
908 & trim(ncname)
909 exit_flag=3
910 RETURN
911 END IF
912 IF (.not.got_var(idvbar).and.(ifile.le.2)) THEN
913 IF (master) WRITE (stdout,50) trim(vname(1,idvbar)), &
914 & trim(ncname)
915 exit_flag=3
916 RETURN
917 END IF
918# ifdef ADJUST_BOUNDARY
919 IF (.not.got_var(idsbry(isfsur)).and.(ifile.eq.3)) THEN
920 IF (master) WRITE (stdout,50) trim(vname(1,idsbry(isfsur))), &
921 & trim(ncname)
922 exit_flag=3
923 RETURN
924 END IF
925 IF (.not.got_var(idsbry(isubar)).and.(ifile.eq.3)) THEN
926 IF (master) WRITE (stdout,50) trim(vname(1,idsbry(isubar))), &
927 & trim(ncname)
928 exit_flag=3
929 RETURN
930 END IF
931 IF (.not.got_var(idsbry(isvbar)).and.(ifile.eq.3)) THEN
932 IF (master) WRITE (stdout,50) trim(vname(1,idsbry(isvbar))), &
933 & trim(ncname)
934 exit_flag=3
935 RETURN
936 END IF
937# endif
938# ifdef SOLVE3D
939 IF (.not.got_var(iduvel).and.(ifile.le.2)) THEN
940 IF (master) WRITE (stdout,50) trim(vname(1,iduvel)), &
941 & trim(ncname)
942 exit_flag=3
943 RETURN
944 END IF
945 IF (.not.got_var(idvvel).and.(ifile.le.2)) THEN
946 IF (master) WRITE (stdout,50) trim(vname(1,idvvel)), &
947 & trim(ncname)
948 exit_flag=3
949 RETURN
950 END IF
951# ifdef ADJUST_BOUNDARY
952 IF (.not.got_var(idsbry(isuvel)).and.(ifile.eq.3)) THEN
953 IF (master) WRITE (stdout,50) trim(vname(1,idsvar(isuvel))), &
954 & trim(ncname)
955 exit_flag=3
956 RETURN
957 END IF
958 IF (.not.got_var(idsbry(isvvel)).and.(ifile.eq.3)) THEN
959 IF (master) WRITE (stdout,50) trim(vname(1,idsvar(isvvel))), &
960 & trim(ncname)
961 exit_flag=3
962 RETURN
963 END IF
964# endif
965# endif
966# ifdef ADJUST_WSTRESS
967 IF (.not.got_var(idusms).and.(ifile.eq.4)) THEN
968 IF (master) WRITE (stdout,50) trim(vname(1,idusms)), &
969 & trim(ncname)
970 exit_flag=3
971 RETURN
972 END IF
973 IF (.not.got_var(idvsms).and.(ifile.eq.4)) THEN
974 IF (master) WRITE (stdout,50) trim(vname(1,idvsms)), &
975 & trim(ncname)
976 exit_flag=3
977 RETURN
978 END IF
979# endif
980# ifdef SOLVE3D
981 DO itrc=1,nt(ng)
982 IF (.not.got_var(idtvar(itrc)).and.(ifile.le.2)) THEN
983 IF (master) WRITE (stdout,50) trim(vname(1,idtvar(itrc))), &
984 & trim(ncname)
985 exit_flag=3
986 RETURN
987 END IF
988# ifdef ADJUST_BOUNDARY
989 IF (.not.got_var(idsbry(istvar(itrc))).and.(ifile.eq.3)) THEN
990 IF (master) WRITE (stdout,50) &
991 & trim(vname(1,idsbry(istvar(itrc)))), &
992 & trim(ncname)
993 exit_flag=3
994 RETURN
995 END IF
996# endif
997# ifdef ADJUST_STFLUX
998 IF (.not.got_var(idtsur(itrc)).and.(ifile.eq.4).and. &
999 & lstflux(itrc,ng)) THEN
1000 IF (master) WRITE (stdout,50) trim(vname(1,idtsur(itrc))), &
1001 & trim(ncname)
1002 exit_flag=3
1003 RETURN
1004 END IF
1005# endif
1006 END DO
1007# endif
1008!
1009! Set unlimited time record dimension to zero to allow other programs
1010! to process and write normalization factors for different variable.
1011!
1012 nrm(ifile,ng)%Rindex=0
1013 END IF query
1014!
1015 10 FORMAT (/,' DEF_NORM_NF90 - unable to create norm NetCDF', &
1016 & ' file: ',a)
1017 20 FORMAT (a,', ',a)
1018 30 FORMAT (1pe11.4,1x,'millimeter')
1019 40 FORMAT (/,' DEF_NORM_NF90 - unable to open norm NetCDF file: ',a)
1020 50 FORMAT (/,' DEF_NORM_NF90 - unable to find variable: ',a,2x, &
1021 & ' in norm NetCDF file: ',a)
1022!
1023 RETURN
1024 END SUBROUTINE def_norm_nf90
1025
1026# if defined PIO_LIB && defined DISTRIBUTE
1027!
1028!***********************************************************************
1029 SUBROUTINE def_norm_pio (ng, model, ifile)
1030!***********************************************************************
1031!
1032 USE mod_pio_netcdf
1033!
1034! Imported variable declarations.
1035!
1036 integer, intent(in) :: ng, model, ifile
1037!
1038! Local variable declarations.
1039!
1040 logical :: ldefine, got_var(nv)
1041!
1042 integer, parameter :: natt = 25
1043
1044 integer :: i, j, nvd3, nvd4
1045 integer :: recdim, status
1046# ifdef ADJUST_BOUNDARY
1047 integer :: iorjdim, ifield
1048# endif
1049 integer :: dimids(ndimid)
1050 integer :: t2dgrd(3), u2dgrd(3), v2dgrd(3)
1051# ifdef ADJUST_BOUNDARY
1052 integer :: t2dobc(3)
1053# endif
1054# ifdef SOLVE3D
1055 integer :: itrc
1056
1057 integer :: t3dgrd(4), u3dgrd(4), v3dgrd(4)
1058# ifdef ADJUST_BOUNDARY
1059 integer :: t3dobc(4)
1060# endif
1061# endif
1062!
1063 real(r8) :: aval(6)
1064!
1065 character (len=60 ) :: text
1066 character (len=256) :: ncname
1067 character (len=MaxLen) :: vinfo(natt)
1068
1069 character (len=*), parameter :: myfile = &
1070 & __FILE__//", def_norm_pio"
1071!
1072 sourcefile=myfile
1073!
1074!=======================================================================
1075! Create a new background covariace normalization file.
1076!=======================================================================
1077!
1078 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1079 ncname=nrm(ifile,ng)%name
1080!
1081 define : IF (ldefnrm(ifile,ng)) THEN
1082 CALL pio_netcdf_create (ng, itlm, trim(ncname), &
1083 & nrm(ifile,ng)%pioFile)
1084 IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
1085 IF (master) WRITE (stdout,10) trim(ncname)
1086 RETURN
1087 END IF
1088!
1089!-----------------------------------------------------------------------
1090! Define the dimensions of staggered fields.
1091!-----------------------------------------------------------------------
1092!
1093 dimids=0
1094!
1095 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1096 & 'xi_rho', iobounds(ng)%xi_rho, dimids( 1))
1097 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1098
1099 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1100 & 'xi_u',iobounds(ng)%xi_u, dimids( 2))
1101 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1102
1103 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1104 & 'xi_v', iobounds(ng)%xi_v, dimids( 3))
1105 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1106
1107 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1108 & 'xi_psi', iobounds(ng)%xi_psi, dimids( 4))
1109 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1110
1111 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1112 & 'eta_rho', iobounds(ng)%eta_rho, dimids( 5))
1113 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1114
1115 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1116 & 'eta_u', iobounds(ng)%eta_u, dimids( 6))
1117 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1118
1119 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1120 & 'eta_v', iobounds(ng)%eta_v, dimids( 7))
1121 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1122
1123 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1124 & 'eta_psi', iobounds(ng)%eta_psi, dimids( 8))
1125 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1126
1127# ifdef ADJUST_BOUNDARY
1128 IF (ifile.eq.3) THEN
1129 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1130 & 'IorJ', iobounds(ng)%IorJ, iorjdim)
1131 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1132 END IF
1133# endif
1134
1135# if defined WRITE_WATER && defined MASKING
1136 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1137 & 'xy_rho', iobounds(ng)%xy_rho, dimids(17))
1138 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1139
1140 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1141 & 'xy_u', iobounds(ng)%xy_u, dimids(18))
1142 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1143
1144 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname,
1145 & 'xy_v', iobounds(ng)%xy_v, dimids(19))
1146 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1147# endif
1148
1149# ifdef SOLVE3D
1150# if defined WRITE_WATER && defined MASKING
1151 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1152 & 'xyz_rho', iobounds(ng)%xy_rho*n(ng), dimids(20))
1153 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1154
1155 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1156 & 'xyz_u', iobounds(ng)%xy_u*n(ng), dimids(21))
1157 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1158
1159 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1160 & 'xyz_v', iobounds(ng)%xy_v*n(ng), dimids(22))
1161 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1162
1163 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1164 & 'xyz_w', iobounds(ng)%xy_rho*(n(ng)+1), &
1165 & dimids(23))
1166 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1167# endif
1168
1169 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1170 & 's_rho', n(ng), dimids( 9))
1171 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1172
1173 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1174 & 's_w', n(ng)+1, dimids(10))
1175 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1176
1177 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1178 & 'tracer', nt(ng), dimids(11))
1179 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1180
1181# ifdef SEDIMENT
1182 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1183 & 'NST', nst, dimids(32))
1184 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1185
1186 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1187 & 'Nbed', nbed, dimids(16))
1188 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1189
1190# if defined WRITE_WATER && defined MASKING
1191 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1192 & 'xybed', iobounds(ng)%xy_rho*nbed, dimids(24))
1193 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1194# endif
1195# endif
1196
1197# ifdef ECOSIM
1198 status=def_dim(ng, inlm, nrm(ifile,ng)%pioFile, ncname, &
1199 & 'Nbands', nbands, dimids(33))
1200 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1201
1202 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1203 & 'Nphy', nphy, dimids(25))
1204 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1205
1206 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1207 & 'Nbac', nbac, dimids(26))
1208 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1209
1210 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname,
1211 & 'Ndom', ndom, dimids(27))
1212 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1213
1214 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1215 & 'Nfec', nfec, dimids(28))
1216 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1217# endif
1218# endif
1219
1220 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1221 & 'boundary', 4, dimids(14))
1222 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1223
1224# ifdef FOUR_DVAR
1225 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1226 & 'Nstate', nstatevar(ng), dimids(29))
1227 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1228# endif
1229
1230 status=def_dim(ng, itlm, nrm(ifile,ng)%pioFile, ncname, &
1231 & trim(adjustl(vname(5,idtime))), &
1232 & pio_unlimited, dimids(12))
1233 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1234
1235 recdim=dimids(12)
1236!
1237! Set number of dimensions for output variables.
1238!
1239# if defined WRITE_WATER && defined MASKING
1240 nvd3=2
1241 nvd4=2
1242# else
1243 nvd3=3
1244 nvd4=4
1245# endif
1246!
1247! Define dimension vectors for staggered tracer type variables.
1248!
1249# if defined WRITE_WATER && defined MASKING
1250 t2dgrd(1)=dimids(17)
1251 t2dgrd(2)=dimids(12)
1252# ifdef SOLVE3D
1253 t3dgrd(1)=dimids(20)
1254 t3dgrd(2)=dimids(12)
1255# endif
1256# else
1257 t2dgrd(1)=dimids( 1)
1258 t2dgrd(2)=dimids( 5)
1259 t2dgrd(3)=dimids(12)
1260# ifdef SOLVE3D
1261 t3dgrd(1)=dimids( 1)
1262 t3dgrd(2)=dimids( 5)
1263 t3dgrd(3)=dimids( 9)
1264 t3dgrd(4)=dimids(12)
1265# endif
1266# endif
1267# ifdef ADJUST_BOUNDARY
1268 IF (ifile.eq.3) THEN
1269 t2dobc(1)=iorjdim
1270 t2dobc(2)=dimids(14)
1271 t2dobc(3)=dimids(12)
1272# ifdef SOLVE3D
1273 t3dobc(1)=iorjdim
1274 t3dobc(2)=dimids( 9)
1275 t3dobc(3)=dimids(14)
1276 t3dobc(4)=dimids(12)
1277# endif
1278 END IF
1279# endif
1280!
1281! Define dimension vectors for staggered u-momentum type variables.
1282!
1283# if defined WRITE_WATER && defined MASKING
1284 u2dgrd(1)=dimids(18)
1285 u2dgrd(2)=dimids(12)
1286# ifdef SOLVE3D
1287 u3dgrd(1)=dimids(21)
1288 u3dgrd(2)=dimids(12)
1289# endif
1290# else
1291 u2dgrd(1)=dimids( 2)
1292 u2dgrd(2)=dimids( 6)
1293 u2dgrd(3)=dimids(12)
1294# ifdef SOLVE3D
1295 u3dgrd(1)=dimids( 2)
1296 u3dgrd(2)=dimids( 6)
1297 u3dgrd(3)=dimids( 9)
1298 u3dgrd(4)=dimids(12)
1299# endif
1300# endif
1301!
1302! Define dimension vectors for staggered v-momentum type variables.
1303!
1304# if defined WRITE_WATER && defined MASKING
1305 v2dgrd(1)=dimids(19)
1306 v2dgrd(2)=dimids(12)
1307# ifdef SOLVE3D
1308 v3dgrd(1)=dimids(22)
1309 v3dgrd(2)=dimids(12)
1310# endif
1311# else
1312 v2dgrd(1)=dimids( 3)
1313 v2dgrd(2)=dimids( 7)
1314 v2dgrd(3)=dimids(12)
1315# ifdef SOLVE3D
1316 v3dgrd(1)=dimids( 3)
1317 v3dgrd(2)=dimids( 7)
1318 v3dgrd(3)=dimids( 9)
1319 v3dgrd(4)=dimids(12)
1320# endif
1321# endif
1322!
1323! Initialize unlimited time record dimension.
1324!
1325 nrm(ifile,ng)%Rindex=0
1326!
1327! Initialize local information variable arrays.
1328!
1329 DO i=1,natt
1330 DO j=1,len(vinfo(1))
1331 vinfo(i)(j:j)=' '
1332 END DO
1333 END DO
1334 DO i=1,6
1335 aval(i)=0.0_r8
1336 END DO
1337!
1338!-----------------------------------------------------------------------
1339! Define time-recordless information variables.
1340!-----------------------------------------------------------------------
1341!
1342 CALL def_info (ng, itlm, nrm(ifile,ng)%pioFile, ncname, dimids)
1343 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1344!
1345!-----------------------------------------------------------------------
1346! Define background covariance normalization variables.
1347!-----------------------------------------------------------------------
1348!
1349! Define model time.
1350!
1351 vinfo( 1)=vname(1,idtime)
1352 vinfo( 2)=vname(2,idtime)
1353 WRITE (vinfo( 3),'(a,a)') 'seconds since ', trim(rclock%string)
1354 vinfo( 4)=trim(rclock%calendar)
1355 vinfo(14)=vname(4,idtime)
1356 vinfo(21)=vname(6,idtime)
1357 nrm(ifile,ng)%pioVar(idtime)%dkind=pio_tout
1358 nrm(ifile,ng)%pioVar(idtime)%gtype=0
1359!
1360 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1361 & nrm(ifile,ng)%pioVar(idtime)%vd, &
1362 & pio_tout, 1, (/recdim/), aval, vinfo, ncname, &
1363 & setparaccess = .true.)
1364 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1365!
1366!-----------------------------------------------------------------------
1367! Initial conditions or model error covariance normalization factors.
1368!-----------------------------------------------------------------------
1369!
1370 IF ((ifile.eq.1).or.(ifile.eq.2)) THEN
1371 IF (ifile.eq.1) THEN
1372 text='initial conditions error covariance normalization'
1373 ELSE IF (ifile.eq.2) THEN
1374 text='model error covariance normalization'
1375 END IF
1376!
1377! Define free-surface normalization factor.
1378!
1379 vinfo( 1)=vname(1,idfsur)
1380 WRITE (vinfo( 2),20) trim(vname(2,idfsur)), trim(text)
1381 vinfo( 3)='nondimensional'
1382 vinfo(14)=vname(4,idfsur)
1383 vinfo(16)=vname(1,idtime)
1384# if defined WRITE_WATER && defined MASKING
1385 vinfo(20)='mask_rho'
1386# endif
1387 vinfo(21)=vname(6,idfsur)
1388 vinfo(22)='coordinates'
1389 aval(5)=real(iinfo(1,idfsur,ng),r8)
1390 nrm(ifile,ng)%pioVar(idfsur)%dkind=pio_fout
1391 nrm(ifile,ng)%pioVar(idfsur)%gtype=r2dvar
1392!
1393 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1394 & nrm(ifile,ng)%pioVar(idfsur)%vd, &
1395 & pio_fout, nvd3, t2dgrd, aval, vinfo, ncname)
1396 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1397!
1398! Define 2D U-momentum normalization factor.
1399!
1400 vinfo( 1)=vname(1,idubar)
1401 WRITE (vinfo( 2),20) trim(vname(2,idubar)), trim(text)
1402 vinfo( 3)='nondimensional'
1403 vinfo(14)=vname(4,idubar)
1404 vinfo(16)=vname(1,idtime)
1405# if defined WRITE_WATER && defined MASKING
1406 vinfo(20)='mask_u'
1407# endif
1408 vinfo(21)=vname(6,idubar)
1409 vinfo(22)='coordinates'
1410 aval(5)=real(iinfo(1,idubar,ng),r8)
1411 nrm(ifile,ng)%pioVar(idubar)%dkind=pio_fout
1412 nrm(ifile,ng)%pioVar(idubar)%gtype=u2dvar
1413!
1414 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1415 & nrm(ifile,ng)%pioVar(idubar)%vd, &
1416 & pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
1417 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1418!
1419! Define 2D V-momentum normalization factor.
1420!
1421 vinfo( 1)=vname(1,idvbar)
1422 WRITE (vinfo( 2),20) trim(vname(2,idvbar)), trim(text)
1423 vinfo( 3)='nondimensional'
1424 vinfo(14)=vname(4,idvbar)
1425 vinfo(16)=vname(1,idtime)
1426# if defined WRITE_WATER && defined MASKING
1427 vinfo(20)='mask_v'
1428# endif
1429 vinfo(21)=vname(6,idvbar)
1430 vinfo(22)='coordinates'
1431 aval(5)=real(iinfo(1,idvbar,ng),r8)
1432 nrm(ifile,ng)%pioVar(idvbar)%dkind=pio_fout
1433 nrm(ifile,ng)%pioVar(idvbar)%gtype=v2dvar
1434!
1435 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1436 & nrm(ifile,ng)%pioVar(idvbar)%vd, &
1437 & pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
1438 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1439
1440# ifdef SOLVE3D
1441!
1442! Define 3D U-momentum normalization factor.
1443!
1444 vinfo( 1)=vname(1,iduvel)
1445 WRITE (vinfo( 2),20) trim(vname(2,iduvel)), trim(text)
1446 vinfo( 3)='nondimensional'
1447 vinfo(14)=vname(4,iduvel)
1448 vinfo(16)=vname(1,idtime)
1449# if defined WRITE_WATER && defined MASKING
1450 vinfo(20)='mask_u'
1451# endif
1452 vinfo(21)=vname(6,iduvel)
1453 vinfo(22)='coordinates'
1454 aval(5)=real(iinfo(1,iduvel,ng),r8)
1455 nrm(ifile,ng)%pioVar(iduvel)%dkind=pio_fout
1456 nrm(ifile,ng)%pioVar(iduvel)%gtype=u3dvar
1457!
1458 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1459 & nrm(ifile,ng)%pioVar(iduvel)%vd, &
1460 & pio_fout, nvd4, u3dgrd, aval, vinfo, ncname)
1461 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1462!
1463! Define 3D V-momentum normalization factor.
1464!
1465 vinfo( 1)=vname(1,idvvel)
1466 WRITE (vinfo( 2),20) trim(vname(2,idvvel)), trim(text)
1467 vinfo( 3)='nondimensional'
1468 vinfo(14)=vname(4,idvvel)
1469 vinfo(16)=vname(1,idtime)
1470# if defined WRITE_WATER && defined MASKING
1471 vinfo(20)='mask_v'
1472# endif
1473 vinfo(21)=vname(6,idvvel)
1474 vinfo(22)='coordinates'
1475 aval(5)=real(iinfo(1,idvvel,ng),r8)
1476 nrm(ifile,ng)%pioVar(idvvel)%dkind=pio_fout
1477 nrm(ifile,ng)%pioVar(idvvel)%gtype=v3dvar
1478!
1479 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1480 & nrm(ifile,ng)%pioVar(idvvel)%vd, &
1481 & pio_fout, nvd4, v3dgrd, aval, vinfo, ncname)
1482 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1483!
1484! Define tracer type normalization factors.
1485!
1486 DO itrc=1,nt(ng)
1487 vinfo( 1)=vname(1,idtvar(itrc))
1488 WRITE (vinfo( 2),20) trim(vname(2,idtvar(itrc))), trim(text)
1489 vinfo( 3)='nondimensional'
1490 vinfo(14)=vname(4,idtvar(itrc))
1491 vinfo(16)=vname(1,idtime)
1492# ifdef SEDIMENT
1493 DO i=1,nst
1494 IF (itrc.eq.idsed(i)) THEN
1495 WRITE (vinfo(19),30) 1000.0_r8*sd50(i,ng)
1496 END IF
1497 END DO
1498# endif
1499# if defined WRITE_WATER && defined MASKING
1500 vinfo(20)='mask_rho'
1501# endif
1502 vinfo(21)=vname(6,idtvar(itrc))
1503 vinfo(22)='coordinates'
1504 aval(5)=real(r3dvar,r8)
1505 nrm(ifile,ng)%pioTrc(itrc)%dkind=pio_fout
1506 nrm(ifile,ng)%pioTrc(itrc)%gtype=r3dvar
1507!
1508 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1509 & nrm(ifile,ng)%pioTrc(itrc)%vd, &
1510 & pio_fout, nvd4, t3dgrd, aval, vinfo, ncname)
1511 END DO
1512 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1513# endif
1514
1515# ifdef ADJUST_BOUNDARY
1516!
1517!-----------------------------------------------------------------------
1518! Boundary conditions error covariance normalization factors.
1519!-----------------------------------------------------------------------
1520!
1521 ELSE IF (ifile.eq.3) THEN
1522 text='error covariance normalization'
1523!
1524! Define free-surface open boundaries.
1525!
1526 IF (any(lobc(:,isfsur,ng))) THEN
1527 ifield=idsbry(isfsur)
1528 vinfo( 1)=vname(1,ifield)
1529 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
1530 vinfo( 3)=vname(3,ifield)
1531 vinfo(14)=vname(4,ifield)
1532 vinfo(16)=vname(1,idtime)
1533 vinfo(21)=vname(6,ifield)
1534 aval(5)=real(iinfo(1,ifield,ng),r8)
1535 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1536 nrm(ifile,ng)%pioVar(ifield)%gtype=r2dobc
1537!
1538 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1539 & nrm(ifile,ng)%pioVar(ifield)%vd, &
1540 & pio_fout, 3, t2dobc, aval, vinfo, ncname, &
1541 & setfillval = .false.)
1542 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1543 END IF
1544!
1545! Define 2D U-momentum component open boundaries.
1546!
1547 IF (any(lobc(:,isubar,ng))) THEN
1548 ifield=idsbry(isubar)
1549 vinfo( 1)=vname(1,ifield)
1550 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
1551 vinfo( 3)=vname(3,ifield)
1552 vinfo(14)=vname(4,ifield)
1553 vinfo(16)=vname(1,idtime)
1554 vinfo(21)=vname(6,ifield)
1555 aval(5)=real(iinfo(1,ifield,ng),r8)
1556 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1557 nrm(ifile,ng)%pioVar(ifield)%gtype=u2dobc
1558!
1559 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1560 & nrm(ifile,ng)%pioVar(ifield)%vd, &
1561 & pio_fout, 3, t2dobc, aval, vinfo, ncname, &
1562 & setfillval = .false.)
1563 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1564 END IF
1565!
1566! Define 2D V-momentum component open boundaries.
1567!
1568 IF (any(lobc(:,isvbar,ng))) THEN
1569 ifield=idsbry(isvbar)
1570 vinfo( 1)=vname(1,ifield)
1571 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
1572 vinfo( 3)=vname(3,ifield)
1573 vinfo(14)=vname(4,ifield)
1574 vinfo(16)=vname(1,idtime)
1575 vinfo(21)=vname(6,ifield)
1576 aval(5)=real(iinfo(1,ifield,ng),r8)
1577 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1578 nrm(ifile,ng)%pioVar(ifield)%gtype=v2dobc
1579!
1580 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1581 & nrm(ifile,ng)%pioVar(ifield)%vd, &
1582 & pio_fout, 3, t2dobc, aval, vinfo, ncname, &
1583 & setfillval = .false.)
1584 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1585 END IF
1586
1587# ifdef SOLVE3D
1588!
1589! Define 3D U-momentum component open boundaries.
1590!
1591 IF (any(lobc(:,isuvel,ng))) THEN
1592 ifield=idsbry(isuvel)
1593 vinfo( 1)=vname(1,ifield)
1594 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
1595 vinfo( 3)=vname(3,ifield)
1596 vinfo(14)=vname(4,ifield)
1597 vinfo(16)=vname(1,idtime)
1598 vinfo(21)=vname(6,ifield)
1599 aval(5)=real(iinfo(1,ifield,ng),r8)
1600 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1601 nrm(ifile,ng)%pioVar(ifield)%gtype=u3dobc
1602!
1603 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1604 & nrm(ifile,ng)%pioVar(ifield)%vd, &
1605 & pio_fout, 4, t3dobc, aval, vinfo, ncname, &
1606 & setfillval = .false.)
1607 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1608 END IF
1609!
1610! Define 3D V-momentum component open boundaries.
1611!
1612 IF (any(lobc(:,isvvel,ng))) THEN
1613 ifield=idsbry(isvvel)
1614 vinfo( 1)=vname(1,ifield)
1615 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
1616 vinfo( 3)=vname(3,ifield)
1617 vinfo(14)=vname(4,ifield)
1618 vinfo(16)=vname(1,idtime)
1619 vinfo(21)=vname(6,ifield)
1620 aval(5)=real(iinfo(1,ifield,ng),r8)
1621 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1622 nrm(ifile,ng)%pioVar(ifield)%gtype=v3dobc
1623!
1624 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1625 & nrm(ifile,ng)%pioVar(ifield)%vd, &
1626 & pio_fout, 4, t3dobc, aval, vinfo, ncname, &
1627 & setfillval = .false.)
1628 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1629 END IF
1630!
1631! Define tracer type variables open boundaries.
1632!
1633 DO itrc=1,nt(ng)
1634 IF (any(lobc(:,istvar(itrc),ng))) THEN
1635 ifield=idsbry(istvar(itrc))
1636 vinfo( 1)=vname(1,ifield)
1637 WRITE (vinfo( 2),20) trim(vname(2,ifield)), trim(text)
1638 vinfo( 3)=vname(3,ifield)
1639 vinfo(14)=vname(4,ifield)
1640 vinfo(16)=vname(1,idtime)
1641# ifdef SEDIMENT
1642 DO i=1,nst
1643 IF (itrc.eq.idsed(i)) THEN
1644 WRITE (vinfo(19),30) 1000.0_r8*sd50(i,ng)
1645 END IF
1646 END DO
1647# endif
1648 vinfo(21)=vname(6,ifield)
1649 aval(5)=real(iinfo(1,ifield,ng),r8)
1650 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1651 nrm(ifile,ng)%pioVar(ifield)%gtype=r3dobc
1652!
1653 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1654 & nrm(ifile,ng)%pioVar(ifield)%vd, &
1655 & pio_fout, 4, t3dobc, aval, vinfo, ncname, &
1656 & setfillval = .false.)
1657 IF (founderror(exit_flag, noerror, &
1658 & __line__, myfile)) RETURN
1659 END IF
1660 END DO
1661# endif
1662# endif
1663
1664# if defined ADJUST_WSTRESS || defined ADJUST_STFLUX
1665!
1666!-----------------------------------------------------------------------
1667! Surface forcing error covariance normalization factors.
1668!-----------------------------------------------------------------------
1669!
1670 ELSE IF (ifile.eq.4) THEN
1671 text='error covariance normalization'
1672
1673# ifdef ADJUST_WSTRESS
1674!
1675! Define surface U-momentum stress normalization factors.
1676!
1677 vinfo( 1)=vname(1,idusms)
1678 WRITE (vinfo( 2),20) trim(vname(2,idusms)), trim(text)
1679 vinfo( 3)='nondimensional'
1680 vinfo(14)=vname(4,idusms)
1681 vinfo(16)=vname(1,idtime)
1682# if defined WRITE_WATER && defined MASKING
1683 vinfo(20)='mask_u'
1684# endif
1685 vinfo(21)=vname(6,idusms)
1686 vinfo(22)='coordinates'
1687 aval(5)=real(u2dvar,r8)
1688 nrm(ifile,ng)%pioVar(idusms)%dkind=pio_fout
1689 nrm(ifile,ng)%pioVar(idusms)%gtype=u2dvar
1690!
1691 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1692 & nrm(ifile,ng)%pioVar(idusms)%vd, &
1693 & pio_fout, nvd3, u2dgrd, aval, vinfo, ncname)
1694 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1695!
1696! Define surface V-momentum stress normalization factors.
1697!
1698 vinfo( 1)=vname(1,idvsms)
1699 WRITE (vinfo( 2),20) trim(vname(2,idvsms)), trim(text)
1700 vinfo( 2)=vname(2,idvsms)
1701 vinfo( 3)='nondimensional'
1702 vinfo(14)=vname(4,idvsms)
1703 vinfo(16)=vname(1,idtime)
1704# if defined WRITE_WATER && defined MASKING
1705 vinfo(20)='mask_v'
1706# endif
1707 vinfo(21)=vname(6,idvsms)
1708 vinfo(22)='coordinates'
1709 aval(5)=real(v2dvar,r8)
1710 nrm(ifile,ng)%pioVar(idvsms)%dkind=pio_fout
1711 nrm(ifile,ng)%pioVar(idvsms)%gtype=v2dvar
1712!
1713 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1714 & nrm(ifile,ng)%pioVar(idvsms)%vd, &
1715 & pio_fout, nvd3, v2dgrd, aval, vinfo, ncname)
1716 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1717# endif
1718
1719# if defined ADJUST_STFLUX && defined SOLVE3D
1720!
1721! Define surface tracer fluxes.
1722!
1723 DO itrc=1,nt(ng)
1724 IF (lstflux(itrc,ng)) THEN
1725 vinfo( 1)=vname(1,idtsur(itrc))
1726 WRITE (vinfo( 2),20) trim(vname(2,idtsur(itrc))), &
1727 & trim(text)
1728 vinfo( 3)='nondimensional'
1729 vinfo(14)=vname(4,idtsur(itrc))
1730 vinfo(16)=vname(1,idtime)
1731# if defined WRITE_WATER && defined MASKING
1732 vinfo(20)='mask_rho'
1733# endif
1734 vinfo(21)=vname(6,idtsur(itrc))
1735 vinfo(22)='coordinates'
1736 aval(5)=real(r2dvar,r8)
1737 nrm(ifile,ng)%pioVar(idtsur(itrc))%dkind=pio_fout
1738 nrm(ifile,ng)%pioVar(idtsur(itrc))%gtype=r2dvar
1739!
1740 status=def_var(ng, itlm, nrm(ifile,ng)%pioFile, &
1741 & nrm(ifile,ng)%pioVar(idtsur(itrc))%vd, &
1742 & pio_fout, nvd3, t2dgrd, aval, vinfo, &
1743 & ncname)
1744 IF (founderror(exit_flag, noerror, &
1745 & __line__, myfile)) RETURN
1746 END IF
1747 END DO
1748# endif
1749# endif
1750 END IF
1751!
1752!-----------------------------------------------------------------------
1753! Leave definition mode.
1754!-----------------------------------------------------------------------
1755!
1756 CALL pio_netcdf_enddef (ng, model, ncname, &
1757 & nrm(ifile,ng)%pioFile)
1758 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1759!
1760!-----------------------------------------------------------------------
1761! Write out time-recordless, information variables. Deactive file
1762! creation switch.
1763!-----------------------------------------------------------------------
1764!
1765 CALL wrt_info (ng, model, nrm(ifile,ng)%pioFile, ncname)
1766 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1767 ldefnrm(ifile,ng)=.false.
1768 END IF define
1769!
1770!=======================================================================
1771! Open an existing normalization file, check its contents, and
1772! prepare for appending data.
1773!=======================================================================
1774!
1775 query : IF (.not.ldefnrm(ifile,ng)) THEN
1776 ncname=nrm(ifile,ng)%name
1777!
1778! Open normalization file for read/write.
1779!
1780 CALL pio_netcdf_open (ng, model, ncname, 1, &
1781 & nrm(ifile,ng)%pioFile)
1782 IF (founderror(exit_flag, noerror, __line__, myfile)) THEN
1783 WRITE (stdout,40) trim(ncname)
1784 RETURN
1785 END IF
1786!
1787! Inquire about the dimensions and check for consistency.
1788!
1789 CALL pio_netcdf_check_dim (ng, model, ncname, &
1790 & piofile = nrm(ifile,ng)%pioFile)
1791 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1792!
1793! Inquire about the variables.
1794!
1795 CALL pio_netcdf_inq_var (ng, model, ncname, &
1796 & piofile = nrm(ifile,ng)%pioFile)
1797 IF (founderror(exit_flag, noerror, __line__, myfile)) RETURN
1798!
1799! Initialize logical switches.
1800!
1801 DO i=1,nv
1802 got_var(i)=.false.
1803 END DO
1804!
1805! Scan variable list from input NetCDF and activate switches for
1806! normalization variables. Get variable IDs.
1807!
1808 DO i=1,n_var
1809 IF (trim(var_name(i)).eq.trim(vname(1,idtime))) THEN
1810 got_var(idtime)=.true.
1811 nrm(ifile,ng)%pioVar(idtime)%vd=var_desc(i)
1812 nrm(ifile,ng)%pioVar(idtime)%dkind=pio_tout
1813 nrm(ifile,ng)%pioVar(idtime)%gtype=0
1814 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idfsur))) THEN
1815 got_var(idfsur)=.true.
1816 nrm(ifile,ng)%pioVar(idfsur)%vd=var_desc(i)
1817 nrm(ifile,ng)%pioVar(idfsur)%dkind=pio_fout
1818 nrm(ifile,ng)%pioVar(idfsur)%gtype=r2dvar
1819 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idubar))) THEN
1820 got_var(idubar)=.true.
1821 nrm(ifile,ng)%pioVar(idubar)%vd=var_desc(i)
1822 nrm(ifile,ng)%pioVar(idubar)%dkind=pio_fout
1823 nrm(ifile,ng)%pioVar(idubar)%gtype=u2dvar
1824 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvbar))) THEN
1825 got_var(idvbar)=.true.
1826 nrm(ifile,ng)%pioVar(idvbar)%vd=var_desc(i)
1827 nrm(ifile,ng)%pioVar(idvbar)%dkind=pio_fout
1828 nrm(ifile,ng)%pioVar(idvbar)%gtype=v2dvar
1829# ifdef ADJUST_BOUNDARY
1830 ELSE IF (trim(var_name(i)).eq. &
1831 & trim(vname(1,idsbry(isfsur)))) THEN
1832 ifield=idsbry(isfsur)
1833 got_var(ifield)=.true.
1834 nrm(ifile,ng)%pioVar(ifield)%vd=var_desc(i)
1835 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1836 nrm(ifile,ng)%pioVar(ifield)%gtype=r2dobc
1837 ELSE IF (trim(var_name(i)).eq. &
1838 & trim(vname(1,idsbry(isubar)))) THEN
1839 ifield=idsbry(isubar)
1840 got_var(ifield)=.true.
1841 nrm(ifile,ng)%pioVar(ifield)%vd=var_desc(i)
1842 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1843 nrm(ifile,ng)%pioVar(ifield)%gtype=u2dobc
1844 ELSE IF (trim(var_name(i)).eq. &
1845 & trim(vname(1,idsbry(isvbar)))) THEN
1846 ifield=idsbry(isvbar)
1847 got_var(ifield)=.true.
1848 nrm(ifile,ng)%pioVar(ifield)%vd=var_desc(i)
1849 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1850 nrm(ifile,ng)%pioVar(ifield)%gtype=v2dobc
1851# endif
1852# ifdef SOLVE3D
1853 ELSE IF (trim(var_name(i)).eq.trim(vname(1,iduvel))) THEN
1854 got_var(iduvel)=.true.
1855 nrm(ifile,ng)%pioVar(iduvel)%vd=var_desc(i)
1856 nrm(ifile,ng)%pioVar(iduvel)%dkind=pio_fout
1857 nrm(ifile,ng)%pioVar(iduvel)%gtype=u3dvar
1858 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvvel))) THEN
1859 got_var(idvvel)=.true.
1860 nrm(ifile,ng)%pioVar(idvvel)%vd=var_desc(i)
1861 nrm(ifile,ng)%pioVar(idvvel)%dkind=pio_fout
1862 nrm(ifile,ng)%pioVar(idvvel)%gtype=v3dvar
1863# ifdef ADJUST_BOUNDARY
1864 ELSE IF (trim(var_name(i)).eq. &
1865 & trim(vname(1,idsbry(isuvel)))) THEN
1866 ifield=idsbry(isuvel)
1867 got_var(ifield)=.true.
1868 nrm(ifile,ng)%pioVar(ifield)%vd=var_desc(i)
1869 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1870 nrm(ifile,ng)%pioVar(ifield)%gtype=u3dobc
1871 ELSE IF (trim(var_name(i)).eq. &
1872 & trim(vname(1,idsbry(isvvel)))) THEN
1873 ifield=idsbry(isvvel)
1874 got_var(ifield)=.true.
1875 nrm(ifile,ng)%pioVar(ifield)%vd=var_desc(i)
1876 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1877 nrm(ifile,ng)%pioVar(ifield)%gtype=v3dobc
1878# endif
1879# endif
1880# ifdef ADJUST_WSTRESS
1881 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idusms))) THEN
1882 got_var(idusms)=.true.
1883 nrm(ifile,ng)%pioVar(idusms)%vd=var_desc(i)
1884 nrm(ifile,ng)%pioVar(idusms)%dkind=pio_fout
1885 nrm(ifile,ng)%pioVar(idusms)%gtype=u2dvar
1886 ELSE IF (trim(var_name(i)).eq.trim(vname(1,idvsms))) THEN
1887 got_var(idvsms)=.true.
1888 nrm(ifile,ng)%pioVar(idvsms)%vd=var_desc(i)
1889 nrm(ifile,ng)%pioVar(idvsms)%dkind=pio_fout
1890 nrm(ifile,ng)%pioVar(idvsms)%gtype=v2dvar
1891# endif
1892 END IF
1893# ifdef SOLVE3D
1894 DO itrc=1,nt(ng)
1895 IF (trim(var_name(i)).eq.trim(vname(1,idtvar(itrc)))) THEN
1896 got_var(idtvar(itrc))=.true.
1897 nrm(ifile,ng)%pioTrc(itrc)%vd=var_desc(i)
1898 nrm(ifile,ng)%pioTrc(itrc)%dkind=pio_fout
1899 nrm(ifile,ng)%pioTrc(itrc)%gtype=r3dvar
1900# ifdef ADJUST_BOUNDARY
1901 ELSE IF (trim(var_name(i)).eq. &
1902 & trim(vname(1,idsbry(istvar(itrc))))) THEN
1903 ifield=idsbry(istvar(itrc))
1904 got_var(ifield)=.true.
1905 nrm(ifile,ng)%pioVar(ifield)%vd=var_desc(i)
1906 nrm(ifile,ng)%pioVar(ifield)%dkind=pio_fout
1907 nrm(ifile,ng)%pioVar(ifield)%gtype=r3dobc
1908# endif
1909# ifdef ADJUST_STFLUX
1910 ELSE IF (trim(var_name(i)).eq. &
1911 & trim(vname(1,idtsur(itrc)))) THEN
1912 got_var(idtsur(itrc))=.true.
1913 nrm(ifile,ng)%pioVar(idtsur(itrc))%vd=var_desc(i)
1914 nrm(ifile,ng)%pioVar(idtsur(itrc))%dkind=pio_fout
1915 nrm(ifile,ng)%pioVar(idtsur(itrc))%gtype=r2dvar
1916# endif
1917 END IF
1918 END DO
1919# endif
1920 END DO
1921!
1922! Check if initialization variables are available in input NetCDF
1923! file.
1924!
1925 IF (.not.got_var(idtime)) THEN
1926 IF (master) WRITE (stdout,50) trim(vname(1,idtime)), &
1927 & trim(ncname)
1928 exit_flag=3
1929 RETURN
1930 END IF
1931 IF (.not.got_var(idfsur).and.(ifile.le.2)) THEN
1932 IF (master) WRITE (stdout,50) trim(vname(1,idfsur)), &
1933 & trim(ncname)
1934 exit_flag=3
1935 RETURN
1936 END IF
1937 IF (.not.got_var(idubar).and.(ifile.le.2)) THEN
1938 IF (master) WRITE (stdout,50) trim(vname(1,idubar)), &
1939 & trim(ncname)
1940 exit_flag=3
1941 RETURN
1942 END IF
1943 IF (.not.got_var(idvbar).and.(ifile.le.2)) THEN
1944 IF (master) WRITE (stdout,50) trim(vname(1,idvbar)), &
1945 & trim(ncname)
1946 exit_flag=3
1947 RETURN
1948 END IF
1949# ifdef ADJUST_BOUNDARY
1950 IF (.not.got_var(idsbry(isfsur)).and.(ifile.eq.3)) THEN
1951 IF (master) WRITE (stdout,50) trim(vname(1,idsbry(isfsur))), &
1952 & trim(ncname)
1953 exit_flag=3
1954 RETURN
1955 END IF
1956 IF (.not.got_var(idsbry(isubar)).and.(ifile.eq.3)) THEN
1957 IF (master) WRITE (stdout,50) trim(vname(1,idsbry(isubar))), &
1958 & trim(ncname)
1959 exit_flag=3
1960 RETURN
1961 END IF
1962 IF (.not.got_var(idsbry(isvbar)).and.(ifile.eq.3)) THEN
1963 IF (master) WRITE (stdout,50) trim(vname(1,idsbry(isvbar))), &
1964 & trim(ncname)
1965 exit_flag=3
1966 RETURN
1967 END IF
1968# endif
1969# ifdef SOLVE3D
1970 IF (.not.got_var(iduvel).and.(ifile.le.2)) THEN
1971 IF (master) WRITE (stdout,50) trim(vname(1,iduvel)), &
1972 & trim(ncname)
1973 exit_flag=3
1974 RETURN
1975 END IF
1976 IF (.not.got_var(idvvel).and.(ifile.le.2)) THEN
1977 IF (master) WRITE (stdout,50) trim(vname(1,idvvel)), &
1978 & trim(ncname)
1979 exit_flag=3
1980 RETURN
1981 END IF
1982# ifdef ADJUST_BOUNDARY
1983 IF (.not.got_var(idsbry(isuvel)).and.(ifile.eq.3)) THEN
1984 IF (master) WRITE (stdout,50) trim(vname(1,idsvar(isuvel))), &
1985 & trim(ncname)
1986 exit_flag=3
1987 RETURN
1988 END IF
1989 IF (.not.got_var(idsbry(isvvel)).and.(ifile.eq.3)) THEN
1990 IF (master) WRITE (stdout,50) trim(vname(1,idsvar(isvvel))), &
1991 & trim(ncname)
1992 exit_flag=3
1993 RETURN
1994 END IF
1995# endif
1996# endif
1997# ifdef ADJUST_WSTRESS
1998 IF (.not.got_var(idusms).and.(ifile.eq.4)) THEN
1999 IF (master) WRITE (stdout,50) trim(vname(1,idusms)), &
2000 & trim(ncname)
2001 exit_flag=3
2002 RETURN
2003 END IF
2004 IF (.not.got_var(idvsms).and.(ifile.eq.4)) THEN
2005 IF (master) WRITE (stdout,50) trim(vname(1,idvsms)), &
2006 & trim(ncname)
2007 exit_flag=3
2008 RETURN
2009 END IF
2010# endif
2011# ifdef SOLVE3D
2012 DO itrc=1,nt(ng)
2013 IF (.not.got_var(idtvar(itrc)).and.(ifile.le.2)) THEN
2014 IF (master) WRITE (stdout,50) trim(vname(1,idtvar(itrc))), &
2015 & trim(ncname)
2016 exit_flag=3
2017 RETURN
2018 END IF
2019# ifdef ADJUST_BOUNDARY
2020 IF (.not.got_var(idsbry(istvar(itrc))).and.(ifile.eq.3)) THEN
2021 IF (master) WRITE (stdout,50) &
2022 & trim(vname(1,idsbry(istvar(itrc)))), &
2023 & trim(ncname)
2024 exit_flag=3
2025 RETURN
2026 END IF
2027# endif
2028# ifdef ADJUST_STFLUX
2029 IF (.not.got_var(idtsur(itrc)).and.(ifile.eq.4).and. &
2030 & lstflux(itrc,ng)) THEN
2031 IF (master) WRITE (stdout,50) trim(vname(1,idtsur(itrc))), &
2032 & trim(ncname)
2033 exit_flag=3
2034 RETURN
2035 END IF
2036# endif
2037 END DO
2038# endif
2039!
2040! Set unlimited time record dimension to zero to allow other programs
2041! to process and write normalization factors for different variable.
2042!
2043 nrm(ifile,ng)%Rindex=0
2044 END IF query
2045!
2046 10 FORMAT (/,' DEF_NORM_PIO - unable to create norm NetCDF', &
2047 & ' file: ',a)
2048 20 FORMAT (a,', ',a)
2049 30 FORMAT (1pe11.4,1x,'millimeter')
2050 40 FORMAT (/,' DEF_NORM_PIO - unable to open norm NetCDF file: ',a)
2051 50 FORMAT (/,' DEF_NORM_PIO - unable to find variable: ',a,2x, &
2052 & ' in norm NetCDF file: ',a)
2053!
2054 RETURN
2055 END SUBROUTINE def_norm_pio
2056# endif
2057#endif
2058 END MODULE def_norm_mod
def_dim_mod::def_dim
Definition def_dim.F:24
def_info_mod::def_info
Definition def_info.F:72
def_var_mod::def_var
Definition def_var.F:88
wrt_info_mod::wrt_info
Definition wrt_info.F:54
def_dim_mod
Definition def_dim.F:2
def_info_mod
Definition def_info.F:2
def_norm_mod
Definition def_norm.F:2
def_norm_mod::def_norm_pio
subroutine, private def_norm_pio(ng, model, ifile)
Definition def_norm.F:1030
def_norm_mod::def_norm_nf90
subroutine, private def_norm_nf90(ng, model, ifile)
Definition def_norm.F:90
def_norm_mod::def_norm
subroutine, public def_norm(ng, model, ifile)
Definition def_norm.F:53
def_var_mod
Definition def_var.F:2
mod_biology
Definition ecosim_mod.h:1
mod_biology::nfec
integer, parameter nfec
Definition ecosim_mod.h:204
mod_biology::nbac
integer, parameter nbac
Definition ecosim_mod.h:202
mod_biology::ndom
integer, parameter ndom
Definition ecosim_mod.h:203
mod_biology::nbands
integer, parameter nbands
Definition ecosim_mod.h:201
mod_biology::nphy
integer, parameter nphy
Definition ecosim_mod.h:205
mod_fourdvar
Definition mod_fourdvar.F:2
mod_fourdvar::nstatevar
integer, dimension(:), allocatable nstatevar
Definition mod_fourdvar.F:362
mod_iounits
Definition mod_iounits.F:2
mod_iounits::nrm
type(t_io), dimension(:,:), allocatable nrm
Definition mod_iounits.F:230
mod_iounits::stdout
integer stdout
Definition mod_iounits.F:195
mod_iounits::sourcefile
character(len=256) sourcefile
Definition mod_iounits.F:317
mod_ncparam
Definition mod_ncparam.F:2
mod_ncparam::io_nf90
integer, parameter io_nf90
Definition mod_ncparam.F:95
mod_ncparam::idubar
integer idubar
Definition mod_ncparam.F:325
mod_ncparam::idvvel
integer idvvel
Definition mod_ncparam.F:381
mod_ncparam::idvsms
integer idvsms
Definition mod_ncparam.F:377
mod_ncparam::nv
integer, parameter nv
Definition mod_ncparam.F:105
mod_ncparam::isvvel
integer isvvel
Definition mod_ncparam.F:507
mod_ncparam::idsbry
integer, dimension(:), allocatable idsbry
Definition mod_ncparam.F:523
mod_ncparam::io_pio
integer, parameter io_pio
Definition mod_ncparam.F:96
mod_ncparam::isvbar
integer isvbar
Definition mod_ncparam.F:505
mod_ncparam::idtsur
integer, dimension(:), allocatable idtsur
Definition mod_ncparam.F:453
mod_ncparam::idfsur
integer idfsur
Definition mod_ncparam.F:206
mod_ncparam::idtvar
integer, dimension(:), allocatable idtvar
Definition mod_ncparam.F:568
mod_ncparam::istvar
integer, dimension(:), allocatable istvar
Definition mod_ncparam.F:520
mod_ncparam::isuvel
integer isuvel
Definition mod_ncparam.F:506
mod_ncparam::isfsur
integer isfsur
Definition mod_ncparam.F:503
mod_ncparam::iduvel
integer iduvel
Definition mod_ncparam.F:343
mod_ncparam::vname
character(len=maxlen), dimension(6, 0:nv) vname
Definition mod_ncparam.F:591
mod_ncparam::idsvar
integer, dimension(:), allocatable idsvar
Definition mod_ncparam.F:522
mod_ncparam::idtime
integer idtime
Definition mod_ncparam.F:307
mod_ncparam::isubar
integer isubar
Definition mod_ncparam.F:504
mod_ncparam::iinfo
integer, dimension(:,:,:), allocatable iinfo
Definition mod_ncparam.F:578
mod_ncparam::idusms
integer idusms
Definition mod_ncparam.F:339
mod_ncparam::ndimid
integer, parameter ndimid
Definition mod_ncparam.F:110
mod_ncparam::idvbar
integer idvbar
Definition mod_ncparam.F:360
mod_netcdf
Definition mod_netcdf.F:2
mod_netcdf::nf_tout
integer, parameter nf_tout
Definition mod_netcdf.F:207
mod_netcdf::netcdf_check_dim
subroutine, public netcdf_check_dim(ng, model, ncname, ncid)
Definition mod_netcdf.F:538
mod_netcdf::netcdf_open
subroutine, public netcdf_open(ng, model, ncname, omode, ncid)
Definition mod_netcdf.F:9293
mod_netcdf::nf_fout
integer, parameter nf_fout
Definition mod_netcdf.F:188
mod_netcdf::netcdf_enddef
subroutine, public netcdf_enddef(ng, model, ncname, ncid)
Definition mod_netcdf.F:9220
mod_netcdf::var_name
character(len=100), dimension(mvars) var_name
Definition mod_netcdf.F:169
mod_netcdf::var_id
integer, dimension(mvars) var_id
Definition mod_netcdf.F:160
mod_netcdf::n_var
integer n_var
Definition mod_netcdf.F:152
mod_netcdf::netcdf_create
subroutine, public netcdf_create(ng, model, ncname, ncid)
Definition mod_netcdf.F:9089
mod_netcdf::netcdf_inq_var
subroutine, public netcdf_inq_var(ng, model, ncname, ncid, myvarname, searchvar, varid, nvardim, nvaratt)
Definition mod_netcdf.F:1371
mod_parallel
Definition mod_parallel.F:2
mod_parallel::master
logical master
Definition mod_parallel.F:40
mod_param
Definition mod_param.F:2
mod_param::u2dobc
integer, parameter u2dobc
Definition mod_param.F:729
mod_param::v3dobc
integer, parameter v3dobc
Definition mod_param.F:733
mod_param::inlm
integer, parameter inlm
Definition mod_param.F:662
mod_param::nbed
integer nbed
Definition mod_param.F:517
mod_param::r2dobc
integer, parameter r2dobc
Definition mod_param.F:728
mod_param::n
integer, dimension(:), allocatable n
Definition mod_param.F:479
mod_param::r3dvar
integer, parameter r3dvar
Definition mod_param.F:721
mod_param::r3dobc
integer, parameter r3dobc
Definition mod_param.F:731
mod_param::v2dobc
integer, parameter v2dobc
Definition mod_param.F:730
mod_param::iobounds
type(t_iobounds), dimension(:), allocatable iobounds
Definition mod_param.F:282
mod_param::u3dobc
integer, parameter u3dobc
Definition mod_param.F:732
mod_param::u3dvar
integer, parameter u3dvar
Definition mod_param.F:722
mod_param::u2dvar
integer, parameter u2dvar
Definition mod_param.F:718
mod_param::itlm
integer, parameter itlm
Definition mod_param.F:663
mod_param::nt
integer, dimension(:), allocatable nt
Definition mod_param.F:489
mod_param::r2dvar
integer, parameter r2dvar
Definition mod_param.F:717
mod_param::v2dvar
integer, parameter v2dvar
Definition mod_param.F:719
mod_param::nst
integer nst
Definition mod_param.F:521
mod_param::v3dvar
integer, parameter v3dvar
Definition mod_param.F:723
mod_pio_netcdf
Definition mod_pio_netcdf.F:2
mod_pio_netcdf::pio_fout
integer, parameter pio_fout
Definition mod_pio_netcdf.F:569
mod_pio_netcdf::var_desc
type(var_desc_t), dimension(:), pointer var_desc
Definition mod_pio_netcdf.F:437
mod_pio_netcdf::pio_netcdf_create
subroutine, public pio_netcdf_create(ng, model, ncname, piofile)
Definition mod_pio_netcdf.F:8997
mod_pio_netcdf::pio_netcdf_inq_var
subroutine, public pio_netcdf_inq_var(ng, model, ncname, piofile, myvarname, searchvar, piovar, nvardim, nvaratt)
Definition mod_pio_netcdf.F:1923
mod_pio_netcdf::pio_netcdf_open
subroutine, public pio_netcdf_open(ng, model, ncname, omode, piofile)
Definition mod_pio_netcdf.F:9148
mod_pio_netcdf::pio_netcdf_check_dim
subroutine, public pio_netcdf_check_dim(ng, model, ncname, piofile)
Definition mod_pio_netcdf.F:1080
mod_pio_netcdf::pio_tout
integer, parameter pio_tout
Definition mod_pio_netcdf.F:588
mod_pio_netcdf::pio_netcdf_enddef
subroutine, public pio_netcdf_enddef(ng, model, ncname, piofile)
Definition mod_pio_netcdf.F:9096
mod_scalars
Definition mod_scalars.F:2
mod_scalars::lobc
logical, dimension(:,:,:), allocatable lobc
Definition mod_scalars.F:1127
mod_scalars::lstflux
logical, dimension(:,:), allocatable lstflux
Definition mod_scalars.F:1156
mod_scalars::rclock
type(t_clock) rclock
Definition mod_scalars.F:120
mod_scalars::ldefnrm
logical, dimension(:,:), allocatable ldefnrm
Definition mod_scalars.F:1007
mod_scalars::exit_flag
integer exit_flag
Definition mod_scalars.F:559
mod_scalars::noerror
integer noerror
Definition mod_scalars.F:561
mod_sediment
Definition sediment_mod.h:1
mod_sediment::idsed
integer, dimension(:), allocatable idsed
Definition sediment_mod.h:88
mod_sediment::sd50
real(r8), dimension(:,:), allocatable sd50
Definition sediment_mod.h:173
strings_mod
Definition strings.F:2
strings_mod::founderror
logical function, public founderror(flag, noerr, line, routine)
Definition strings.F:52
wrt_info_mod
Definition wrt_info.F:2