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:             read3DXYSlice = 1
  
           case (HDF)
              print*,'Not implemented'
  
        end select
  
  
negative.code: 
        1001 format('[ingestcm1',3(':',A),']: Error: ',A)
        1002 format('[ingestcm1',3(':',A),']: Warning: ',A)
        1004 format('[ingestcm1',3(':',A),']: ',A)
  
        select case (self%dtype)
  