Changeset 8557 for branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

Timestamp:

2017-09-22T13:06:31+02:00 (7 years ago)

Author:

jpalmier

Message:

JPALM -- trdtrc -- fix SBC trend (2D not 3D) -- remove DMP from XIOS group to avoid NaNs?

File:

• branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90 (modified) (5 diffs)

branches/NERC/dev_r5518_GO6_package_trdtrc/NEMOGCM/NEMO/TOP_SRC/TRP/trcsbc.F90

@@ -68 +68 @@
       CHARACTER (len=22) :: charout
       REAL(wp), POINTER, DIMENSION(:,:  ) :: zsfx
-      REAL(wp), POINTER, DIMENSION(:,:,:) :: ztrtrd
+      REAL(wp), POINTER, DIMENSION(:,:  ) :: ztrtrd
 
       !!---------------------------------------------------------------------
@@ -76 +76 @@
       ! Allocate temporary workspace
                       CALL wrk_alloc( jpi, jpj,      zsfx   )
-      IF( l_trdtrc )  CALL wrk_alloc( jpi, jpj, jpk, ztrtrd )
+      IF( l_trdtrc )  CALL wrk_alloc( jpi, jpj,      ztrtrd )
       !
       zrtrn = 1.e-15_wp
@@ -140 +140 @@
       DO jn = 1, jptra
          !
-         IF( l_trdtrc ) ztrtrd(:,:,:) = tra(:,:,:,jn)  ! save trends
+         IF( l_trdtrc ) ztrtrd(:,:) = tra(:,:,1,jn)  ! save trends
          !                                             ! add the trend to the general tracer trend
 
@@ -184 +184 @@
          !
          IF( l_trdtrc ) THEN
-            ztrtrd(:,:,:) = tra(:,:,:,jn) - ztrtrd(:,:,:)
+            ztrtrd(:,:) = tra(:,:,1,jn) - ztrtrd(:,:)
             CALL trd_tra( kt, 'TRC', jn, jptra_nsr, ztrtrd )
          END IF
@@ -207 +207 @@
                                            CALL prt_ctl_trc( tab4d=tra, mask=tmask, clinfo=ctrcnm, clinfo2='trd' )
       ENDIF
-                      CALL wrk_dealloc( jpi, jpj,      zsfx   )
-      IF( l_trdtrc )  CALL wrk_dealloc( jpi, jpj, jpk, ztrtrd )
+                      CALL wrk_dealloc( jpi, jpj, zsfx   )
+      IF( l_trdtrc )  CALL wrk_dealloc( jpi, jpj, ztrtrd )
       !
       IF( nn_timing == 1 )  CALL timing_stop('trc_sbc')