anchor.diff:              allocate(nodeslice(ni,nj))
               do nf=0, nodex*nodey-1
   
  +              nodeslice = 0.0
                 read(unit_base+nf, rec=idx) nodeslice
                 sj = nf / nodex + 1
                 si = nf - (sj-1)*nodex + 1
  
positive.code:             do fu=0, nodex*nodey-1
                 close(unit_base+fu)
              end do
              readMultStop = 1
              ismult = 0
  
  
negative.code:                      slice((si-1)*ni+i,(sj-1)*(nj-1)+j) = nodeslice(i,j)
                     end do
                     end do
                   case default
                     print *,'[ingest_cm1::read3DXYSlice]: Unsupported MPI grid, something bad happened.'
                end select
  
              end do
  