anchor.diff:                                    grid%jms:grid%jme), stat=err)
               if (err /= 0) stop "variable:grid:2d: Allocation request denied"
   
  +            this%data_2d = 0
  +
               if (trim(this%forcing_var) /= "") then
                   allocate(this%dqdt_2d(grid%ims:grid%ime,    &
                                         grid%jms:grid%jme), stat=err)
                   if (err /= 0) stop "variable:grid:dqdt_2d: Allocation request denied"
  +                
  +                this%dqdt_2d = 0
               endif
   
           endif
  
positive.code:             allocate(this%data_3d(dims(1), dims(2), dims(3)), stat=err)
              if (err /= 0) stop "variable:dims:3d: Allocation request denied"
  
              if (trim(this%forcing_var) /= "") then
                  if (allocated(this%dqdt_3d)) deallocate(this%dqdt_3d)
                  allocate(this%dqdt_3d(dims(1), dims(2), dims(3)), stat=err)
                  if (err /= 0) stop "variable:dims:dqdt_3d: Allocation request denied"
              endif
          endif
  
  
negative.code:                                   grid%jms:grid%jme), stat=err)
              if (err /= 0) stop "variable:grid:3d: Allocation request denied"
  
              this%data_3d = 0
  
              if (trim(this%forcing_var) /= "") then
                  allocate(this%dqdt_3d(grid%ims:grid%ime,    &
                                        grid%kms:grid%kme,    &
                                        grid%jms:grid%jme), stat=err)
                  if (err /= 0) stop "variable:grid:dqdt_3d: Allocation request denied"
  
                  this%dqdt_3d = 0
              endif
  
          endif
  