anchor.diff:          type(variable_t) :: input_data
           ! number of layers has to be used when subsetting for update_pressure (for now)
           integer :: nz
  +        logical :: var_is_u, var_is_v
   
           update_only = .False.
           if (present(update)) update_only = update
  
positive.code:                        temporary_data, this%grid)
          this%longitude%data_2d = temporary_data(this%grid%ims:this%grid%ime, this%grid%jms:this%grid%jme)
  
          ! Read the u-grid longitude data if specified, other wise interpolate from mass grid
          if (options%parameters%ulon_hi /= "") then
              call load_data(options%parameters%init_conditions_file,   &
  
negative.code:             this%v_latitude%data_2d = temporary_data(this%v_grid%ims:this%v_grid%ime, this%v_grid%jms:this%v_grid%jme)
  
          else
              ! load the mass grid data again to get the full grid
              call load_data(options%parameters%init_conditions_file,   &
                             options%parameters%lat_hi,                 &
                             temporary_data, this%grid)
  
              associate(vlat => this%v_latitude%data_2d, grid => this%v_grid)
  
                  call array_offset_y(temporary_data, temp_offset)
                  vlat = temp_offset(grid%ims:grid%ime, grid%jms:grid%jme)
              end associate
          endif
  
  