anchor.diff:  ! do transpose among all  processes with same 1-coord (me1)
   !---------------------------------------------------------------------
         ierr = geopm_prof_region(c_char_'AT_PHASE_MPI_ATA_2'//c_null_char, &
  -                               GEOPM_REGION_HINT_UNKNOWN, region_id)
  +                               GEOPM_REGION_HINT_NETWORK, region_id)
         ierr = geopm_prof_enter(region_id)
         if (timers_enabled)call timer_start(T_transxzglo)
         call mpi_alltoall(xin, d1*d2*d3/np2, dc_type, &
  
positive.code:              chunk_size = 1
           end if
           ierr = geopm_prof_region(c_char_'AT_PHASE_INDEXMAP_1'//c_null_char, &
                                   GEOPM_REGION_HINT_UNKNOWN, region_id)
           ierr = geopm_prof_enter(region_id)
  !$omp parallel default(shared) private(i,j,k,ii,ii2,jj,ij2,kk,tid)
           tid = omp_get_thread_num()
  
negative.code:              chunk_size = 1
           end if
           ierr = geopm_prof_region(c_char_'AT_PHASE_INDEXMAP_2'//c_null_char, &
                                    GEOPM_REGION_HINT_COMPUTE, region_id)
           ierr = geopm_prof_enter(region_id)
  !$omp parallel default(shared) private(i,j,k,ii,ii2,jj,ij2,kk,tid)
           tid = omp_get_thread_num()
  