Changeset 7701 for branches/UKMO/dev_r5518_MEDUSA_optim_RH/NEMOGCM/NEMO/TOP_SRC/trcstp.F90

Timestamp:

2017-02-20T11:07:50+01:00 (7 years ago)

Author:

frrh

Message:

Improve control of numstr unit.

In addition to only producing output, and doing the
related global sums when we really need them, we
need to restrict it to one instance on the master PE
in all circumstances and to explicitly close it at the
end of the run. (Currently if lwp = true you get a separate
file for every PE containing identical information and none
of the tracer.stat files are explicitly closed.)

File:

• branches/UKMO/dev_r5518_MEDUSA_optim_RH/NEMOGCM/NEMO/TOP_SRC/trcstp.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_MEDUSA_optim_RH/NEMOGCM/NEMO/TOP_SRC/trcstp.F90

@@ -55 +55 @@
       !!              Update the passive tracers
       !!-------------------------------------------------------------------
+
+      USE dom_oce, ONLY: narea 
+
       INTEGER, INTENT( in ) ::  kt      ! ocean time-step index
       INTEGER               ::  jk, jn  ! dummy loop indices
@@ -136 +139 @@
          ! Hence we restrict its use by protecting it with the ln_ctl RTL
          ! which should normally only be used under debugging conditions
-         ! and not in operational runs. 
-         ztrai = 0._wp                                                   !  content of all tracers
+         ! and not in operational runs. We also need to restrict output 
+         ! to the master PE since there's no point duplicating the same results
+         ! on all processors.    
+         ztrai = 0._wp                                !  content of all tracers
          DO jn = 1, jptra
             ztrai = ztrai + glob_sum( trn(:,:,:,jn) * cvol(:,:,:)   )
          END DO
-         IF( lwp ) WRITE(numstr,9300) kt,  ztrai / areatot
+         IF (narea == 1) WRITE(numstr,9300) kt,  ztrai / areatot
 9300     FORMAT(i10,e18.10)
       ENDIF