[Met_help] [rt.rap.ucar.edu #78380] History for Mask files for Europe
John Halley Gotway via RT
met_help at ucar.edu
Fri Oct 28 11:05:10 MDT 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello!
I've succesfully installed and tested the MET software, but I've come to
a point with no solution - how to obtain a NetCDF mask file for a region
over Europe? What kind of GRIB data is needed? I've tried with standard
GDAS observation data and GFS GRIB2 data, but the only thing which
appeared was Segmentation fault (core dumped).
MAybe there is some issue with .poly file (I've attached it to email)
Many thanks
---
ppor. Łukasz Pajurek
Szefostwo Służby Hydrometeorologicznej SZ RP
Wydział Zabezpieczenia Informatycznego
2LT Łukasz Pajurek
Hydrometeorological Service of the Polish Armed Forces
IT Support Branch
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Mask files for Europe
From: John Halley Gotway
Time: Tue Oct 25 14:05:04 2016
Hello,
I'm see that you have compiled MET version 5.2 and would now like to
define
some new masking regions. Thanks for sending your sample POLAND.poly
polyline file. When I opened it on my Linux machine, I noticed that
it
contains "^M" characters at the end of each line. Those show up in
files
generated on Windows machines. When I ran your polyline file through
the
gen_vx_mask, the "^M" caused a runtime error since NetCDF doesn't like
including that character in the NetCDF variable name.
One easy way to strip them off is the "dos2unix" utility:
dos2unix POLAND.poly
Next, I grabbed a sample GFS file:
wget http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/
gfs.2016102512/gfs.t12z.pgrb2.0p50.f000
Then I used the MET regrid_data_plane tool to put some data on a
lat/lon
grid over Europe:
met-5.2/bin/regrid_data_plane gfs.t12z.pgrb2.0p50.f000 'latlon 100
100
25 0 0.5 0.5' gfs_euro.nc -field 'name="TMP"; level="Z2";'
Next, I ran the MET gen_vx_mask tool to apply your polyline to the
European
domain:
met-5.2/bin/gen_vx_mask gfs_euro.nc POLAND.poly POLAND_mask.nc
And then I ran the MET plot_data_plane tool to display the resulting
mask
field:
met-5.2/bin/plot_data_plane POLAND_mask.nc POLAND_mask.ps
'name="POLAND"; level="(*,*)";'
And I attached a PNG version of the resulting output.
As you can see, the mask is in roughly the right spot, but there are
obvious problems with the latitude and longitude values you've used to
define that mask for Poland.
Hope that helps.
Thanks,
John Halley Gotway
met_help at ucar.edu
On Tue, Oct 25, 2016 at 5:25 AM, l.pajurek via RT <met_help at ucar.edu>
wrote:
>
> Tue Oct 25 05:25:32 2016: Request 78380 was acted upon.
> Transaction: Ticket created by l.pajurek at meteo.mil.pl
> Queue: met_help
> Subject: Mask files for Europe
> Owner: Nobody
> Requestors: l.pajurek at meteo.mil.pl
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=78380 >
>
>
> Hello!
>
> I've succesfully installed and tested the MET software, but I've
come to
> a point with no solution - how to obtain a NetCDF mask file for a
region
> over Europe? What kind of GRIB data is needed? I've tried with
standard
> GDAS observation data and GFS GRIB2 data, but the only thing which
> appeared was Segmentation fault (core dumped).
> MAybe there is some issue with .poly file (I've attached it to
email)
>
> Many thanks
> ---
> ppor. Łukasz Pajurek
> Szefostwo Służby Hydrometeorologicznej SZ RP
> Wydział Zabezpieczenia Informatycznego
>
> 2LT Łukasz Pajurek
> Hydrometeorological Service of the Polish Armed Forces
> IT Support Branch
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #78380] Mask files for Europe
From: l.pajurek
Time: Fri Oct 28 02:18:49 2016
W dniu 2016-10-25 20:05, John Halley Gotway via RT napisał(a):
> Hello,
>
> I'm see that you have compiled MET version 5.2 and would now like to
> define
> some new masking regions. Thanks for sending your sample
POLAND.poly
> polyline file. When I opened it on my Linux machine, I noticed that
it
> contains "^M" characters at the end of each line. Those show up in
> files
> generated on Windows machines. When I ran your polyline file
through
> the
> gen_vx_mask, the "^M" caused a runtime error since NetCDF doesn't
like
> including that character in the NetCDF variable name.
>
> One easy way to strip them off is the "dos2unix" utility:
> dos2unix POLAND.poly
>
> Next, I grabbed a sample GFS file:
> wget http://www.ftp.ncep.noaa.gov/data/nccf/com/gfs/prod/
> gfs.2016102512/gfs.t12z.pgrb2.0p50.f000
>
> Then I used the MET regrid_data_plane tool to put some data on a
> lat/lon
> grid over Europe:
>
> met-5.2/bin/regrid_data_plane gfs.t12z.pgrb2.0p50.f000 'latlon
100
> 100
> 25 0 0.5 0.5' gfs_euro.nc -field 'name="TMP"; level="Z2";'
>
> Next, I ran the MET gen_vx_mask tool to apply your polyline to the
> European
> domain:
>
> met-5.2/bin/gen_vx_mask gfs_euro.nc POLAND.poly POLAND_mask.nc
>
> And then I ran the MET plot_data_plane tool to display the resulting
> mask
> field:
>
> met-5.2/bin/plot_data_plane POLAND_mask.nc POLAND_mask.ps
> 'name="POLAND"; level="(*,*)";'
>
> And I attached a PNG version of the resulting output.
>
> As you can see, the mask is in roughly the right spot, but there are
> obvious problems with the latitude and longitude values you've used
to
> define that mask for Poland.
>
> Hope that helps.
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
> On Tue, Oct 25, 2016 at 5:25 AM, l.pajurek via RT
<met_help at ucar.edu>
> wrote:
>
>>
>> Tue Oct 25 05:25:32 2016: Request 78380 was acted upon.
>> Transaction: Ticket created by l.pajurek at meteo.mil.pl
>> Queue: met_help
>> Subject: Mask files for Europe
>> Owner: Nobody
>> Requestors: l.pajurek at meteo.mil.pl
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=78380
>> >
>>
>>
>> Hello!
>>
>> I've succesfully installed and tested the MET software, but I've
come
>> to
>> a point with no solution - how to obtain a NetCDF mask file for a
>> region
>> over Europe? What kind of GRIB data is needed? I've tried with
>> standard
>> GDAS observation data and GFS GRIB2 data, but the only thing which
>> appeared was Segmentation fault (core dumped).
>> MAybe there is some issue with .poly file (I've attached it to
email)
>>
>> Many thanks
>> ---
>> ppor. Łukasz Pajurek
>> Szefostwo Służby Hydrometeorologicznej SZ RP
>> Wydział Zabezpieczenia Informatycznego
>>
>> 2LT Łukasz Pajurek
>> Hydrometeorological Service of the Polish Armed Forces
>> IT Support Branch
>>
Thank You very much for Your help! I've been using MET 5.0 vwersion,
but
now I've got this 5.2 version and it looks like it works.
I will try to make this mask file more suitable for my domain.
Thank You again fro Your cooperation.
Many thanks
---
ppor. Łukasz Pajurek
Szefostwo Służby Hydrometeorologicznej SZ RP
Wydział Zabezpieczenia Informatycznego
2LT Łukasz Pajurek
Hydrometeorological Service of the Polish Armed Forces
IT Support Branch
------------------------------------------------
More information about the Met_help
mailing list