[Met_help] [rt.rap.ucar.edu #96981] History for GFS/regrid_data_plane call
Julie Prestopnik via RT
met_help at ucar.edu
Wed Oct 7 08:52:31 MDT 2020
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
To MET help,
We are having issues with regrid_data_plane to regrid the 1-hr QPF for the
GFS using MET V9.1. We have set up a test call to regrid_data_plane in a
script with the I/O filenames hard coded into the script.
The test script is:
/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/save/qpfvs/processing/regrid_test/test.sh
Note that there is a call to ncdump prior to regrid_data_plane, which
gives:
APCP_01:name = "APCP_01" ;
APCP_01:long_name = "Total Precipitation" ;
APCP_01:level = "A1" ;
with the name/level indicated in the call to regrid_data_plane:
-field "name="APCP_01"; level="A1";"
It should be noted that the input netcdf file is the output from
pcp_combine to compute a one hour QPF for the GFS.
The error from the call to regrid_data_plane is:
ERROR :
ERROR : do_assign_id() -> identifier "APCP_01" not defined in this scope!
ERROR :
Note that we tried many variations of name/level with no success, hence
this request to MET help. We have been using regrid_data_plane
successfully for years now, but cannot get past this hurdle this time.
Lastly, as -name is not required for regrid_data_plane , I left it out of
the call, but we tried many variations with no success.
Thanks in advance for your help.......Brad Diehl
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: GFS/regrid_data_plane call
From: Julie Prestopnik
Time: Tue Oct 06 15:37:37 2020
Hi Brad.
I see that you are having trouble getting regrid_data_plane to work.
Can
you please tell me what type of file
"/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co"
is so that I can take a look at it? Unfortunately, the "file" command
is
not very helpful:
[Julie.Prestopnik at m71a1 diehl_20201006]$ file
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co
>
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co:
> data
>
I tried ncdump and h5dump, but neither were successful.
Thanks!
Julie
On Tue, Oct 6, 2020 at 3:13 PM Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> Tue Oct 06 15:12:31 2020: Request 96981 was acted upon.
> Transaction: Ticket created by brad.diehl at noaa.gov
> Queue: met_help
> Subject: GFS/regrid_data_plane call
> Owner: Nobody
> Requestors: brad.diehl at noaa.gov
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96981 >
>
>
> To MET help,
>
> We are having issues with regrid_data_plane to regrid the 1-hr QPF
for the
> GFS using MET V9.1. We have set up a test call to regrid_data_plane
in a
> script with the I/O filenames hard coded into the script.
>
> The test script is:
>
>
>
/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/save/qpfvs/processing/regrid_test/test.sh
>
> Note that there is a call to ncdump prior to regrid_data_plane,
which
> gives:
>
> APCP_01:name = "APCP_01" ;
> APCP_01:long_name = "Total Precipitation" ;
> APCP_01:level = "A1" ;
>
> with the name/level indicated in the call to regrid_data_plane:
>
> -field "name="APCP_01"; level="A1";"
>
> It should be noted that the input netcdf file is the output from
> pcp_combine to compute a one hour QPF for the GFS.
>
> The error from the call to regrid_data_plane is:
>
>
> ERROR :
> ERROR : do_assign_id() -> identifier "APCP_01" not defined in this
scope!
> ERROR :
>
> Note that we tried many variations of name/level with no success,
hence
> this request to MET help. We have been using regrid_data_plane
> successfully for years now, but cannot get past this hurdle this
time.
> Lastly, as -name is not required for regrid_data_plane , I left it
out of
> the call, but we tried many variations with no success.
>
> Thanks in advance for your help.......Brad Diehl
>
>
--
Julie Prestopnik (she/her/hers)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: GFS/regrid_data_plane call
From: Julie Prestopnik
Time: Tue Oct 06 15:56:21 2020
Ah! I see that it is a GRIB2 file (from using "strings blend_co"). I
will
follow up once I know more.
Julie
On Tue, Oct 6, 2020 at 3:37 PM Julie Prestopnik <jpresto at ucar.edu>
wrote:
> Hi Brad.
>
> I see that you are having trouble getting regrid_data_plane to work.
Can
> you please tell me what type of file
> "/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co"
> is so that I can take a look at it? Unfortunately, the "file"
command is
> not very helpful:
>
> [Julie.Prestopnik at m71a1 diehl_20201006]$ file
>> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co
>>
>> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co:
>> data
>>
>
> I tried ncdump and h5dump, but neither were successful.
>
> Thanks!
>
> Julie
>
> On Tue, Oct 6, 2020 at 3:13 PM Brad Diehl - NOAA Federal via RT <
> met_help at ucar.edu> wrote:
>
>>
>> Tue Oct 06 15:12:31 2020: Request 96981 was acted upon.
>> Transaction: Ticket created by brad.diehl at noaa.gov
>> Queue: met_help
>> Subject: GFS/regrid_data_plane call
>> Owner: Nobody
>> Requestors: brad.diehl at noaa.gov
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96981 >
>>
>>
>> To MET help,
>>
>> We are having issues with regrid_data_plane to regrid the 1-hr QPF
for
>> the
>> GFS using MET V9.1. We have set up a test call to
regrid_data_plane in a
>> script with the I/O filenames hard coded into the script.
>>
>> The test script is:
>>
>>
>>
/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/save/qpfvs/processing/regrid_test/test.sh
>>
>> Note that there is a call to ncdump prior to regrid_data_plane,
which
>> gives:
>>
>> APCP_01:name = "APCP_01" ;
>> APCP_01:long_name = "Total Precipitation" ;
>> APCP_01:level = "A1" ;
>>
>> with the name/level indicated in the call to regrid_data_plane:
>>
>> -field "name="APCP_01"; level="A1";"
>>
>> It should be noted that the input netcdf file is the output from
>> pcp_combine to compute a one hour QPF for the GFS.
>>
>> The error from the call to regrid_data_plane is:
>>
>>
>> ERROR :
>> ERROR : do_assign_id() -> identifier "APCP_01" not defined in this
scope!
>> ERROR :
>>
>> Note that we tried many variations of name/level with no success,
hence
>> this request to MET help. We have been using regrid_data_plane
>> successfully for years now, but cannot get past this hurdle this
time.
>> Lastly, as -name is not required for regrid_data_plane , I left it
out of
>> the call, but we tried many variations with no success.
>>
>> Thanks in advance for your help.......Brad Diehl
>>
>>
>
> --
> Julie Prestopnik (she/her/hers)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
--
Julie Prestopnik (she/her/hers)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: GFS/regrid_data_plane call
From: Julie Prestopnik
Time: Tue Oct 06 16:07:12 2020
Hi Brad.
Ok. I see the following in your test.sh file:
$METexec/regrid_data_plane
>
/gpfs/dell2/ptmp/Brad.Diehl/matching_urma_2020092312_co_qpf01.75056/match_co_qpf01_2020092312.75056/012/gfscon
> c2p520092300qpf01f012_combine
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co
> /gpfs/dell2/ptmp/Br
>
ad.Diehl/matching_urma_2020092312_co_qpf01.75056/match_co_qpf01_2020092312.75056/012/gfsconc2p520092300qpf01f012
> -method BUDGET -width 1
> -v 4 -field "name="APCP_01"; level="A1";"
>
I copied over gfsconc2p520092300qpf01f012_combine and blend_co to my
local
area and ran the following command successfully:
[Julie.Prestopnik at m71a1 diehl_20201006]$ regrid_data_plane
> gfsconc2p520092300qpf01f012_combine blend_co output.nc -field
> 'name="APCP_01"; level="A1";' -method BUDGET -width 1
> WARNING:
> WARNING: RegridInfo::validate() -> Resetting the regridding method
from
> "BUDGET" to "NEAREST" since the regridding width is 1.
> WARNING:
> DEBUG 1: Reading data file: gfsconc2p520092300qpf01f012_combine
> DEBUG 2: Input grid: Projection: Lat/Lon Nx: 1440 Ny: 721 lat_ll:
-90.000
> lon_ll: -0.000 delta_lat: 0.250 delta_lon: 0.250
> DEBUG 2: Output grid: Projection: Lambert Conformal Nx: 2345 Ny:
1597
> Lat_LL: 19.229 Lon_LL: 126.277 Lon_orient: 95.000 Alpha: 6509.082
Cone:
> 0.423 Bx: 1288.0039 By: 5483.6648
> DEBUG 2: Interpolation options: method = NEAREST, width = 1, shape =
> SQUARE, vld_thresh = 0.5
> DEBUG 2: Range of input data (name="APCP_01"; level="A1";) is 0 to
81.625.
> DEBUG 2: Range of regridded data (name="APCP_01"; level="A1";) is 0
to
> 14.4375.
> DEBUG 1: Writing output file: output.nc
Note that I changed the order of the arguments to regrid_data_plane
and
also changed the double quote that surrounded:
> name="APCP_01"; level="A1";
to single quotes. But, that is all I changed.
I was not able to reproduce the error you described:
> ERROR :
> ERROR : do_assign_id() -> identifier "APCP_01" not defined in this
scope!
> ERROR :
>
When I ran the command the way you had it listed in your test.sh file,
I
simply received the usage statement for regrid_data_plane, which is
why I
made the changes I described above:
> regrid_data_plane gfsconc2p520092300qpf01f012_combine blend_co
output.nc
> -method BUDGET -width 1-v 4 -field "name="APCP"; level="A1";"
Give that a try and let me know how it goes please. I hope it helps!
Thanks,
Julie
>
On Tue, Oct 6, 2020 at 3:56 PM Julie Prestopnik <jpresto at ucar.edu>
wrote:
> Ah! I see that it is a GRIB2 file (from using "strings blend_co").
I will
> follow up once I know more.
>
> Julie
>
> On Tue, Oct 6, 2020 at 3:37 PM Julie Prestopnik <jpresto at ucar.edu>
wrote:
>
>> Hi Brad.
>>
>> I see that you are having trouble getting regrid_data_plane to
work. Can
>> you please tell me what type of file
>> "/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co"
>> is so that I can take a look at it? Unfortunately, the "file"
command is
>> not very helpful:
>>
>> [Julie.Prestopnik at m71a1 diehl_20201006]$ file
>>> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co
>>>
>>> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co:
>>> data
>>>
>>
>> I tried ncdump and h5dump, but neither were successful.
>>
>> Thanks!
>>
>> Julie
>>
>> On Tue, Oct 6, 2020 at 3:13 PM Brad Diehl - NOAA Federal via RT <
>> met_help at ucar.edu> wrote:
>>
>>>
>>> Tue Oct 06 15:12:31 2020: Request 96981 was acted upon.
>>> Transaction: Ticket created by brad.diehl at noaa.gov
>>> Queue: met_help
>>> Subject: GFS/regrid_data_plane call
>>> Owner: Nobody
>>> Requestors: brad.diehl at noaa.gov
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96981 >
>>>
>>>
>>> To MET help,
>>>
>>> We are having issues with regrid_data_plane to regrid the 1-hr
QPF for
>>> the
>>> GFS using MET V9.1. We have set up a test call to
regrid_data_plane in a
>>> script with the I/O filenames hard coded into the script.
>>>
>>> The test script is:
>>>
>>>
>>>
/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/save/qpfvs/processing/regrid_test/test.sh
>>>
>>> Note that there is a call to ncdump prior to regrid_data_plane,
which
>>> gives:
>>>
>>> APCP_01:name = "APCP_01" ;
>>> APCP_01:long_name = "Total Precipitation" ;
>>> APCP_01:level = "A1" ;
>>>
>>> with the name/level indicated in the call to regrid_data_plane:
>>>
>>> -field "name="APCP_01"; level="A1";"
>>>
>>> It should be noted that the input netcdf file is the output from
>>> pcp_combine to compute a one hour QPF for the GFS.
>>>
>>> The error from the call to regrid_data_plane is:
>>>
>>>
>>> ERROR :
>>> ERROR : do_assign_id() -> identifier "APCP_01" not defined in
this
>>> scope!
>>> ERROR :
>>>
>>> Note that we tried many variations of name/level with no success,
hence
>>> this request to MET help. We have been using regrid_data_plane
>>> successfully for years now, but cannot get past this hurdle this
time.
>>> Lastly, as -name is not required for regrid_data_plane , I left it
out of
>>> the call, but we tried many variations with no success.
>>>
>>> Thanks in advance for your help.......Brad Diehl
>>>
>>>
>>
>> --
>> Julie Prestopnik (she/her/hers)
>> Software Engineer
>> National Center for Atmospheric Research
>> Research Applications Laboratory
>> Email: jpresto at ucar.edu
>>
>> My working day may not be your working day. Please do not feel
obliged
>> to reply to this email outside of your normal working hours.
>>
>
>
> --
> Julie Prestopnik (she/her/hers)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
--
Julie Prestopnik (she/her/hers)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: GFS/regrid_data_plane call
From: Brad Diehl - NOAA Federal
Time: Wed Oct 07 08:41:20 2020
Julie,
I made the simple change you recommended for the single quotes as
follows:
'name="APCP_01"; level="A1";'
and that did the trick, and when I looked back at the original call I
had
in my earlier script, the single quotes were there, don't know how I
did'nt
notice. I guess you do this work long enough you see what you want to
see
and can't see what's right in front of you, same is true in life too I
guess. At any rate, many thanks for your help on this matter and
consider
it closed........Brad
On Tue, Oct 6, 2020 at 6:07 PM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Brad.
>
> Ok. I see the following in your test.sh file:
>
> $METexec/regrid_data_plane
> >
>
/gpfs/dell2/ptmp/Brad.Diehl/matching_urma_2020092312_co_qpf01.75056/match_co_qpf01_2020092312.75056/012/gfscon
> > c2p520092300qpf01f012_combine
> >
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co
> > /gpfs/dell2/ptmp/Br
> >
>
ad.Diehl/matching_urma_2020092312_co_qpf01.75056/match_co_qpf01_2020092312.75056/012/gfsconc2p520092300qpf01f012
> > -method BUDGET -width 1
> > -v 4 -field "name="APCP_01"; level="A1";"
> >
>
> I copied over gfsconc2p520092300qpf01f012_combine and blend_co to my
local
> area and ran the following command successfully:
>
> [Julie.Prestopnik at m71a1 diehl_20201006]$ regrid_data_plane
> > gfsconc2p520092300qpf01f012_combine blend_co output.nc -field
> > 'name="APCP_01"; level="A1";' -method BUDGET -width 1
> > WARNING:
> > WARNING: RegridInfo::validate() -> Resetting the regridding method
from
> > "BUDGET" to "NEAREST" since the regridding width is 1.
> > WARNING:
> > DEBUG 1: Reading data file: gfsconc2p520092300qpf01f012_combine
> > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 1440 Ny: 721 lat_ll:
-90.000
> > lon_ll: -0.000 delta_lat: 0.250 delta_lon: 0.250
> > DEBUG 2: Output grid: Projection: Lambert Conformal Nx: 2345 Ny:
1597
> > Lat_LL: 19.229 Lon_LL: 126.277 Lon_orient: 95.000 Alpha: 6509.082
Cone:
> > 0.423 Bx: 1288.0039 By: 5483.6648
> > DEBUG 2: Interpolation options: method = NEAREST, width = 1, shape
=
> > SQUARE, vld_thresh = 0.5
> > DEBUG 2: Range of input data (name="APCP_01"; level="A1";) is 0 to
> 81.625.
> > DEBUG 2: Range of regridded data (name="APCP_01"; level="A1";) is
0 to
> > 14.4375.
> > DEBUG 1: Writing output file: output.nc
>
>
> Note that I changed the order of the arguments to regrid_data_plane
and
> also changed the double quote that surrounded:
>
> > name="APCP_01"; level="A1";
>
> to single quotes. But, that is all I changed.
>
> I was not able to reproduce the error you described:
>
> > ERROR :
> > ERROR : do_assign_id() -> identifier "APCP_01" not defined in
this
> scope!
> > ERROR :
> >
>
> When I ran the command the way you had it listed in your test.sh
file, I
> simply received the usage statement for regrid_data_plane, which is
why I
> made the changes I described above:
>
> > regrid_data_plane gfsconc2p520092300qpf01f012_combine blend_co
output.nc
> > -method BUDGET -width 1-v 4 -field "name="APCP"; level="A1";"
>
>
> Give that a try and let me know how it goes please. I hope it
helps!
>
> Thanks,
> Julie
>
> >
>
>
> On Tue, Oct 6, 2020 at 3:56 PM Julie Prestopnik <jpresto at ucar.edu>
wrote:
>
> > Ah! I see that it is a GRIB2 file (from using "strings
blend_co"). I
> will
> > follow up once I know more.
> >
> > Julie
> >
> > On Tue, Oct 6, 2020 at 3:37 PM Julie Prestopnik <jpresto at ucar.edu>
> wrote:
> >
> >> Hi Brad.
> >>
> >> I see that you are having trouble getting regrid_data_plane to
work.
> Can
> >> you please tell me what type of file
> >>
> "/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co"
> >> is so that I can take a look at it? Unfortunately, the "file"
command
> is
> >> not very helpful:
> >>
> >> [Julie.Prestopnik at m71a1 diehl_20201006]$ file
> >>>
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co
> >>>
> >>>
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co:
> >>> data
> >>>
> >>
> >> I tried ncdump and h5dump, but neither were successful.
> >>
> >> Thanks!
> >>
> >> Julie
> >>
> >> On Tue, Oct 6, 2020 at 3:13 PM Brad Diehl - NOAA Federal via RT <
> >> met_help at ucar.edu> wrote:
> >>
> >>>
> >>> Tue Oct 06 15:12:31 2020: Request 96981 was acted upon.
> >>> Transaction: Ticket created by brad.diehl at noaa.gov
> >>> Queue: met_help
> >>> Subject: GFS/regrid_data_plane call
> >>> Owner: Nobody
> >>> Requestors: brad.diehl at noaa.gov
> >>> Status: new
> >>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96981
> >
> >>>
> >>>
> >>> To MET help,
> >>>
> >>> We are having issues with regrid_data_plane to regrid the 1-hr
QPF for
> >>> the
> >>> GFS using MET V9.1. We have set up a test call to
regrid_data_plane
> in a
> >>> script with the I/O filenames hard coded into the script.
> >>>
> >>> The test script is:
> >>>
> >>>
> >>>
>
/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/save/qpfvs/processing/regrid_test/test.sh
> >>>
> >>> Note that there is a call to ncdump prior to regrid_data_plane,
which
> >>> gives:
> >>>
> >>> APCP_01:name = "APCP_01" ;
> >>> APCP_01:long_name = "Total Precipitation" ;
> >>> APCP_01:level = "A1" ;
> >>>
> >>> with the name/level indicated in the call to regrid_data_plane:
> >>>
> >>> -field "name="APCP_01"; level="A1";"
> >>>
> >>> It should be noted that the input netcdf file is the output from
> >>> pcp_combine to compute a one hour QPF for the GFS.
> >>>
> >>> The error from the call to regrid_data_plane is:
> >>>
> >>>
> >>> ERROR :
> >>> ERROR : do_assign_id() -> identifier "APCP_01" not defined in
this
> >>> scope!
> >>> ERROR :
> >>>
> >>> Note that we tried many variations of name/level with no
success, hence
> >>> this request to MET help. We have been using regrid_data_plane
> >>> successfully for years now, but cannot get past this hurdle this
time.
> >>> Lastly, as -name is not required for regrid_data_plane , I left
it out
> of
> >>> the call, but we tried many variations with no success.
> >>>
> >>> Thanks in advance for your help.......Brad Diehl
> >>>
> >>>
> >>
> >> --
> >> Julie Prestopnik (she/her/hers)
> >> Software Engineer
> >> National Center for Atmospheric Research
> >> Research Applications Laboratory
> >> Email: jpresto at ucar.edu
> >>
> >> My working day may not be your working day. Please do not feel
obliged
> >> to reply to this email outside of your normal working hours.
> >>
> >
> >
> > --
> > Julie Prestopnik (she/her/hers)
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged
> to
> > reply to this email outside of your normal working hours.
> >
>
>
> --
> Julie Prestopnik (she/her/hers)
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
------------------------------------------------
Subject: GFS/regrid_data_plane call
From: Julie Prestopnik
Time: Wed Oct 07 08:51:28 2020
Great! Sometimes those are the hardest problems to see. I've done
that
myself several times. :) I'm glad the problem was easy to solve.
I'll go
ahead and close this ticket.
Julie
On Wed, Oct 7, 2020 at 8:41 AM Brad Diehl - NOAA Federal via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96981 >
>
> Julie,
>
> I made the simple change you recommended for the single quotes as
follows:
> 'name="APCP_01"; level="A1";'
>
> and that did the trick, and when I looked back at the original call
I had
> in my earlier script, the single quotes were there, don't know how I
did'nt
> notice. I guess you do this work long enough you see what you want
to see
> and can't see what's right in front of you, same is true in life too
I
> guess. At any rate, many thanks for your help on this matter and
consider
> it closed........Brad
>
>
>
>
>
> On Tue, Oct 6, 2020 at 6:07 PM Julie Prestopnik via RT
<met_help at ucar.edu>
> wrote:
>
> > Hi Brad.
> >
> > Ok. I see the following in your test.sh file:
> >
> > $METexec/regrid_data_plane
> > >
> >
>
/gpfs/dell2/ptmp/Brad.Diehl/matching_urma_2020092312_co_qpf01.75056/match_co_qpf01_2020092312.75056/012/gfscon
> > > c2p520092300qpf01f012_combine
> > >
> >
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co
> > > /gpfs/dell2/ptmp/Br
> > >
> >
>
ad.Diehl/matching_urma_2020092312_co_qpf01.75056/match_co_qpf01_2020092312.75056/012/gfsconc2p520092300qpf01f012
> > > -method BUDGET -width 1
> > > -v 4 -field "name="APCP_01"; level="A1";"
> > >
> >
> > I copied over gfsconc2p520092300qpf01f012_combine and blend_co to
my
> local
> > area and ran the following command successfully:
> >
> > [Julie.Prestopnik at m71a1 diehl_20201006]$ regrid_data_plane
> > > gfsconc2p520092300qpf01f012_combine blend_co output.nc -field
> > > 'name="APCP_01"; level="A1";' -method BUDGET -width 1
> > > WARNING:
> > > WARNING: RegridInfo::validate() -> Resetting the regridding
method from
> > > "BUDGET" to "NEAREST" since the regridding width is 1.
> > > WARNING:
> > > DEBUG 1: Reading data file: gfsconc2p520092300qpf01f012_combine
> > > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 1440 Ny: 721
lat_ll:
> -90.000
> > > lon_ll: -0.000 delta_lat: 0.250 delta_lon: 0.250
> > > DEBUG 2: Output grid: Projection: Lambert Conformal Nx: 2345 Ny:
1597
> > > Lat_LL: 19.229 Lon_LL: 126.277 Lon_orient: 95.000 Alpha:
6509.082 Cone:
> > > 0.423 Bx: 1288.0039 By: 5483.6648
> > > DEBUG 2: Interpolation options: method = NEAREST, width = 1,
shape =
> > > SQUARE, vld_thresh = 0.5
> > > DEBUG 2: Range of input data (name="APCP_01"; level="A1";) is 0
to
> > 81.625.
> > > DEBUG 2: Range of regridded data (name="APCP_01"; level="A1";)
is 0 to
> > > 14.4375.
> > > DEBUG 1: Writing output file: output.nc
> >
> >
> > Note that I changed the order of the arguments to
regrid_data_plane and
> > also changed the double quote that surrounded:
> >
> > > name="APCP_01"; level="A1";
> >
> > to single quotes. But, that is all I changed.
> >
> > I was not able to reproduce the error you described:
> >
> > > ERROR :
> > > ERROR : do_assign_id() -> identifier "APCP_01" not defined in
this
> > scope!
> > > ERROR :
> > >
> >
> > When I ran the command the way you had it listed in your test.sh
file, I
> > simply received the usage statement for regrid_data_plane, which
is why I
> > made the changes I described above:
> >
> > > regrid_data_plane gfsconc2p520092300qpf01f012_combine blend_co
> output.nc
> > > -method BUDGET -width 1-v 4 -field "name="APCP"; level="A1";"
> >
> >
> > Give that a try and let me know how it goes please. I hope it
helps!
> >
> > Thanks,
> > Julie
> >
> > >
> >
> >
> > On Tue, Oct 6, 2020 at 3:56 PM Julie Prestopnik <jpresto at ucar.edu>
> wrote:
> >
> > > Ah! I see that it is a GRIB2 file (from using "strings
blend_co"). I
> > will
> > > follow up once I know more.
> > >
> > > Julie
> > >
> > > On Tue, Oct 6, 2020 at 3:37 PM Julie Prestopnik
<jpresto at ucar.edu>
> > wrote:
> > >
> > >> Hi Brad.
> > >>
> > >> I see that you are having trouble getting regrid_data_plane to
work.
> > Can
> > >> you please tell me what type of file
> > >>
> >
> "/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co"
> > >> is so that I can take a look at it? Unfortunately, the "file"
command
> > is
> > >> not very helpful:
> > >>
> > >> [Julie.Prestopnik at m71a1 diehl_20201006]$ file
> > >>>
> >
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co
> > >>>
> > >>>
> >
> /gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/mdl-
verification/fix/urma_mask/blend_co:
> > >>> data
> > >>>
> > >>
> > >> I tried ncdump and h5dump, but neither were successful.
> > >>
> > >> Thanks!
> > >>
> > >> Julie
> > >>
> > >> On Tue, Oct 6, 2020 at 3:13 PM Brad Diehl - NOAA Federal via RT
<
> > >> met_help at ucar.edu> wrote:
> > >>
> > >>>
> > >>> Tue Oct 06 15:12:31 2020: Request 96981 was acted upon.
> > >>> Transaction: Ticket created by brad.diehl at noaa.gov
> > >>> Queue: met_help
> > >>> Subject: GFS/regrid_data_plane call
> > >>> Owner: Nobody
> > >>> Requestors: brad.diehl at noaa.gov
> > >>> Status: new
> > >>> Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=96981
> > >
> > >>>
> > >>>
> > >>> To MET help,
> > >>>
> > >>> We are having issues with regrid_data_plane to regrid the 1-
hr QPF
> for
> > >>> the
> > >>> GFS using MET V9.1. We have set up a test call to
regrid_data_plane
> > in a
> > >>> script with the I/O filenames hard coded into the script.
> > >>>
> > >>> The test script is:
> > >>>
> > >>>
> > >>>
> >
>
/gpfs/dell3/mdl/mdlverif/noscrub/usr/Brad.Diehl/save/qpfvs/processing/regrid_test/test.sh
> > >>>
> > >>> Note that there is a call to ncdump prior to
regrid_data_plane, which
> > >>> gives:
> > >>>
> > >>> APCP_01:name = "APCP_01" ;
> > >>> APCP_01:long_name = "Total Precipitation" ;
> > >>> APCP_01:level = "A1" ;
> > >>>
> > >>> with the name/level indicated in the call to
regrid_data_plane:
> > >>>
> > >>> -field "name="APCP_01"; level="A1";"
> > >>>
> > >>> It should be noted that the input netcdf file is the output
from
> > >>> pcp_combine to compute a one hour QPF for the GFS.
> > >>>
> > >>> The error from the call to regrid_data_plane is:
> > >>>
> > >>>
> > >>> ERROR :
> > >>> ERROR : do_assign_id() -> identifier "APCP_01" not defined in
this
> > >>> scope!
> > >>> ERROR :
> > >>>
> > >>> Note that we tried many variations of name/level with no
success,
> hence
> > >>> this request to MET help. We have been using
regrid_data_plane
> > >>> successfully for years now, but cannot get past this hurdle
this
> time.
> > >>> Lastly, as -name is not required for regrid_data_plane , I
left it
> out
> > of
> > >>> the call, but we tried many variations with no success.
> > >>>
> > >>> Thanks in advance for your help.......Brad Diehl
> > >>>
> > >>>
> > >>
> > >> --
> > >> Julie Prestopnik (she/her/hers)
> > >> Software Engineer
> > >> National Center for Atmospheric Research
> > >> Research Applications Laboratory
> > >> Email: jpresto at ucar.edu
> > >>
> > >> My working day may not be your working day. Please do not feel
> obliged
> > >> to reply to this email outside of your normal working hours.
> > >>
> > >
> > >
> > > --
> > > Julie Prestopnik (she/her/hers)
> > > Software Engineer
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Email: jpresto at ucar.edu
> > >
> > > My working day may not be your working day. Please do not feel
obliged
> > to
> > > reply to this email outside of your normal working hours.
> > >
> >
> >
> > --
> > Julie Prestopnik (she/her/hers)
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged
> to
> > reply to this email outside of your normal working hours.
> >
> >
>
>
--
Julie Prestopnik (she/her/hers)
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
More information about the Met_help
mailing list