Changeset 8372 for branches/NERC/dev_r5518_GO6_CleanMedusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

Timestamp:

2017-07-25T18:58:24+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- debug switch TOP module call order ; remove trc2d and trc3d from MEDUSA (obsolete); rm trcrst_medusa (obsolete)

File:

• branches/NERC/dev_r5518_GO6_CleanMedusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90 (modified) (4 diffs)

branches/NERC/dev_r5518_GO6_CleanMedusa/NEMOGCM/NEMO/TOP_SRC/MEDUSA/bio_medusa_fin.F90

@@ -51 +51 @@
                                    zn_sed_c, zn_sed_ca, zn_sed_fe,      &
                                    zn_sed_n, zn_sed_si
-      USE trc,               ONLY: ctrc2d, ctrc3d, ln_diatrc,           &
-                                   med_diag, nittrc000,                 &
-                                   trbio, trc2d, trc3d
+      USE trc,               ONLY: med_diag, nittrc000 
       USE trcnam_trp,        ONLY: ln_trcadv_cen2, ln_trcadv_tvd
  
@@ -100 +98 @@
          zn_sed_ca(:,:) = za_sed_ca(:,:)
       endif
-      IF( ln_diatrc ) THEN
-         DO jj = 2,jpjm1
-            DO ji = 2,jpim1
-               trc2d(ji,jj,131) = za_sed_n(ji,jj)
-               trc2d(ji,jj,132) = za_sed_fe(ji,jj)
-               trc2d(ji,jj,133) = za_sed_c(ji,jj)
-               trc2d(ji,jj,134) = za_sed_si(ji,jj)
-               trc2d(ji,jj,135) = za_sed_ca(ji,jj)
-            ENDDO
-         ENDDO
-      ENDIF 
       !!
       if (ibenthic.eq.2) then
@@ -314 +301 @@
 #  endif
 
-      IF( ln_diatrc ) THEN
-         !!-----------------------------------------------------------------
-         !! Output several accumulated diagnostics
-         !!   - biomass-average phytoplankton limitation terms
-         !!   - integrated tendency terms
-         !!-----------------------------------------------------------------
-         !! 
-         DO jj = 2,jpjm1
-            DO ji = 2,jpim1
-               !! non-diatom phytoplankton limitations
-               trc2d(ji,jj,25)  = trc2d(ji,jj,25) / MAX(ftot_pn(ji,jj), rsmall)
-               trc2d(ji,jj,26)  = trc2d(ji,jj,26) / MAX(ftot_pn(ji,jj), rsmall)
-               trc2d(ji,jj,27)  = trc2d(ji,jj,27) / MAX(ftot_pn(ji,jj), rsmall)
-               !! diatom phytoplankton limitations
-               trc2d(ji,jj,28)  = trc2d(ji,jj,28) / MAX(ftot_pd(ji,jj), rsmall)
-               trc2d(ji,jj,29)  = trc2d(ji,jj,29) / MAX(ftot_pd(ji,jj), rsmall)
-               trc2d(ji,jj,30)  = trc2d(ji,jj,30) / MAX(ftot_pd(ji,jj), rsmall)
-               trc2d(ji,jj,31)  = trc2d(ji,jj,31) / MAX(ftot_pd(ji,jj), rsmall)
-               trc2d(ji,jj,32)  = trc2d(ji,jj,32) / MAX(ftot_pd(ji,jj), rsmall)
-               !! tendency terms
-               trc2d(ji,jj,76)  = fflx_n(ji,jj)
-               trc2d(ji,jj,77)  = fflx_si(ji,jj)
-               trc2d(ji,jj,78)  = fflx_fe(ji,jj)
-               !! Integrated biomass.
-               !! integrated non-diatom phytoplankton
-               trc2d(ji,jj,79)  = ftot_pn(ji,jj)
-               !! integrated diatom phytoplankton
-               trc2d(ji,jj,80)  = ftot_pd(ji,jj)
-               !! Integrated microzooplankton
-               trc2d(ji,jj,217) = ftot_zmi(ji,jj)
-               !! Integrated mesozooplankton
-               trc2d(ji,jj,218) = ftot_zme(ji,jj)
-               !! Integrated slow detritus, nitrogen
-               trc2d(ji,jj,219) = ftot_det(ji,jj)
-               !! Integrated slow detritus, carbon
-               trc2d(ji,jj,220) = ftot_dtc(ji,jj)
-# if defined key_roam
-               !! The balance of nitrogen production/consumption.
-               !! integrated nitrogen production
-               trc2d(ji,jj,111) = fnit_prod(ji,jj)
-               !! integrated nitrogen consumption
-               trc2d(ji,jj,112) = fnit_cons(ji,jj)
-               !! The balance of carbon production/consumption.
-               !! integrated carbon production
-               trc2d(ji,jj,113) = fcar_prod(ji,jj)
-               !! integrated carbon consumption
-               trc2d(ji,jj,114) = fcar_cons(ji,jj)
-               !! The balance of oxygen production/consumption.
-               !! integrated oxygen production
-               trc2d(ji,jj,115) = foxy_prod(ji,jj)
-               !! integrated oxygen consumption
-               trc2d(ji,jj,116) = foxy_cons(ji,jj)
-               !! integrated unrealised oxygen consumption
-               trc2d(ji,jj,117) = foxy_anox(ji,jj)
-# endif
-            ENDDO
-         ENDDO
-         
-# if defined key_roam
-#  if defined key_axy_nancheck
-         !!---------------------------------------------------------------
-         !! Check for NaNs in diagnostic outputs
-         !!---------------------------------------------------------------
-         !! 
-         !! 2D diagnostics
-         DO jn = 1,150
-            fq0 = SUM(trc2d(:,:,jn))
-            !! AXY (30/01/14): "isnan" problem on HECTOR
-            !! if (fq0 /= fq0 ) then
-            if ( ieee_is_nan( fq0 ) ) then
-               !! there's a NaN here
-               if (lwp) write(numout,*)                                     &
-                  'NAN detected in 2D diagnostic field', jn, 'at time',     &
-                   kt, 'at position:'
-               DO jj = 1,jpj
-                  DO ji = 1,jpi
-                     if ( ieee_is_nan( trc2d(ji,jj,jn) ) ) then
-                        if (lwp) write (numout,'(a,3i6)') 'NAN-CHECK',      &
-                                 ji, jj, jn
-                     endif
-                  ENDDO
-               ENDDO
-          CALL ctl_stop( 'trcbio_medusa, NAN in 2D diagnostic field' )
-            endif
-         ENDDO
-         !!
-         !! 3D diagnostics
-         DO jn = 1,5
-            fq0 = SUM(trc3d(:,:,:,jn))
-            !! AXY (30/01/14): "isnan" problem on HECTOR
-            !! if (fq0 /= fq0 ) then
-            if ( ieee_is_nan( fq0 ) ) then
-               !! there's a NaN here
-               if (lwp) write(numout,*)                                     &
-                  'NAN detected in 3D diagnostic field', jn, 'at time',     &
-                  kt, 'at position:'
-               DO jk = 1,jpk
-                  DO jj = 1,jpj
-                     DO ji = 1,jpi
-                        if ( ieee_is_nan( trc3d(ji,jj,jk,jn) ) ) then
-                           if (lwp) write (numout,'(a,4i6)') 'NAN-CHECK',   &
-                                    ji, jj, jk, jn
-                        endif
-                     ENDDO
-                  ENDDO
-               ENDDO
-          CALL ctl_stop( 'trcbio_medusa, NAN in 3D diagnostic field' )
-            endif
-         ENDDO
-    CALL flush(numout)
-#  endif
-# endif
-
-         !!----------------------------------------------------------------
-         !! Don't know what this does; belongs to someone else ...
-         !!----------------------------------------------------------------
-         !! 
-         !! Lateral boundary conditions on trc2d
-         DO jn=1,jp_medusa_2d
-             CALL lbc_lnk(trc2d(:,:,jn),'T',1. )
-         ENDDO 
-
-         !! Lateral boundary conditions on trc3d
-         DO jn=1,jp_medusa_3d
-             CALL lbc_lnk(trc3d(:,:,1,jn),'T',1. )
-         ENDDO 
-
-
-# if defined key_axy_nodiag
-         !!----------------------------------------------------------------
-         !! Blank diagnostics as a NaN-trap
-         !!----------------------------------------------------------------
-         !! 
-         !! blank 2D diagnostic array
-         trc2d(:,:,:) = 0.e0
-         !!
-         !! blank 3D diagnostic array
-         trc3d(:,:,:,:) = 0.e0
-# endif
-
-
-         !!----------------------------------------------------------------
-         !! Add in XML diagnostics stuff
-         !!----------------------------------------------------------------
-         !!
-         !! ** 2D diagnostics
-         DO jn=1,jp_medusa_2d
-            CALL iom_put(TRIM(ctrc2d(jn)), trc2d(:,:,jn))
-         END DO
-!! AXY (17/02/14): don't think I need this if I modify the above for 
-!! all diagnostics
-!! #  if defined key_roam
-!!          DO jn=91,jp_medusa_2d
-!!             CALL iom_put(TRIM(ctrc2d(jn)), trc2d(:,:,jn))
-!!          END DO      
-!! #  endif
-         !!
-         !! ** 3D diagnostics
-         DO jn=1,jp_medusa_3d
-            CALL iom_put(TRIM(ctrc3d(jn)), trc3d(:,:,:,jn))
-         END DO
-!! AXY (17/02/14): don't think I need this if I modify the above for 
-!! all diagnostics
-!! #  if defined key_roam
-!!          CALL iom_put(TRIM(ctrc3d(5)), trc3d(:,:,:,5))
-!! #  endif
-
-
-      ELSE IF( lk_iomput .AND. .NOT. ln_diatrc ) THEN
          !!!---------------------------------------------------------------
          !! Add very last diag calculations 
@@ -1173 +991 @@
           DEALLOCATE( zw2d )
 
-       ENDIF                    ! end of ln_diatrc option
-
    END SUBROUTINE bio_medusa_fin