[gradsusr] concatenating files

Mon Jul 11 16:03:53 EDT 2016

Hi everybody

I am a new Grads user.  Right now I am running simulations using NHRCM and GCM models to get downscaling data (temperature, precipitation, humidity, etc.) most them as a ctl or binary output files.

I should do bias corrections.  For instance for the temperature case I should compare or correlate values of observations versus simulation values.  So, I would like to know if you could help me with a script which allow  get the data or information (temperature, precipation, etc.)  from these ctl files, allows to sort this data (small value to the largest one), and graph them to get the correlation coefficent....and so on.

Thanks so much in advance

Kind regards

Rei

________________________________
From: gradsusr-bounces at gradsusr.org <gradsusr-bounces at gradsusr.org> on behalf of Frank Colby <Frank_Colby at uml.edu>
Sent: Monday, July 11, 2016 6:51:46 PM
To: GrADS Users Forum
Subject: Re: [gradsusr] concatenating files

Jennifer,

Thanks for the help  I didn't realize how much more meta data was included in the grib2 files.
I had found something online that seemed to suggest grib1 files worked the same way as grib2.
Now it makes sense that it didn't work out.

Frank

On 7/8/2016 9:18 AM, Jennifer M Adams wrote:
Hi, Frank —
You can’t concatenate grib1 ensemble files and expect gribmap to distinguish between them — there isn’t any ensemble metadata in a grib1 header. This is discussed in the documentation for ensembles under ‘how to organize your data’: http://cola.gmu.edu/grads/gadoc/ensembles.html#data<https://urldefense.proofpoint.com/v2/url?u=http-3A__cola.gmu.edu_grads_gadoc_ensembles.html-23data&d=CwMF-g&c=lqHimbpwJeF7VTDNof4ddl8H-RbXeAdbMI2MFE1TXqA&r=6tSfNfW9eok2dvPp1qvMgNtd3pXLhNg0FMEbTXyDrbI&m=QPqA9k9AYEgTmQJHxaUoZ_ddngbq5ay2qvumidlIA3A&s=rfyXUHXXIldohh5hjQltp04HTFH4kkXgAH3khibQUIk&e=>. For GFS ensembles, I use a descriptor file like this:

dset ^gfsens.%e.%iy4%im2%id2%ih2.fh%f3.grib2
dtype grib2
title NCEP 5D Ensemble Forecast Initialized 00z08jul2016 6-hourly out to 192 hours
index ^gfsens.2016070800.map
options template pascals
undef -9.99e33
xdef 360 linear  0.0 1.0
ydef 181 linear -90.0 1.0
zdef  10 levels 100000 92500 85000 70000 50000 25000 20000 10000 5000 1000
tdef  33 linear 00z08jul2016 6hr
edef 22
*name length initialtime <type,pert or deriv>
avg  33  00z08jul2016   0
c00  33  00z08jul2016  1,0
p01  33  00z08jul2016  3,1
p02  33  00z08jul2016  3,2
p03  33  00z08jul2016  3,3
p04  33  00z08jul2016  3,4
p05  33  00z08jul2016  3,5
p06  33  00z08jul2016  3,6
p07  33  00z08jul2016  3,7
p08  33  00z08jul2016  3,8
p09  33  00z08jul2016  3,9
p10  33  00z08jul2016  3,10
p11  33  00z08jul2016  3,11
p12  33  00z08jul2016  3,12
p13  33  00z08jul2016  3,13
p14  33  00z08jul2016  3,14
p15  33  00z08jul2016  3,15
p16  33  00z08jul2016  3,16
p17  33  00z08jul2016  3,17
p18  33  00z08jul2016  3,18
p19  33  00z08jul2016  3,19
p20  33  00z08jul2016  3,20
endedef
vars
....
endvars

On Jul 7, 2016, at 5:21 PM, Frank Colby <Frank_Colby at uml.edu<mailto:Frank_Colby at uml.edu>> wrote:

Hi GradsUsers,

Does anyone have an easy way to handle GFS Ensemble files?  I have
output from ensemble runs that have separate files for each ensemble
member at each forecast hour.  I would like to put the ensemble members
together into a single file.  I used cat to put two files together as a
test, and wgrib (these are grib1 files) shows that both ensemble members
are indeed there.  However, when I use grib2ctl and gribmap to make the
ctl and index files, and open the file, it only has the first ensemble
member.   It's as if grib2ctl never read that there was a second
ensemble member.

What am I missing here?

Thanks,
Frank Colby
UMass Lowell

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

--
Jennifer Miletta Adams
Center for Ocean-Land-Atmosphere Studies (COLA)
George Mason University

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

***************Este mensaje (y sus adjuntos), en adelante "mensaje", ha sido enviado exclusivamente a su(s) destinatario(s) y es confidencial. Si usted recibe este mensaje por error, por favor bórrelo y comunique inmediatamente al remitente. Toda utilización o publicación, total o parcial, queda prohibida salvo autorización expresa. La UTP no podrá ser considerada responsable si el mensaje ha sido modificado y/o utilizado sin autorización. This message (and any attachments), are confidential intended solely for the people whose addresses appear. If you have received this message by error, please delete it and immediately notify the sender. Any use, dissemination or disclosure, either whole or partial, without formal approval is prohibited. The UTP will not therefore be liable for the message if modified and/or used without approval.**************
***************Este mensaje (y sus adjuntos), en adelante "mensaje", ha sido enviado exclusivamente a su(s) destinatario(s) y es confidencial. Si usted recibe este mensaje por error, por favor bórrelo y comunique inmediatamente al remitente. Toda utilización o publicación, total o parcial, queda prohibida salvo autorización expresa. La UTP no podrá ser considerada responsable si el mensaje ha sido modificado y/o utilizado sin autorización. This message (and any attachments), are confidential intended solely for the people whose addresses appear. If you have received this message by error, please delete it and immediately notify the sender. Any use, dissemination or disclosure, either whole or partial, without formal approval is prohibited. The UTP will not therefore be liable for the message if modified and/or used without approval.**************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20160711/3e63cb8a/attachment-0001.html 