[gradsusr] Templating 3B42_daily.*nc TRMM files

[Matei Georgescu]

Thanks for your reply John.

The tdef number is appropriate for 31 days (number of days for the month), but i could not even plot anything for t = 1, so that wasn't the main problem.

It was indeed the control file. I've modified it to:
DSET ^3B42_daily.2005.07.%d2.6.nc
DTYPE netcdf
OPTIONS template
UNDEF -9999.9
xdef 1440 linear -179.875 0.25
ydef 400  linear -49.875 0.25
ZDEF 1 linear 1 1
TDEF 31 linear 00Z1JUL2005 1dy
VARS 1
hrf=>prate 1 t,y,x daily precipitation [mm]
ENDVARS

and it works fine now. Primary parameter that I changed was the variable name definition and ordering of dimensions.

thank you,
matt

----- 
Matei ("Matt") Georgescu 
Post-doctoral Fellow 
Program on Food Security and the Environment 
Stanford University 
Y2E2 Bldg - MC4205 
473 Via Ortega 
Stanford, CA 94305 
homepage - http://stanford.edu/~mateig 

----- Original Message -----
From: "John Huddleston" <Huddleston at cira.colostate.edu>
To: "GrADS Users Forum" <gradsusr at gradsusr.org>
Sent: Thursday, May 6, 2010 11:02:06 AM
Subject: Re: [gradsusr] Templating 3B42_daily.*nc TRMM files

Matt,

See http://grads.iges.org/grads/gadoc/templates.html

Your TDEF number looks to be too small.

How many total days are in all the files in your template? Change 31 to
that number.

John Huddleston, PhD


-----Original Message-----
From: gradsusr-bounces at gradsusr.org
[mailto:gradsusr-bounces at gradsusr.org] On Behalf Of Matei Georgescu
Sent: Thursday, May 06, 2010 11:43 AM
To: GrADS Users Forum
Subject: [gradsusr] Templating 3B42_daily.*nc TRMM files

Greetings all,

My question relates to the use and templating of netcdf daily TRMM
precipitation data (3B42), with opengrads 2.0.a7.oga.3. I've noticed a
few threads on this subject but haven't been able to locate anything
that helps me.

I can successfully open individual files. For example:
ga-> sdfopen 3B42_daily.2005.07.01.6.nc
Scanning self-describing file: 3B42_daily.2005.07.01.6.nc
SDF file 3B42_daily.2005.07.01.6.nc is open as file 1
LON set to 0 360
LAT set to -49.875 49.875
LEV set to 0 0
Time values set: 2005:7:1:0 2005:7:1:0
E set to 1 1
ga-> d hrf
Contouring: 0 to 330 interval 30

Everything has the proper appearance.

I can do this for several files ... but I need to open tens of files at
once, and need to template.

I created the following control file (control.ctl) based on users'
suggestions, to template all my netcdf files (to start, for just 1
month):

DSET ^3B42_daily.2005.07.%d2.6.nc
DTYPE netcdf
OPTIONS template
UNDEF -9999.0
xdef 1440 linear -179.875 0.25
ydef 400 linear -49.875 0.25
ZDEF 1 linear 1 1
TDEF 31 linear 00Z1JUL2005 1dy
VARS 1
hrf 1 z,x,y daily precipitation [mm]
ENDVARS

Upon opening:
ga-> sdfopen control.ctl
Scanning self-describing file: control.ctl
gaopnc error: nc_open failed to open file control.ctl
read_metadata: gaopnc failed
gadsdf: Couldn't ingest SDF metadata.

Or else:
ga-> open control.ctl
Scanning description file: control.ctl
Data file 3B42_daily.2005.07.%d2.6.nc is open as file 1
LON set to 0 360
LAT set to -50 49.75
LEV set to 1 1
Time values set: 2005:7:1:0 2005:7:1:0
E set to 1 1
ga-> d hrf
NetCDF Error (gancrow, nc_get_vara_double): NetCDF: Start+count exceeds
dimension bound
Data Request Error: Error for variable 'hrf'
Error ocurred at column 1
DISPLAY error: Invalid expression
Expression = hrf

No doubt I'm doing something silly ... I would appreciate any
suggestions.

thank you,
matt

----- Matei ("Matt") Georgescu
Post-doctoral Fellow
Program on Food Security and the Environment
Stanford University
Y2E2 Bldg - MC4205
473 Via Ortega
Stanford, CA 94305
homepage - http://stanford.edu/~mateig
_______________________________________________ gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr

_______________________________________________ gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr