add code for dataPtr_area when creating the grid and attaching area t…

…o grid
NOAA-GFDL · Feb 26, 2019 · 76b0f4d · 76b0f4d
1 parent 8447a7d
commit 76b0f4d
Showing 1 changed file with 5 additions and 2 deletions.
diff --git a/config_src/nuopc_driver/mom_cap.F90 b/config_src/nuopc_driver/mom_cap.F90
@@ -1513,11 +1513,14 @@ subroutine InitializeRealize(gcomp, importState, exportState, clock, rc)
            dataPtr_mask(i1,j1)  = ocean_grid%mask2dT(ig,jg)
            dataPtr_xcen(i1,j1)  = ocean_grid%geolonT(ig,jg)
            dataPtr_ycen(i1,j1)  = ocean_grid%geolatT(ig,jg)
+           if(grid_attach_area) then
+            dataPtr_area(i1,j1)  = ocean_grid%areaT(ig,jg)
+           end if
          end do
        end do
 
-                       jlast = jec
-       if(jec .eq. nyg)jlast = jec+1
+                     jlast = jec
+       if(jec == nyg)jlast = jec+1
 
        do j = jsc, jlast
          j1 = j + lbnd4 - jsc