Changeset 6223 for branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO

Timestamp:

2016-01-07T16:48:36+01:00 (8 years ago)

Author:

jpalmier

Message:

JPALM -- 07-01-2016 -- add debugg print statements

Location:

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC

Files:

• trcrst.F90 (modified) (6 diffs)
• trcstp.F90 (modified) (3 diffs)

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcrst.F90

@@ -47 +47 @@
    PUBLIC   trc_rst_wri       ! called by ???
    PUBLIC   trc_rst_cal
+   PUBLIC   trc_rst_stat
 
    !! * Substitutions
@@ -296 +297 @@
       !!
       !! calculate stats on these fields
+      call trc_rst_dia_stat(qint_idtra(:,:,1), 'qint_IDTRA')
       IF(lwp) WRITE(numout,*) ' IDTRA averaged properties stats (min, max, sum) ...'
       fq0 = MINVAL(qint_idtra(:,:,1))
@@ -319 +321 @@
          !!
          !! calculate stats on these fields
+         call trc_rst_dia_stat(qint_cfc(:,:,jl), 'qint_'//ctrcnm(jn))
          IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...'
          fq0 = MINVAL(qint_cfc(:,:,jl))
@@ -457 +460 @@
       !! calculate stats on these fields
       IF(lwp) WRITE(numout,*) ' IDTRA averaged properties stats (min, max, sum) ...'
+      call trc_rst_dia_stat(qint_idtra(:,:,1), 'qint_IDTRA')
       fq0 = MINVAL(qint_idtra(:,:,1))
       fq1 = MAXVAL(qint_idtra(:,:,1))
@@ -470 +474 @@
          !!
          !! calculate stats on these fields
+         call trc_rst_dia_stat(qint_cfc(:,:,jl), 'qint_'//ctrcnm(jn))
          IF(lwp) WRITE(numout,*) ' CFC averaged properties stats (min, max, sum) ...'
          fq0 = MINVAL(qint_cfc(:,:,jl))
@@ -640 +645 @@
          IF(lwp) WRITE(numout,9000) jn, TRIM( ctrcnm(jn) ), zmean, zmin, zmax, zdrift
       END DO
-      WRITE(numout,*) 
+      IF(lwp) WRITE(numout,*) 
 9000  FORMAT(' tracer nb :',i2,'    name :',a10,'    mean :',e18.10,'    min :',e18.10, &
       &      '    max :',e18.10,'    drift :',e18.10, ' %')
       !
    END SUBROUTINE trc_rst_stat
+
+
+   SUBROUTINE trc_rst_dia_stat
+      !!----------------------------------------------------------------------
+      !!                    ***  trc_rst_dia_stat  ***
+      !!
+      !! ** purpose  :   Compute tracers statistics
+      !!----------------------------------------------------------------------
+      REAL(wp), DIMENSION(jpi,jpj) , INTENT(in) ::   dgtr      ! ocean time-step
+      CHARACTER(len=*)             , INTENT(in) ::   names       ! "READ"/"WRITE" flag
+      !!---------------------------------------------------------------------
+      INTEGER  :: jk, jn
+      REAL(wp) :: ztraf, zmin, zmax, zmean, areasf
+      REAL(wp), DIMENSION(jpi,jpj) :: zvol
+      !!----------------------------------------------------------------------
+
+      IF( lwp )  WRITE(numout,*) 'STAT- ', names
+      !
+      zvol(:,:) = e1e2t(:,:) * fse3t_a(:,:,1) * tmask(:,:,1)
+      ztraf = glob_sum( dgtr(:,:) * zvol(:,:) )
+      areasf = glob_sum(e1e2t(:,:) * tmask(:,:,1) )
+      zmin  = MINVAL( dgtr(:,:), mask= ((tmask(:,:,1)*SPREAD(tmask_i,DIM=2,NCOPIES=1).NE.0.)) )
+      zmax  = MAXVAL( dgtr(:,:), mask= ((tmask(:,:,1)*SPREAD(tmask_i,DIM=2,NCOPIES=1).NE.0.)) )
+      IF( lk_mpp ) THEN
+         CALL mpp_min( zmin )      ! min over the global domain
+         CALL mpp_max( zmax )      ! max over the global domain
+      END IF
+      zmean  = ztraf / areatot
+      IF(lwp) WRITE(numout,9001) TRIM( names ), zmean, zmin, zmax
+      !
+      IF(lwp) WRITE(numout,*)
+9001  FORMAT(' tracer name :',a10,'    mean :',e18.10,'    min :',e18.10, &
+      &      '    max :',e18.10 )
+      !
+   END SUBROUTINE trc_rst_dia_stat
+
 
 #else

branches/NERC/dev_r5518_NOC_MEDUSA_Stable/NEMOGCM/NEMO/TOP_SRC/trcstp.F90

@@ -92 +92 @@
 # endif
                                    CALL trc_rst_opn  ( kt )       ! Open tracer restart file 
+# if defined key_debug_medusa
+                                   CALL trc_rst_stat 
+# endif
          IF( lrst_trc )            CALL trc_rst_cal  ( kt, 'WRITE' )   ! calendar
          IF( lk_iomput ) THEN  ;   CALL trc_wri      ( kt )       ! output of passive tracers with iom I/O manager
@@ -99 +102 @@
 # if defined key_debug_medusa
          IF(lwp) WRITE(numout,*) ' MEDUSA trc_stp SMS complete at kt =', kt
+         CALL trc_rst_stat
          CALL flush(numout)
 # endif
@@ -104 +108 @@
 # if defined key_debug_medusa
          IF(lwp) WRITE(numout,*) ' MEDUSA trc_stp transport complete at kt =', kt
+         CALL trc_rst_stat
          CALL flush(numout)
 # endif