gen.scrip.grids.ncl

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

;; create SCRIP gridfiles for each model + UDEL


rcms = (/"CRCM","ECP2","HRM3","MM5I","RCM3","WRFG","UDEL"/)
nm = dimsizes(rcms)

do m=0,nm-1
  model = rcms(m)

  print(""+model)

  infile = "landmask."+model+".nc"
  
  fin = addfile(infile,"r")
  lat = fin->lat
  lon = fin->lon
  land = fin->landmask
  delete(fin)  

  outfile = "scrip.grid."+model+".nc"

  opt = True
  opt@ForceOverwrite = True

  if (model .eq. "UDEL") then
    rectilinear_to_SCRIP(outfile,lat,lon,opt)
  else
    curvilinear_to_SCRIP(outfile,lat,lon,opt)
  end if

  outfile = "scrip.mask."+model+".nc"
  opt@Mask2D = land

  if (model .eq. "UDEL") then
    rectilinear_to_SCRIP(outfile,lat,lon,opt)
  else
    curvilinear_to_SCRIP(outfile,lat,lon,opt)
  end if

  delete([/lat,lon,land,infile,outfile,opt/])

end do


;; Copyright 2013 Univ. Corp. for Atmos. Research
;; Author: Seth McGinnis, mcginnis@ucar.edu