R/sidfex.osisaf.nc2rda.R

In helgegoessling/SIDFEx: Tools for the Sea Ice Drift Forecast Experiment (SIDFEx)

sidfex.osisaf.nc2rda <- function(filename, do.noSIC2NA = T){
  if(any(grepl("package:SVGAnnotation", search()))) detach("package:SVGAnnotation") 
  sic_fl = nc_open(filename)
  sic = ncvar_get(sic_fl, "ice_conc")
  sic_lon = ncvar_get(sic_fl, "lon")
  sic_lat = ncvar_get(sic_fl, "lat")
  nc_close(sic_fl)
  lonlati = sl.lonlat2D.c2i(lon.c=sic_lon,lat.c=sic_lat)
  if (do.noSIC2NA) {sic[sic==0]=NA}
  out = list(sic, lonlati)
  return(out)
}