Changeset 10810 for NEMO/branches/UKMO/dev_r10037_vorticity_trends/src/OCE/IOM/iom.F90

Timestamp:

2019-03-29T09:43:38+01:00 (5 years ago)

Author:

davestorkey

Message:

UKMO/dev_r10037_vorticity_trends branch: commit code changes.

File:

• NEMO/branches/UKMO/dev_r10037_vorticity_trends/src/OCE/IOM/iom.F90 (modified) (4 diffs)

NEMO/branches/UKMO/dev_r10037_vorticity_trends/src/OCE/IOM/iom.F90

@@ -153 +153 @@
          CALL set_grid( "U", glamu, gphiu, .FALSE., .FALSE. )
          CALL set_grid( "V", glamv, gphiv, .FALSE., .FALSE. )
+         CALL set_grid( "F", glamf, gphif, .FALSE., .FALSE. )
          CALL set_grid( "W", glamt, gphit, .FALSE., .FALSE. )
          CALL set_grid_znl( gphit )
@@ -160 +161 @@
             CALL iom_set_domain_attr("grid_U", area = e1e2u(nldi:nlei, nldj:nlej))
             CALL iom_set_domain_attr("grid_V", area = e1e2v(nldi:nlei, nldj:nlej))
+            CALL iom_set_domain_attr("grid_F", area = e1e2f(nldi:nlei, nldj:nlej))
             CALL iom_set_domain_attr("grid_W", area = e1e2t(nldi:nlei, nldj:nlej))
             CALL set_grid_bounds( "T", glamf, gphif, glamt, gphit )
             CALL set_grid_bounds( "U", glamv, gphiv, glamu, gphiu )
             CALL set_grid_bounds( "V", glamu, gphiu, glamv, gphiv )
+            CALL set_grid_bounds( "F", glamt, gphit, glamf, gphif )
             CALL set_grid_bounds( "W", glamf, gphif, glamt, gphit )
          ENDIF
@@ -174 +177 @@
          CALL set_grid( "U", glamu_crs, gphiu_crs, .FALSE., .FALSE. ) 
          CALL set_grid( "V", glamv_crs, gphiv_crs, .FALSE., .FALSE. ) 
+         CALL set_grid( "F", glamf_crs, gphif_crs, .FALSE., .FALSE. ) 
          CALL set_grid( "W", glamt_crs, gphit_crs, .FALSE., .FALSE. ) 
          CALL set_grid_znl( gphit_crs )
@@ -183 +187 @@
             CALL iom_set_domain_attr("grid_U", area = e1u_crs(nldi:nlei, nldj:nlej) * e2u_crs(nldi:nlei, nldj:nlej))
             CALL iom_set_domain_attr("grid_V", area = e1v_crs(nldi:nlei, nldj:nlej) * e2v_crs(nldi:nlei, nldj:nlej))
+            CALL iom_set_domain_attr("grid_F", area = e1f_crs(nldi:nlei, nldj:nlej) * e2f_crs(nldi:nlei, nldj:nlej))
             CALL iom_set_domain_attr("grid_W", area = e1e2t_crs(nldi:nlei, nldj:nlej))
             CALL set_grid_bounds( "T", glamf_crs, gphif_crs, glamt_crs, gphit_crs )
             CALL set_grid_bounds( "U", glamv_crs, gphiv_crs, glamu_crs, gphiu_crs )
             CALL set_grid_bounds( "V", glamu_crs, gphiu_crs, glamv_crs, gphiv_crs )
+            CALL set_grid_bounds( "F", glamt_crs, gphit_crs, glamf_crs, gphif_crs )
             CALL set_grid_bounds( "W", glamf_crs, gphif_crs, glamt_crs, gphit_crs )
          ENDIF