Changeset 10205 for branches/UKMO/AMM15_v3_6_STABLE_package_bgc_updates/NEMOGCM/NEMO/OPA_SRC/DIA/diatmb.F90

Timestamp:

2018-10-19T13:36:08+02:00 (5 years ago)

Author:

dford

Message:

Add 25h mean and TMB diagnostics for FABM-ERSEM variables.

File:

• branches/UKMO/AMM15_v3_6_STABLE_package_bgc_updates/NEMOGCM/NEMO/OPA_SRC/DIA/diatmb.F90 (modified) (3 diffs)

branches/UKMO/AMM15_v3_6_STABLE_package_bgc_updates/NEMOGCM/NEMO/OPA_SRC/DIA/diatmb.F90

@@ -11 +11 @@
    USE iom             ! I/0 library
    USE wrk_nemo        ! working arrays
+#if defined key_fabm
+   USE trc, ONLY: trn
+   USE par_fabm
+   USE fabm, ONLY: fabm_get_bulk_diagnostic_data
+#endif
 
 
@@ -133 +138 @@
       REAL(wp), POINTER, DIMENSION(:,:,:) :: zwtmb    ! temporary workspace 
       REAL(wp)                         ::   zmdi      ! set masked values
+      INTEGER                          ::   jn        ! loop counter
 
       zmdi=1.e+20 !missing data indicator for maskin
@@ -162 +168 @@
          CALL iom_put( "bot_v" , zwtmb(:,:,3) )    ! tmb  V Velocity
 !Called in  dynspg_ts.F90       CALL iom_put( "baro_v" , vn_b )    ! Barotropic  V Velocity
+
+#if defined key_fabm
+         DO jn = 1, jp_fabm
+            CALL dia_calctmb( trn(:,:,:,jp_fabm_m1+jn), zwtmb )
+            CALL iom_put( "top_"//TRIM(model%state_variables(jn)%name) , zwtmb(:,:,1) )
+            CALL iom_put( "mid_"//TRIM(model%state_variables(jn)%name) , zwtmb(:,:,2) )
+            CALL iom_put( "bot_"//TRIM(model%state_variables(jn)%name) , zwtmb(:,:,3) )
+         END DO
+         DO jn = 1, jp_fabm_3d
+            CALL dia_calctmb( fabm_get_bulk_diagnostic_data(model, jn), zwtmb )
+            CALL iom_put( "top_"//TRIM(model%diagnostic_variables(jn)%name) , zwtmb(:,:,1) )
+            CALL iom_put( "mid_"//TRIM(model%diagnostic_variables(jn)%name) , zwtmb(:,:,2) )
+            CALL iom_put( "bot_"//TRIM(model%diagnostic_variables(jn)%name) , zwtmb(:,:,3) )
+         END DO
+#endif
       ELSE
          CALL ctl_warn('dia_tmb: tmb diagnostic is set to false you should not have seen this')