[Met_help] [rt.rap.ucar.edu #83497] History for regrid_data_plane question
John Halley Gotway via RT
met_help at ucar.edu
Fri Jan 5 15:47:44 MST 2018
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
To whom it may concern,
I was trying to run a simple script that takes the netcdf output file of
pcp_combine and then uses it as input to regrid_data_plane, with the script
pasted below. The results of the script is pasted below that.
It seems to get getting hung up on getting field "name="APCP"; level="A6";
I know I should have ftp'd the input files but could not find the
instructions from earlier threads, if you could point me to the anonymous
ftp instructions again I'll be sure to include them if need be.
Thanks for your help,
Brad Diehl
#!/bin/sh
###
############
source /u/Brad.Diehl/met_env/tut_env
in_file1='/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/hrrrcogb3p018010415qpfff12'
in_file2='/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/hrrrcogb3p018010415qpfff06'
pcp_combine -subtract $in_file1 12 $in_file2 06
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
regrid_data_plane /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc \
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/urma_co \
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/re_grid_out.nc \
-method BUDGET \
-width 2 \
-v 4 \
-field 'name="APCP"; level="A6";'
echo after regrid
Brad.Diehl at g14a2: /mdlstat/save/usr/Brad.Diehl/pcp_hrrr$ ./pcp_hrrr_test.sh
DEBUG 1: Reading input file:
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/hrrrcogb3p018010415qpfff12
DEBUG 1: Reading input file:
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/hrrrcogb3p018010415qpfff06
DEBUG 1: Writing output file: /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
DEBUG 1: Reading data file: /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_NcMet".
DEBUG 4: Lambert Conformal Grid Data:
DEBUG 4: scale_lat_1 = 38.5
DEBUG 4: scale_lat_2 = 38.5
DEBUG 4: lat_pin = 21.1381
DEBUG 4: lon_pin = 122.72
DEBUG 4: x_pin = 0
DEBUG 4: y_pin = 0
DEBUG 4: lon_orient = 97.5
DEBUG 4: d_km = 3
DEBUG 4: r_km = 6371.2
DEBUG 4: nx = 1799
DEBUG 4: ny = 1059
DEBUG 2: Input grid: Projection: Lambert Conformal Nx: 1799 Ny: 1059
Lat_LL: 21.138 Lon_LL: 122.720 Lon_orient: 97.500 Alpha: 4203.494 Cone:
0.623 Bx: 899.1693 By: 3198.9984
DEBUG 3: Use the grid defined by file
"/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/urma_co".
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb2".
DEBUG 4:
DEBUG 4: Lambert Conformal Grid Data:
DEBUG 4: scale_lat_1: 25
DEBUG 4: scale_lat_2: 25
DEBUG 4: lat_pin: 20.192
DEBUG 4: lon_pin: 121.554
DEBUG 4: x_pin: 0
DEBUG 4: y_pin: 0
DEBUG 4: lon_orient: 95
DEBUG 4: d_km: 2.5397
DEBUG 4: r_km: 6371.2
DEBUG 4: nx: 2145
DEBUG 4: ny: 1377
DEBUG 4:
DEBUG 3: Grid Definition: Projection: Lambert Conformal Nx: 2145 Ny: 1377
Lat_LL: 20.192 Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone:
0.423 Bx: 1088.0030 By: 5483.6617
DEBUG 2: Output grid: Projection: Lambert Conformal Nx: 2145 Ny: 1377
Lat_LL: 20.192 Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082 Cone:
0.423 Bx: 1088.0030 By: 5483.6617
DEBUG 2: Interpolation options: method = BUDGET, width = 2, vld_thresh = 0.5
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object of
type "FileType_NcMet".
ERROR :
ERROR : process_data_file() -> trouble getting field "name="APCP";
level="A6";" from file "/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc"
ERROR :
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: regrid_data_plane question
From: John Halley Gotway
Time: Thu Jan 04 15:25:55 2018
Hi Brad,
I see that you're trying to run pcp_combine followed by
regrid_data_plane.
I suspect you just need to modify your script slightly. pcp_combine
writes
out a NetCDF file and you should take a look at the output variable
that it
writes:
ncdump -h /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
I'm guessing the data variable is named "APCP_06".
When you pass this file to regrid_data_plane, you need to tell it the
name
of NetCDF variable to process:
regrid_data_plane out.nc regrid_out.nc -field 'name="APCP_06";
level="(*,*)";'
Unfortunately, there's a difference in how you specify the name and
level
for GRIB versus NetCDF files and I'm guessing that's the source of the
confusion here.
FYI, instructions for posting to our anonymous ftp site can be found
here:
https://dtcenter.org/met/users/support/met_help.php#ftp
Thanks,
John Halley Gotway
On Thu, Jan 4, 2018 at 3:17 PM, Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> Thu Jan 04 15:17:58 2018: Request 83497 was acted upon.
> Transaction: Ticket created by brad.diehl at noaa.gov
> Queue: met_help
> Subject: regrid_data_plane question
> Owner: Nobody
> Requestors: brad.diehl at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83497 >
>
>
> To whom it may concern,
>
> I was trying to run a simple script that takes the netcdf output
file of
> pcp_combine and then uses it as input to regrid_data_plane, with the
script
> pasted below. The results of the script is pasted below that.
> It seems to get getting hung up on getting field "name="APCP";
level="A6";
>
> I know I should have ftp'd the input files but could not find the
> instructions from earlier threads, if you could point me to the
anonymous
> ftp instructions again I'll be sure to include them if need be.
>
> Thanks for your help,
>
> Brad Diehl
>
>
> #!/bin/sh
> ###
> ############
> source /u/Brad.Diehl/met_env/tut_env
>
>
> in_file1='/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/
> hrrrcogb3p018010415qpfff12'
> in_file2='/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/
> hrrrcogb3p018010415qpfff06'
>
>
> pcp_combine -subtract $in_file1 12 $in_file2 06
> /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
>
> regrid_data_plane /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc \
> /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/urma_co \
> /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/re_grid_out.nc \
> -method BUDGET \
> -width 2 \
> -v 4 \
> -field 'name="APCP"; level="A6";'
>
> echo after regrid
>
>
>
>
>
>
> Brad.Diehl at g14a2: /mdlstat/save/usr/Brad.Diehl/pcp_hrrr$
> ./pcp_hrrr_test.sh
> DEBUG 1: Reading input file:
> /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/hrrrcogb3p018010415qpfff12
> DEBUG 1: Reading input file:
> /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/hrrrcogb3p018010415qpfff06
> DEBUG 1: Writing output file:
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
> DEBUG 1: Reading data file:
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_NcMet".
> DEBUG 4: Lambert Conformal Grid Data:
> DEBUG 4: scale_lat_1 = 38.5
> DEBUG 4: scale_lat_2 = 38.5
> DEBUG 4: lat_pin = 21.1381
> DEBUG 4: lon_pin = 122.72
> DEBUG 4: x_pin = 0
> DEBUG 4: y_pin = 0
> DEBUG 4: lon_orient = 97.5
> DEBUG 4: d_km = 3
> DEBUG 4: r_km = 6371.2
> DEBUG 4: nx = 1799
> DEBUG 4: ny = 1059
> DEBUG 2: Input grid: Projection: Lambert Conformal Nx: 1799 Ny: 1059
> Lat_LL: 21.138 Lon_LL: 122.720 Lon_orient: 97.500 Alpha: 4203.494
Cone:
> 0.623 Bx: 899.1693 By: 3198.9984
> DEBUG 3: Use the grid defined by file
> "/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/urma_co".
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_Gb2".
> DEBUG 4:
> DEBUG 4: Lambert Conformal Grid Data:
> DEBUG 4: scale_lat_1: 25
> DEBUG 4: scale_lat_2: 25
> DEBUG 4: lat_pin: 20.192
> DEBUG 4: lon_pin: 121.554
> DEBUG 4: x_pin: 0
> DEBUG 4: y_pin: 0
> DEBUG 4: lon_orient: 95
> DEBUG 4: d_km: 2.5397
> DEBUG 4: r_km: 6371.2
> DEBUG 4: nx: 2145
> DEBUG 4: ny: 1377
> DEBUG 4:
> DEBUG 3: Grid Definition: Projection: Lambert Conformal Nx: 2145 Ny:
1377
> Lat_LL: 20.192 Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082
Cone:
> 0.423 Bx: 1088.0030 By: 5483.6617
> DEBUG 2: Output grid: Projection: Lambert Conformal Nx: 2145 Ny:
1377
> Lat_LL: 20.192 Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082
Cone:
> 0.423 Bx: 1088.0030 By: 5483.6617
> DEBUG 2: Interpolation options: method = BUDGET, width = 2,
vld_thresh =
> 0.5
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> type "FileType_NcMet".
> ERROR :
> ERROR : process_data_file() -> trouble getting field "name="APCP";
> level="A6";" from file
"/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc"
> ERROR :
>
>
------------------------------------------------
Subject: regrid_data_plane question
From: Brad Diehl - NOAA Federal
Time: Fri Jan 05 10:32:50 2018
As always, prompt service and the right answers. It was successful
after
making the field name and level changes. You can close this ticket.
Many thanks again, Brad Diehl
On Thu, Jan 4, 2018 at 10:25 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hi Brad,
>
> I see that you're trying to run pcp_combine followed by
regrid_data_plane.
>
> I suspect you just need to modify your script slightly. pcp_combine
writes
> out a NetCDF file and you should take a look at the output variable
that it
> writes:
>
> ncdump -h /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
>
> I'm guessing the data variable is named "APCP_06".
>
> When you pass this file to regrid_data_plane, you need to tell it
the name
> of NetCDF variable to process:
>
> regrid_data_plane out.nc regrid_out.nc -field 'name="APCP_06";
> level="(*,*)";'
>
> Unfortunately, there's a difference in how you specify the name and
level
> for GRIB versus NetCDF files and I'm guessing that's the source of
the
> confusion here.
>
> FYI, instructions for posting to our anonymous ftp site can be found
here:
> https://dtcenter.org/met/users/support/met_help.php#ftp
>
> Thanks,
> John Halley Gotway
>
> On Thu, Jan 4, 2018 at 3:17 PM, Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Thu Jan 04 15:17:58 2018: Request 83497 was acted upon.
> > Transaction: Ticket created by brad.diehl at noaa.gov
> > Queue: met_help
> > Subject: regrid_data_plane question
> > Owner: Nobody
> > Requestors: brad.diehl at noaa.gov
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=83497 >
> >
> >
> > To whom it may concern,
> >
> > I was trying to run a simple script that takes the netcdf output
file of
> > pcp_combine and then uses it as input to regrid_data_plane, with
the
> script
> > pasted below. The results of the script is pasted below that.
> > It seems to get getting hung up on getting field "name="APCP";
> level="A6";
> >
> > I know I should have ftp'd the input files but could not find the
> > instructions from earlier threads, if you could point me to the
anonymous
> > ftp instructions again I'll be sure to include them if need be.
> >
> > Thanks for your help,
> >
> > Brad Diehl
> >
> >
> > #!/bin/sh
> > ###
> > ############
> > source /u/Brad.Diehl/met_env/tut_env
> >
> >
> > in_file1='/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/
> > hrrrcogb3p018010415qpfff12'
> > in_file2='/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/
> > hrrrcogb3p018010415qpfff06'
> >
> >
> > pcp_combine -subtract $in_file1 12 $in_file2 06
> > /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
> >
> > regrid_data_plane /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc \
> > /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/urma_co \
> > /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/re_grid_out.nc \
> > -method BUDGET \
> > -width 2 \
> > -v 4 \
> > -field 'name="APCP"; level="A6";'
> >
> > echo after regrid
> >
> >
> >
> >
> >
> >
> > Brad.Diehl at g14a2: /mdlstat/save/usr/Brad.Diehl/pcp_hrrr$
> > ./pcp_hrrr_test.sh
> > DEBUG 1: Reading input file:
> > /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/hrrrcogb3p018010415qpfff12
> > DEBUG 1: Reading input file:
> > /mdlstat/save/usr/Brad.Diehl/pcp_hrrr/hrrrcogb3p018010415qpfff06
> > DEBUG 1: Writing output file:
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/
> out.nc
> > DEBUG 1: Reading data file:
/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc
> > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > Met2dDataFile object of type "FileType_NcMet".
> > DEBUG 4: Lambert Conformal Grid Data:
> > DEBUG 4: scale_lat_1 = 38.5
> > DEBUG 4: scale_lat_2 = 38.5
> > DEBUG 4: lat_pin = 21.1381
> > DEBUG 4: lon_pin = 122.72
> > DEBUG 4: x_pin = 0
> > DEBUG 4: y_pin = 0
> > DEBUG 4: lon_orient = 97.5
> > DEBUG 4: d_km = 3
> > DEBUG 4: r_km = 6371.2
> > DEBUG 4: nx = 1799
> > DEBUG 4: ny = 1059
> > DEBUG 2: Input grid: Projection: Lambert Conformal Nx: 1799 Ny:
1059
> > Lat_LL: 21.138 Lon_LL: 122.720 Lon_orient: 97.500 Alpha: 4203.494
Cone:
> > 0.623 Bx: 899.1693 By: 3198.9984
> > DEBUG 3: Use the grid defined by file
> > "/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/urma_co".
> > DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
> > Met2dDataFile object of type "FileType_Gb2".
> > DEBUG 4:
> > DEBUG 4: Lambert Conformal Grid Data:
> > DEBUG 4: scale_lat_1: 25
> > DEBUG 4: scale_lat_2: 25
> > DEBUG 4: lat_pin: 20.192
> > DEBUG 4: lon_pin: 121.554
> > DEBUG 4: x_pin: 0
> > DEBUG 4: y_pin: 0
> > DEBUG 4: lon_orient: 95
> > DEBUG 4: d_km: 2.5397
> > DEBUG 4: r_km: 6371.2
> > DEBUG 4: nx: 2145
> > DEBUG 4: ny: 1377
> > DEBUG 4:
> > DEBUG 3: Grid Definition: Projection: Lambert Conformal Nx: 2145
Ny: 1377
> > Lat_LL: 20.192 Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082
Cone:
> > 0.423 Bx: 1088.0030 By: 5483.6617
> > DEBUG 2: Output grid: Projection: Lambert Conformal Nx: 2145 Ny:
1377
> > Lat_LL: 20.192 Lon_LL: 121.554 Lon_orient: 95.000 Alpha: 6509.082
Cone:
> > 0.423 Bx: 1088.0030 By: 5483.6617
> > DEBUG 2: Interpolation options: method = BUDGET, width = 2,
vld_thresh =
> > 0.5
> > DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object of
> > type "FileType_NcMet".
> > ERROR :
> > ERROR : process_data_file() -> trouble getting field
"name="APCP";
> > level="A6";" from file
"/mdlstat/save/usr/Brad.Diehl/pcp_hrrr/out.nc"
> > ERROR :
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list