anchor.diff:  
     call roms_read_hgrid(r%fid,r%hgrid,err)
   
  -  nullify(r%vgrid%s_r)
  -  nullify(r%vgrid%s_w)
  -  nullify(r%vgrid%Cs_r)
  -  nullify(r%vgrid%Cs_w)
  +  nullify(r%vgrid%S%sc_r)
  +  nullify(r%vgrid%S%sc_w)
  +  nullify(r%vgrid%S%Cs_r)
  +  nullify(r%vgrid%S%Cs_w)
   
     call roms_read_vgrid(r%fid,r%vgrid,err)
   
  
positive.code: use cdf
  
  implicit none
  
  TYPE MHgrid
    integer                             :: xi_rho  = -1
  
negative.code:     return
    endif
    err = NF90_INQUIRE_DIMENSION(r%fid,ii,len=r%vgrid%N)
    r%vgrid%S%N = r%vgrid%N
  
    err = NF90_INQ_DIMID(r%fid,'s_w',ii)
    if (err.ne.NF90_NOERR) then
  