Subsetting NLDAS variables

From: Noel Aloysius <noel.aloysius_at_nyahnyahspammersnyahnyah>
Date: Tue Mar 12 2013 - 10:11:45 MDT

Hello NCL-talk,

I want to extract a subset (lat 39.5 to 41.5 and lon -83.5 to -81.5) of
hourly precipitation for the period 1979-2012 from the NLDAS-2 OPeNDAP
archives.

I am not sure how to define the subset.

f1 = "http://hydro1.sci.gsfc.nasa.gov:80/dods/NLDAS_FORA0125_H.002"
exists = isfilepresent(f1)
print("File: <" + f1 + ">, existence: " + exists) ; returns the following

  (0) File: <
http://hydro1.sci.gsfc.nasa.gov:80/dods/NLDAS_FORA0125_H.002>, existence:
True

fin = addfile(f1, "r")

At this point, I am not sure how to define the subset of the variable
apcpsfc (precipitation hourly total (kg/m^2)

Thanks in advance for the help.

I am using NCL 6.1.2 DAP-enabled version.

Thanks,
Noel

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Mar 12 08:11:57 2013

This archive was generated by hypermail 2.1.8 : Wed Mar 13 2013 - 14:19:38 MDT