anchor.diff:        do i = 1, self%nt
            self%times(i) = self%t0 + ((i-1)*self%dt)
         end do
  -      print 500,'T',self%nt
  +      write (output,500) 'T',self%nt
  +      call cm1log(self, LOG_MSG, 'read_ctl', trim(output))
   
         ! variables
         read(self%ctl_unit,*) tmp, self%nv
  
positive.code:       ! Does the variable exist in this dataset?
        varid = self%getVarByName(varname)
        if (varid.eq.0) then
           print *,'[ingest_cm1::read2DMult]: variable not found: ',trim(varname)
           read2DMult = 0
           return
        end if
  
        print *,'[ingest_cm1::read2DMult]: Reading: ',trim(varname)
        ! Read the variable from the dataset
        status = self%read3DXYSlice(varid, 0, Field2D(:,:))
  
  
negative.code: 
        self%isopen = 0
        close_cm1 = 1
        call cm1log(self, LOG_MSG, 'close_cm1', 'Dataset closed:')
  
     end function
  
  