[gradsusr] using EDEF to create an ensemble in Grads

Jennifer Adams jma at cola.iges.org
Tue Oct 5 12:58:46 EDT 2010


Well, with templated data sets, if GrADS can't find the file for your  
requested time, it just gives missing data. So you should make sure  
that all the files exist where they are expected to and that the names  
match up correctly with your template substitution strings in your  
DSET entry. Another check is to write a descriptor for one member,  
without the edef and substituting YSU for %e and make sure that works.
--Jennifer

On Oct 5, 2010, at 12:43 PM, Jeffrey Duda wrote:

> Jennifer,
> Here's what my actual control file looks like (not just the output  
> from 'q ctlinfo'):
>
> dset ^summer_%e_%y4-%m2-%d2_%h2:%n2.dat
> options  byteswapped template
> undef 1.e30
> title  OUTPUT FROM WRF V3.1.1 MODEL
> pdef 325 249 nps  163.030  375.118     -97.1640   23.0459518
> xdef 1021 linear -150.3   0.10381060
> ydef  489 linear   12.4   0.10381060
> zdef   26 linear 1 1
> tdef   61 linear 00Z14JUL2010      60MN
> edef    7 names YSU MYJ QNSE MYNN25 MYNN3 Pleim BouLac
> VARS   20
> U             26  0  x-wind component (m s-1)
> V             26  0  y-wind component (m s-1)
> W             26  0  z-wind component (m s-1)
> WW            26  0  omega (Pa s-1)
> PH            26  0  perturbation geopotential (m2 s-2)
> PHB           26  0  base-state geopotential (m2 s-2)
> P             26  0  perturbation pressure (Pa)
> PB            26  0  BASE STATE PRESSURE (Pa)
> Q2             1  0  QV at 2 M (kg kg-1)
> T2             1  0  TEMP at 2 M (K)
> PSFC           1  0  SFC PRESSURE (Pa)
> U10            1  0  U at 10 M (m s-1)
> V10            1  0  V at 10 M (m s-1)
> QVAPOR        26  0  Water vapor mixing ratio (kg kg-1)
> HGT            1  0  Terrain Height (m)
> RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
> RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
> PBLH           1  0  PBL HEIGHT (m)
> tc            26  0  Temperature (C)
> slp            1  0  Sea Levelp Pressure (hPa)
> ENDVARS
>
> You can see that 'options template' and the 'edef' line is in there.
>
> Jeff
>
> On Tue, Oct 5, 2010 at 11:34 AM, Jennifer Adams <jma at cola.iges.org>  
> wrote:
> Hi, Jeff --
> You appear to have neglected to add 'options template' in your  
> descriptor file. Also, if all your ensembles are the same length and  
> have the same start time, you can simplify your edef as  'edef 7  
> names ysu myj ... boulac '
> --Jennifer
>
>
>
> On Oct 5, 2010, at 12:15 PM, Jeffrey Duda wrote:
>
>> I wish to concatenate several separate sets of WRF output in which  
>> only the PBL scheme is changed into an ensemble.  I edited a  
>> control file accordingly, and made sure to place all binary files  
>> in the same directory.  However, when I open the control file and  
>> attempt to display any variables (whether the e dimension is  
>> varying or not), I get the "cannot contour grid - all undefined  
>> values" error.  What am I doing wrong?
>>
>> Here is the cltinfo:
>>
>> dset summer_%e_%y4-%m2-%d2_%h2:%n2.dat
>> title OUTPUT FROM WRF V3.1.1 MODEL
>> undef 1e+30
>> xdef 1021 linear -150.3 0.103811
>> ydef 489 linear 12.4 0.103811
>> zdef 26 linear 1 1
>> tdef 61 linear 00Z14JUL2010 60mn
>> edef 7
>> ysu 61 00Z14JUL2010
>> myj 61 00Z14JUL2010
>> qnse 61 00Z14JUL2010
>> mynn25 61 00Z14JUL2010
>> mynn3 61 00Z14JUL2010
>> pleim 61 00Z14JUL2010
>> boulac 61 00Z14JUL2010
>> endedef
>> vars 20
>> u  26  0  x-wind component (m s-1)
>> v  26  0  y-wind component (m s-1)
>> w  26  0  z-wind component (m s-1)
>> ww  26  0  omega (Pa s-1)
>> ph  26  0  perturbation geopotential (m2 s-2)
>> phb  26  0  base-state geopotential (m2 s-2)
>> p  26  0  perturbation pressure (Pa)
>> pb  26  0  BASE STATE PRESSURE (Pa)
>> q2  1  0  QV at 2 M (kg kg-1)
>> t2  1  0  TEMP at 2 M (K)
>> psfc  1  0  SFC PRESSURE (Pa)
>> u10  1  0  U at 10 M (m s-1)
>> v10  1  0  V at 10 M (m s-1)
>> qvapor  26  0  Water vapor mixing ratio (kg kg-1)
>> hgt  1  0  Terrain Height (m)
>> rainc  1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
>> rainnc  1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
>> pblh  1  0  PBL HEIGHT (m)
>> tc  26  0  Temperature (C)
>> slp  1  0  Sea Levelp Pressure (hPa)
>> endvars
>>
>> Jeff Duda
>>
>> -- 
>> Jeff Duda
>> Iowa State University
>> Meteorology Graduate Student
>> 3134 Agronomy Hall
>> www.meteor.iastate.edu/~jdduda
>> _______________________________________________
>> gradsusr mailing list
>> gradsusr at gradsusr.org
>> http://gradsusr.org/mailman/listinfo/gradsusr
>
> --
> Jennifer M. Adams
> IGES/COLA
> 4041 Powder Mill Road, Suite 302
> Calverton, MD 20705
> jma at cola.iges.org
>
>
>
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
>
>
>
> -- 
> Jeff Duda
> Iowa State University
> Meteorology Graduate Student
> 3134 Agronomy Hall
> www.meteor.iastate.edu/~jdduda
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr

--
Jennifer M. Adams
IGES/COLA
4041 Powder Mill Road, Suite 302
Calverton, MD 20705
jma at cola.iges.org



-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20101005/ae7f221b/attachment-0003.html 


More information about the gradsusr mailing list