Changeset 7599 for branches/UKMO/dev_r5518_GO6_package_for_static_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90

Timestamp:

2017-01-23T16:12:59+01:00 (8 years ago)

Author:

timgraham

Message:

Fully working version

File:

• branches/UKMO/dev_r5518_GO6_package_for_static_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90 (modified) (1 diff)

branches/UKMO/dev_r5518_GO6_package_for_static_diagnostics/NEMOGCM/NEMO/OPA_SRC/TRA/trabbc.F90

@@ -91 +91 @@
       IF( kt == nit000 ) THEN
       !Output the geothermal heat flux once for CMIP6 diagnostics
-         wrk_alloc(jpi, jpj, ztmp)
-         ztmp(:,:) = qgh_trd0(:,:)*rcp_rau0
+         CALL wrk_alloc(jpi, jpj, ztmp)
+         DO jj=1,jpj
+            DO ji=1,jpi
+               ztmp(ji,jj) = qgh_trd0(ji,jj)*rau0_rcp
+            ENDDO
+         ENDDO
          CALL iom_put( 'hfgeou', ztmp )
-         wrk_dealloc(jpi, jpj, ztmp)
+         CALL wrk_dealloc(jpi, jpj, ztmp)
       ENDIF