[Met_help] [rt.rap.ucar.edu #92553] History for MODE error
John Halley Gotway via RT
met_help at ucar.edu
Fri Oct 11 09:47:53 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
I'm trying to run MODE on some hail data, and am getting an error I haven't seen before:
ERROR :
ERROR :
ERROR : ShapeData::shrink(const int) -> value too large ... new grid is empty
ERROR :
HDF5-DIAG: Error detected in HDF5 (1.8.14) thread 0:
#000: H5T.c line 1763 in H5Tclose(): not a datatype
major: Invalid arguments to routine
minor: Inappropriate type
terminate called after throwing an instance of 'netCDF::exceptions::NcHdfErr'
what(): NcHdfErr: An error was reported by the HDF5 layer.
file: ncFile.cpp line:15
Abort
Any idea what this means? It seems like it doesn't like something in the netCDF file, though I'm not sure what. I have some files in /scratch/WEEKLY/Tina
Thanks,
Tina
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: MODE error
From: John Halley Gotway
Time: Thu Oct 10 11:41:25 2019
Tina,
I moved this over to met-help. Whenever debugging MET's ability to
read gridded data, I always test using plot_data_plane:
/usr/local/met-8.1.1/bin/plot_data_plane \
news-e_MET_12_20190518_0000_0100.nc \
news-e_MET_12_20190518_0000_0100.ps \
'name="HAILCAST"; level="(0,*,*)";'
That also fails, but with a different error message:
ERROR : Box::set_lrbt() -> width and/or height negative!
Based on that, I suspect there's a problem with the definition of the
grid:
:Projection = "Lambert Conformal" ;
:scale_lat_1 = "60.00000" ;
:scale_lat_2 = "30.00000" ;
:lat_pin = 35.07053f ;
:lon_pin = -105.6247f ;
:x_pin = "0.000000" ;
:y_pin = "0.000000" ;
:lon_orient = -100.5997f ;
:d_km = "3.000000" ;
:r_km = "6371.200000" ;
:nx = 300LL ;
:ny = 300LL ;
Those lines for nx and ny = 300LL look a little suspicious to me.
They're usually listed as strings, like this:
:nx = "185" ;
:ny = "129 grid_points" ;
So I tried changing them to strings:
ncatted -a nx,global,o,c,"300" \
-a ny,global,o,c,"300" \
news-e_MET_12_20190518_0000_0100.nc \
-o news-e_MET_12_20190518_0000_0100_JHG.nc
And now plot_data_plane is happy again:
/usr/local/met-8.1.1/bin/plot_data_plane \
news-e_MET_12_20190518_0000_0100_JHG.nc \
news-e_MET_12_20190518_0000_0100.ps \
'name="HAILCAST"; level="(0,*,*)";'
And produces the attached image. Presumably, that's the underlying
problem in MODE as well.
Thanks,
John
------------------------------------------------
Subject: MODE error
From: Minna Win
Time: Thu Oct 10 18:39:14 2019
Hi Tina,
I'm going to assign this ticket to Randy, who handles our MODE and MDT
issues.
Thanks,
Minna
---------------
Minna Win
National Center for Atmospheric Research
Developmental Testbed Center
Phone: 303-497-8423
Fax: 303-497-8401
On Thu, Oct 10, 2019 at 11:32 AM John Halley Gotway via RT <
met_help at ucar.edu> wrote:
>
> Thu Oct 10 11:32:13 2019: Request 92553 was acted upon.
> Transaction: Ticket created by johnhg
> Queue: met_help
> Subject: MODE error
> Owner: johnhg
> Requestors: kalb at ucar.edu
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=92553 >
>
>
> I'm trying to run MODE on some hail data, and am getting an error I
> haven't seen before:
>
> ERROR :
> ERROR :
> ERROR : ShapeData::shrink(const int) -> value too large ... new
grid is
> empty
> ERROR :
> HDF5-DIAG: Error detected in HDF5 (1.8.14) thread 0:
> #000: H5T.c line 1763 in H5Tclose(): not a datatype
> major: Invalid arguments to routine
> minor: Inappropriate type
> terminate called after throwing an instance of
> 'netCDF::exceptions::NcHdfErr'
> what(): NcHdfErr: An error was reported by the HDF5 layer.
> file: ncFile.cpp line:15
> Abort
>
> Any idea what this means? It seems like it doesn't like something
in the
> netCDF file, though I'm not sure what. I have some files in
> /scratch/WEEKLY/Tina
>
> Thanks,
> Tina
>
------------------------------------------------
Subject: MODE error
From: Christina Kalb
Time: Thu Oct 10 20:06:18 2019
Thanks! I made that change, and MODE ran fine now :-D
Tina
On Thu, Oct 10, 2019 at 11:41 AM John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Tina,
>
> I moved this over to met-help. Whenever debugging MET's ability to
read
> gridded data, I always test using plot_data_plane:
>
> /usr/local/met-8.1.1/bin/plot_data_plane \
> news-e_MET_12_20190518_0000_0100.nc \
> news-e_MET_12_20190518_0000_0100.ps \
> 'name="HAILCAST"; level="(0,*,*)";'
>
> That also fails, but with a different error message:
>
> ERROR : Box::set_lrbt() -> width and/or height negative!
>
> Based on that, I suspect there's a problem with the definition of
the grid:
>
> :Projection = "Lambert Conformal" ;
> :scale_lat_1 = "60.00000" ;
> :scale_lat_2 = "30.00000" ;
> :lat_pin = 35.07053f ;
> :lon_pin = -105.6247f ;
> :x_pin = "0.000000" ;
> :y_pin = "0.000000" ;
> :lon_orient = -100.5997f ;
> :d_km = "3.000000" ;
> :r_km = "6371.200000" ;
> :nx = 300LL ;
> :ny = 300LL ;
>
> Those lines for nx and ny = 300LL look a little suspicious to me.
They're
> usually listed as strings, like this:
> :nx = "185" ;
> :ny = "129 grid_points" ;
> So I tried changing them to strings:
> ncatted -a nx,global,o,c,"300" \
> -a ny,global,o,c,"300" \
> news-e_MET_12_20190518_0000_0100.nc \
> -o news-e_MET_12_20190518_0000_0100_JHG.nc
>
> And now plot_data_plane is happy again:
> /usr/local/met-8.1.1/bin/plot_data_plane \
> news-e_MET_12_20190518_0000_0100_JHG.nc \
> news-e_MET_12_20190518_0000_0100.ps \
> 'name="HAILCAST"; level="(0,*,*)";'
>
> And produces the attached image. Presumably, that's the underlying
> problem in MODE as well.
>
> Thanks,
> John
>
>
>
>
>
>
--
*************************************************
Christina Kalb
Associate Scientist
National Center for Atmospheric Research
3450 Mitchell Lane
Boulder, CO 80301
kalb at ucar.edu
303-497-2819
*************************************************
------------------------------------------------
More information about the Met_help
mailing list