Changeset 7829 for branches/UKMO/dev_r5518_nemo_fabm_ukmo/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcwri_my_trc.F90

Timestamp:

2017-03-23T15:36:14+01:00 (7 years ago)

Author:

dford

Message:

Add a version of the NEMO-FABM coupling code. In theory, this should give equivalent results to PML gitlab commit 2e51db55.

File:

• branches/UKMO/dev_r5518_nemo_fabm_ukmo/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcwri_my_trc.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_nemo_fabm_ukmo/NEMOGCM/NEMO/TOP_SRC/MY_TRC/trcwri_my_trc.F90

@@ -19 +19 @@
 
    PUBLIC trc_wri_my_trc 
+#if defined key_tracer_budget
+   REAL(wp), ALLOCATABLE, DIMENSION(:,:,:,:), SAVE :: trb_temp ! slwa
+#endif
+
 
 #  include "top_substitute.h90"
 CONTAINS
 
+#if defined key_tracer_budget
+   SUBROUTINE trc_wri_my_trc (kt, fl) ! slwa
+#else
    SUBROUTINE trc_wri_my_trc
+#endif
       !!---------------------------------------------------------------------
       !!                     ***  ROUTINE trc_wri_trc  ***
@@ -29 +37 @@
       !! ** Purpose :   output passive tracers fields 
       !!---------------------------------------------------------------------
+#if defined key_tracer_budget
+      INTEGER, INTENT( in ), OPTIONAL     :: fl 
+      INTEGER, INTENT( in )               :: kt
+      REAL(wp), DIMENSION(jpi,jpj,jpk)    :: trpool !tracer pool temporary output
+#endif
       CHARACTER (len=20)   :: cltra
-      INTEGER              :: jn
+      INTEGER              :: jn,jk
       !!---------------------------------------------------------------------
  
       ! write the tracer concentrations in the file
       ! ---------------------------------------
+
+
+#if defined key_tracer_budget
+      IF( PRESENT(fl)) THEN
+! depth integrated
+! for strict budgetting write this out at end of timestep as an average between 'now' and 'after' at kt
+         DO jn = jp_myt0, jp_myt1 
+            trpool(:,:,:) = 0.5 * ( trn(:,:,:,jn) * fse3t_a(:,:,:) +  &
+                                        trb_temp(:,:,:,jn) * fse3t(:,:,:) )
+!
+            cltra = TRIM( ctrcnm(jn) )                  ! output of tracer density 
+            CALL iom_put( cltra, trpool(:,:,:) / (0.5* (fse3t(:,:,:) + fse3t_a(:,:,:) ) ) )
+!
+            cltra = TRIM( ctrcnm(jn) )//"_pool"     ! volume integrated output
+            DO jk = 1, jpk
+               trpool(:,:,jk) = trpool(:,:,jk) * e1t(:,:) * e2t(:,:)
+            END DO
+            CALL iom_put( cltra, trpool)
+
+!           cltra = TRIM( ctrcnm(jn) )//"_pool"     ! volume integrated output
+!           DO jk = 1, jpk
+!              trpool(:,:,jk) = 0.5 * ( trn(:,:,jk,jn) * fse3t_a(:,:,jk) +  & 
+!                                       trb_temp(:,:,jk,jn) * fse3t(:,:,jk) ) * & 
+!                                       e1t(:,:) * e2t(:,:)
+!           END DO
+!           CALL iom_put( cltra, trpool)
+!           cltra = TRIM( ctrcnm(jn) )                  ! output of tracer density 
+!           CALL iom_put( cltra, trpool(:,:,:) / (0.5* (fse3t(:,:,:) + fse3t_a(:,:,:) ) ) )
+         END DO
+         CALL iom_put( "DEPTH" , 0.5* (fse3t(:,:,:) + fse3t_a(:,:,:) ) )  !  equivalent 'depth' at same time as tracer pool output
+      ELSE
+
+         IF( kt == nittrc000 ) THEN
+           ALLOCATE(trb_temp(jpi,jpj,jpk,jptra))  ! slwa
+         ENDIF
+         trb_temp(:,:,:,:)=trn(:,:,:,:) ! slwa save for tracer budget (unfiltered trn)
+
+!        DO jn = jp_myt0, jp_myt1
+!           cltra = TRIM( ctrcnm(jn) )                  ! short title for tracer
+!           CALL iom_put( cltra, trn(:,:,:,jn) ) 
+!        END DO
+! write out depths and areas in double precision for tracer budget calculations
+         CALL iom_put( "AREA" , e1t(:,:) * e2t(:,:))
+!        CALL iom_put( "DEPTH" , fse3t(:,:,:) )  ! need depth at same time as tracer output
+
+      END IF
+#else
       DO jn = jp_myt0, jp_myt1
          cltra = TRIM( ctrcnm(jn) )                  ! short title for tracer
          CALL iom_put( cltra, trn(:,:,:,jn) )
       END DO
+#endif
       !
    END SUBROUTINE trc_wri_my_trc