New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
trcdia.F90 in trunk/NEMO/TOP_SRC – NEMO

source: trunk/NEMO/TOP_SRC/trcdia.F90 @ 1353

Last change on this file since 1353 was 1353, checked in by cetlod, 15 years ago

correction of calendar in TOP outputs, see ticket:368
  • Property svn:eol-style set to native
  • Property svn:keywords set to Id
File size: 29.5 KB
Line 
1MODULE trcdia
2   !!======================================================================
3   !!                       *** MODULE trcdia ***
4   !! TOP :   Output of passive tracers
5   !!======================================================================
6   !! History :    -   !  1995-01 (M. Levy)  Original code
7   !!              -   !  1998-01 (C. Levy) NETCDF format using ioipsl interface
8   !!              -   !  1999-01 (M.A. Foujols) adapted for passive tracer
9   !!              -   !  1999-09 (M.A. Foujols) split into three parts
10   !!             1.0  !  2005-03 (O. Aumont, A. El Moussaoui) F90
11   !!                  !  2008-05 (C. Ethe re-organization)
12   !!----------------------------------------------------------------------
13#if defined key_top
14   !!----------------------------------------------------------------------
15   !!   'key_top'                                                TOP models
16   !!----------------------------------------------------------------------
17   !! trc_dia     : main routine of output passive tracer
18   !! trcdit_wr   : outputs of concentration fields
19   !! trcdid_wr   : outputs of dvection-diffusion trends
20   !! trcdii_wr   : outputs of additional 2D/3D diagnostics
21   !! trcdib_wr   : outputs of biological fields
22   !!----------------------------------------------------------------------
23   USE oce_trc
24   USE trc
25   USE trp_trc
26   USE trdmld_trc_oce, ONLY : luttrd
27   USE dianam    ! build name of file (routine)
28   USE in_out_manager  ! I/O manager
29   USE lib_mpp
30   USE ioipsl
31
32   IMPLICIT NONE
33   PRIVATE
34
35   PUBLIC trc_dia     
36
37   INTEGER  ::   nit5      !: id for tracer output file
38   INTEGER  ::   ndepit5   !: id for depth mesh
39   INTEGER  ::   nhorit5   !: id for horizontal mesh
40   INTEGER  ::   ndimt50   !: number of ocean points in index array
41   INTEGER  ::   ndimt51   !: number of ocean points in index array
42   REAL(wp) ::   xjulian   !: ????   not DOCTOR !
43   INTEGER , DIMENSION (jpij*jpk) ::   ndext50   !: integer arrays for ocean 3D index
44   INTEGER , DIMENSION (jpij)     ::   ndext51   !: integer arrays for ocean surface index
45# if defined key_trc_diaadd
46   INTEGER  ::   nitd      !: id for additional array output file
47   INTEGER  ::   ndepitd   !: id for depth mesh
48   INTEGER  ::   nhoritd   !: id for horizontal mesh
49# endif
50# if defined key_trc_diatrd
51   INTEGER , DIMENSION (jptra) ::   nit6      !: id for additional array output file
52   INTEGER , DIMENSION (jptra) ::   ndepit6   !: id for depth mesh
53   INTEGER , DIMENSION (jptra) ::   nhorit6   !: id for horizontal mesh
54# endif
55# if defined key_trc_diabio
56   INTEGER  ::   nitb        !:         id.         for additional array output file
57   INTEGER  ::   ndepitb   !:  id for depth mesh
58   INTEGER  ::   nhoritb   !:  id for horizontal mesh
59# endif
60
61   !! * Substitutions
62#  include "top_substitute.h90"
63   !!----------------------------------------------------------------------
64   !! NEMO/TOP 1.0 , LOCEAN-IPSL (2005)
65   !! $Id$
66   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
67   !!----------------------------------------------------------------------
68
69CONTAINS
70
71   SUBROUTINE trc_dia( kt, kindic ) 
72      !!---------------------------------------------------------------------
73      !!                     ***  ROUTINE trc_dia  ***
74      !!
75      !! ** Purpose :   output passive tracers fields
76      !!---------------------------------------------------------------------
77      INTEGER, INTENT( in ) :: kt, kindic
78      !!---------------------------------------------------------------------
79     
80      CALL trcdit_wr( kt, kindic )      ! outputs for tracer concentration
81      CALL trcdid_wr( kt, kindic )      ! outputs for dynamical trends
82      CALL trcdii_wr( kt, kindic )      ! outputs for additional arrays
83      CALL trcdib_wr( kt, kindic )      ! outputs for biological trends
84
85      !
86   END SUBROUTINE trc_dia
87
88   SUBROUTINE trcdit_wr( kt, kindic )
89      !!----------------------------------------------------------------------
90      !!                     ***  ROUTINE trcdit_wr  ***
91      !!
92      !! ** Purpose :   Standard output of passive tracer : concentration fields
93      !!
94      !! ** Method  :   At the beginning of the first time step (nit000), define all
95      !!             the NETCDF files and fields for concentration of passive tracer
96      !!
97      !!        At each time step call histdef to compute the mean if necessary
98      !!        Each nwritetrc time step, output the instantaneous or mean fields
99      !!
100      !!        IF kindic <0, output of fields before the model interruption.
101      !!        IF kindic =0, time step loop
102      !!        IF kindic >0, output of fields before the time step loop
103      !!----------------------------------------------------------------------
104      INTEGER, INTENT( in ) ::   kt          ! ocean time-step
105      INTEGER, INTENT( in ) ::   kindic      ! indicator of abnormal termination
106      !!
107      INTEGER ::   jn
108      LOGICAL ::   ll_print = .FALSE.
109      CHARACTER (len=40) :: clhstnam, clop
110      CHARACTER (len=20) :: cltra, cltrau
111      CHARACTER (len=80) :: cltral
112      REAL(wp) :: zsto, zout, zdt
113      INTEGER  :: iimi, iima, ijmi, ijma, ipk, it, itmod
114      !!----------------------------------------------------------------------
115
116      ! Initialisation
117      ! --------------
118
119      ! local variable for debugging
120      ll_print = .FALSE.                  ! change it to true for more control print
121      ll_print = ll_print .AND. lwp
122
123      ! Define frequency of output and means
124      zdt = rdt
125      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!)
126      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time)
127      ENDIF
128# if defined key_diainstant
129      zsto = nwritetrc * rdt
130      clop = "inst("//TRIM(clop)//")"
131# else
132      zsto = zdt
133      clop = "ave("//TRIM(clop)//")"
134# endif
135      zout = nwritetrc * zdt
136
137      ! Define indices of the horizontal output zoom and vertical limit storage
138      iimi = 1      ;      iima = jpi
139      ijmi = 1      ;      ijma = jpj
140      ipk = jpk
141
142      ! define time axis
143      itmod = kt - nittrc000 + 1
144      it    = kt
145
146      ! Define NETCDF files and fields at beginning of first time step
147      ! --------------------------------------------------------------
148
149      IF(ll_print)WRITE(numout,*)'trcdit_wr kt=',kt,' kindic ',kindic
150     
151      IF( kt == nittrc000 ) THEN
152
153         ! Compute julian date from starting date of the run
154         CALL ymds2ju( nyear, nmonth, nday, rdt, xjulian )
155         xjulian = xjulian - adatrj   !   set calendar origin to the beginning of the experiment
156         IF(lwp)WRITE(numout,*)' ' 
157         IF(lwp)WRITE(numout,*)' Date 0 used :', nittrc000                         &
158            &                 ,' YEAR ', nyear, ' MONTH ', nmonth, ' DAY ', nday   &
159            &                 ,'Julian day : ', xjulian 
160 
161         IF(lwp) WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma,  &
162            &                    ' limit storage in depth = ', ipk
163
164
165         ! Define the NETCDF files for passive tracer concentration
166         CALL dia_nam( clhstnam, nwritetrc, 'ptrc_T' )
167         IF(lwp)WRITE(numout,*)" Name of NETCDF file ", clhstnam
168
169         ! Horizontal grid : glamt and gphit
170         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,     &
171            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,         & 
172            &          nittrc000-ndttrc, xjulian, zdt, nhorit5, nit5 , domain_id=nidom)
173
174         ! Vertical grid for tracer : gdept
175         CALL histvert( nit5, 'deptht', 'Vertical T levels', 'm', ipk, gdept_0, ndepit5)
176
177         ! Index of ocean points in 3D and 2D (surface)
178         CALL wheneq( jpi*jpj*ipk, tmask, 1, 1., ndext50, ndimt50 )
179         CALL wheneq( jpi*jpj    , tmask, 1, 1., ndext51, ndimt51 )
180
181         ! Declare all the output fields as NETCDF variables
182         DO jn = 1, jptra
183            IF( lutsav(jn) ) THEN
184               cltra  = ctrcnm(jn)   ! short title for tracer
185               cltral = ctrcnl(jn)   ! long title for tracer
186               cltrau = ctrcun(jn)   ! UNIT for tracer
187               CALL histdef( nit5, cltra, cltral, cltrau, jpi, jpj, nhorit5,  &
188                  &          ipk, 1, ipk,  ndepit5, 32, clop, zsto, zout ) 
189            ENDIF
190         END DO
191
192         ! end netcdf files header
193         CALL histend( nit5 )
194         IF(lwp) WRITE(numout,*)
195         IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization in trcdit_wr'
196         IF( ll_print )   CALL FLUSH(numout )
197
198      ENDIF
199
200      ! Start writing the tracer concentrations
201      ! ---------------------------------------
202
203      IF( lwp .AND. MOD( itmod, nwritetrc ) == 0 ) THEN
204         WRITE(numout,*) 'trcdit_wr : write NetCDF passive tracer concentrations at ', kt, 'time-step'
205         WRITE(numout,*) '~~~~~~~~~ '
206      ENDIF
207
208      DO jn = 1, jptra
209         IF( lutsav(jn) ) THEN
210            cltra = ctrcnm(jn)      ! short title for tracer
211            CALL histwrite( nit5, cltra, it, trn(:,:,:,jn), ndimt50, ndext50 )
212         ENDIF
213      END DO
214
215      ! close the file
216      ! --------------
217      IF( kt == nitend .OR. kindic < 0 )   CALL histclo( nit5 )
218      !
219   END SUBROUTINE trcdit_wr
220
221# if defined key_trc_diatrd
222
223   SUBROUTINE trcdid_wr( kt, kindic )
224      !!----------------------------------------------------------------------
225      !!                     ***  ROUTINE trcdid_wr  ***
226      !!
227      !! ** Purpose :   output of passive tracer : advection-diffusion trends
228      !!
229      !! ** Method  :   At the beginning of the first time step (nit000), define all
230      !!             the NETCDF files and fields for concentration of passive tracer
231      !!
232      !!        At each time step call histdef to compute the mean if necessary
233      !!        Each nwritetrc time step, output the instantaneous or mean fields
234      !!
235      !!        IF kindic <0, output of fields before the model interruption.
236      !!        IF kindic =0, time step loop
237      !!        IF kindic >0, output of fields before the time step loop
238      !!----------------------------------------------------------------------
239      INTEGER, INTENT( in ) ::   kt          ! ocean time-step
240      INTEGER, INTENT( in ) ::   kindic      ! indicator of abnormal termination
241      !!
242      LOGICAL ::   ll_print = .FALSE.
243      CHARACTER (len=40) ::   clhstnam, clop
244      CHARACTER (len=20) ::   cltra, cltrau
245      CHARACTER (len=80) ::   cltral
246      CHARACTER (len=10) ::   csuff
247      INTEGER  ::   jn, jl
248      INTEGER  ::   iimi, iima, ijmi, ijma, ipk, it, itmod
249      REAL(wp) ::   zsto, zout, zdt
250      !!----------------------------------------------------------------------
251
252      ! 0. Initialisation
253      ! -----------------
254
255      ! local variable for debugging
256      ll_print = .FALSE.
257      ll_print = ll_print .AND. lwp
258      !
259      ! Define frequency of output and means
260      zdt = rdt
261      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!)
262      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time)
263      ENDIF
264#  if defined key_diainstant
265      zsto = nwritetrd * rdt
266      clop = "inst("//TRIM(clop)//")"
267#  else
268      zsto = zdt
269      clop = "ave("//TRIM(clop)//")"
270#  endif
271      zout = nwritetrd * zdt
272
273      ! Define indices of the horizontal output zoom and vertical limit storage
274      iimi = 1      ;      iima = jpi
275      ijmi = 1      ;      ijma = jpj
276      ipk = jpk
277
278      ! define time axis
279      itmod = kt - nittrc000 + 1
280      it    = kt
281
282      ! Define the NETCDF files (one per tracer)
283      IF( ll_print ) WRITE(numout,*) 'trcdid kt=', kt, ' kindic ', kindic
284     
285     
286      IF( kt == nittrc000 ) THEN
287
288         DO jn = 1, jptra
289            !
290            IF( luttrd(jn) ) THEN      ! Define the file for dynamical trends - one per each tracer IF required
291
292               IF(lwp)WRITE(numout,*) ' indexes of zoom = ', iimi, iima, ijmi, ijma,  &
293                   &                   ' limit storage in depth = ', ipk
294               csuff='DY_'//ctrcnm(jn)
295               CALL dia_nam( clhstnam, nwritetrd, csuff )
296               IF(lwp)WRITE(numout,*)   " Name of NETCDF file for dynamical trends",   &
297                  &                     " of tracer number : ",clhstnam
298
299               CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,       &
300                  &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,   &
301                  &          nittrc000-ndttrc, xjulian, zdt, nhorit6(jn),  &
302                  &          nit6(jn) , domain_id=nidom )
303
304               ! Vertical grid for tracer trend - one per each tracer IF needed
305               CALL histvert( nit6(jn), 'deptht', 'Vertical T levels', 'm', ipk, gdept_0, ndepit6(jn) ) 
306             END IF
307          END DO
308
309          ! Declare all the output fields as NETCDF variables
310
311          ! trends for tracer concentrations
312          DO jn = 1, jptra
313            IF( luttrd(jn) ) THEN
314                DO jl = 1, jpdiatrc
315                  IF( jl == 1 ) THEN
316                      ! short and long title for x advection for tracer
317                      WRITE (cltra,'("XAD_",16a)') ctrcnm(jn)
318                      WRITE (cltral,'("X advective trend for ",58a)')  &
319                         &      ctrcnl(jn)(1:58)
320                  END IF
321                  IF( jl == 2 ) THEN
322                      ! short and long title for y advection for tracer
323                      WRITE (cltra,'("YAD_",16a)') ctrcnm(jn)
324                      WRITE (cltral,'("Y advective trend for ",58a)')  &
325                         &      ctrcnl(jn)(1:58)
326                  END IF
327                  IF( jl == 3 ) THEN
328                      ! short and long title for Z advection for tracer
329                      WRITE (cltra,'("ZAD_",16a)') ctrcnm(jn)
330                      WRITE (cltral,'("Z advective trend for ",58a)')  &
331                         &      ctrcnl(jn)(1:58)
332                  END IF
333                  IF( jl == 4 ) THEN
334                      ! short and long title for X diffusion for tracer
335                      WRITE (cltra,'("XDF_",16a)') ctrcnm(jn)
336                      WRITE (cltral,'("X diffusion trend for ",58a)')  &
337                         &      ctrcnl(jn)(1:58)
338                  END IF
339                  IF( jl == 5 ) THEN
340                      ! short and long title for Y diffusion for tracer
341                      WRITE (cltra,'("YDF_",16a)') ctrcnm(jn)
342                      WRITE (cltral,'("Y diffusion trend for ",58a)')  &
343                         &      ctrcnl(jn)(1:58)
344                  END IF
345                  IF( jl == 6 ) THEN
346                      ! short and long title for Z diffusion for tracer
347                      WRITE (cltra,'("ZDF_",16a)') ctrcnm(jn)
348                      WRITE (cltral,'("Z diffusion trend for ",58a)')  &
349                         &      ctrcnl(jn)(1:58)
350                  END IF
351# if defined key_trcldf_eiv
352                  IF( jl == 7 ) THEN
353                      ! short and long title for x gent velocity for tracer
354                      WRITE (cltra,'("XGV_",16a)') ctrcnm(jn)
355                      WRITE (cltral,'("X gent velocity trend for ",53a)')  &
356                         &      ctrcnl(jn)(1:53)
357                  END IF
358                  IF( jl == 8 ) THEN
359                      ! short and long title for y gent velocity for tracer
360                      WRITE (cltra,'("YGV_",16a)') ctrcnm(jn)
361                      WRITE (cltral,'("Y gent velocity trend for ",53a)')  &
362                         &      ctrcnl(jn)(1:53)
363                  END IF
364                  IF( jl == 9 ) THEN
365                      ! short and long title for Z gent velocity for tracer
366                      WRITE (cltra,'("ZGV_",16a)') ctrcnm(jn)
367                      WRITE (cltral,'("Z gent velocity trend for ",53a)')  &
368                         &      ctrcnl(jn)(1:53)
369                  END IF
370# endif
371# if defined key_trcdmp
372                  IF( jl == jpdiatrc - 1 ) THEN
373                      ! last trends for tracer damping : short and long title
374                      WRITE (cltra,'("TDM_",16a)') ctrcnm(jn)
375                      WRITE (cltral,'("Tracer damping trend for ",55a)')  &
376                         &      ctrcnl(jn)(1:55)
377                  END IF
378# endif
379                  IF( jl == jpdiatrc ) THEN
380                      ! last trends for tracer damping : short and long title
381                      WRITE (cltra,'("SBC_",16a)') ctrcnm(jn)
382                      WRITE (cltral,'("Surface boundary flux ",58a)')  &
383                      &      ctrcnl(jn)(1:58)
384                  END IF
385
386                  CALL FLUSH( numout )
387                  cltrau = ctrcun(jn)      ! UNIT for tracer /trends
388                  CALL histdef( nit6(jn), cltra, cltral, cltrau, jpi,jpj,  &
389                     &          nhorit6(jn), ipk, 1, ipk,  ndepit6(jn), 32, clop ,  &
390                     &          zsto,zout )
391               END DO
392            END IF
393         END DO
394
395         ! CLOSE netcdf Files
396          DO jn = 1, jptra
397             IF( luttrd(jn) )   CALL histend( nit6(jn) )
398          END DO
399
400         IF(lwp) WRITE(numout,*)
401         IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization in trcdid'
402         IF(ll_print) CALL FLUSH(numout )
403         !
404      ENDIF
405
406      ! SOME diagnostics to DO first time
407
408      ! Start writing data
409      ! ---------------------
410
411      ! trends for tracer concentrations
412
413      IF( lwp .AND. MOD( itmod, nwritetrd ) == 0 ) THEN
414         WRITE(numout,*) 'trcdid_wr : write NetCDF dynamical trends at ', kt, 'time-step'
415         WRITE(numout,*) '~~~~~~ '
416      ENDIF
417
418      DO jn = 1, jptra
419         IF( luttrd(jn) ) THEN
420            DO jl = 1, jpdiatrc
421               ! short titles
422               IF( jl == 1)   WRITE (cltra,'("XAD_",16a)') ctrcnm(jn)      ! x advection for tracer
423               IF( jl == 2)   WRITE (cltra,'("YAD_",16a)') ctrcnm(jn)      ! z advection for tracer
424               IF( jl == 3)   WRITE (cltra,'("ZAD_",16a)') ctrcnm(jn)      ! z advection for tracer
425               IF( jl == 4)   WRITE (cltra,'("XDF_",16a)') ctrcnm(jn)      ! x diffusion for tracer
426               IF( jl == 5)   WRITE (cltra,'("YDF_",16a)') ctrcnm(jn)      ! y diffusion for tracer
427               IF( jl == 6)   WRITE (cltra,'("ZDF_",16a)') ctrcnm(jn)      ! z diffusion for tracer
428# if defined key_trcldf_eiv
429               IF( jl == 7)   WRITE (cltra,'("XGV_",16a)') ctrcnm(jn)      ! x gent velocity for tracer
430               IF( jl == 8)   WRITE (cltra,'("YGV_",16a)') ctrcnm(jn)      ! y gent velocity for tracer
431               IF( jl == 9)   WRITE (cltra,'("ZGV_",16a)') ctrcnm(jn)      ! z gent velocity for tracer
432# endif
433# if defined key_trcdmp
434               IF( jl == jpdiatrc - 1 )   WRITE (cltra,'("TDM_",16a)') ctrcnm(jn)      ! damping
435# endif
436               IF( jl == jpdiatrc )   WRITE (cltra,'("SBC_",a)') ctrcnm(jn)      ! surface boundary conditions
437               !
438               CALL histwrite(nit6(jn), cltra, it, trtrd(:,:,:,ikeep(jn),jl),ndimt50, ndext50)
439            END DO
440         END IF
441      END DO
442
443      ! Closing all files
444      ! -----------------
445      IF( kt == nitend .OR. kindic < 0 ) THEN
446         DO jn = 1, jptra
447            IF( luttrd(jn) )   CALL histclo( nit6(jn) )
448         END DO
449      ENDIF
450      !
451   END SUBROUTINE trcdid_wr
452
453# else
454
455   SUBROUTINE trcdid_wr( kt, kindic )                      ! Dummy routine
456      INTEGER, INTENT ( in ) ::   kt, kindic
457   END SUBROUTINE trcdid_wr
458
459# endif
460
461#if defined key_trc_diaadd
462
463   SUBROUTINE trcdii_wr( kt, kindic )
464      !!----------------------------------------------------------------------
465      !!                     ***  ROUTINE trcdii_wr  ***
466      !!
467      !! ** Purpose :   output of passive tracer : additional 2D and 3D arrays
468      !!
469      !! ** Method  :   At the beginning of the first time step (nit000), define all
470      !!             the NETCDF files and fields for concentration of passive tracer
471      !!
472      !!        At each time step call histdef to compute the mean if necessary
473      !!        Each nwritetrc time step, output the instantaneous or mean fields
474      !!
475      !!        IF kindic <0, output of fields before the model interruption.
476      !!        IF kindic =0, time step loop
477      !!        IF kindic >0, output of fields before the time step loop
478      !!----------------------------------------------------------------------
479      INTEGER, INTENT( in ) ::   kt          ! ocean time-step
480      INTEGER, INTENT( in ) ::   kindic      ! indicator of abnormal termination
481      !!
482      LOGICAL ::   ll_print = .FALSE.
483      CHARACTER (len=40) ::   clhstnam, clop
484      CHARACTER (len=20) ::   cltra, cltrau
485      CHARACTER (len=80) ::   cltral
486      INTEGER  ::   jn
487      INTEGER  ::   iimi, iima, ijmi, ijma, ipk, it, itmod
488      REAL(wp) ::   zsto, zout, zdt
489      !!----------------------------------------------------------------------
490
491      ! Initialisation
492      ! --------------
493
494      ! local variable for debugging
495      ll_print = .FALSE.
496      ll_print = ll_print .AND. lwp
497      !
498      ! Define frequency of output and means
499      zdt = rdt
500      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!)
501      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time)
502      ENDIF
503#  if defined key_diainstant
504      zsto = nwritedia * zdt
505      clop = "inst("//TRIM(clop)//")"
506#  else
507      zsto = zdt
508      clop = "ave("//TRIM(clop)//")"
509#  endif
510      zout = nwritedia * zdt
511
512      ! Define indices of the horizontal output zoom and vertical limit storage
513      iimi = 1      ;      iima = jpi
514      ijmi = 1      ;      ijma = jpj
515      ipk = jpk
516
517      ! define time axis
518      itmod = kt - nittrc000 + 1
519      it    = kt
520
521      ! 1. Define NETCDF files and fields at beginning of first time step
522      ! -----------------------------------------------------------------
523
524      IF( ll_print ) WRITE(numout,*) 'trcdii_wr kt=', kt, ' kindic ', kindic
525
526      IF( kt == nittrc000 ) THEN
527
528         ! Define the NETCDF files for additional arrays : 2D or 3D
529
530         ! Define the T grid file for tracer auxiliary files
531
532         CALL dia_nam( clhstnam, nwrite, 'diad_T' )
533         IF(lwp) WRITE(numout,*) " Name of NETCDF file ", clhstnam
534
535         ! Define a netcdf FILE for 2d and 3d arrays
536
537         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,             &
538            &          iimi, iima-iimi+1, ijmi, ijma-ijmi+1,         &
539            &          nittrc000-ndttrc, xjulian, zdt, nhoritd, nitd , domain_id=nidom )
540
541         ! Vertical grid for 2d and 3d arrays
542
543         CALL histvert( nitd, 'deptht', 'Vertical T levels','m', ipk, gdept_0, ndepitd)
544
545         ! Declare all the output fields as NETCDF variables
546
547         ! more 3D horizontal arrays
548         DO jn = 1, jpdia3d
549            cltra  = ctrc3d(jn)   ! short title for 3D diagnostic
550            cltral = ctrc3l(jn)   ! long title for 3D diagnostic
551            cltrau = ctrc3u(jn)   ! UNIT for 3D diagnostic
552            CALL histdef( nitd, cltra, cltral, cltrau, jpi, jpj, nhoritd,   &
553               &          ipk, 1, ipk,  ndepitd, 32, clop, zsto, zout )
554         END DO
555
556         ! more 2D horizontal arrays
557         DO jn = 1, jpdia2d
558            cltra  = ctrc2d(jn)    ! short title for 2D diagnostic
559            cltral = ctrc2l(jn)   ! long title for 2D diagnostic
560            cltrau = ctrc2u(jn)   ! UNIT for 2D diagnostic
561            CALL histdef( nitd, cltra, cltral, cltrau, jpi, jpj, nhoritd,  &
562               &          1, 1, 1,  -99, 32, clop, zsto, zout )
563         END DO
564
565         ! TODO: more 2D vertical sections arrays : I or J indice fixed
566
567         ! CLOSE netcdf Files
568         CALL histend( nitd )
569
570         IF(lwp) WRITE(numout,*)
571         IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization in trcdii_wr'
572         IF( ll_print )   CALL FLUSH(numout )
573         !
574      ENDIF
575
576      ! 2. Start writing data
577      ! ---------------------
578
579      IF( lwp .AND. MOD( itmod, nwritedia ) == 0 ) THEN
580         WRITE(numout,*) 'trcdii_wr : write NetCDF additional arrays at ', kt, 'time-step'
581         WRITE(numout,*) '~~~~~~ '
582      ENDIF
583
584      ! more 3D horizontal arrays
585      DO jn = 1, jpdia3d
586         cltra = ctrc3d(jn)   ! short title for 3D diagnostic
587         CALL histwrite( nitd, cltra, it, trc3d(:,:,:,jn), ndimt50 ,ndext50)
588      END DO
589
590      ! more 2D horizontal arrays
591      DO jn = 1, jpdia2d
592         cltra = ctrc2d(jn)   ! short title for 2D diagnostic
593         CALL histwrite(nitd, cltra, it, trc2d(:,:,jn), ndimt51  ,ndext51)
594      END DO
595
596      ! Closing all files
597      ! -----------------
598      IF( kt == nitend .OR. kindic < 0 )   CALL histclo(nitd)
599      !
600   END SUBROUTINE trcdii_wr
601
602# else
603
604   SUBROUTINE trcdii_wr( kt, kindic )                      ! Dummy routine
605      INTEGER, INTENT ( in ) :: kt, kindic
606   END SUBROUTINE trcdii_wr
607
608# endif
609
610# if defined key_trc_diabio
611
612   SUBROUTINE trcdib_wr( kt, kindic )
613      !!----------------------------------------------------------------------
614      !!                     ***  ROUTINE trcdib_wr  ***
615      !!
616      !! ** Purpose :   output of passive tracer : biological fields
617      !!
618      !! ** Method  :   At the beginning of the first time step (nit000), define all
619      !!             the NETCDF files and fields for concentration of passive tracer
620      !!
621      !!        At each time step call histdef to compute the mean if necessary
622      !!        Each nwritetrc time step, output the instantaneous or mean fields
623      !!
624      !!        IF kindic <0, output of fields before the model interruption.
625      !!        IF kindic =0, time step loop
626      !!        IF kindic >0, output of fields before the time step loop
627      !!----------------------------------------------------------------------
628      !!
629      INTEGER, INTENT( in ) ::   kt          ! ocean time-step
630      INTEGER, INTENT( in ) ::   kindic      ! indicator of abnormal termination
631      !!
632      LOGICAL ::   ll_print = .FALSE.
633      CHARACTER (len=40) ::   clhstnam, clop
634      CHARACTER (len=20) ::   cltra, cltrau
635      CHARACTER (len=80) ::   cltral
636      INTEGER  ::   ji, jj, jk, jn
637      INTEGER  ::   iimi, iima, ijmi, ijma, ipk, it, itmod
638      REAL(wp) ::   zsto, zout, zdt
639      !!----------------------------------------------------------------------
640
641      ! Initialisation
642      ! --------------
643
644      ! local variable for debugging
645      ll_print = .FALSE.
646      ll_print = ll_print .AND. lwp
647
648      ! Define frequency of output and means
649      zdt = rdt
650      IF( ln_mskland )   THEN   ;   clop = "only(x)"   ! put 1.e+20 on land (very expensive!!)
651      ELSE                      ;   clop = "x"         ! no use of the mask value (require less cpu time)
652      ENDIF
653#        if defined key_diainstant
654      zsto = nwritebio * zdt
655      clop = "inst("//TRIM(clop)//")"
656#        else
657      zsto = zdt
658      clop = "ave("//TRIM(clop)//")"
659#        endif
660      zout = nwritebio * zdt
661
662      ! Define indices of the horizontal output zoom and vertical limit storage
663      iimi = 1      ;      iima = jpi
664      ijmi = 1      ;      ijma = jpj
665      ipk = jpk
666
667      ! define time axis
668      itmod = kt - nittrc000 + 1
669      it    = kt
670
671      ! Define NETCDF files and fields at beginning of first time step
672      ! --------------------------------------------------------------
673
674      IF(ll_print) WRITE(numout,*)'trcdib_wr kt=',kt,' kindic ',kindic
675
676      IF( kt == nittrc000 ) THEN
677
678         ! Define the NETCDF files for biological trends
679
680         CALL dia_nam(clhstnam,nwrite,'biolog')
681         IF(lwp)WRITE(numout,*) " Name of NETCDF file for biological trends ", clhstnam
682         ! Horizontal grid : glamt and gphit
683         CALL histbeg( clhstnam, jpi, glamt, jpj, gphit,      &
684            &    iimi, iima-iimi+1, ijmi, ijma-ijmi+1,          &
685            &    nittrc000-ndttrc, xjulian, zdt, nhoritb, nitb , domain_id=nidom )
686         ! Vertical grid for biological trends
687         CALL histvert(nitb, 'deptht', 'Vertical T levels', 'm', ipk, gdept_0, ndepitb)
688
689         ! Declare all the output fields as NETCDF variables
690         ! biological trends
691         DO jn = 1, jpdiabio
692            cltra  = ctrbio(jn)   ! short title for biological diagnostic
693            cltral = ctrbil(jn)   ! long title for biological diagnostic
694            cltrau = ctrbiu(jn)   ! UNIT for biological diagnostic
695            CALL histdef( nitb, cltra, cltral, cltrau, jpi, jpj, nhoritb,  &
696               &         ipk, 1, ipk,  ndepitb, 32, clop, zsto, zout)
697         END DO
698
699         ! CLOSE netcdf Files
700          CALL histend( nitb )
701
702         IF(lwp) WRITE(numout,*)
703         IF(lwp) WRITE(numout,*) 'End of NetCDF Initialization in trcdib_wr'
704         IF(ll_print) CALL FLUSH(numout )
705         !
706      ENDIF
707
708      ! Start writing data
709      ! ------------------
710
711      ! biological trends
712      IF( lwp .AND. MOD( itmod, nwritebio ) == 0 ) THEN
713         WRITE(numout,*) 'trcdit_wr : write NetCDF biological trends at ', kt, 'time-step'
714         WRITE(numout,*) '~~~~~~ '
715      ENDIF
716
717      DO jn = 1, jpdiabio
718         cltra = ctrbio(jn)  ! short title for biological diagnostic
719         CALL histwrite(nitb, cltra, it, trbio(:,:,:,jn), ndimt50,ndext50)
720      END DO
721
722      ! Closing all files
723      ! -----------------
724      IF( kt == nitend .OR. kindic < 0 )   CALL histclo( nitb )
725      !
726   END SUBROUTINE trcdib_wr
727
728# else
729
730   SUBROUTINE trcdib_wr( kt, kindic )                      ! Dummy routine
731      INTEGER, INTENT ( in ) ::   kt, kindic
732   END SUBROUTINE trcdib_wr
733
734# endif 
735
736#else
737   !!----------------------------------------------------------------------
738   !!  Dummy module :                                     No passive tracer
739   !!----------------------------------------------------------------------
740CONTAINS
741   SUBROUTINE trc_dia                      ! Empty routine   
742   END SUBROUTINE trc_dia   
743
744#endif
745
746   !!======================================================================
747END MODULE trcdia
Note: See TracBrowser for help on using the repository browser.