anchor.diff:    endif
   
   
  -  ! ... Read h
  +  ! ... Read h and hraw
     ! ...
     err = NF90_INQ_VARID(r%fid,'h',ii)
     if (err.NE.NF90_NOERR) then
  
positive.code:   allocate(r%hgrid%h(r%hgrid%xi_rho,r%hgrid%eta_rho))
    err = NF90_GET_VAR(r%fid,ii,r%hgrid%h)
    if (err.NE.NF90_NOERR) stop 'ERROR while reading h'
  
    ! ... masks
    ! ...
  
negative.code: 
    err = NF90_INQ_VARID(r%fid,'lon_psi',ii)
    if (err.EQ.NF90_NOERR) then
      err = NF90_GET_VAR(r%fid,ii,r%hgrid%lon_psi)
      if (err.NE.NF90_NOERR) stop 'ERROR: Reading variable lon_psi'
    endif
  
    err = NF90_INQ_VARID(r%fid,'lat_psi',ii)
    if (err.EQ.NF90_NOERR) then
      err = NF90_GET_VAR(r%fid,ii,r%hgrid%lat_psi)
      if (err.NE.NF90_NOERR) stop 'ERROR: Reading variable lat_psi'
    endif
  