map mesh migrate should happen only for read map scenario

also, remove #ifdef HAVE_MOAB #endif in this file

Author: iulian787

driver-moab/main/prep_ocn_mod.F90

@@ -420,7 +420,6 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
                'seq_maps.rc','atm2ocn_fmapname:','atm2ocn_fmaptype:',samegrid_ao, &
                'mapper_Fa2o initialization', esmf_map_flag, no_match)
           call shr_sys_flush(logunit)
-#ifdef HAVE_MOAB
           ! Call moab intx only if atm and ocn are init in moab
           if ((mbaxid .ge. 0) .and.  (mboxid .ge. 0)) then
             if (iamroot_CPLID) then
@@ -605,8 +604,7 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
             endif ! if (.not. samegrid_ao)
 
          endif ! if ((mbaxid .ge. 0) .and. (mboxid .ge. 0))
-! endif HAVE_MOAB
-#endif
+
        end if ! if (atm_present)
 
        ! atm_c2_ice flag is here because ICE and OCN are constrained to be on the same
@@ -694,7 +692,6 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
              write(logunit,F00) 'Initializing mapper_SFi2o'
           end if
           call seq_map_init_rearrolap(mapper_SFi2o, ice(1), ocn(1), 'mapper_SFi2o')
-#ifdef HAVE_MOAB
           if ( (mbixid .ge. 0) .and. (mboxid .ge. 0)) then
              if (iamroot_CPLID) then
                write(logunit,*) ' '
@@ -737,7 +734,6 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
             mapper_SFi2o%intx_context = ocn(1)%cplcompid
             mapper_SFi2o%mbname = 'mapper_SFi2o'
          endif
-#endif
        endif ! if (ice_present)
        call shr_sys_flush(logunit)
 
@@ -754,7 +750,6 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
           call seq_comm_getData(CPLID, mpicom=mpicom_CPLID, iamroot=iamroot_CPLID)
           call seq_comm_getData(CPLID, mpigrp=mpigrp_CPLID)   !  second group, the coupler group CPLID is global variable
 
-#ifdef HAVE_MOAB
           appname = "ROF_OCN_COU"//CHAR(0)
             ! rmapid  is a unique external number of MOAB app that takes care of map between rof and ocn mesh
           rmapid = 100*rof(1)%cplcompid + ocn(1)%cplcompid ! something different, to differentiate it
@@ -835,7 +830,6 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
          mapper_Rr2o_liq%weight_identifier = wgtIdFr2ol
          mapper_Rr2o_liq%mbname = 'mapper_Rr2o_liq'
 
-#endif
 
          if (iamroot_CPLID) then
             write(logunit,*) ' '
@@ -846,7 +840,6 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
             'seq_maps.rc', 'rof2ocn_ice_rmapname:', 'rof2ocn_ice_rmaptype:',samegrid_ro, &
             'mapper_Rr2o_ice  initialization', esmf_map_flag, no_match )
 ! us the same one for mapper_Rr2o_ice and mapper_Fr2o
-#ifdef HAVE_MOAB
 ! now take care of the mapper for MOAB mapper_Rr2o_ice
          if (iamroot_CPLID) then
             write(logunit,*) ' '
@@ -868,7 +861,6 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
          mapper_Rr2o_ice%intx_context = rmapid ! read map is the same context as intersection now
          mapper_Rr2o_ice%weight_identifier = wgtIdFr2oi
          mapper_Rr2o_ice%mbname = 'mapper_Rr2o_ice'
-#endif
          if (flood_present) then
             if (iamroot_CPLID) then
                write(logunit,*) ' '
@@ -878,8 +870,6 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
             call seq_map_init_rcfile( mapper_Fr2o, rof(1), ocn(1), &
                   'seq_maps.rc', 'rof2ocn_fmapname:', 'rof2ocn_fmaptype:',samegrid_ro, &
                   string='mapper_Fr2o initialization', esmf_map=esmf_map_flag, no_match=no_match )
-
-#ifdef HAVE_MOAB
             ! now take care of the mapper for MOAB mapper_Fr2o
             if (iamroot_CPLID) then
                write(logunit,*) ' '
@@ -892,30 +882,29 @@ subroutine prep_ocn_init(infodata, atm_c2_ocn, atm_c2_ice, ice_c2_ocn, rof_c2_oc
                call moab_map_init_rcfile( mbrxid, mboxid, mbintxro, type_grid, &
                      'seq_maps.rc', 'rof2ocn_fmapname:', 'rof2ocn_fmaptype:', samegrid_ro, &
                      arearead, wgtIdFr2o, 'mapper_Fr2o MOAB initialization', esmf_map_flag, wgtIdFr2ol )
-
             end if
-
             mapper_Fr2o%src_mbid = mbrxid
             mapper_Fr2o%tgt_mbid = mboxid ! special
             mapper_Fr2o%intx_mbid = mbintxro
             mapper_Fr2o%src_context = rof(1)%cplcompid
             mapper_Fr2o%intx_context = rmapid ! read map is the same context as intersection now
             mapper_Fr2o%weight_identifier = wgtIdFr2o
             mapper_Fr2o%mbname = 'mapper_Fr2o'
-#endif
          endif
 
-         context_id = rmapid ! ocn(1)%cplcompid
+         context_id = rmapid ! ocn(1)%cplcompid*100+rof(1)%cplcompid
+         type_grid = 3 ! this is type of grid
+         if (.not. compute_maps_online_r2o) then
+            ierr = iMOAB_MigrateMapMesh ( mbrxid, mbintxro, mpicom_CPLID, mpigrp_CPLID, &
+                                          mpigrp_CPLID, type_grid, rof(1)%cplcompid, context_id )
+            if (ierr .ne. 0) then
+               write(logunit,*) subname,' error in migrating rof mesh for map rof c2 ocn '
+               call shr_sys_abort(subname//' ERROR in migrating rof mesh for map rof c2 ocn ')
+            endif
+         endif
          ! this creates a parallel communication graph between mbrxid and mbintxro,
          ! with ids rof(1)%cplcompid, rmapid (rmapid is 100*src+tgt)
          ! this will be used in send/receive mappers
-         type_grid = 3 ! this is type of grid
-         ierr = iMOAB_MigrateMapMesh ( mbrxid, mbintxro, mpicom_CPLID, mpigrp_CPLID, &
-                                       mpigrp_CPLID, type_grid, rof(1)%cplcompid, context_id )
-         if (ierr .ne. 0) then
-            write(logunit,*) subname,' error in migrating rof mesh for map rof c2 ocn '
-            call shr_sys_abort(subname//' ERROR in migrating rof mesh for map rof c2 ocn ')
-         endif
          ierr = iMOAB_ComputeCommGraph( mbrxid, mbintxro, mpicom_CPLID, mpigrp_CPLID, mpigrp_CPLID, &
                                           type_grid, type_grid, rof(1)%cplcompid, rmapid )
          if (ierr .ne. 0) then