[gradsusr] gradsusr Digest, Vol 6, Issue 29
jl sa
esjluis at yahoo.com.mx
Mon Aug 23 18:48:36 EDT 2010
--- El lun 23-ago-10, gradsusr-request at gradsusr.org <gradsusr-request at gradsusr.org> escribió:
De: gradsusr-request at gradsusr.org <gradsusr-request at gradsusr.org>
Asunto: gradsusr Digest, Vol 6, Issue 29
A: gradsusr at gradsusr.org
Fecha: lunes, 23 de agosto de 2010, 13:16
Send gradsusr mailing list submissions to
gradsusr at gradsusr.org
To subscribe or unsubscribe via the World Wide Web, visit
http://gradsusr.org/mailman/listinfo/gradsusr
or, via email, send a message with subject or body 'help' to
gradsusr-request at gradsusr.org
You can reach the person managing the list at
gradsusr-owner at gradsusr.org
When replying, please edit your Subject line so it is more specific
than "Re: Contents of gradsusr digest..."
Today's Topics:
1. Re: ensemble reforecast runs (Wesley Ebisuzaki)
2. Re: ensemble reforecast runs (Daniele Mastrangelo)
3.RE: lats4d: 2nd try (Matt Brewer)
4. Re: Isobars: draw "H" and "L" (Michael Fiorino)
----------------------------------------------------------------------
Message: 1
Date: Mon, 23 Aug 2010 10:43:08 -0400
From: Wesley Ebisuzaki <Wesley.Ebisuzaki at noaa.gov>
Subject: Re: [gradsusr] ensemble reforecast runs
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID: <4C7288FC.7010800 at noaa.gov>
Content-Type: text/plain; charset=UTF-8; format=flowed
Daniele,
g2ctl assumes that the length of the forecast is the same for all
ensemble members.
If not, you can do either
(1) You could try using the longest ensemble as the 1st ensemble.
(It may work.)
(2) You have to set the TDEF by hand.
Wesley Ebisuzaki
Daniele Mastrangelo wrote:
> Dear GrADS users,
>
> I'd like to use GrADS v 2.0.a8 to evaluate a set of 21 reforecast
> ensemble runs. These runs are obtained running the same global model
> initialized at the same hour and date of 21 different years.
>
> I obtained a single ctl file through "g2ctl.pl -ens" (g2ctl.pl v
> 0.0.4o) but, when I used gribmap to produce the corresponding idx
> file, I got the following error message:
>
> Open Error: ensemble 1 (e2) extends beyond the time axis limits
> The data file was not opened.
>
> I see that the problem is that these runs cover different periods;
> moreover, I get the same error when I try to open two of the 21
> different reforecast runs in a normal GrADS session (as file 1 and
> file 2). Therefore, my question is: can GrADS treat a set reforecast
> runs used as ensemble runs?
> If it cannot directly do that, is there a way to bypass the problem?
>
> Thank you,
> Daniele
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> gradsusr mailing list
> gradsusr at gradsusr.org
> http://gradsusr.org/mailman/listinfo/gradsusr
>
------------------------------
Message: 2
Date: Mon, 23 Aug 2010 16:02:59 +0000 (GMT)
From: Daniele Mastrangelo <daniele_mastrangelo at yahoo.it>
Subject: Re: [gradsusr] ensemble reforecast runs
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID: <878904.28876.qm at web28509.mail.ukl.yahoo.com>
Content-Type: text/plain; charset="utf-8"
Hi Wesley,
thanks for your reply.
My previous email was not totally clear. Actually, I have 21 ensemble members
with the same forecast length; these members are produced initializing the model
at the same hour, day, and month (but) of 21 different years (12 UTC 29 aug
1989, 12 UTC 29 aug 1990, ... , 12 UTC 29 aug 2009).
I tried, however, to apply something similar to your point 2): I aggregated 2
members (e1, e2) in a single ctl file and changed by hand the year of the second
member in the EDEF section, making it the same as the first member
...
tdef 67 linear 12Z29aug1989 12hr
zdef 1 linear 1 1
vars 8
...
ENDVARS
EDEF 2
e1 67 12Z29aug1989
e2 67 12Z29aug1989 ---1989 instead of 1990---
ENDEDEF
Gribmap ran with no error messages but then, in GrADS, I could not use the e2
fields: "Cannot contour grid - all undefined values"
Is there some other way I could try?
Best regards,
Daniele
________________________________
Da: Wesley Ebisuzaki <Wesley.Ebisuzaki at noaa.gov>
A: GrADS Users Forum <gradsusr at gradsusr.org>
Inviato: Lun 23 agosto 2010, 16:43:08
Oggetto: Re: [gradsusr] ensemble reforecast runs
Daniele,
g2ctl assumes that the length of the forecast is the same for all
ensemble members.
If not, you can do either
(1) You could try using the longest ensemble as the 1st ensemble.
(It may work.)
(2) You have to set the TDEF by hand.
Wesley Ebisuzaki
Daniele Mastrangelo wrote:
> Dear GrADS users,
>
> I'd like to use GrADS v 2.0.a8 to evaluate a set of 21 reforecast
> ensemble runs. These runs are obtained running the same global model
> initialized at the same hour and date of 21 different years.
>
> I obtained a single ctl file through "g2ctl.pl -ens" (g2ctl.pl v
> 0.0.4o) but, when I used gribmap to produce the corresponding idx
> file, I got the following error message:
>
> Open Error: ensemble 1 (e2) extends beyond the time axis limits
> The data file was not opened.
>
> I see that the problem is that these runs cover different periods;
> moreover, I get the same error when I try to open two of the 21
> different reforecast runs in a normal GrADS session (as file 1 and
> file 2). Therefore, my question is: can GrADS treat a set reforecast
> runs used as ensemble runs?
> If it cannot directly do that, is there a way to bypass the problem?
>
> Thank you,
> Daniele
>
> ------------------------------------------------------------------------
>
> _______________________________________________
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100823/60d66e97/attachment-0001.html
------------------------------
Message: 3
Date: Mon, 23 Aug 2010 11:01:31 -0700 (PDT)
From: Matt Brewer <mcbrewer83 at yahoo.com>
Subject: [gradsusr] lats4d: 2nd try
To: GrADS Users Forum <gradsusr at gradsusr.org>
Message-ID: <734216.86942.qm at web56302.mail.re3.yahoo.com>
Content-Type: text/plain; charset="us-ascii"
Folks
Didn't get any help the first time I tried the gradsusr list with this problem,
so I thought I would give it another try. Currently, I have been using grads to
convert GRIB to netCDF files using lats4d. Basically, when I convert from GRIB
to netCDF, it skews the data. I input something like this in the command line:
lats4d -i narr-a_221_19790101_0000_000.sub -o narr_2005 -ftype ctl -lon -160 -90
-lat 25 65 -time 0z1jan2005 21z31dec2005 -mxtimes 2920 -levs 925mb -xvars
ugrdprs vgrdprs ugrd10m vgrd10m msletmsl prmslmsl tmp10m tmpprs
That will leaves me with geo height data at 925mb. Notice how I ask for data
from 160 to 90W and from 25 to 65N. However, when I run the above code, it ends
up giving me data from 160 to 89.6667W and from 25 to 65.333N (The data in the 3
hourly NCEP NARR dataset is every third degree, so it seems to just add on an
extra line of gridpoints, so I get a 122x212 matrix rather than 121x211). So on
the first latitude line of 25N, all the points are correct, until the very last
point at 25N 89.667W. Here, the program takes the value of the point that should
be at 25.333N 160W, and sticks it at 25N 89.667W. An thus, on the 25.333N array,
all the values are shifted one place to the left from where they should be. And
the last 2 grid points of that line are taken from the 2 most left points on the
next line up...and so on. I hope that makes sense. Please let me know if you
have any ideas or if you know anyone that might know the answer or if there is a
better way to convert GRIB to netCDF, though I would really like to make this
work with GrADS! Thank you so much for any help. By they way, I am using
version 2.0.a8.oga.1 on a linux machine and I got lats4d from downloading it off
surgeforce with the 2.0.a8.oga.1 bundle Here is my ctl file, as I have
suspicions that the answer lies in the PDEF part....I just don't know how. I
also included below the ctl file what happens when I run lats4d.sh in verbose
mode as it shows that I truely am just asking for 160-90W and 25-65N. Thank
you!!!
Regards,
Matt Brewer
dset ^narr-a_221_%y4%m2%d2_%h200_000.sub.grb
index ^narr-a_221_19790101_0000_000.sub.grb.idx
undef 9.999E+20
title narr-a_221_19790101_0000_000.sub.grb
* produced by grib2ctl v0.9.12.5p46
dtype grib 221
options template
pdef 349 277 lcc 1 -145.5 1 1 50 50 -107 32463 32463
xdef 1111 linear -250 0.333333
ydef 247 linear 8 0.333333
tdef 90584 linear 00Z01jan1979 3hr
* z has 3 levels, for prs
zdef 3 levels
925 850 500
vars 9
HGTprs 3 7,100,0 ** (profile) Geopotential height [gpm]
MSLETmsl 0 130,102,0 ** mean-sea level Mean sea level pressure (ETA model) [Pa]
PRMSLmsl 0 2,102,0 ** mean-sea level Pressure reduced to MSL [Pa]
TMPprs 3 11,100,0 ** (profile) Temp. [K]
TMP10m 0 11,105,10 ** 10 m above ground Temp. [K]
UGRDprs 3 33,100,0 ** (profile) u wind [m/s]
UGRD10m 0 33,105,10 ** 10 m above ground u wind [m/s]
VGRDprs 3 34,100,0 ** (profile) v wind [m/s]
VGRD10m 0 34,105,10 ** 10 m above ground v wind [m/s]
ENDVARS
lats4d in verbose mode
3 drizzle% lats4d.sh -i narr-a_221_19790101_0000_000.sub -o narrtest -ftype ctl
-v -lon -160 -90 -lat 25 65 -time 0z1jan2005 21z1jan2005 -levs 925mb -xvars
ugrdprs vgrdprs ugrd10m vgrd10m msletmsl prmslmsl tmp10m tmpprs
/bin/ls: cannot access 1.*: No such file or directory
/opt/opengrads/grads -blc 'run lats4d.gs -q -i narr-a_221_19790101_0000_000.sub
-o narrtest -ftype ctl -v -lon -160 -90 -lat 25 65 -time 0z1jan2005 21z1jan2005
-levs 925mb -xvars ugrdprs vgrdprs ugrd10m vgrd10m msletmsl prmslmsl tmp10m
tmpprs '
Welcome to the OpenGrADS Bundle Distribution
--------------------------------------------
For additional information enter "grads -h".
Starting "/opt/opengrads/Linux/Versions/2.0.a8.oga.1/i686/grads -blc 'run
lats4d.gs -q -i narr-a_221_19790101_0000_000.sub -o narrtest -ftype ctl -v -lon
-160 -90 -lat 25 65 -time 0z1jan2005 21z1jan2005 -levs 925mb -xvars ugrdprs
vgrdprs ugrd10m vgrd10m msletmsl prmslmsl tmp10m tmpprs ' " ...
Grid Analysis and Display System (GrADS) Version 2.0.a8.oga.1
Copyright (c) 1988-2010 by Brian Doty and the
Institute for Global Environment and Society (IGES)
GrADS comes with ABSOLUTELY NO WARRANTY
See file COPYRIGHT for more information
Config: v2.0.a8.oga.1 little-endian readline printim grib2 netcdf hdf4-sds hdf5
opendap-grids,stn athena geotiff shapefile
Issue 'q config' command for more information.
Loading User Defined Extensions table
</opt/opengrads/Linux/Versions/2.0.a8.oga.1/i686/gex/udxt> ... ok.
GX Package Initialization: Size = 11 8.5
Running in Batch mode
lats4d: Version 2.1.5 of 28 June 2010
No hardcopy metafile open
All files closed; all defined objects released;
All GrADS attributes have been reinitialized
Notice: Implied interpolation for file narr-a_221_19790101_0000_000.sub
Interpolation will be performed on any data displayed from this file
lats4d: Data file is
File 1 : narr-a_221_19790101_0000_000.sub.grb
Descriptor: narr-a_221_19790101_0000_000.sub.ctl
Binary: narr-a_221_%y4%m2%d2_%h200_000.sub.grb
Type = Gridded
Xsize = 1111 Ysize = 247 Zsize = 3 Tsize = 90584 Esize = 1
Number of Variables = 9
hgtprs 3 7 ** (profile) Geopotential height [gpm]
msletmsl 0 130 ** mean-sea level Mean sea level pressure (ETA model)
[Pa]
prmslmsl 0 2 ** mean-sea level Pressure reduced to MSL [Pa]
tmpprs 3 11 ** (profile) Temp. [K]
tmp10m 0 11 ** 10 m above ground Temp. [K]
ugrdprs 3 33 ** (profile) u wind [m/s]
ugrd10m 0 33 ** 10 m above ground u wind [m/s]
vgrdprs 3 34 ** (profile) v wind [m/s]
vgrd10m 0 34 ** 10 m above ground v wind [m/s]
lats4d: Dimension environment file same as data file
lats4d: excluding variable msletmsl
lats4d: excluding variable prmslmsl
lats4d: excluding variable tmpprs
lats4d: excluding variable tmp10m
lats4d: excluding variable ugrdprs
lats4d: excluding variable ugrd10m
lats4d: excluding variable vgrdprs
lats4d: excluding variable vgrd10m
lats4d: creating LATS PARAMETER TABLE file .grads.lats.table
lats4d: time range: 00Z01JAN2005 21Z01JAN2005 by 1, delta t: 3 hourly
lats4d: vertical levels: 925mb
lats4d: using PRESSURE for vertical coordinate
lats4d: latitudinal range: 25 65
lats4d: longitudinal range: -160 -90
lats4d: upper air variables: hgtprs
lats4d: writing to coards file narrtest.nc on 00Z01JAN2005
lats4d: writing to coards file narrtest.nc on 03Z01JAN2005
lats4d: writing to coards file narrtest.nc on 06Z01JAN2005
lats4d: writing to coards file narrtest.nc on 09Z01JAN2005
lats4d: writing to coards file narrtest.nc on 12Z01JAN2005
lats4d: writing to coards file narrtest.nc on 15Z01JAN2005
lats4d: writing to coards file narrtest.nc on 18Z01JAN2005
lats4d: writing to coards file narrtest.nc on 21Z01JAN2005
lats4d: created coards file narrtest.nc
lats4d: exiting from GrADS...
No hardcopy metafile open
GX package terminated
otra forma de crear archivos netcdf a partir de grib es con el programa VCDAT el cual es para graficar archivos netcdf ,pero puedes abrir el ctl desde el vcdat y de ahi definir la variable a desplegar (o incluso mas variables) y de ahi guardas como archivo netcdf, listo la instalacion del VCDAT es sencilla. disculapa es traduccion.
another
way to create NetCDF files from the program grib is VCDAT which is to
plot NetCDF files, but you can open the ctl from vcdat and from there
define the variable to display (or even more variables) and then save
it as NetCDF file, ready installation of the VCDAT is simple. apology is translation
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100823/8b0d7bf8/attachment-0001.html
------------------------------
Message: 4
Date: Mon, 23 Aug 2010 22:18:17 +0000
From: Michael Fiorino <Michael.Fiorino at noaa.gov>
Subject: Re: [gradsusr] Isobars: draw "H" and "L"
To: GrADS Users Forum <gradsusr at gradsusr.org>
Cc: Daniele Gandini <daniele.gandini at arpa.piemonte.it>
Message-ID: <4C72F3A9.4030507 at noaa.gov>
Content-Type: text/plain; charset="iso-8859-1"
Daniele et al., my mfhilo function is still not working correctly, but I am working on
it...hope to have a better version in the next couple of weeks that i'll put on
opengrads.org. sorry I can't be more helpful. best /r mike
On 8/23/10 12:22 PM, Daniele Gandini wrote:
> Dear Jonas,
>
> in GrADS 1.9 you can use clhilo.f and philo2.gs created by Mike Fiorino as
> User Defined Functions.
>
> In GrADS 2.* there is mfhilo; I haven't still used it but here is
> the link:
>
> http://opengrads.org/doc/udxt/libmf/
>
> Good luck
>
> Daniele
>
> ----- Messaggio originale -----
> Da: "Jonas Kaufmann"<me at j0nes.de>
> A: "GrADS Users Forum"<gradsusr at gradsusr.org>
> Inviato: Luned?, 23 agosto 2010 13:56:01
> Oggetto: [gradsusr] Isobars: draw "H" and "L"
>
> Hello,
>
> I am trying to generate isobar maps from GFS output. So far everything
> is looking fine, but I have one questions: I have not found a way to
> automatically add "H" and "L" marks in the center of a high- or low
> pressure system.
>
> Is there a way to do this?
>
> Thank you in advance!
>
> Best regards,
> Jonas Kaufmann
> _______________________________________________
> 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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100823/efbac109/attachment.html
------------------------------
_______________________________________________
gradsusr mailing list
gradsusr at gradsusr.org
http://gradsusr.org/mailman/listinfo/gradsusr
End of gradsusr Digest, Vol 6, Issue 29
***************************************
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://gradsusr.org/pipermail/gradsusr/attachments/20100823/acc9e502/attachment-0003.html
More information about the gradsusr
mailing list