[Met_help] [rt.rap.ucar.edu #85886] History for Mask collapsing in a pixel
John Halley Gotway via RT
met_help at ucar.edu
Tue Jul 9 12:03:26 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear Met Helpdesk,
I'm trying to create a mask using the command:
gen_vx_mask -type poly
/scratch4/NCEPDEV/ohd/Francesca.Viterbo/Ellicott_NWM/STAGEIV_cut_1h_2728052018/
ST4.2018052723.01h_cut.nc ELLICOTT_unix.poly prova_poly_ellicott.nc -v 10
The script runs but the result is a mask with only one pixel.
The watershed I want to add is 9.5 km2 on a grid of 3 km2 of regridded
stage IV data so this might be possible but I would expect at least two
pixels. I have created the mask using arcgis and I can attach the polyline
mask to this email thread if necessary and the netcdf result.
Let me know if you need more informations about that.
Thank you very much,
Francesca Viterbo
--
--
Francesca Viterbo, Ph.D.
Physical Sciences Division
NOAA Earth System Research Laboratory
R/PSD2, 325 Broadway
Boulder, CO 80305-3337
303-497-6702
francesca.viterbo at noaa.gov
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Mask collapsing in a pixel
From: John Halley Gotway
Time: Tue Jun 26 14:28:47 2018
Hello Francesca,
I see you have a question about defining a masking region for MET
using the
gen_vx_mask tool. I logged on to theia and took a look at your data.
I
reran your case and confirmed that the mask contains a single point.
I did
run the plot_data_plane tool to generate a picture and at least see
that
it's in the right spot (see attached mask.png).
In the beta version for met-7.1, we've added support for defining
masking
regions using shapefiles. So I went hunting around for a shapefile
corresponding to this region.
I pulled the Mid-Atlantic RFC from this site:
https://www.nohrsc.noaa.gov/gisdatasets/
Then I ran the gis_dump_dbf utility in MET and see this list of
descriptions:
gis_dump_dbf b_marfc.dbf | grep NAME | cut -d'=' -f2 | tr -d '\n'
Do any of these ring a bell?
Thanks,
John
"MTGN4" "NEVN6" "SHBN6" "CINN6" "CRTN6" "RCKN6" "UNDN6" "ESIN6"
"BAIN6"
"GNEN6" "CNON6" "WSRN6" "CMPN6" "CKLN6" "WVYN6" "CRNN6" "WCRN6"
"LDYN6"
"CMGN6" "TOWP1" "LWVP1" "TGRP1" "TGJP1" "TNKP1" "MHPP1" "MFDP1"
"MONP1"
"OFRP1" "CEDP1" "WBRP1" "LOYP1" "LYCP1" "SFFP1" "WESP1" "WILP1"
"BMBP1"
"LHVP1" "ROVP1" "SNNP1" "SBYP1" "BECP1" "BCHP1" "JRSP1" "PNCP1"
"MSBP1"
"KARP1" "HTVP1" "HERP1" "NPTP1" "CRVP1" "DMLP1" "LNCP1" "HARP1"
"LWSP1"
"MPLP1" "SMDP1" "PORP1" "HUJP1" "RTBP1" "CPHP1" "HGSP1" "SLYP1"
"MRTP1"
"SXTP1" "FAVM2" "FDKM2" "CWCM2" "SACM2" "HNKM2" "LUKM2" "SRVM2"
"DAWM2"
"BNMW2" "PORM2" "SHEW2" "CBEM2" "MBGW2" "PAWW2" "BRKM2" "KITM2"
"LEEV2"
"SPRW2" "MILW2" "FROV2" "RENV2" "COOV2" "STGV2" "FDBV2" "PETW2"
"MOFW2"
"FKLW2" "SFSV2" "CLPV2" "LYWV2" "PYAV2" "WAYV2" "JCKV2" "BVSV2"
"COWV2"
"RMDV2" "SVEV2" "CARV2" "CVGV2" "HRKV2" "BECV2" "LIRV2" "MTCV2"
"FRMV2"
"MTXV2" "BNNV2" "WALN6" "PEPN6" "CNNN6" "HVDN6" "FSHN6" "NVRN6"
"HLEN6"
"CKFN6" "CCRN6" "BRGN6" "BRYN6" "HWYP1" "MAWN4" "WNQN4" "MTMP1"
"PPTN4"
"MCPN4" "PPPN4" "SHOP1" "LODN4" "TKSN4" "BONN4" "LTFN4" "WHTP1"
"BOTN4"
"MNSP1" "PINN4" "BVDN4" "MILN4" "PRRP1" "CAMN4" "LHTP1" "RRTN4"
"STTN4"
"WNTP1" "BDKN4" "ALNP1" "MNVN4" "RGLN4" "BETP1" "BREP1" "BKWN4"
"TACN4"
"TREN4" "GRAP1" "LNGP1" "RDRP1" "BMRP1" "PTNP1" "PEBN4" "PADP1"
"CDFP1"
"VSTN6" "GCPW2" "WIBP1" "SPKP1" "LWBP1
On Mon, Jun 25, 2018 at 6:29 PM francesca.viterbo at noaa.gov via RT <
met_help at ucar.edu> wrote:
>
> Mon Jun 25 18:29:16 2018: Request 85886 was acted upon.
> Transaction: Ticket created by francesca.viterbo at noaa.gov
> Queue: met_help
> Subject: Mask collapsing in a pixel
> Owner: Nobody
> Requestors: francesca.viterbo at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85886 >
>
>
> Dear Met Helpdesk,
> I'm trying to create a mask using the command:
>
> gen_vx_mask -type poly
>
>
/scratch4/NCEPDEV/ohd/Francesca.Viterbo/Ellicott_NWM/STAGEIV_cut_1h_2728052018/
> ST4.2018052723.01h_cut.nc ELLICOTT_unix.poly prova_poly_ellicott.nc
-v 10
>
> The script runs but the result is a mask with only one pixel.
> The watershed I want to add is 9.5 km2 on a grid of 3 km2 of
regridded
> stage IV data so this might be possible but I would expect at least
two
> pixels. I have created the mask using arcgis and I can attach the
polyline
> mask to this email thread if necessary and the netcdf result.
> Let me know if you need more informations about that.
>
> Thank you very much,
>
> Francesca Viterbo
> --
> --
> Francesca Viterbo, Ph.D.
> Physical Sciences Division
> NOAA Earth System Research Laboratory
> R/PSD2, 325 Broadway
> Boulder, CO 80305-3337
> 303-497-6702
> francesca.viterbo at noaa.gov
>
>
------------------------------------------------
Subject: Mask collapsing in a pixel
From: francesca.viterbo at noaa.gov
Time: Wed Jun 27 11:10:44 2018
Dear John,
thank you very much for your help. I really appreciate you attempt to
reproduce my experiment in order to double check the results.
Unfortunately none of the names rings me a bell, but I'm generally
having a
very bad memory. Nevertheless I can attach the shapefile of the
watershed
attached to this email, if this can help you some way and try that on
the
7.1 MET version.
Thank you,
Francesca
On Tue, Jun 26, 2018 at 2:28 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hello Francesca,
>
> I see you have a question about defining a masking region for MET
using the
> gen_vx_mask tool. I logged on to theia and took a look at your
data. I
> reran your case and confirmed that the mask contains a single point.
I did
> run the plot_data_plane tool to generate a picture and at least see
that
> it's in the right spot (see attached mask.png).
>
> In the beta version for met-7.1, we've added support for defining
masking
> regions using shapefiles. So I went hunting around for a shapefile
> corresponding to this region.
>
> I pulled the Mid-Atlantic RFC from this site:
> https://www.nohrsc.noaa.gov/gisdatasets/
>
> Then I ran the gis_dump_dbf utility in MET and see this list of
> descriptions:
>
> gis_dump_dbf b_marfc.dbf | grep NAME | cut -d'=' -f2 | tr -d '\n'
>
>
> Do any of these ring a bell?
>
>
> Thanks,
> John
>
>
> "MTGN4" "NEVN6" "SHBN6" "CINN6" "CRTN6" "RCKN6" "UNDN6" "ESIN6"
"BAIN6"
> "GNEN6" "CNON6" "WSRN6" "CMPN6" "CKLN6" "WVYN6" "CRNN6" "WCRN6"
"LDYN6"
> "CMGN6" "TOWP1" "LWVP1" "TGRP1" "TGJP1" "TNKP1" "MHPP1" "MFDP1"
"MONP1"
> "OFRP1" "CEDP1" "WBRP1" "LOYP1" "LYCP1" "SFFP1" "WESP1" "WILP1"
"BMBP1"
> "LHVP1" "ROVP1" "SNNP1" "SBYP1" "BECP1" "BCHP1" "JRSP1" "PNCP1"
"MSBP1"
> "KARP1" "HTVP1" "HERP1" "NPTP1" "CRVP1" "DMLP1" "LNCP1" "HARP1"
"LWSP1"
> "MPLP1" "SMDP1" "PORP1" "HUJP1" "RTBP1" "CPHP1" "HGSP1" "SLYP1"
"MRTP1"
> "SXTP1" "FAVM2" "FDKM2" "CWCM2" "SACM2" "HNKM2" "LUKM2" "SRVM2"
"DAWM2"
> "BNMW2" "PORM2" "SHEW2" "CBEM2" "MBGW2" "PAWW2" "BRKM2" "KITM2"
"LEEV2"
> "SPRW2" "MILW2" "FROV2" "RENV2" "COOV2" "STGV2" "FDBV2" "PETW2"
"MOFW2"
> "FKLW2" "SFSV2" "CLPV2" "LYWV2" "PYAV2" "WAYV2" "JCKV2" "BVSV2"
"COWV2"
> "RMDV2" "SVEV2" "CARV2" "CVGV2" "HRKV2" "BECV2" "LIRV2" "MTCV2"
"FRMV2"
> "MTXV2" "BNNV2" "WALN6" "PEPN6" "CNNN6" "HVDN6" "FSHN6" "NVRN6"
"HLEN6"
> "CKFN6" "CCRN6" "BRGN6" "BRYN6" "HWYP1" "MAWN4" "WNQN4" "MTMP1"
"PPTN4"
> "MCPN4" "PPPN4" "SHOP1" "LODN4" "TKSN4" "BONN4" "LTFN4" "WHTP1"
"BOTN4"
> "MNSP1" "PINN4" "BVDN4" "MILN4" "PRRP1" "CAMN4" "LHTP1" "RRTN4"
"STTN4"
> "WNTP1" "BDKN4" "ALNP1" "MNVN4" "RGLN4" "BETP1" "BREP1" "BKWN4"
"TACN4"
> "TREN4" "GRAP1" "LNGP1" "RDRP1" "BMRP1" "PTNP1" "PEBN4" "PADP1"
"CDFP1"
> "VSTN6" "GCPW2" "WIBP1" "SPKP1" "LWBP1
>
>
>
> On Mon, Jun 25, 2018 at 6:29 PM francesca.viterbo at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Mon Jun 25 18:29:16 2018: Request 85886 was acted upon.
> > Transaction: Ticket created by francesca.viterbo at noaa.gov
> > Queue: met_help
> > Subject: Mask collapsing in a pixel
> > Owner: Nobody
> > Requestors: francesca.viterbo at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85886 >
> >
> >
> > Dear Met Helpdesk,
> > I'm trying to create a mask using the command:
> >
> > gen_vx_mask -type poly
> >
> > /scratch4/NCEPDEV/ohd/Francesca.Viterbo/Ellicott_
> NWM/STAGEIV_cut_1h_2728052018/
> > ST4.2018052723.01h_cut.nc ELLICOTT_unix.poly
prova_poly_ellicott.nc -v
> 10
> >
> > The script runs but the result is a mask with only one pixel.
> > The watershed I want to add is 9.5 km2 on a grid of 3 km2 of
regridded
> > stage IV data so this might be possible but I would expect at
least two
> > pixels. I have created the mask using arcgis and I can attach the
> polyline
> > mask to this email thread if necessary and the netcdf result.
> > Let me know if you need more informations about that.
> >
> > Thank you very much,
> >
> > Francesca Viterbo
> > --
> > --
> > Francesca Viterbo, Ph.D.
> > Physical Sciences Division
> > NOAA Earth System Research Laboratory
> > R/PSD2, 325 Broadway
> > Boulder, CO 80305-3337
> > 303-497-6702
> > francesca.viterbo at noaa.gov
> >
> >
>
>
--
--
Francesca Viterbo, Ph.D.
Physical Sciences Division
NOAA Earth System Research Laboratory
R/PSD2, 325 Broadway
Boulder, CO 80305-3337
303-497-6702
francesca.viterbo at noaa.gov
------------------------------------------------
Subject: Mask collapsing in a pixel
From: John Halley Gotway
Time: Wed Jun 27 13:39:49 2018
Randy,
Can you please take a look at this shapefile the Francesca sent. I
tried
running it through gen_vx_mask, but all the grid points were always
turned
on. Running it through the debugger reveals that we're reading
lat/lon
values which include... (1967090, -1630850).
Running gis_dump_shp produces the output listed below.
Any idea where the problem may lie?
Thanks,
John
file = "GlobalWatershedPoint.shp"
File Header ...
| file_code = 9994
| file_length = 128
| version = 1000
| shape_type = shape_type_point
|
| (x_min, x_max) = (1.63064e+06, 1.63064e+06)
| (y_min, y_max) = (1.9702e+06, 1.9702e+06)
| (z_min, z_max) = (0, 0)
| (m_min, m_max) = (0, 0)
Record Header ...
| record_number = 0 (0-based)
| content_length = 20
Record Data ...
| shape_type = shape_type_point
|
| (x, y) = (1.63064e+06, 1.9702e+06)
On Wed, Jun 27, 2018 at 11:11 AM francesca.viterbo at noaa.gov via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85886 >
>
> Dear John,
> thank you very much for your help. I really appreciate you attempt
to
> reproduce my experiment in order to double check the results.
>
> Unfortunately none of the names rings me a bell, but I'm generally
having a
> very bad memory. Nevertheless I can attach the shapefile of the
watershed
> attached to this email, if this can help you some way and try that
on the
> 7.1 MET version.
>
> Thank you,
> Francesca
>
>
> On Tue, Jun 26, 2018 at 2:28 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Hello Francesca,
> >
> > I see you have a question about defining a masking region for MET
using
> the
> > gen_vx_mask tool. I logged on to theia and took a look at your
data. I
> > reran your case and confirmed that the mask contains a single
point. I
> did
> > run the plot_data_plane tool to generate a picture and at least
see that
> > it's in the right spot (see attached mask.png).
> >
> > In the beta version for met-7.1, we've added support for defining
masking
> > regions using shapefiles. So I went hunting around for a
shapefile
> > corresponding to this region.
> >
> > I pulled the Mid-Atlantic RFC from this site:
> > https://www.nohrsc.noaa.gov/gisdatasets/
> >
> > Then I ran the gis_dump_dbf utility in MET and see this list of
> > descriptions:
> >
> > gis_dump_dbf b_marfc.dbf | grep NAME | cut -d'=' -f2 | tr -d '\n'
> >
> >
> > Do any of these ring a bell?
> >
> >
> > Thanks,
> > John
> >
> >
> > "MTGN4" "NEVN6" "SHBN6" "CINN6" "CRTN6" "RCKN6" "UNDN6" "ESIN6"
"BAIN6"
> > "GNEN6" "CNON6" "WSRN6" "CMPN6" "CKLN6" "WVYN6" "CRNN6" "WCRN6"
"LDYN6"
> > "CMGN6" "TOWP1" "LWVP1" "TGRP1" "TGJP1" "TNKP1" "MHPP1" "MFDP1"
"MONP1"
> > "OFRP1" "CEDP1" "WBRP1" "LOYP1" "LYCP1" "SFFP1" "WESP1" "WILP1"
"BMBP1"
> > "LHVP1" "ROVP1" "SNNP1" "SBYP1" "BECP1" "BCHP1" "JRSP1" "PNCP1"
"MSBP1"
> > "KARP1" "HTVP1" "HERP1" "NPTP1" "CRVP1" "DMLP1" "LNCP1" "HARP1"
"LWSP1"
> > "MPLP1" "SMDP1" "PORP1" "HUJP1" "RTBP1" "CPHP1" "HGSP1" "SLYP1"
"MRTP1"
> > "SXTP1" "FAVM2" "FDKM2" "CWCM2" "SACM2" "HNKM2" "LUKM2" "SRVM2"
"DAWM2"
> > "BNMW2" "PORM2" "SHEW2" "CBEM2" "MBGW2" "PAWW2" "BRKM2" "KITM2"
"LEEV2"
> > "SPRW2" "MILW2" "FROV2" "RENV2" "COOV2" "STGV2" "FDBV2" "PETW2"
"MOFW2"
> > "FKLW2" "SFSV2" "CLPV2" "LYWV2" "PYAV2" "WAYV2" "JCKV2" "BVSV2"
"COWV2"
> > "RMDV2" "SVEV2" "CARV2" "CVGV2" "HRKV2" "BECV2" "LIRV2" "MTCV2"
"FRMV2"
> > "MTXV2" "BNNV2" "WALN6" "PEPN6" "CNNN6" "HVDN6" "FSHN6" "NVRN6"
"HLEN6"
> > "CKFN6" "CCRN6" "BRGN6" "BRYN6" "HWYP1" "MAWN4" "WNQN4" "MTMP1"
"PPTN4"
> > "MCPN4" "PPPN4" "SHOP1" "LODN4" "TKSN4" "BONN4" "LTFN4" "WHTP1"
"BOTN4"
> > "MNSP1" "PINN4" "BVDN4" "MILN4" "PRRP1" "CAMN4" "LHTP1" "RRTN4"
"STTN4"
> > "WNTP1" "BDKN4" "ALNP1" "MNVN4" "RGLN4" "BETP1" "BREP1" "BKWN4"
"TACN4"
> > "TREN4" "GRAP1" "LNGP1" "RDRP1" "BMRP1" "PTNP1" "PEBN4" "PADP1"
"CDFP1"
> > "VSTN6" "GCPW2" "WIBP1" "SPKP1" "LWBP1
> >
> >
> >
> > On Mon, Jun 25, 2018 at 6:29 PM francesca.viterbo at noaa.gov via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Mon Jun 25 18:29:16 2018: Request 85886 was acted upon.
> > > Transaction: Ticket created by francesca.viterbo at noaa.gov
> > > Queue: met_help
> > > Subject: Mask collapsing in a pixel
> > > Owner: Nobody
> > > Requestors: francesca.viterbo at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85886
> >
> > >
> > >
> > > Dear Met Helpdesk,
> > > I'm trying to create a mask using the command:
> > >
> > > gen_vx_mask -type poly
> > >
> > > /scratch4/NCEPDEV/ohd/Francesca.Viterbo/Ellicott_
> > NWM/STAGEIV_cut_1h_2728052018/
> > > ST4.2018052723.01h_cut.nc ELLICOTT_unix.poly
prova_poly_ellicott.nc -v
> > 10
> > >
> > > The script runs but the result is a mask with only one pixel.
> > > The watershed I want to add is 9.5 km2 on a grid of 3 km2 of
regridded
> > > stage IV data so this might be possible but I would expect at
least two
> > > pixels. I have created the mask using arcgis and I can attach
the
> > polyline
> > > mask to this email thread if necessary and the netcdf result.
> > > Let me know if you need more informations about that.
> > >
> > > Thank you very much,
> > >
> > > Francesca Viterbo
> > > --
> > > --
> > > Francesca Viterbo, Ph.D.
> > > Physical Sciences Division
> > > NOAA Earth System Research Laboratory
> > > R/PSD2, 325 Broadway
> > > Boulder, CO 80305-3337
> > > 303-497-6702
> > > francesca.viterbo at noaa.gov
> > >
> > >
> >
> >
>
>
> --
> --
> Francesca Viterbo, Ph.D.
> Physical Sciences Division
> NOAA Earth System Research Laboratory
> R/PSD2, 325 Broadway
> Boulder, CO 80305-3337
> 303-497-6702
> francesca.viterbo at noaa.gov
>
>
------------------------------------------------
Subject: Mask collapsing in a pixel
From: John Halley Gotway
Time: Mon Jul 09 17:10:41 2018
Randy,
Can you please take a look at this issue from Francesca?
Thanks,
John
On Wed, Jun 27, 2018 at 1:39 PM John Halley Gotway <johnhg at ucar.edu>
wrote:
> Randy,
>
> Can you please take a look at this shapefile the Francesca sent. I
tried
> running it through gen_vx_mask, but all the grid points were always
turned
> on. Running it through the debugger reveals that we're reading
lat/lon
> values which include... (1967090, -1630850).
>
> Running gis_dump_shp produces the output listed below.
>
> Any idea where the problem may lie?
>
> Thanks,
> John
>
> file = "GlobalWatershedPoint.shp"
>
> File Header ...
> | file_code = 9994
> | file_length = 128
> | version = 1000
> | shape_type = shape_type_point
> |
> | (x_min, x_max) = (1.63064e+06, 1.63064e+06)
> | (y_min, y_max) = (1.9702e+06, 1.9702e+06)
> | (z_min, z_max) = (0, 0)
> | (m_min, m_max) = (0, 0)
>
> Record Header ...
> | record_number = 0 (0-based)
> | content_length = 20
>
> Record Data ...
> | shape_type = shape_type_point
> |
> | (x, y) = (1.63064e+06, 1.9702e+06)
>
>
>
> On Wed, Jun 27, 2018 at 11:11 AM francesca.viterbo at noaa.gov via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85886 >
>>
>> Dear John,
>> thank you very much for your help. I really appreciate you attempt
to
>> reproduce my experiment in order to double check the results.
>>
>> Unfortunately none of the names rings me a bell, but I'm generally
having
>> a
>> very bad memory. Nevertheless I can attach the shapefile of the
watershed
>> attached to this email, if this can help you some way and try that
on the
>> 7.1 MET version.
>>
>> Thank you,
>> Francesca
>>
>>
>> On Tue, Jun 26, 2018 at 2:28 PM, John Halley Gotway via RT <
>> met_help at ucar.edu> wrote:
>>
>> > Hello Francesca,
>> >
>> > I see you have a question about defining a masking region for MET
using
>> the
>> > gen_vx_mask tool. I logged on to theia and took a look at your
data. I
>> > reran your case and confirmed that the mask contains a single
point. I
>> did
>> > run the plot_data_plane tool to generate a picture and at least
see that
>> > it's in the right spot (see attached mask.png).
>> >
>> > In the beta version for met-7.1, we've added support for defining
>> masking
>> > regions using shapefiles. So I went hunting around for a
shapefile
>> > corresponding to this region.
>> >
>> > I pulled the Mid-Atlantic RFC from this site:
>> > https://www.nohrsc.noaa.gov/gisdatasets/
>> >
>> > Then I ran the gis_dump_dbf utility in MET and see this list of
>> > descriptions:
>> >
>> > gis_dump_dbf b_marfc.dbf | grep NAME | cut -d'=' -f2 | tr -d '\n'
>> >
>> >
>> > Do any of these ring a bell?
>> >
>> >
>> > Thanks,
>> > John
>> >
>> >
>> > "MTGN4" "NEVN6" "SHBN6" "CINN6" "CRTN6" "RCKN6" "UNDN6" "ESIN6"
"BAIN6"
>> > "GNEN6" "CNON6" "WSRN6" "CMPN6" "CKLN6" "WVYN6" "CRNN6" "WCRN6"
"LDYN6"
>> > "CMGN6" "TOWP1" "LWVP1" "TGRP1" "TGJP1" "TNKP1" "MHPP1" "MFDP1"
"MONP1"
>> > "OFRP1" "CEDP1" "WBRP1" "LOYP1" "LYCP1" "SFFP1" "WESP1" "WILP1"
"BMBP1"
>> > "LHVP1" "ROVP1" "SNNP1" "SBYP1" "BECP1" "BCHP1" "JRSP1" "PNCP1"
"MSBP1"
>> > "KARP1" "HTVP1" "HERP1" "NPTP1" "CRVP1" "DMLP1" "LNCP1" "HARP1"
"LWSP1"
>> > "MPLP1" "SMDP1" "PORP1" "HUJP1" "RTBP1" "CPHP1" "HGSP1" "SLYP1"
"MRTP1"
>> > "SXTP1" "FAVM2" "FDKM2" "CWCM2" "SACM2" "HNKM2" "LUKM2" "SRVM2"
"DAWM2"
>> > "BNMW2" "PORM2" "SHEW2" "CBEM2" "MBGW2" "PAWW2" "BRKM2" "KITM2"
"LEEV2"
>> > "SPRW2" "MILW2" "FROV2" "RENV2" "COOV2" "STGV2" "FDBV2" "PETW2"
"MOFW2"
>> > "FKLW2" "SFSV2" "CLPV2" "LYWV2" "PYAV2" "WAYV2" "JCKV2" "BVSV2"
"COWV2"
>> > "RMDV2" "SVEV2" "CARV2" "CVGV2" "HRKV2" "BECV2" "LIRV2" "MTCV2"
"FRMV2"
>> > "MTXV2" "BNNV2" "WALN6" "PEPN6" "CNNN6" "HVDN6" "FSHN6" "NVRN6"
"HLEN6"
>> > "CKFN6" "CCRN6" "BRGN6" "BRYN6" "HWYP1" "MAWN4" "WNQN4" "MTMP1"
"PPTN4"
>> > "MCPN4" "PPPN4" "SHOP1" "LODN4" "TKSN4" "BONN4" "LTFN4" "WHTP1"
"BOTN4"
>> > "MNSP1" "PINN4" "BVDN4" "MILN4" "PRRP1" "CAMN4" "LHTP1" "RRTN4"
"STTN4"
>> > "WNTP1" "BDKN4" "ALNP1" "MNVN4" "RGLN4" "BETP1" "BREP1" "BKWN4"
"TACN4"
>> > "TREN4" "GRAP1" "LNGP1" "RDRP1" "BMRP1" "PTNP1" "PEBN4" "PADP1"
"CDFP1"
>> > "VSTN6" "GCPW2" "WIBP1" "SPKP1" "LWBP1
>> >
>> >
>> >
>> > On Mon, Jun 25, 2018 at 6:29 PM francesca.viterbo at noaa.gov via RT
<
>> > met_help at ucar.edu> wrote:
>> >
>> > >
>> > > Mon Jun 25 18:29:16 2018: Request 85886 was acted upon.
>> > > Transaction: Ticket created by francesca.viterbo at noaa.gov
>> > > Queue: met_help
>> > > Subject: Mask collapsing in a pixel
>> > > Owner: Nobody
>> > > Requestors: francesca.viterbo at noaa.gov
>> > > Status: new
>> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85886
>> >
>> > >
>> > >
>> > > Dear Met Helpdesk,
>> > > I'm trying to create a mask using the command:
>> > >
>> > > gen_vx_mask -type poly
>> > >
>> > > /scratch4/NCEPDEV/ohd/Francesca.Viterbo/Ellicott_
>> > NWM/STAGEIV_cut_1h_2728052018/
>> > > ST4.2018052723.01h_cut.nc ELLICOTT_unix.poly
prova_poly_ellicott.nc
>> -v
>> > 10
>> > >
>> > > The script runs but the result is a mask with only one pixel.
>> > > The watershed I want to add is 9.5 km2 on a grid of 3 km2 of
regridded
>> > > stage IV data so this might be possible but I would expect at
least
>> two
>> > > pixels. I have created the mask using arcgis and I can attach
the
>> > polyline
>> > > mask to this email thread if necessary and the netcdf result.
>> > > Let me know if you need more informations about that.
>> > >
>> > > Thank you very much,
>> > >
>> > > Francesca Viterbo
>> > > --
>> > > --
>> > > Francesca Viterbo, Ph.D.
>> > > Physical Sciences Division
>> > > NOAA Earth System Research Laboratory
>> > > R/PSD2, 325 Broadway
>> > > Boulder, CO 80305-3337
>> > > 303-497-6702
>> > > francesca.viterbo at noaa.gov
>> > >
>> > >
>> >
>> >
>>
>>
>> --
>> --
>> Francesca Viterbo, Ph.D.
>> Physical Sciences Division
>> NOAA Earth System Research Laboratory
>> R/PSD2, 325 Broadway
>> Boulder, CO 80305-3337
>> 303-497-6702
>> francesca.viterbo at noaa.gov
>>
>>
------------------------------------------------
Subject: Mask collapsing in a pixel
From: Randy Bullock
Time: Wed Jul 11 13:09:55 2018
Hello Francesca -
Looking at the GlobalWatershed.shp file you sent us, I see numbers
like
these:
| | (1.63079e+06, 1.96716e+06)
| | (1.63077e+06, 1.96716e+06)
| | (1.63077e+06, 1.96717e+06)
| | (1.63074e+06, 1.96717e+06)
| | (1.63074e+06, 1.96718e+06)
| | (1.6307e+06, 1.96718e+06)
| | (1.6307e+06, 1.96719e+06)
| | (1.63067e+06, 1.96719e+06)
| | (1.63067e+06, 1.9672e+06)
| | (1.63065e+06, 1.9672e+06)
| | (1.63065e+06, 1.96721e+06)
| | (1.63064e+06, 1.96721e+06)
| | (1.63064e+06, 1.96722e+06)
| | (1.63063e+06, 1.96722e+06)
| | (1.63063e+06, 1.96723e+06)
| | (1.63062e+06, 1.96723e+06)
| | (1.63062e+06, 1.96724e+06)
| | (1.63061e+06, 1.96724e+06)
| | (1.63061e+06, 1.96725e+06)
| | (1.6306e+06, 1.96725e+06)
These are from the gis_dump_shp utility.
These are not latitudes and longitudes. Our GIS code asssumes that
the
numbers in a shp file are lats and lons. Could these numbers be
universal
transverse mercator coordinates?
In any case, our code is not set up to handle data that is not lat/lon
pairs. Is there any way you could obtain (or make) an equivalent file
with
the latitudes and longitudes?
Hope this helps.
Randy Bullock
On Wed, Jun 27, 2018 at 11:10 AM, francesca.viterbo at noaa.gov via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85886 >
>
> Dear John,
> thank you very much for your help. I really appreciate you attempt
to
> reproduce my experiment in order to double check the results.
>
> Unfortunately none of the names rings me a bell, but I'm generally
having a
> very bad memory. Nevertheless I can attach the shapefile of the
watershed
> attached to this email, if this can help you some way and try that
on the
> 7.1 MET version.
>
> Thank you,
> Francesca
>
>
> On Tue, Jun 26, 2018 at 2:28 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
> > Hello Francesca,
> >
> > I see you have a question about defining a masking region for MET
using
> the
> > gen_vx_mask tool. I logged on to theia and took a look at your
data. I
> > reran your case and confirmed that the mask contains a single
point. I
> did
> > run the plot_data_plane tool to generate a picture and at least
see that
> > it's in the right spot (see attached mask.png).
> >
> > In the beta version for met-7.1, we've added support for defining
masking
> > regions using shapefiles. So I went hunting around for a
shapefile
> > corresponding to this region.
> >
> > I pulled the Mid-Atlantic RFC from this site:
> > https://www.nohrsc.noaa.gov/gisdatasets/
> >
> > Then I ran the gis_dump_dbf utility in MET and see this list of
> > descriptions:
> >
> > gis_dump_dbf b_marfc.dbf | grep NAME | cut -d'=' -f2 | tr -d '\n'
> >
> >
> > Do any of these ring a bell?
> >
> >
> > Thanks,
> > John
> >
> >
> > "MTGN4" "NEVN6" "SHBN6" "CINN6" "CRTN6" "RCKN6" "UNDN6" "ESIN6"
"BAIN6"
> > "GNEN6" "CNON6" "WSRN6" "CMPN6" "CKLN6" "WVYN6" "CRNN6" "WCRN6"
"LDYN6"
> > "CMGN6" "TOWP1" "LWVP1" "TGRP1" "TGJP1" "TNKP1" "MHPP1" "MFDP1"
"MONP1"
> > "OFRP1" "CEDP1" "WBRP1" "LOYP1" "LYCP1" "SFFP1" "WESP1" "WILP1"
"BMBP1"
> > "LHVP1" "ROVP1" "SNNP1" "SBYP1" "BECP1" "BCHP1" "JRSP1" "PNCP1"
"MSBP1"
> > "KARP1" "HTVP1" "HERP1" "NPTP1" "CRVP1" "DMLP1" "LNCP1" "HARP1"
"LWSP1"
> > "MPLP1" "SMDP1" "PORP1" "HUJP1" "RTBP1" "CPHP1" "HGSP1" "SLYP1"
"MRTP1"
> > "SXTP1" "FAVM2" "FDKM2" "CWCM2" "SACM2" "HNKM2" "LUKM2" "SRVM2"
"DAWM2"
> > "BNMW2" "PORM2" "SHEW2" "CBEM2" "MBGW2" "PAWW2" "BRKM2" "KITM2"
"LEEV2"
> > "SPRW2" "MILW2" "FROV2" "RENV2" "COOV2" "STGV2" "FDBV2" "PETW2"
"MOFW2"
> > "FKLW2" "SFSV2" "CLPV2" "LYWV2" "PYAV2" "WAYV2" "JCKV2" "BVSV2"
"COWV2"
> > "RMDV2" "SVEV2" "CARV2" "CVGV2" "HRKV2" "BECV2" "LIRV2" "MTCV2"
"FRMV2"
> > "MTXV2" "BNNV2" "WALN6" "PEPN6" "CNNN6" "HVDN6" "FSHN6" "NVRN6"
"HLEN6"
> > "CKFN6" "CCRN6" "BRGN6" "BRYN6" "HWYP1" "MAWN4" "WNQN4" "MTMP1"
"PPTN4"
> > "MCPN4" "PPPN4" "SHOP1" "LODN4" "TKSN4" "BONN4" "LTFN4" "WHTP1"
"BOTN4"
> > "MNSP1" "PINN4" "BVDN4" "MILN4" "PRRP1" "CAMN4" "LHTP1" "RRTN4"
"STTN4"
> > "WNTP1" "BDKN4" "ALNP1" "MNVN4" "RGLN4" "BETP1" "BREP1" "BKWN4"
"TACN4"
> > "TREN4" "GRAP1" "LNGP1" "RDRP1" "BMRP1" "PTNP1" "PEBN4" "PADP1"
"CDFP1"
> > "VSTN6" "GCPW2" "WIBP1" "SPKP1" "LWBP1
> >
> >
> >
> > On Mon, Jun 25, 2018 at 6:29 PM francesca.viterbo at noaa.gov via RT
<
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Mon Jun 25 18:29:16 2018: Request 85886 was acted upon.
> > > Transaction: Ticket created by francesca.viterbo at noaa.gov
> > > Queue: met_help
> > > Subject: Mask collapsing in a pixel
> > > Owner: Nobody
> > > Requestors: francesca.viterbo at noaa.gov
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=85886
> >
> > >
> > >
> > > Dear Met Helpdesk,
> > > I'm trying to create a mask using the command:
> > >
> > > gen_vx_mask -type poly
> > >
> > > /scratch4/NCEPDEV/ohd/Francesca.Viterbo/Ellicott_
> > NWM/STAGEIV_cut_1h_2728052018/
> > > ST4.2018052723.01h_cut.nc ELLICOTT_unix.poly
prova_poly_ellicott.nc -v
> > 10
> > >
> > > The script runs but the result is a mask with only one pixel.
> > > The watershed I want to add is 9.5 km2 on a grid of 3 km2 of
regridded
> > > stage IV data so this might be possible but I would expect at
least two
> > > pixels. I have created the mask using arcgis and I can attach
the
> > polyline
> > > mask to this email thread if necessary and the netcdf result.
> > > Let me know if you need more informations about that.
> > >
> > > Thank you very much,
> > >
> > > Francesca Viterbo
> > > --
> > > --
> > > Francesca Viterbo, Ph.D.
> > > Physical Sciences Division
> > > NOAA Earth System Research Laboratory
> > > R/PSD2, 325 Broadway
> > > Boulder, CO 80305-3337
> > > 303-497-6702
> > > francesca.viterbo at noaa.gov
> > >
> > >
> >
> >
>
>
> --
> --
> Francesca Viterbo, Ph.D.
> Physical Sciences Division
> NOAA Earth System Research Laboratory
> R/PSD2, 325 Broadway
> Boulder, CO 80305-3337
> 303-497-6702
> francesca.viterbo at noaa.gov
>
>
------------------------------------------------
More information about the Met_help
mailing list