anchor.diff:             self%ni = self%nx/self%nodex
              self%nj = 1+(self%ny-1)/self%nodey
           case default
  -           print *,'Unsupported grid'
  +           call cm1log(self, LOG_ERROR, 'cm1_set_nodes', 'Unsupported grid')
              stop
         end select
  -      print *,'[ingest_cm1::cm1_set_nodes]: ni=',self%ni,' nj=',self%nj
  +      write(output,102) self%ni, self%nj
  +      call cm1log(self, LOG_INFO, 'cm1_set_nodes', trim(output))
   
         allocate(recltest(self%ni,self%nj))
         inquire(iolength=self%mpireclen) recltest
         deallocate(recltest)
  -      print *,'[ingest_cm1::cm1_set_nodes]: mpi recl = ',self%mpireclen
  +      write(output,103) recltest
  +      call cm1log(self, LOG_INFO, 'cm1_set_nodes', trim(output))
   
         cm1_set_nodes = self%nodex*self%nodey
      end function cm1_set_nodes
  
positive.code: 
        status = self%readMultStart(time)
        if (status.eq.0) then
           print *,'[ingest_cm1:read3D]: data open failure, aborting'
           read3D = 0
           return
        endif
  
        status = self%read3DMult(varname, Field3D)
        if (status.eq.0) then
           print *,'[ingest_cm1:read3D]: data read failure, aborting'
           read3D = 0
           status = self%readMultStop()
           return
  
negative.code:                ! open dat file
                 open(newunit=self%dat_units(fu+1),file=datfile,form='unformatted',access='direct',recl=self%mpireclen,status='old')
              end do
              call cm1log(self, LOG_INFO, 'read3DMultStart', 'Multiread started for time: '//trim(dtime))
              readMultStart = 1
              self%ismult = .true.
  
           case (HDF)
              call cm1log(self, LOG_ERROR, 'read3DMultStart', 'Not implemented.')
  
        end select
  
  