[Met_help] [rt.rap.ucar.edu #71310] History for interpolate stage iv precip to wrf grid
John Halley Gotway via RT
met_help at ucar.edu
Fri Apr 3 15:01:35 MDT 2015
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
I’m trying to interpolate stage IV rainfall to my wrf grid for ensemble cases to see how results compare when I interpolated all ensemble member wrfout rainfall to the stage IV domain. I had success interpolating the wrfout to the stage IV grid, but am getting the error copygb: unsupported output grid ”255 when I try to interpolate stage IV to my wrf grid like so,
/home/brburgha/UPPV2.2/bin/copygb.exe -xg”255 3 327 339 25751 -106189 8 -99698 4000 4000 0 64 30000 60000” ST4.2013052001.01h ST4.2013052001.01h_lambert.grb
I’ve followed the example from the website, but i’m not sure why i’m getting the unsupported grid type 255. Any ideas?
doing a wgrib -V on my UPP converted wrfout gives:
rec 1:0:date 2013051500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0) grid=255 sfc 0-0hr acc:
APCP=Total precipitation [kg/m^2]
timerange 4 P1 0 P2 0 TimeU 1 nx 327 ny 339 GDS grid 3 num_in_ave 0 missing 0
center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
Lambert Conf: Lat1 25.751000 Lon1 -106.189000 Lov -101.000000
Latin1 60.000000 Latin2 30.000000 LatSP 0.000000 LonSP 0.000000
North Pole (327 x 339) Dx 4.000000 Dy 4.000000 scan 64 mode 136
min/max data 0 0 num bits 0 BDS_Ref 0 DecScale 0 BinScale 0
Thanks,
Brock
Brock Burghardt
Atmospheric Science PhD Student
Dept. of Geosciences
Texas Tech University
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: interpolate stage iv precip to wrf grid
From: John Halley Gotway
Time: Fri Apr 03 13:02:48 2015
Hello Brock,
I see that you're getting a runtime error from copygb when converting
to a
lambert grid. When I run similar commands, I get good results...
# Pull a sample ST4 file
wget
ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/hourly/prod/nam_pcpn_anal.20150401/ST4.2015040112.01h.gz
gunzip ST4.2015040112.01h.gz
# Run copygb to regrid
copygb -xg"255 3 327 339 25751 -106189 8 -99698 4000 4000 0 64 30000
60000" ST4.2015040112.01h ST4.2015040112.01h_lambert.grb
# Run plot_data_plane to plot copygb input/output
plot_data_plane ST4.2015040112.01h ST4.2015040112.01h.ps 'name="APCP";
level="A1";'
plot_data_plane ST4.2015040112.01h_lambert.grb
ST4.2015040112.01h_lambert.ps
'name="APCP"; level="A1";'
The resulting images are attached.
Now why are you getting that error? I'm able to reproduce your error
by
running that same copygb command *WITHOUT* the quotes around the grid
definition:
copygb -xg255 3 327 339 25751 -106189 8 -99698 4000 4000 0 64 30000
60000
ST4.2015040112.01h ST4.2015040112.01h_lambert.grb
copygb: unsupported output grid 255
My guess is that's what's happening here, and that this is really a
scripting problem rather than a copygb problem. If you're calling
copygb
from a script, look at the command to see if you have quotes around
the
grid definition. If not, add them and try again. Either single or
double
quotes will work fine.
Hope that helps.
Thanks,
John Halley Gotway
met_help at ucar.edu
On Fri, Apr 3, 2015 at 10:26 AM, Burghardt, Brock via RT
<met_help at ucar.edu>
wrote:
>
> Fri Apr 03 10:26:15 2015: Request 71310 was acted upon.
> Transaction: Ticket created by brock.burghardt at ttu.edu
> Queue: met_help
> Subject: interpolate stage iv precip to wrf grid
> Owner: Nobody
> Requestors: brock.burghardt at ttu.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71310 >
>
>
> Hello,
>
> I’m trying to interpolate stage IV rainfall to my wrf grid for
ensemble
> cases to see how results compare when I interpolated all ensemble
member
> wrfout rainfall to the stage IV domain. I had success interpolating
the
> wrfout to the stage IV grid, but am getting the error copygb:
unsupported
> output grid ”255 when I try to interpolate stage IV to my wrf
grid like
> so,
>
> /home/brburgha/UPPV2.2/bin/copygb.exe -xg”255 3 327 339 25751
-106189 8
> -99698 4000 4000 0 64 30000 60000” ST4.2013052001.01h
> ST4.2013052001.01h_lambert.grb
>
> I’ve followed the example from the website, but i’m not sure why i’m
> getting the unsupported grid type 255. Any ideas?
>
> doing a wgrib -V on my UPP converted wrfout gives:
>
> rec 1:0:date 2013051500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
> grid=255 sfc 0-0hr acc:
> APCP=Total precipitation [kg/m^2]
> timerange 4 P1 0 P2 0 TimeU 1 nx 327 ny 339 GDS grid 3 num_in_ave
0
> missing 0
> center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
> Lambert Conf: Lat1 25.751000 Lon1 -106.189000 Lov -101.000000
> Latin1 60.000000 Latin2 30.000000 LatSP 0.000000 LonSP
0.000000
> North Pole (327 x 339) Dx 4.000000 Dy 4.000000 scan 64 mode
136
> min/max data 0 0 num bits 0 BDS_Ref 0 DecScale 0 BinScale 0
>
> Thanks,
> Brock
>
>
> Brock Burghardt
> Atmospheric Science PhD Student
> Dept. of Geosciences
> Texas Tech University
>
>
>
>
>
------------------------------------------------
Subject: interpolate stage iv precip to wrf grid
From: John Halley Gotway
Time: Fri Apr 03 13:08:57 2015
Brock,
Also, I wanted to mention that we have a great new enhancement that
will be
included in the next release of MET, version 5.1, due out in late
August.
Rather than requiring the user to regrid their data to a common grid,
we've
added that capability into MET.
In version 5.1, you'll can pass your model GRIB file and StageIV file
directly to the MET grid-to-grid tools, such as Grid-Stat and MODE.
To
verify on your model domain, you'd set "to_grid = FCST" in the config
file. To verify on the StageIV domain, you'd set "to_grid = OBS" in
the
config file. And you can select what regridding interpolation method
should be used.
Hopefully that'll make the MET tools much easier to use.
Thanks,
John
On Fri, Apr 3, 2015 at 1:02 PM, John Halley Gotway <johnhg at ucar.edu>
wrote:
> Hello Brock,
>
> I see that you're getting a runtime error from copygb when
converting to a
> lambert grid. When I run similar commands, I get good results...
>
> # Pull a sample ST4 file
> wget
>
ftp://ftpprd.ncep.noaa.gov/pub/data/nccf/com/hourly/prod/nam_pcpn_anal.20150401/ST4.2015040112.01h.gz
> gunzip ST4.2015040112.01h.gz
>
> # Run copygb to regrid
> copygb -xg"255 3 327 339 25751 -106189 8 -99698 4000 4000 0 64
30000
> 60000" ST4.2015040112.01h ST4.2015040112.01h_lambert.grb
>
> # Run plot_data_plane to plot copygb input/output
> plot_data_plane ST4.2015040112.01h ST4.2015040112.01h.ps
'name="APCP";
> level="A1";'
> plot_data_plane ST4.2015040112.01h_lambert.grb
> ST4.2015040112.01h_lambert.ps 'name="APCP"; level="A1";'
>
> The resulting images are attached.
>
> Now why are you getting that error? I'm able to reproduce your
error by
> running that same copygb command *WITHOUT* the quotes around the
grid
> definition:
> copygb -xg255 3 327 339 25751 -106189 8 -99698 4000 4000 0 64 30000
60000
> ST4.2015040112.01h ST4.2015040112.01h_lambert.grb
> copygb: unsupported output grid 255
>
> My guess is that's what's happening here, and that this is really a
> scripting problem rather than a copygb problem. If you're calling
copygb
> from a script, look at the command to see if you have quotes around
the
> grid definition. If not, add them and try again. Either single or
double
> quotes will work fine.
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
>
> On Fri, Apr 3, 2015 at 10:26 AM, Burghardt, Brock via RT <
> met_help at ucar.edu> wrote:
>
>>
>> Fri Apr 03 10:26:15 2015: Request 71310 was acted upon.
>> Transaction: Ticket created by brock.burghardt at ttu.edu
>> Queue: met_help
>> Subject: interpolate stage iv precip to wrf grid
>> Owner: Nobody
>> Requestors: brock.burghardt at ttu.edu
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=71310 >
>>
>>
>> Hello,
>>
>> I’m trying to interpolate stage IV rainfall to my wrf grid for
ensemble
>> cases to see how results compare when I interpolated all ensemble
member
>> wrfout rainfall to the stage IV domain. I had success interpolating
the
>> wrfout to the stage IV grid, but am getting the error copygb:
unsupported
>> output grid ”255 when I try to interpolate stage IV to my wrf
grid like
>> so,
>>
>> /home/brburgha/UPPV2.2/bin/copygb.exe -xg”255 3 327 339 25751
-106189 8
>> -99698 4000 4000 0 64 30000 60000” ST4.2013052001.01h
>> ST4.2013052001.01h_lambert.grb
>>
>> I’ve followed the example from the website, but i’m not sure why
i’m
>> getting the unsupported grid type 255. Any ideas?
>>
>> doing a wgrib -V on my UPP converted wrfout gives:
>>
>> rec 1:0:date 2013051500 APCP kpds5=61 kpds6=1 kpds7=0 levels=(0,0)
>> grid=255 sfc 0-0hr acc:
>> APCP=Total precipitation [kg/m^2]
>> timerange 4 P1 0 P2 0 TimeU 1 nx 327 ny 339 GDS grid 3
num_in_ave 0
>> missing 0
>> center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
>> Lambert Conf: Lat1 25.751000 Lon1 -106.189000 Lov -101.000000
>> Latin1 60.000000 Latin2 30.000000 LatSP 0.000000 LonSP
0.000000
>> North Pole (327 x 339) Dx 4.000000 Dy 4.000000 scan 64 mode
136
>> min/max data 0 0 num bits 0 BDS_Ref 0 DecScale 0 BinScale 0
>>
>> Thanks,
>> Brock
>>
>>
>> Brock Burghardt
>> Atmospheric Science PhD Student
>> Dept. of Geosciences
>> Texas Tech University
>>
>>
>>
>>
>>
>
------------------------------------------------
Subject: interpolate stage iv precip to wrf grid
From: Burghardt, Brock
Time: Fri Apr 03 13:58:15 2015
Hi John,
Thank you for pointing that out, I figured it was something very
simple. In fact I realized the real issue, I had typed out the command
in TextEdit on mac and forgot it now defaults to rich text formatting
which automatically converts quotes to curly quotes, which i’m
guessing the unix platform just ignores. When I copied your example
with the normal quotes it works fine and when I forced textedit to
default to .txt instead of .rft it kept the normal quotes.
copygb -xg"255 3 327 339 25751 -106189 8 -99698 4000 4000 0 64 30000
60000" ST4.2013052001.01h ST4.2013052001.01h_lambert.grb #this works
copygb -xg”255 3 327 339 25751 -106189 8 -99698 4000 4000 0 64 30000
60000” ST4.2013052001.01h ST4.2013052001.01h_lambert.grb #this does
not
That’ll be a nice feature in the next release, i’ll keep an eye out
later this summer for the new download.
Thanks again, i’ll likely be sending another email to the help mail
later today or monday, but you can close this issue.
Brock
On Apr 3, 2015, at 2:08 PM, John Halley Gotway via RT
<met_help at ucar.edu<mailto:met_help at ucar.edu>> wrote:
home/brburgha/UPPV2.2/bin/copygb.exe -xg”255 3 327 339 25751 -106189 8
-99698 4000 4000 0 64 30000 60000” ST4.2013052001.01h
ST4.2013052001.01h_lambert.grb
------------------------------------------------
More information about the Met_help
mailing list