SUBROUTINE GET_GRID_NETCDF()

  ! UW Land Surface Hydrology Group implementation of SAC/SNOW17 model
  ! modified from NLDAS implementation
  ! author: Ted Bohn, tbohn@hydro.washington.edu

  ! OPENS NETCDF-FORMAT LSC FILE AND GETS DATA DIMENSIONS AND LANDMASK

  ! driverMod contains definitions for all global driver variables
  USE driverMod

  IMPLICIT NONE
      

  ! Define local variables
  INTEGER ndims,nvars,ngatts,natts,unlimited
  INTEGER xdimid,ydimid,landid,zdimid,landcoverdimid,tstepdimid
  CHARACTER*20 name
  INTEGER xtype,dimids(2)
  INTEGER varid

  ! Open LSC file and get data dimensions
  ! Assume that LSC file is NOT compressed
  status = NF_OPEN(LSC, 0, LSC_NCID)
  status = NF_INQ(LSC_NCID, ndims, nvars, ngatts, unlimited)
  status = NF_INQ_DIMID(LSC_NCID,'z',zdimid)
  status = NF_INQ_DIMID(LSC_NCID,'Longitude',ydimid)
  status = NF_INQ_DIMID(LSC_NCID,'Latitude',xdimid)
  status = NF_INQ_DIMID(LSC_NCID,'t',tstepdimid)
  status = NF_INQ_DIMLEN(LSC_NCID,zdimid,MAXNSOIL)
  status = NF_INQ_DIMLEN(LSC_NCID,ydimid,ylen)
  status = NF_INQ_DIMLEN(LSC_NCID,xdimid,xlen)
  status = NF_INQ_DIMLEN(LSC_NCID,tstepdimid,NMONTHS)

  ! Get LANDMASK
  ALLOCATE (LANDMASK(xlen,ylen))
  status = NF_INQ_VARID(LSC_NCID,'land',varid)
  status = NF_GET_VAR_INT(LSC_NCID,varid,LANDMASK)
  ALLOCATE (CID(xlen,ylen))
  status = NF_INQ_VARID(LSC_NCID,'cid',varid)
  status = NF_GET_VAR_INT(LSC_NCID,varid,CID)
  ALLOCATE (LAT(xlen,ylen))
  status = NF_INQ_VARID(LSC_NCID,'lat',varid)
  status = NF_GET_VAR_REAL(LSC_NCID,varid,LAT)
  ALLOCATE (LON(xlen,ylen))
  status = NF_INQ_VARID(LSC_NCID,'lon',varid)
  status = NF_GET_VAR_REAL(LSC_NCID,varid,LON)

END SUBROUTINE GET_GRID_NETCDF