Changeset 7582 for branches/UKMO/dev_r5518_GO6_package_for_static_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

Timestamp:

2017-01-20T14:53:20+01:00 (7 years ago)

Author:

timgraham

Message:

Added code to allow output of some diagnostics which are constant in Met Office CMIP6 models
Changes described in UKESM ticket 324 (​https://code.metoffice.gov.uk/trac/UKESM/ticket/324)

File:

• branches/UKMO/dev_r5518_GO6_package_for_static_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90 (modified) (2 diffs)

branches/UKMO/dev_r5518_GO6_package_for_static_diagnostics/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

@@ -146 +146 @@
       ENDIF
 
+      !CMIP6 output of areacello and sftof (In NEMO this is currently just the T
+      !mask converted to a percentage)
+      CALL iom_put("areacello", e12t)
+      CALL iom_put("sftof", tmask(:,:,1)*100 )
+
+      !CMIP6 at timestep 1 output depth of bottom W level
+      !This will account for any modifications made to the bathy (e.g.
+      !suppression of isolated grid points
+      IF(kt == nit000) THEN
+         z2d(:,:)=0._wp
+         DO jj=1,jpj
+            DO ji=1,jpi
+               jk=mbathy(ji,jj)
+               IF(jk > 0) z2d(ji,jj) = gdepw_0(ji,jj,jk+1)
+            ENDDO
+         ENDDO
+         CALL iom_put( 'deptho', z2d )
+      ENDIF
+      
+
       ! Output of initial vertical scale factor
       CALL iom_put("e3t_0", e3t_0(:,:,:) )
@@ -257 +277 @@
       IF( iom_use('logavt') )   CALL iom_put( "logavt", LOG( MAX( 1.e-20_wp, avt  (:,:,:) ) ) )
       IF( iom_use('logavs') )   CALL iom_put( "logavs", LOG( MAX( 1.e-20_wp, fsavs(:,:,:) ) ) )
+
+      IF(kt == nit000) THEN
+         !Output 3D fields of background vertical diffusion (avtb and avmb)
+         DO jj=1,jpj
+            DO ji=1,jpi
+               z3d(ji,jj,:) = avtb(:)
+            ENDDO
+         ENDDO
+         CALL iom_put('avtb', z3d)
+         DO jj=1,jpj
+            DO ji=1,jpi
+               z3d(ji,jj,:) = avmb(:)
+            ENDDO
+         ENDDO
+         CALL iom_put('avmb', z3d)
+      ENDIF
 
       IF ( iom_use("sstgrad") .OR. iom_use("sstgrad2") ) THEN