anchor.diff:                  
                   do i = 1, v_ndims
                       var_dims(i) = &
  -                        ncdr_files(file_ncid)%dims( &
  -                            ncdr_files(file_ncid)%vars(var_index)%var_dim_inds(i) &
  +                        ncdr_files(nc_diag_read_get_index_from_ncid(file_ncid))%dims( &
  +                            ncdr_files(nc_diag_read_get_index_from_ncid(file_ncid))%vars(var_index)%var_dim_inds(i) &
                               )%dim_size
                   end do
               end if
  
positive.code:             
              var_index = nc_diag_read_id_assert_var(file_ncid, var_name)
              
              var_type = ncdr_files(file_ncid)%vars(var_index)%var_type
          end function nc_diag_read_id_get_var_type
          
          function nc_diag_read_noid_get_var_type(var_name) result(var_type)
  
negative.code:             
              var_index = nc_diag_read_id_assert_var(file_ncid, var_name)
              
              var_ndims = ncdr_files(nc_diag_read_get_index_from_ncid(file_ncid))%vars(var_index)%var_ndims
          end function nc_diag_read_id_get_var_ndims
          
          function nc_diag_read_noid_get_var_ndims(var_name) result(var_ndims)
  