Changeset 9330 for branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90

Timestamp:

2018-02-16T18:29:04+01:00 (6 years ago)

Author:

dford

Message:

Add required diagnostics.

File:

• branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90 (modified) (3 diffs)

branches/UKMO/CO6_KD490_amm7_oper_fabm_hem08/NEMOGCM/NEMO/TOP_SRC/FABM/trcsms_fabm.F90

@@ -33 +33 @@
    USE inputs_fabm
    USE vertical_movement_fabm
+   USE zdfmxl
 
    !USE fldread         !  time interpolation
@@ -113 +114 @@
 
       CALL st2d_fabm_nxt( kt )
+      
+      CALL asmdiags_fabm
 
       IF( l_trdtrc )  CALL wrk_alloc( jpi, jpj, jpk, ztrfabm )
@@ -130 +133 @@
 
    END SUBROUTINE trc_sms_fabm
+   
+   SUBROUTINE asmdiags_fabm
+      INTEGER :: ji,jj,jk,jkmax
+   
+      PGROW_AVG(:,:) = 1.0
+      PLOSS_AVG(:,:) = 1.0
+      PHYT_AVG(:,:)  = 0.0
+
+      DO jj = 1, jpj
+         DO ji = 1, jpi
+            !
+            jkmax = jpk-1
+            DO jk = jpk-1, 1, -1
+               IF ( ( hmld_tref(ji,jj) >  gdepw_n(ji,jj,jk)   ) .AND. &
+                  & ( hmld_tref(ji,jj) <= gdepw_n(ji,jj,jk+1) ) ) THEN
+                  hmld_tref(ji,jj) = gdepw_n(ji,jj,jk+1)
+                  jkmax = jk
+               ENDIF
+            END DO
+            !
+            DO jk = 1, jkmax
+               PHYT_AVG(ji,jj) = PHYT_AVG(ji,jj) + &
+                  &              trn(ji,jj,jk,jp_fabm_p1n) + &
+                  &              trn(ji,jj,jk,jp_fabm_p2n) + &
+                  &              trn(ji,jj,jk,jp_fabm_p3n) + &
+                  &              trn(ji,jj,jk,jp_fabm_p4n)
+            END DO
+            
+            PHYT_AVG(ji,jj) = PHYT_AVG(ji,jj) / REAL(jkmax)
+
+            IF ( hmld_tref(ji,jj) .GT. MLD_MAX(ji,jj) ) THEN
+               MLD_MAX(ji,jj) = hmld_tref(ji,jj)
+            ENDIF
+            !
+         END DO
+      END DO
+   
+   END SUBROUTINE asmdiags_fabm
 
    SUBROUTINE compute_fabm()