anchor.diff:               chunk_size = 1
            end if
            ierr = geopm_prof_region(c_char_'AT_PHASE_INDEXMAP_2'//c_null_char, &
  -                                  GEOPM_REGION_HINT_UNKNOWN, region_id)
  +                                  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()
  
positive.code:       end if
  
        ierr = geopm_prof_region(c_char_'AT_PHASE_TRANS_1'//c_null_char, &
                                 GEOPM_REGION_HINT_UNKNOWN, region_id)
           ierr = geopm_prof_enter(region_id)
  !$omp parallel default(shared) private(i,j,ii,jj,z,tid)
           tid = omp_get_thread_num()
  
negative.code:           chunk_size = 1
        end if
        ierr = geopm_prof_region(c_char_'AT_PHASE_FFT_3'//c_null_char, &
                                GEOPM_REGION_HINT_COMPUTE, region_id)
        ierr = geopm_prof_enter(region_id)
  !$omp parallel default(shared) private(i,j,k,ii,y, tid)
        tid = omp_get_thread_num()
  