[Met_help] [rt.rap.ucar.edu #87383] History for NcCfFile getData error
Julie Prestopnik via RT
met_help at ucar.edu
Mon Oct 22 11:58:53 MDT 2018
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear Sir/Madam,
I am Lintao, trying to use regrid_data_plane tool to regrid wrf output data
based on StageIV data set. However, I got the error message as:
ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane &) ->
needed 3 arguments for variable accu_precip, got 2
I run MET on Cheyenne, using the follow code as an example:
regrid_data_plane \
/glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-31_23:0 \
/glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h \
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-08-31_
23:00 \
-field 'name="accu_precip"; level="surf";' \
-method UW_MEAN -width 4
May I ask what is the problem please?
Please feel free to go to my working directory and perhaps check the data
files.
Thank you in advance.
Best regards,
Lintao
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: NcCfFile getData error
From: Julie Prestopnik
Time: Mon Oct 15 14:29:43 2018
Thank you for providing the command that you ran. It is useful in
reproducing the problem. I changed your value of level from "surf" to
"(0,*,*)". For more information on this change and specifying the
level
for a NetCDF file please see page 59 in our met-8.0 User's Guide:
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v8.0.pdf
Running the following command:
> regrid_data_plane hrly_precip_d01_2008-08-31_23:0 ST4.2008083123.01h
> regridded_precip_2008-08-31_23:00 -field 'name="accu_precip";
> level="(0,*,*)";' -method UW_MEAN -width 4
>
I did receive an output file, however the timing information will be
incorrect. Please note the warning in the output:
> DEBUG 1: Reading data file: hrly_precip_d01_2008-08-31_23:0
> WARNING:
> WARNING: NcCfFile::open() -> could not determine the valid time,
using 0.
> WARNING:
> DEBUG 2: Input grid: Projection: Lat/Lon Nx: 615 Ny: 559 lat_ll:
-112.815
> lon_ll: -27.730 delta_lat: 0.005 delta_lon: 0.005
> DEBUG 2: Output grid: Projection: Stereographic Nx: 1121 Ny: 881
> IsNorthHemisphere: true Lon_orient: 105.000 Bx: 399.440 By: 1599.329
Alpha:
> 2496.0783
> DEBUG 2: Interpolation options: method = UW_MEAN, width = 4, shape =
> SQUARE, vld_thresh = 0.5
> DEBUG 2: NcCfFile::getData -> data was flipped to north.
> DEBUG 2: Range of input data (name="accu_precip"; level="(0,*,*)";)
is 0
> to 38.4619.
> DEBUG 2: Range of regridded data (name="accu_precip";
level="(0,*,*)";) is
> -9999 to -9999.
> DEBUG 1: Writing output file: regridded_precip_2008-08-31_23:00
>
Using ncdump to see the file attributes, I see this timing
information:
> XTIME:standard_name = "time" ;
> XTIME:units = "minutes since 2008-4-1 00:00:00" ;
>
However, that is not the CF-compliant way of specifying the timing
info.
So, while MET will be able to read the data, the timing info won't be
accurate in the MET output files. Instead, you'll see "19700101" which
is
unixtime = 0.
I would suggest making your file more CF-compliant. Information for
the CF
coordinate time convention can be found here:
http://cfconventions.org/cf-conventions/v1.6.0/cf-
conventions.html#time-coordinate
I hope that helps!
Julie
On Mon, Oct 15, 2018 at 9:09 AM Lintao Li via RT <met_help at ucar.edu>
wrote:
>
> Mon Oct 15 09:09:34 2018: Request 87383 was acted upon.
> Transaction: Ticket created by lintaoli168 at gmail.com
> Queue: met_help
> Subject: NcCfFile getData error
> Owner: Nobody
> Requestors: lintaoli168 at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
>
>
> Dear Sir/Madam,
> I am Lintao, trying to use regrid_data_plane tool to regrid wrf
output data
> based on StageIV data set. However, I got the error message as:
> ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane &)
->
> needed 3 arguments for variable accu_precip, got 2
>
> I run MET on Cheyenne, using the follow code as an example:
>
> regrid_data_plane \
>
> /glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-
31_23:0 \
>
> /glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h \
>
>
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-
08-31_
> 23:00 \
>
> -field 'name="accu_precip"; level="surf";' \
>
> -method UW_MEAN -width 4
>
> May I ask what is the problem please?
> Please feel free to go to my working directory and perhaps check the
data
> files.
>
> Thank you in advance.
>
> Best regards,
> Lintao
>
>
------------------------------------------------
Subject: NcCfFile getData error
From: Lintao Li
Time: Mon Oct 15 15:49:19 2018
Hi Julie,
Many thanks for your very detailed explanation. It is very helpful. I
appreciate it very much.
Best regards,
Lintao
On Mon, Oct 15, 2018 at 2:29 PM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Thank you for providing the command that you ran. It is useful in
> reproducing the problem. I changed your value of level from "surf"
to
> "(0,*,*)". For more information on this change and specifying the
level
> for a NetCDF file please see page 59 in our met-8.0 User's Guide:
>
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v8.0.pdf
>
> Running the following command:
>
> > regrid_data_plane hrly_precip_d01_2008-08-31_23:0
ST4.2008083123.01h
> > regridded_precip_2008-08-31_23:00 -field 'name="accu_precip";
> > level="(0,*,*)";' -method UW_MEAN -width 4
> >
>
> I did receive an output file, however the timing information will be
> incorrect. Please note the warning in the output:
>
> > DEBUG 1: Reading data file: hrly_precip_d01_2008-08-31_23:0
> > WARNING:
> > WARNING: NcCfFile::open() -> could not determine the valid time,
using 0.
> > WARNING:
> > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 615 Ny: 559 lat_ll:
-112.815
> > lon_ll: -27.730 delta_lat: 0.005 delta_lon: 0.005
> > DEBUG 2: Output grid: Projection: Stereographic Nx: 1121 Ny: 881
> > IsNorthHemisphere: true Lon_orient: 105.000 Bx: 399.440 By:
1599.329
> Alpha:
> > 2496.0783
> > DEBUG 2: Interpolation options: method = UW_MEAN, width = 4, shape
=
> > SQUARE, vld_thresh = 0.5
> > DEBUG 2: NcCfFile::getData -> data was flipped to north.
> > DEBUG 2: Range of input data (name="accu_precip";
level="(0,*,*)";) is 0
> > to 38.4619.
> > DEBUG 2: Range of regridded data (name="accu_precip";
level="(0,*,*)";)
> is
> > -9999 to -9999.
> > DEBUG 1: Writing output file: regridded_precip_2008-08-31_23:00
> >
>
> Using ncdump to see the file attributes, I see this timing
information:
>
> > XTIME:standard_name = "time" ;
> > XTIME:units = "minutes since 2008-4-1 00:00:00" ;
> >
>
>
> However, that is not the CF-compliant way of specifying the timing
info.
> So, while MET will be able to read the data, the timing info won't
be
> accurate in the MET output files. Instead, you'll see "19700101"
which is
> unixtime = 0.
>
> I would suggest making your file more CF-compliant. Information for
the CF
> coordinate time convention can be found here:
>
> http://cfconventions.org/cf-conventions/v1.6.0/cf-
conventions.html#time-coordinate
>
> I hope that helps!
>
> Julie
>
>
> On Mon, Oct 15, 2018 at 9:09 AM Lintao Li via RT <met_help at ucar.edu>
> wrote:
>
> >
> > Mon Oct 15 09:09:34 2018: Request 87383 was acted upon.
> > Transaction: Ticket created by lintaoli168 at gmail.com
> > Queue: met_help
> > Subject: NcCfFile getData error
> > Owner: Nobody
> > Requestors: lintaoli168 at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
> >
> >
> > Dear Sir/Madam,
> > I am Lintao, trying to use regrid_data_plane tool to regrid wrf
output
> data
> > based on StageIV data set. However, I got the error message as:
> > ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane
&) ->
> > needed 3 arguments for variable accu_precip, got 2
> >
> > I run MET on Cheyenne, using the follow code as an example:
> >
> > regrid_data_plane \
> >
> > /glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-
31_23:0 \
> >
> > /glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h
\
> >
> >
>
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-
08-31_
> > 23:00 \
> >
> > -field 'name="accu_precip"; level="surf";' \
> >
> > -method UW_MEAN -width 4
> >
> > May I ask what is the problem please?
> > Please feel free to go to my working directory and perhaps check
the data
> > files.
> >
> > Thank you in advance.
> >
> > Best regards,
> > Lintao
> >
> >
>
>
------------------------------------------------
Subject: NcCfFile getData error
From: Julie Prestopnik
Time: Tue Oct 16 13:42:20 2018
Hi Lintao. One of our developers discovered that our MET code is
swapping
the lat/lon values for accu_precip, so there will be missing values
after
regridding. We are working on fixing this bug with MET. We will
release a
patch for met-8.0, hopefully, later this week, which I will then
install on
cheyenne. Would you like me to follow up with you once the patch has
been
installed on cheyenne?
Thank you!
Julie
On Mon, Oct 15, 2018 at 3:49 PM Lintao Li via RT <met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
>
> Hi Julie,
> Many thanks for your very detailed explanation. It is very helpful.
I
> appreciate it very much.
>
> Best regards,
> Lintao
>
>
> On Mon, Oct 15, 2018 at 2:29 PM Julie Prestopnik via RT
<met_help at ucar.edu
> >
> wrote:
>
> > Thank you for providing the command that you ran. It is useful in
> > reproducing the problem. I changed your value of level from
"surf" to
> > "(0,*,*)". For more information on this change and specifying the
level
> > for a NetCDF file please see page 59 in our met-8.0 User's Guide:
> >
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v8.0.pdf
> >
> > Running the following command:
> >
> > > regrid_data_plane hrly_precip_d01_2008-08-31_23:0
ST4.2008083123.01h
> > > regridded_precip_2008-08-31_23:00 -field 'name="accu_precip";
> > > level="(0,*,*)";' -method UW_MEAN -width 4
> > >
> >
> > I did receive an output file, however the timing information will
be
> > incorrect. Please note the warning in the output:
> >
> > > DEBUG 1: Reading data file: hrly_precip_d01_2008-08-31_23:0
> > > WARNING:
> > > WARNING: NcCfFile::open() -> could not determine the valid time,
using
> 0.
> > > WARNING:
> > > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 615 Ny: 559 lat_ll:
> -112.815
> > > lon_ll: -27.730 delta_lat: 0.005 delta_lon: 0.005
> > > DEBUG 2: Output grid: Projection: Stereographic Nx: 1121 Ny: 881
> > > IsNorthHemisphere: true Lon_orient: 105.000 Bx: 399.440 By:
1599.329
> > Alpha:
> > > 2496.0783
> > > DEBUG 2: Interpolation options: method = UW_MEAN, width = 4,
shape =
> > > SQUARE, vld_thresh = 0.5
> > > DEBUG 2: NcCfFile::getData -> data was flipped to north.
> > > DEBUG 2: Range of input data (name="accu_precip";
level="(0,*,*)";) is
> 0
> > > to 38.4619.
> > > DEBUG 2: Range of regridded data (name="accu_precip";
level="(0,*,*)";)
> > is
> > > -9999 to -9999.
> > > DEBUG 1: Writing output file: regridded_precip_2008-08-31_23:00
> > >
> >
> > Using ncdump to see the file attributes, I see this timing
information:
> >
> > > XTIME:standard_name = "time" ;
> > > XTIME:units = "minutes since 2008-4-1 00:00:00"
;
> > >
> >
> >
> > However, that is not the CF-compliant way of specifying the timing
info.
> > So, while MET will be able to read the data, the timing info won't
be
> > accurate in the MET output files. Instead, you'll see "19700101"
which is
> > unixtime = 0.
> >
> > I would suggest making your file more CF-compliant. Information
for the
> CF
> > coordinate time convention can be found here:
> >
> >
> http://cfconventions.org/cf-conventions/v1.6.0/cf-
conventions.html#time-coordinate
> >
> > I hope that helps!
> >
> > Julie
> >
> >
> > On Mon, Oct 15, 2018 at 9:09 AM Lintao Li via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > Mon Oct 15 09:09:34 2018: Request 87383 was acted upon.
> > > Transaction: Ticket created by lintaoli168 at gmail.com
> > > Queue: met_help
> > > Subject: NcCfFile getData error
> > > Owner: Nobody
> > > Requestors: lintaoli168 at gmail.com
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383
> >
> > >
> > >
> > > Dear Sir/Madam,
> > > I am Lintao, trying to use regrid_data_plane tool to regrid wrf
output
> > data
> > > based on StageIV data set. However, I got the error message as:
> > > ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane
&) ->
> > > needed 3 arguments for variable accu_precip, got 2
> > >
> > > I run MET on Cheyenne, using the follow code as an example:
> > >
> > > regrid_data_plane \
> > >
> > > /glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-
31_23:0
> \
> > >
> > >
/glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h \
> > >
> > >
> >
>
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-
08-31_
> > > 23:00 \
> > >
> > > -field 'name="accu_precip"; level="surf";' \
> > >
> > > -method UW_MEAN -width 4
> > >
> > > May I ask what is the problem please?
> > > Please feel free to go to my working directory and perhaps check
the
> data
> > > files.
> > >
> > > Thank you in advance.
> > >
> > > Best regards,
> > > Lintao
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #87383] NcCfFile getData error
From: Lintao Li
Time: Tue Oct 16 15:00:00 2018
Hi Julie,
Thanks a lot for your email. Yes please. I will be appreciated very
much.
Best regards,
Lintao
> On 16 Oct 2018, at 13:42, Julie Prestopnik via RT
<met_help at ucar.edu> wrote:
>
> Hi Lintao. One of our developers discovered that our MET code is
swapping
> the lat/lon values for accu_precip, so there will be missing values
after
> regridding. We are working on fixing this bug with MET. We will
release a
> patch for met-8.0, hopefully, later this week, which I will then
install on
> cheyenne. Would you like me to follow up with you once the patch
has been
> installed on cheyenne?
>
> Thank you!
>
> Julie
>
>> On Mon, Oct 15, 2018 at 3:49 PM Lintao Li via RT
<met_help at ucar.edu> wrote:
>>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
>>
>> Hi Julie,
>> Many thanks for your very detailed explanation. It is very helpful.
I
>> appreciate it very much.
>>
>> Best regards,
>> Lintao
>>
>>
>> On Mon, Oct 15, 2018 at 2:29 PM Julie Prestopnik via RT
<met_help at ucar.edu
>>>
>> wrote:
>>
>>> Thank you for providing the command that you ran. It is useful in
>>> reproducing the problem. I changed your value of level from
"surf" to
>>> "(0,*,*)". For more information on this change and specifying the
level
>>> for a NetCDF file please see page 59 in our met-8.0 User's Guide:
>>>
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v8.0.pdf
>>>
>>> Running the following command:
>>>
>>>> regrid_data_plane hrly_precip_d01_2008-08-31_23:0
ST4.2008083123.01h
>>>> regridded_precip_2008-08-31_23:00 -field 'name="accu_precip";
>>>> level="(0,*,*)";' -method UW_MEAN -width 4
>>>>
>>>
>>> I did receive an output file, however the timing information will
be
>>> incorrect. Please note the warning in the output:
>>>
>>>> DEBUG 1: Reading data file: hrly_precip_d01_2008-08-31_23:0
>>>> WARNING:
>>>> WARNING: NcCfFile::open() -> could not determine the valid time,
using
>> 0.
>>>> WARNING:
>>>> DEBUG 2: Input grid: Projection: Lat/Lon Nx: 615 Ny: 559 lat_ll:
>> -112.815
>>>> lon_ll: -27.730 delta_lat: 0.005 delta_lon: 0.005
>>>> DEBUG 2: Output grid: Projection: Stereographic Nx: 1121 Ny: 881
>>>> IsNorthHemisphere: true Lon_orient: 105.000 Bx: 399.440 By:
1599.329
>>> Alpha:
>>>> 2496.0783
>>>> DEBUG 2: Interpolation options: method = UW_MEAN, width = 4,
shape =
>>>> SQUARE, vld_thresh = 0.5
>>>> DEBUG 2: NcCfFile::getData -> data was flipped to north.
>>>> DEBUG 2: Range of input data (name="accu_precip";
level="(0,*,*)";) is
>> 0
>>>> to 38.4619.
>>>> DEBUG 2: Range of regridded data (name="accu_precip";
level="(0,*,*)";)
>>> is
>>>> -9999 to -9999.
>>>> DEBUG 1: Writing output file: regridded_precip_2008-08-31_23:00
>>>>
>>>
>>> Using ncdump to see the file attributes, I see this timing
information:
>>>
>>>> XTIME:standard_name = "time" ;
>>>> XTIME:units = "minutes since 2008-4-1 00:00:00" ;
>>>>
>>>
>>>
>>> However, that is not the CF-compliant way of specifying the timing
info.
>>> So, while MET will be able to read the data, the timing info won't
be
>>> accurate in the MET output files. Instead, you'll see "19700101"
which is
>>> unixtime = 0.
>>>
>>> I would suggest making your file more CF-compliant. Information
for the
>> CF
>>> coordinate time convention can be found here:
>>>
>>>
>> http://cfconventions.org/cf-conventions/v1.6.0/cf-
conventions.html#time-coordinate
>>>
>>> I hope that helps!
>>>
>>> Julie
>>>
>>>
>>> On Mon, Oct 15, 2018 at 9:09 AM Lintao Li via RT
<met_help at ucar.edu>
>>> wrote:
>>>
>>>>
>>>> Mon Oct 15 09:09:34 2018: Request 87383 was acted upon.
>>>> Transaction: Ticket created by lintaoli168 at gmail.com
>>>> Queue: met_help
>>>> Subject: NcCfFile getData error
>>>> Owner: Nobody
>>>> Requestors: lintaoli168 at gmail.com
>>>> Status: new
>>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383
>>>
>>>>
>>>>
>>>> Dear Sir/Madam,
>>>> I am Lintao, trying to use regrid_data_plane tool to regrid wrf
output
>>> data
>>>> based on StageIV data set. However, I got the error message as:
>>>> ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane
&) ->
>>>> needed 3 arguments for variable accu_precip, got 2
>>>>
>>>> I run MET on Cheyenne, using the follow code as an example:
>>>>
>>>> regrid_data_plane \
>>>>
>>>> /glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-
31_23:0
>> \
>>>>
>>>> /glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h
\
>>>>
>>>>
>>>
>>
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-
08-31_
>>>> 23:00 \
>>>>
>>>> -field 'name="accu_precip"; level="surf";' \
>>>>
>>>> -method UW_MEAN -width 4
>>>>
>>>> May I ask what is the problem please?
>>>> Please feel free to go to my working directory and perhaps check
the
>> data
>>>> files.
>>>>
>>>> Thank you in advance.
>>>>
>>>> Best regards,
>>>> Lintao
>>>>
>>>>
>>>
>>>
>>
>>
>
------------------------------------------------
Subject: NcCfFile getData error
From: Lintao Li
Time: Fri Oct 19 12:40:30 2018
Hi Julie,
I saw there is a new folder 8.0 at /glade/p/ral/jntp/MET/MET_releases.
And
it seems I can load the module by this command: "module load met/8.0".
Therefore, I am writing to ask if the regridding problem you mentioned
has
been solved already, and whether the MET version 8.0 on Cheyenne is
ready
to use?
Thank you very much.
Best regards,
Lintao
On Tue, Oct 16, 2018 at 1:42 PM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Lintao. One of our developers discovered that our MET code is
swapping
> the lat/lon values for accu_precip, so there will be missing values
after
> regridding. We are working on fixing this bug with MET. We will
release a
> patch for met-8.0, hopefully, later this week, which I will then
install on
> cheyenne. Would you like me to follow up with you once the patch
has been
> installed on cheyenne?
>
> Thank you!
>
> Julie
>
> On Mon, Oct 15, 2018 at 3:49 PM Lintao Li via RT <met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
> >
> > Hi Julie,
> > Many thanks for your very detailed explanation. It is very
helpful. I
> > appreciate it very much.
> >
> > Best regards,
> > Lintao
> >
> >
> > On Mon, Oct 15, 2018 at 2:29 PM Julie Prestopnik via RT <
> met_help at ucar.edu
> > >
> > wrote:
> >
> > > Thank you for providing the command that you ran. It is useful
in
> > > reproducing the problem. I changed your value of level from
"surf" to
> > > "(0,*,*)". For more information on this change and specifying
the
> level
> > > for a NetCDF file please see page 59 in our met-8.0 User's
Guide:
> > >
>
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v8.0.pdf
> > >
> > > Running the following command:
> > >
> > > > regrid_data_plane hrly_precip_d01_2008-08-31_23:0
ST4.2008083123.01h
> > > > regridded_precip_2008-08-31_23:00 -field 'name="accu_precip";
> > > > level="(0,*,*)";' -method UW_MEAN -width 4
> > > >
> > >
> > > I did receive an output file, however the timing information
will be
> > > incorrect. Please note the warning in the output:
> > >
> > > > DEBUG 1: Reading data file: hrly_precip_d01_2008-08-31_23:0
> > > > WARNING:
> > > > WARNING: NcCfFile::open() -> could not determine the valid
time,
> using
> > 0.
> > > > WARNING:
> > > > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 615 Ny: 559
lat_ll:
> > -112.815
> > > > lon_ll: -27.730 delta_lat: 0.005 delta_lon: 0.005
> > > > DEBUG 2: Output grid: Projection: Stereographic Nx: 1121 Ny:
881
> > > > IsNorthHemisphere: true Lon_orient: 105.000 Bx: 399.440 By:
1599.329
> > > Alpha:
> > > > 2496.0783
> > > > DEBUG 2: Interpolation options: method = UW_MEAN, width = 4,
shape =
> > > > SQUARE, vld_thresh = 0.5
> > > > DEBUG 2: NcCfFile::getData -> data was flipped to north.
> > > > DEBUG 2: Range of input data (name="accu_precip";
level="(0,*,*)";)
> is
> > 0
> > > > to 38.4619.
> > > > DEBUG 2: Range of regridded data (name="accu_precip";
> level="(0,*,*)";)
> > > is
> > > > -9999 to -9999.
> > > > DEBUG 1: Writing output file: regridded_precip_2008-08-
31_23:00
> > > >
> > >
> > > Using ncdump to see the file attributes, I see this timing
information:
> > >
> > > > XTIME:standard_name = "time" ;
> > > > XTIME:units = "minutes since 2008-4-1
00:00:00" ;
> > > >
> > >
> > >
> > > However, that is not the CF-compliant way of specifying the
timing
> info.
> > > So, while MET will be able to read the data, the timing info
won't be
> > > accurate in the MET output files. Instead, you'll see "19700101"
which
> is
> > > unixtime = 0.
> > >
> > > I would suggest making your file more CF-compliant. Information
for
> the
> > CF
> > > coordinate time convention can be found here:
> > >
> > >
> >
> http://cfconventions.org/cf-conventions/v1.6.0/cf-
conventions.html#time-coordinate
> > >
> > > I hope that helps!
> > >
> > > Julie
> > >
> > >
> > > On Mon, Oct 15, 2018 at 9:09 AM Lintao Li via RT
<met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > Mon Oct 15 09:09:34 2018: Request 87383 was acted upon.
> > > > Transaction: Ticket created by lintaoli168 at gmail.com
> > > > Queue: met_help
> > > > Subject: NcCfFile getData error
> > > > Owner: Nobody
> > > > Requestors: lintaoli168 at gmail.com
> > > > Status: new
> > > > Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383
> > >
> > > >
> > > >
> > > > Dear Sir/Madam,
> > > > I am Lintao, trying to use regrid_data_plane tool to regrid
wrf
> output
> > > data
> > > > based on StageIV data set. However, I got the error message
as:
> > > > ERROR : NcCfFile::getData(NcVar *, const LongArray &,
DataPlane &)
> ->
> > > > needed 3 arguments for variable accu_precip, got 2
> > > >
> > > > I run MET on Cheyenne, using the follow code as an example:
> > > >
> > > > regrid_data_plane \
> > > >
> > > >
> /glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-
31_23:0
> > \
> > > >
> > > >
/glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h \
> > > >
> > > >
> > >
> >
>
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-
08-31_
> > > > 23:00 \
> > > >
> > > > -field 'name="accu_precip"; level="surf";' \
> > > >
> > > > -method UW_MEAN -width 4
> > > >
> > > > May I ask what is the problem please?
> > > > Please feel free to go to my working directory and perhaps
check the
> > data
> > > > files.
> > > >
> > > > Thank you in advance.
> > > >
> > > > Best regards,
> > > > Lintao
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: NcCfFile getData error
From: John Halley Gotway
Time: Fri Oct 19 20:10:40 2018
Hi,
This is John. I work with Julie. Sorry for the delay in getting this
build of MET patched. I just posted the bug fix that should fix your
issue.
Hopefully, Julie will be able to patch this build of met-8.0 on
Monday, or
shortly thereafter.
Thanks
John
On Fri, Oct 19, 2018 at 12:41 PM Lintao Li via RT <met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
>
> Hi Julie,
>
> I saw there is a new folder 8.0 at
/glade/p/ral/jntp/MET/MET_releases. And
> it seems I can load the module by this command: "module load
met/8.0".
>
> Therefore, I am writing to ask if the regridding problem you
mentioned has
> been solved already, and whether the MET version 8.0 on Cheyenne is
ready
> to use?
>
> Thank you very much.
> Best regards,
> Lintao
>
>
> On Tue, Oct 16, 2018 at 1:42 PM Julie Prestopnik via RT
<met_help at ucar.edu
> >
> wrote:
>
> > Hi Lintao. One of our developers discovered that our MET code is
> swapping
> > the lat/lon values for accu_precip, so there will be missing
values after
> > regridding. We are working on fixing this bug with MET. We will
> release a
> > patch for met-8.0, hopefully, later this week, which I will then
install
> on
> > cheyenne. Would you like me to follow up with you once the patch
has
> been
> > installed on cheyenne?
> >
> > Thank you!
> >
> > Julie
> >
> > On Mon, Oct 15, 2018 at 3:49 PM Lintao Li via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
> > >
> > > Hi Julie,
> > > Many thanks for your very detailed explanation. It is very
helpful. I
> > > appreciate it very much.
> > >
> > > Best regards,
> > > Lintao
> > >
> > >
> > > On Mon, Oct 15, 2018 at 2:29 PM Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > >
> > > wrote:
> > >
> > > > Thank you for providing the command that you ran. It is
useful in
> > > > reproducing the problem. I changed your value of level from
"surf"
> to
> > > > "(0,*,*)". For more information on this change and specifying
the
> > level
> > > > for a NetCDF file please see page 59 in our met-8.0 User's
Guide:
> > > >
> >
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v8.0.pdf
> > > >
> > > > Running the following command:
> > > >
> > > > > regrid_data_plane hrly_precip_d01_2008-08-31_23:0
> ST4.2008083123.01h
> > > > > regridded_precip_2008-08-31_23:00 -field
'name="accu_precip";
> > > > > level="(0,*,*)";' -method UW_MEAN -width 4
> > > > >
> > > >
> > > > I did receive an output file, however the timing information
will be
> > > > incorrect. Please note the warning in the output:
> > > >
> > > > > DEBUG 1: Reading data file: hrly_precip_d01_2008-08-31_23:0
> > > > > WARNING:
> > > > > WARNING: NcCfFile::open() -> could not determine the valid
time,
> > using
> > > 0.
> > > > > WARNING:
> > > > > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 615 Ny: 559
lat_ll:
> > > -112.815
> > > > > lon_ll: -27.730 delta_lat: 0.005 delta_lon: 0.005
> > > > > DEBUG 2: Output grid: Projection: Stereographic Nx: 1121 Ny:
881
> > > > > IsNorthHemisphere: true Lon_orient: 105.000 Bx: 399.440 By:
> 1599.329
> > > > Alpha:
> > > > > 2496.0783
> > > > > DEBUG 2: Interpolation options: method = UW_MEAN, width = 4,
shape
> =
> > > > > SQUARE, vld_thresh = 0.5
> > > > > DEBUG 2: NcCfFile::getData -> data was flipped to north.
> > > > > DEBUG 2: Range of input data (name="accu_precip";
level="(0,*,*)";)
> > is
> > > 0
> > > > > to 38.4619.
> > > > > DEBUG 2: Range of regridded data (name="accu_precip";
> > level="(0,*,*)";)
> > > > is
> > > > > -9999 to -9999.
> > > > > DEBUG 1: Writing output file: regridded_precip_2008-08-
31_23:00
> > > > >
> > > >
> > > > Using ncdump to see the file attributes, I see this timing
> information:
> > > >
> > > > > XTIME:standard_name = "time" ;
> > > > > XTIME:units = "minutes since 2008-4-1
00:00:00" ;
> > > > >
> > > >
> > > >
> > > > However, that is not the CF-compliant way of specifying the
timing
> > info.
> > > > So, while MET will be able to read the data, the timing info
won't be
> > > > accurate in the MET output files. Instead, you'll see
"19700101"
> which
> > is
> > > > unixtime = 0.
> > > >
> > > > I would suggest making your file more CF-compliant.
Information for
> > the
> > > CF
> > > > coordinate time convention can be found here:
> > > >
> > > >
> > >
> >
> http://cfconventions.org/cf-conventions/v1.6.0/cf-
conventions.html#time-coordinate
> > > >
> > > > I hope that helps!
> > > >
> > > > Julie
> > > >
> > > >
> > > > On Mon, Oct 15, 2018 at 9:09 AM Lintao Li via RT
<met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > >
> > > > > Mon Oct 15 09:09:34 2018: Request 87383 was acted upon.
> > > > > Transaction: Ticket created by lintaoli168 at gmail.com
> > > > > Queue: met_help
> > > > > Subject: NcCfFile getData error
> > > > > Owner: Nobody
> > > > > Requestors: lintaoli168 at gmail.com
> > > > > Status: new
> > > > > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383
> > > >
> > > > >
> > > > >
> > > > > Dear Sir/Madam,
> > > > > I am Lintao, trying to use regrid_data_plane tool to regrid
wrf
> > output
> > > > data
> > > > > based on StageIV data set. However, I got the error message
as:
> > > > > ERROR : NcCfFile::getData(NcVar *, const LongArray &,
DataPlane &)
> > ->
> > > > > needed 3 arguments for variable accu_precip, got 2
> > > > >
> > > > > I run MET on Cheyenne, using the follow code as an example:
> > > > >
> > > > > regrid_data_plane \
> > > > >
> > > > >
> > /glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-
31_23:0
> > > \
> > > > >
> > > > >
/glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h \
> > > > >
> > > > >
> > > >
> > >
> >
>
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-
08-31_
> > > > > 23:00 \
> > > > >
> > > > > -field 'name="accu_precip"; level="surf";' \
> > > > >
> > > > > -method UW_MEAN -width 4
> > > > >
> > > > > May I ask what is the problem please?
> > > > > Please feel free to go to my working directory and perhaps
check
> the
> > > data
> > > > > files.
> > > > >
> > > > > Thank you in advance.
> > > > >
> > > > > Best regards,
> > > > > Lintao
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: NcCfFile getData error
From: Julie Prestopnik
Time: Mon Oct 22 11:08:32 2018
Hi Lintao. I have patched met-8.0 on Cheyenne. Please let us know if
you
encounter any problems.
Thanks!
Julie
On Fri, Oct 19, 2018 at 8:10 PM John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
>
> Hi,
>
> This is John. I work with Julie. Sorry for the delay in getting
this
> build of MET patched. I just posted the bug fix that should fix
your
> issue.
>
> Hopefully, Julie will be able to patch this build of met-8.0 on
Monday, or
> shortly thereafter.
>
> Thanks
> John
>
> On Fri, Oct 19, 2018 at 12:41 PM Lintao Li via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
> >
> > Hi Julie,
> >
> > I saw there is a new folder 8.0 at
/glade/p/ral/jntp/MET/MET_releases.
> And
> > it seems I can load the module by this command: "module load
met/8.0".
> >
> > Therefore, I am writing to ask if the regridding problem you
mentioned
> has
> > been solved already, and whether the MET version 8.0 on Cheyenne
is ready
> > to use?
> >
> > Thank you very much.
> > Best regards,
> > Lintao
> >
> >
> > On Tue, Oct 16, 2018 at 1:42 PM Julie Prestopnik via RT <
> met_help at ucar.edu
> > >
> > wrote:
> >
> > > Hi Lintao. One of our developers discovered that our MET code
is
> > swapping
> > > the lat/lon values for accu_precip, so there will be missing
values
> after
> > > regridding. We are working on fixing this bug with MET. We
will
> > release a
> > > patch for met-8.0, hopefully, later this week, which I will then
> install
> > on
> > > cheyenne. Would you like me to follow up with you once the
patch has
> > been
> > > installed on cheyenne?
> > >
> > > Thank you!
> > >
> > > Julie
> > >
> > > On Mon, Oct 15, 2018 at 3:49 PM Lintao Li via RT
<met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383
>
> > > >
> > > > Hi Julie,
> > > > Many thanks for your very detailed explanation. It is very
helpful. I
> > > > appreciate it very much.
> > > >
> > > > Best regards,
> > > > Lintao
> > > >
> > > >
> > > > On Mon, Oct 15, 2018 at 2:29 PM Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > >
> > > > wrote:
> > > >
> > > > > Thank you for providing the command that you ran. It is
useful in
> > > > > reproducing the problem. I changed your value of level from
"surf"
> > to
> > > > > "(0,*,*)". For more information on this change and
specifying the
> > > level
> > > > > for a NetCDF file please see page 59 in our met-8.0 User's
Guide:
> > > > >
> > >
>
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v8.0.pdf
> > > > >
> > > > > Running the following command:
> > > > >
> > > > > > regrid_data_plane hrly_precip_d01_2008-08-31_23:0
> > ST4.2008083123.01h
> > > > > > regridded_precip_2008-08-31_23:00 -field
'name="accu_precip";
> > > > > > level="(0,*,*)";' -method UW_MEAN -width 4
> > > > > >
> > > > >
> > > > > I did receive an output file, however the timing information
will
> be
> > > > > incorrect. Please note the warning in the output:
> > > > >
> > > > > > DEBUG 1: Reading data file: hrly_precip_d01_2008-08-
31_23:0
> > > > > > WARNING:
> > > > > > WARNING: NcCfFile::open() -> could not determine the valid
time,
> > > using
> > > > 0.
> > > > > > WARNING:
> > > > > > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 615 Ny: 559
lat_ll:
> > > > -112.815
> > > > > > lon_ll: -27.730 delta_lat: 0.005 delta_lon: 0.005
> > > > > > DEBUG 2: Output grid: Projection: Stereographic Nx: 1121
Ny: 881
> > > > > > IsNorthHemisphere: true Lon_orient: 105.000 Bx: 399.440
By:
> > 1599.329
> > > > > Alpha:
> > > > > > 2496.0783
> > > > > > DEBUG 2: Interpolation options: method = UW_MEAN, width =
4,
> shape
> > =
> > > > > > SQUARE, vld_thresh = 0.5
> > > > > > DEBUG 2: NcCfFile::getData -> data was flipped to north.
> > > > > > DEBUG 2: Range of input data (name="accu_precip";
> level="(0,*,*)";)
> > > is
> > > > 0
> > > > > > to 38.4619.
> > > > > > DEBUG 2: Range of regridded data (name="accu_precip";
> > > level="(0,*,*)";)
> > > > > is
> > > > > > -9999 to -9999.
> > > > > > DEBUG 1: Writing output file: regridded_precip_2008-08-
31_23:00
> > > > > >
> > > > >
> > > > > Using ncdump to see the file attributes, I see this timing
> > information:
> > > > >
> > > > > > XTIME:standard_name = "time" ;
> > > > > > XTIME:units = "minutes since 2008-4-1
00:00:00" ;
> > > > > >
> > > > >
> > > > >
> > > > > However, that is not the CF-compliant way of specifying the
timing
> > > info.
> > > > > So, while MET will be able to read the data, the timing info
won't
> be
> > > > > accurate in the MET output files. Instead, you'll see
"19700101"
> > which
> > > is
> > > > > unixtime = 0.
> > > > >
> > > > > I would suggest making your file more CF-compliant.
Information
> for
> > > the
> > > > CF
> > > > > coordinate time convention can be found here:
> > > > >
> > > > >
> > > >
> > >
> >
> http://cfconventions.org/cf-conventions/v1.6.0/cf-
conventions.html#time-coordinate
> > > > >
> > > > > I hope that helps!
> > > > >
> > > > > Julie
> > > > >
> > > > >
> > > > > On Mon, Oct 15, 2018 at 9:09 AM Lintao Li via RT <
> met_help at ucar.edu>
> > > > > wrote:
> > > > >
> > > > > >
> > > > > > Mon Oct 15 09:09:34 2018: Request 87383 was acted upon.
> > > > > > Transaction: Ticket created by lintaoli168 at gmail.com
> > > > > > Queue: met_help
> > > > > > Subject: NcCfFile getData error
> > > > > > Owner: Nobody
> > > > > > Requestors: lintaoli168 at gmail.com
> > > > > > Status: new
> > > > > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383
> > > > >
> > > > > >
> > > > > >
> > > > > > Dear Sir/Madam,
> > > > > > I am Lintao, trying to use regrid_data_plane tool to
regrid wrf
> > > output
> > > > > data
> > > > > > based on StageIV data set. However, I got the error
message as:
> > > > > > ERROR : NcCfFile::getData(NcVar *, const LongArray &,
DataPlane
> &)
> > > ->
> > > > > > needed 3 arguments for variable accu_precip, got 2
> > > > > >
> > > > > > I run MET on Cheyenne, using the follow code as an
example:
> > > > > >
> > > > > > regrid_data_plane \
> > > > > >
> > > > > >
> > > /glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-
31_23:0
> > > > \
> > > > > >
> > > > > >
> /glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h \
> > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
>
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-
08-31_
> > > > > > 23:00 \
> > > > > >
> > > > > > -field 'name="accu_precip"; level="surf";' \
> > > > > >
> > > > > > -method UW_MEAN -width 4
> > > > > >
> > > > > > May I ask what is the problem please?
> > > > > > Please feel free to go to my working directory and perhaps
check
> > the
> > > > data
> > > > > > files.
> > > > > >
> > > > > > Thank you in advance.
> > > > > >
> > > > > > Best regards,
> > > > > > Lintao
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: NcCfFile getData error
From: Lintao Li
Time: Mon Oct 22 11:54:37 2018
Hi Julie,
Hi John,
That was nice! Thanks a lot for your support! I appreciated very much.
Best regards,
Lintao
On Mon, Oct 22, 2018 at 11:08 AM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Lintao. I have patched met-8.0 on Cheyenne. Please let us know
if you
> encounter any problems.
>
> Thanks!
>
> Julie
>
> On Fri, Oct 19, 2018 at 8:10 PM John Halley Gotway via RT <
> met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
> >
> > Hi,
> >
> > This is John. I work with Julie. Sorry for the delay in getting
this
> > build of MET patched. I just posted the bug fix that should fix
your
> > issue.
> >
> > Hopefully, Julie will be able to patch this build of met-8.0 on
Monday,
> or
> > shortly thereafter.
> >
> > Thanks
> > John
> >
> > On Fri, Oct 19, 2018 at 12:41 PM Lintao Li via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
> > >
> > > Hi Julie,
> > >
> > > I saw there is a new folder 8.0 at
/glade/p/ral/jntp/MET/MET_releases.
> > And
> > > it seems I can load the module by this command: "module load
met/8.0".
> > >
> > > Therefore, I am writing to ask if the regridding problem you
mentioned
> > has
> > > been solved already, and whether the MET version 8.0 on Cheyenne
is
> ready
> > > to use?
> > >
> > > Thank you very much.
> > > Best regards,
> > > Lintao
> > >
> > >
> > > On Tue, Oct 16, 2018 at 1:42 PM Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > >
> > > wrote:
> > >
> > > > Hi Lintao. One of our developers discovered that our MET code
is
> > > swapping
> > > > the lat/lon values for accu_precip, so there will be missing
values
> > after
> > > > regridding. We are working on fixing this bug with MET. We
will
> > > release a
> > > > patch for met-8.0, hopefully, later this week, which I will
then
> > install
> > > on
> > > > cheyenne. Would you like me to follow up with you once the
patch has
> > > been
> > > > installed on cheyenne?
> > > >
> > > > Thank you!
> > > >
> > > > Julie
> > > >
> > > > On Mon, Oct 15, 2018 at 3:49 PM Lintao Li via RT
<met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383 >
> > > > >
> > > > > Hi Julie,
> > > > > Many thanks for your very detailed explanation. It is very
> helpful. I
> > > > > appreciate it very much.
> > > > >
> > > > > Best regards,
> > > > > Lintao
> > > > >
> > > > >
> > > > > On Mon, Oct 15, 2018 at 2:29 PM Julie Prestopnik via RT <
> > > > met_help at ucar.edu
> > > > > >
> > > > > wrote:
> > > > >
> > > > > > Thank you for providing the command that you ran. It is
useful
> in
> > > > > > reproducing the problem. I changed your value of level
from
> "surf"
> > > to
> > > > > > "(0,*,*)". For more information on this change and
specifying
> the
> > > > level
> > > > > > for a NetCDF file please see page 59 in our met-8.0 User's
Guide:
> > > > > >
> > > >
> >
https://dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v8.0.pdf
> > > > > >
> > > > > > Running the following command:
> > > > > >
> > > > > > > regrid_data_plane hrly_precip_d01_2008-08-31_23:0
> > > ST4.2008083123.01h
> > > > > > > regridded_precip_2008-08-31_23:00 -field
'name="accu_precip";
> > > > > > > level="(0,*,*)";' -method UW_MEAN -width 4
> > > > > > >
> > > > > >
> > > > > > I did receive an output file, however the timing
information will
> > be
> > > > > > incorrect. Please note the warning in the output:
> > > > > >
> > > > > > > DEBUG 1: Reading data file: hrly_precip_d01_2008-08-
31_23:0
> > > > > > > WARNING:
> > > > > > > WARNING: NcCfFile::open() -> could not determine the
valid
> time,
> > > > using
> > > > > 0.
> > > > > > > WARNING:
> > > > > > > DEBUG 2: Input grid: Projection: Lat/Lon Nx: 615 Ny: 559
> lat_ll:
> > > > > -112.815
> > > > > > > lon_ll: -27.730 delta_lat: 0.005 delta_lon: 0.005
> > > > > > > DEBUG 2: Output grid: Projection: Stereographic Nx: 1121
Ny:
> 881
> > > > > > > IsNorthHemisphere: true Lon_orient: 105.000 Bx: 399.440
By:
> > > 1599.329
> > > > > > Alpha:
> > > > > > > 2496.0783
> > > > > > > DEBUG 2: Interpolation options: method = UW_MEAN, width
= 4,
> > shape
> > > =
> > > > > > > SQUARE, vld_thresh = 0.5
> > > > > > > DEBUG 2: NcCfFile::getData -> data was flipped to north.
> > > > > > > DEBUG 2: Range of input data (name="accu_precip";
> > level="(0,*,*)";)
> > > > is
> > > > > 0
> > > > > > > to 38.4619.
> > > > > > > DEBUG 2: Range of regridded data (name="accu_precip";
> > > > level="(0,*,*)";)
> > > > > > is
> > > > > > > -9999 to -9999.
> > > > > > > DEBUG 1: Writing output file: regridded_precip_2008-08-
31_23:00
> > > > > > >
> > > > > >
> > > > > > Using ncdump to see the file attributes, I see this timing
> > > information:
> > > > > >
> > > > > > > XTIME:standard_name = "time" ;
> > > > > > > XTIME:units = "minutes since 2008-4-1
> 00:00:00" ;
> > > > > > >
> > > > > >
> > > > > >
> > > > > > However, that is not the CF-compliant way of specifying
the
> timing
> > > > info.
> > > > > > So, while MET will be able to read the data, the timing
info
> won't
> > be
> > > > > > accurate in the MET output files. Instead, you'll see
"19700101"
> > > which
> > > > is
> > > > > > unixtime = 0.
> > > > > >
> > > > > > I would suggest making your file more CF-compliant.
Information
> > for
> > > > the
> > > > > CF
> > > > > > coordinate time convention can be found here:
> > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
> http://cfconventions.org/cf-conventions/v1.6.0/cf-
conventions.html#time-coordinate
> > > > > >
> > > > > > I hope that helps!
> > > > > >
> > > > > > Julie
> > > > > >
> > > > > >
> > > > > > On Mon, Oct 15, 2018 at 9:09 AM Lintao Li via RT <
> > met_help at ucar.edu>
> > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > Mon Oct 15 09:09:34 2018: Request 87383 was acted upon.
> > > > > > > Transaction: Ticket created by lintaoli168 at gmail.com
> > > > > > > Queue: met_help
> > > > > > > Subject: NcCfFile getData error
> > > > > > > Owner: Nobody
> > > > > > > Requestors: lintaoli168 at gmail.com
> > > > > > > Status: new
> > > > > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=87383
> > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Dear Sir/Madam,
> > > > > > > I am Lintao, trying to use regrid_data_plane tool to
regrid wrf
> > > > output
> > > > > > data
> > > > > > > based on StageIV data set. However, I got the error
message as:
> > > > > > > ERROR : NcCfFile::getData(NcVar *, const LongArray &,
> DataPlane
> > &)
> > > > ->
> > > > > > > needed 3 arguments for variable accu_precip, got 2
> > > > > > >
> > > > > > > I run MET on Cheyenne, using the follow code as an
example:
> > > > > > >
> > > > > > > regrid_data_plane \
> > > > > > >
> > > > > > >
> > > >
> /glade/scratch/lintao/mtd_cu/fcst_data/hrly_precip_d01_2008-08-
31_23:0
> > > > > \
> > > > > > >
> > > > > > >
> > /glade/scratch/lintao/mtd_cu/oct08/test/stage4/ST4.2008083123.01h
\
> > > > > > >
> > > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
>
/glade/scratch/lintao/mtd_cu/oct08/test/regrid/regridded_precip_2008-
08-31_
> > > > > > > 23:00 \
> > > > > > >
> > > > > > > -field 'name="accu_precip"; level="surf";' \
> > > > > > >
> > > > > > > -method UW_MEAN -width 4
> > > > > > >
> > > > > > > May I ask what is the problem please?
> > > > > > > Please feel free to go to my working directory and
perhaps
> check
> > > the
> > > > > data
> > > > > > > files.
> > > > > > >
> > > > > > > Thank you in advance.
> > > > > > >
> > > > > > > Best regards,
> > > > > > > Lintao
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > >
> > > >
> > >
> > >
> >
> >
>
>
------------------------------------------------
More information about the Met_help
mailing list