[gradsusr] Reading CMAQ netcdf files in GrADS

[Huddleston, John]

Yes, I do. I am busy right now. I will get back to you. John

________________________________
From: arlindo.dasilva at gmail.com [arlindo.dasilva at gmail.com] On Behalf Of Arlindo da Silva [dasilva at alum.mit.edu]
Sent: Thursday, December 15, 2011 12:54 PM
To: Huddleston, John
Cc: GrADS Users Forum
Subject: Reading CMAQ netcdf files in GrADS

Hi John,

    Something tells me that you might know this. I have these CMAQ netcdf files from the DISCOVER-AQ campaign (see ncdump below). Do you happen to know the PDEF line that goes with these? I have the grid metrics in a separate file, I could always write some code to remap these, but if grads can do it on the fly it would be good enough for the exploratory data analysis that I need to perform.

   Thanks,

       Arlindo

-----

netcdf UMA_AOD {
dimensions:
TSTEP = UNLIMITED ; // (48 currently)
DATE-TIME = 2 ;
LAY = 1 ;
VAR = 1 ;
ROW = 79 ;
COL = 112 ;
variables:
int TFLAG(TSTEP, VAR, DATE-TIME) ;
TFLAG:units = "<YYYYDDD,HHMMSS>" ;
TFLAG:long_name = "TFLAG           " ;
TFLAG:var_desc = "Timestep-valid flags:  (1) YYYYDDD or (2) HHMMSS                                " ;
float AOD(TSTEP, LAY, ROW, COL) ;
AOD:long_name = "AOD             " ;
AOD:units = "unitless        " ;
AOD:var_desc = "AOD by Mie in VIS file                                                          " ;

// global attributes:
:IOAPI_VERSION = "$Id: @(#) ioapi library version 3.0 $                                           " ;
:EXEC_ID = "????????????????                                                                " ;
:FTYPE = 1 ;
:CDATE = 2011182 ;
:CTIME = 164809 ;
:WDATE = 2011182 ;
:WTIME = 164809 ;
:SDATE = 2011182 ;
:STIME = 130000 ;
:TSTEP = 10000 ;
:NTHIK = 1 ;
:NCOLS = 112 ;
:NROWS = 79 ;
:NLAYS = 1 ;
:NVARS = 1 ;
:GDTYP = 2 ;
:P_ALP = 33. ;
:P_BET = 45. ;
:P_GAM = -97. ;
:XCENT = -97. ;
:YCENT = 40. ;
:XORIG = 984000. ;
:YORIG = -384000. ;
:XCELL = 12000. ;
:YCELL = 12000. ;
:VGTYP = 1 ;
:VGTOP = 200.f ;
:VGLVLS = 1.f, 0.9959855f ;
:GDNAM = "AQF_CONUS_5x    " ;
:UPNAM = "setup           " ;
:VAR-LIST = "AOD             " ;
:FILEDESC = "hourly 1-layer cross-point Visibility Information

--
Arlindo da Silva
dasilva at alum.mit.edu<mailto:dasilva at alum.mit.edu>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20111215/01438183/attachment-0003.html