to display NARCCAP data

Juan Matias Mendez Perez matias at TROPOSFERA.ATMOSFCU.UNAM.MX
Mon Jun 2 08:22:03 EDT 2008


Hello GrADS users,

I have NARCCAP data and users suggest NCL to open it, but I'd like to do
it with GrADS. I just tried to display using GrADS but I got this
error:

ga-> sdfopen uas_WRFP_2001010103.nc
Scanning self-describing file:  uas_WRFP_2001010103.nc
SDF file has no discernable X coordinate.

What does that mean? and How can I open it?.

I read GrADS documentation and I tried to construct a netcdf descriptor
file but without successful. I'd thank you so much if someone can help
me with this.

Best regards,
Matias

PS
This is my netcdf file header:

netcdf uas_WRFP_2001010103 {
dimensions:
        time = UNLIMITED ; // (11440 currently)
        yc = 109 ;
        xc = 134 ;
variables:
        char Lambert_Conformal ;
                Lambert_Conformal:grid_mapping_name =
"lambert_conformal_conic" ;
                Lambert_Conformal:standard_parallel = 30.f, 60.f ;
                Lambert_Conformal:longitude_of_central_meridian =
-97.f ;
                Lambert_Conformal:latitude_of_projection_origin =
47.5f ;
        double yc(yc) ;
                yc:long_name = "y-coordinate in Cartesian system" ;
                yc:standard_name = "projection_y_coordinate" ;
                yc:axis = "Y" ;
                yc:units = "m" ;
        double xc(xc) ;
                xc:long_name = "x-coordinate in Cartesian system" ;
                xc:standard_name = "projection_x_coordinate" ;
                xc:axis = "X" ;
                xc:units = "m" ;
        double lon(yc, xc) ;
                lon:long_name = "longitude" ;
                lon:standard_name = "longitude" ;
                lon:units = "degrees_east" ;
                lon:axis = "X" ;
        double lat(yc, xc) ;
                lat:long_name = "latitude" ;
                lat:standard_name = "latitude" ;
                lat:units = "degrees_north" ;
                lat:axis = "Y" ;
        double time(time) ;
                time:long_name = "time" ;
                time:standard_name = "time" ;
                time:units = "days since 1979-09-01 00:00:0.0" ;
                time:axis = "T" ;
                time:calendar = "gregorian" ;
        float uas(time, yc, xc) ;
                uas:long_name = "Zonal Surface Wind Speed" ;
                uas:standard_name = "eastward_wind" ;
                uas:units = "m s-1" ;
                uas:_FillValue = 1.e+20f ;
                uas:missing_values = 1.e+20f ;
                uas:coordinates = "lon lat" ;
                uas:grid_mapping = "Lambert_Conformal" ;
                uas:original_name = "U10 (U at 10 M) in wrfout" ;
                uas:original_units = "m s-1" ;

// global attributes:
                :institution = "PNNL (Pacific Northwest National
Laboratory, Richland, WA, USA)" ;
                :source = "PNNL/NCAR WRF(2006): atmosphere:WRF; sst/sea
ice:AMIP2; Land:Noah;; Convcetion:KF; WSM5; CAM Radiation" ;
                :project_id = "NARCCAP" ;
                :table_id = "Table 2" ;
                :realization = 1 ;
                :experiment_id = "present-day climate using NCEP/DOE
Reanalysis " ;
                :contact = "L. Ruby Leung (Ruby.Leung at pnl.gov)" ;
                :history = "Thu Jan 31 09:41:37 2008: ncrcat uas_2001.nc
uas_2002.nc uas_2003.nc uas_2004.nc uas_WRFP_2000010103.nc\n",
                        "Wed Jan 30 14:12:36 2008: ncrcat
uas_WRFP_2001a.nc uas_WRFP_2001b.nc uas_2001.nc\n",
                        "Wed Jan 30 14:12:34 2008: ncrcat
uas_WRFP_2001010103.nc uas_WRFP_2001020103.nc uas_WRFP_2001030103.nc
uas_WRFP_2001040103.nc uas_WRFP_2001050103.nc uas_WRFP_2001060103.nc
uas_WRFP_2001070103.nc uas_WRFP_2001080103.nc uas_WRFP_2001090103.nc
uas_WRFP_2001a.nc\n",
                        "subsetted from raw WRF output using
pullwrf_rain.f90" ;
                :title = "PNNL model output prepared for NARCCAP
present-day climate using NCEP/DOE Reanalysis" ;
                :nco_openmp_thread_number = 2 ;
                :Conventions = "CF-1.0" ;



More information about the gradsusr mailing list