Back to using MAPL_RemapBounds - gfortran has trouble with bounds remapping

Author: pchakraborty

DynCore_GridCompMod.F90

@@ -27,6 +27,7 @@ module FVdycoreCubed_GridComp
    use MAPL_AbstractRegridderMod, only: AbstractRegridder
    ! pchakrab - TODO: need MAPL3 equivalent
    ! use MAPL_SunMod, only: MAPL_SunOrbit, MAPL_SunGetInsolation
+   use MAPL_BaseMod, only: MAPL_RemapBounds
    use MAPL_GridManagerMod, only: grid_manager
    use MAPL_RegridderManagerMod, only: regridder_manager
    use MAPL_RegridMethods, only: REGRID_METHOD_BILINEAR
@@ -3583,7 +3584,10 @@ subroutine PULL_Q(self, import, QQQ, iNXQ, InFieldName, rc)
          self%vars%TRACER(N)%TNAME = fieldname
          if ( self%vars%TRACER(N)%IS_R4 ) then
             call ESMF_ArrayGet(array, localDE=0, farrayptr=ptr_r4, _RC)
-            self%vars%tracer(n)%content_r4(i1:in, j1:jn, 1:km) => ptr_r4
+            ! gfortran has trouble with pointer bounds remapping
+            ! (Error: Rank remapping target must be rank 1 or simply contiguous)
+            ! self%vars%tracer(n)%content_r4(i1:in, j1:jn, 1:km) => ptr_r4
+            self%vars%tracer(n)%content_r4 => MAPL_RemapBounds(ptr_r4, i1, in, j1, jn, 1, km)
             if (fieldname == QFieldName) then
                qqq%is_r4 = .true.
                qqq%content_r4 => self%vars%tracer(n)%content_r4