anchor.diff:              call setup_varids(ncid,options)
           else
               ! otherwise, create a new file
  -            call create_file(filename,options)
  +            call create_file(filename, options, domain%model_time)
               ! and write constant (in time) variables
               call check( nf90_put_var(ncid, lat_id,    domain%lat), trim(filename)//":Latitude" )
               call check( nf90_put_var(ncid, lon_id,    domain%lon), trim(filename)//":Longitude" )
  
positive.code:     !! @param options   model wide options so the correct output can be created
      !!
      !!------------------------------------------------------------
      subroutine create_file(filename,options)
          character(len=255), intent(in) :: filename
          type(options_type), intent(in) :: options
  
          ! store real (not model) timestamp
          character(len=19) :: todays_date_time
  
negative.code: 
          call check( nf90_def_var(ncid, "time", NF90_DOUBLE, t_id, time_id), trim(err)//"time" )
          call check( nf90_put_att(ncid,time_id,"standard_name","time"))
          call check( nf90_put_att(ncid,time_id,"units",time%units()))
          call check( nf90_put_att(ncid,time_id,"UTCoffset","0"))
  
          call check( nf90_def_var(ncid, "lev", NF90_REAL, dimids(3), lev_id), trim(err)//"lev" )
  