[Met_help] [rt.rap.ucar.edu #62507] History for pcp_combine output creates error in MODE
John Halley Gotway via RT
met_help at ucar.edu
Mon Aug 5 15:23:16 MDT 2013
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hi,
I'm having a problem with pcp_combine output when trying to use it as
MODE input.
At first I thought it was due to the undefined gridcells in the files
I'm using as pcp_combine input. In these files (created by copygb),
there are cells containing "NANf" values. After passing it through
pcp_combine, the NANf receives the value "_".
APCP_06 =
_, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
As sanity check, I dumped the netcdf created by pcp_combine to an
ascii file and replaced all "_" by "0". Then I used ncgen to create a
new netcdf. I used this new netcdf with zeros as input for MODE but it
crashed just as before:
ERROR : process_fcst_obs_files() -> can't get data from file
"2004-08-20_06_wrf6hAcc.nc"
I know the error is coming from pcp_combine. I tested the copygb
output directly in MODE and it works.
I'm using the following command line:
***** Copygb *****
frediani at yslogin4:/glade/p/work/frediani/Validation>
$UPP_Dir/bin/copygb.exe -g"255 5 1121 881 23117 -119023 8 -105000 4763
4763 0 64" -x -k'4*-1 61' -i3
/glade/p/work/frediani/WRF_Sim/USA_NE_hornet/2004-08-20/postprd/WRFPRS_d03.006
2004-08-20_06-06_gridst4
frediani at yslogin4:/glade/p/work/frediani/Validation>
$UPP_Dir/bin/copygb.exe -g"255 5 1121 881 23117 -119023 8 -105000 4763
4763 0 64" -x -k'4*-1 61' -i3
/glade/p/work/frediani/WRF_Sim/USA_NE_hornet/2004-08-20/postprd/WRFPRS_d03.000
2004-08-20_06-00_gridst4
*****pcp_combine *****
pcp_combine -subtract 2004-08-20_06-06_gridst4 06
2004-08-20_06-00_gridst4 00 2004-08-20_06_wrf6hAcc.nc -v 10
DEBUG 1: Reading input file: 2004-08-20_06-06_gridst4
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb1".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "APCP/A6" in GRIB record 1 of GRIB file
"2004-08-20_06-06_gridst4".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "APCP/A6" in GRIB file "2004-08-20_06-06_gridst4".
DEBUG 1: Reading input file: 2004-08-20_06-00_gridst4
DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
Met2dDataFile object of type "FileType_Gb1".
DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo object
of type "FileType_Gb1".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for
VarInfo "APCP/A0" in GRIB record 1 of GRIB file
"2004-08-20_06-00_gridst4".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records
matching VarInfo "APCP/A0" in GRIB file "2004-08-20_06-00_gridst4".
DEBUG 2: Performing subtraction command.
DEBUG 1: Writing output file: 2004-08-20_06_wrf6hAcc.nc
***** Mode *****
frediani at yslogin4:/glade/p/work/frediani/Validation> mode
2004-08-20_06_wrf6hAcc.nc
GCIP-EOP-Pcp-NCEP-EMC-4KmStageIV/EOL/2004/ST4.2004082006.06h
MODEConfig
DEBUG 1: Default Config File:
/glade/u/apps/opt/METv4.0/4.0/intel/12.1.5/data/config/MODEConfig_default
DEBUG 1: Match Config File: MODEConfig
DEBUG 1: Merge Config File: MODEConfig
DEBUG 1: Forecast File: 2004-08-20_06_wrf6hAcc.nc
DEBUG 1: Observation File:
GCIP-EOP-Pcp-NCEP-EMC-4KmStageIV/EOL/2004/ST4.2004082006.06h
ERROR :
ERROR : process_fcst_obs_files() -> can't get data from file
"2004-08-20_06_wrf6hAcc.nc"
ERROR :
pcp_combine -subtract WRFPRS_d03.006_obsgrid 06 WRFPRS_d03.000_obsgrid
00 2004-08-20_06_wrf6hAcc.nc -v 10
Do you have any idea what could be going wrong?
Thanks
Maria Eugenia B. Frediani
-------------------------------------------------------------------------------------
PhD Candidate
University of Connecticut
School of Engineering
261 Glenbrook Rd
Storrs, CT 06269
maria.frediani at uconn.edu
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62507] pcp_combine output creates error in MODE
From: John Halley Gotway
Time: Mon Aug 05 12:09:01 2013
Maria,
Looking at the information you sent, it sounds like the problem is
that MODE isn't able to read the data from the file named "2004-08-
20_06_wrf6hAcc.nc". So the questions are... what does that file
contain and what does the MODEConfig configuration file contain?
I assume that 2004-08-20_06_wrf6hAcc.nc will contain a variable named
"APCP_06". But try running "ncdump -h 2004-08-20_06_wrf6hAcc.nc" to
make sure.
And for MODE since your forecast file is NetCDF and your observation
file is GRIB, the config file should contain something like this:
fcst = {
field = {
name = "APCP_06";
level = "(*,*)";
};
raw_thresh = >=0.0;
conv_radius = 5;
conv_thresh = >=2.5;
vld_thresh = 0.5;
area_thresh = >=0.0;
inten_perc_value = 100;
inten_perc_thresh = >=0.0;
merge_thresh = >=1.25;
merge_flag = THRESH;
};
obs = {
field = {
name = "APCP";
level = "A6";
};
raw_thresh = >=0.0;
conv_radius = 5;
conv_thresh = >=2.5;
vld_thresh = 0.5;
area_thresh = >=0.0;
inten_perc_value = 100;
inten_perc_thresh = >=0.0;
merge_thresh = >=1.25;
merge_flag = THRESH;
};
The "name" and "level" settings differ depending on the type for file
format being used.
Please compare this to what you have in the MODE config file and
adjust as needed. If you're still stuck, you could send me the data
you're passing to MODE, and I could figure it out here. To send
us data, you can post it on our anonymous ftp site following these
instructions:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
Lastly, your call to pcp_combine seems a little odd. Why are you
subtracting 0-hours of precip from 6-hours of precip? The precip at
forecast hour zero should be 0's everywhere anyway. You could
just use the 6-hour accumulation of precip you already have in the
GRIB file named "WRFPRS_d03.006_obsgrid".
Thanks,
John Halley Gotway
met_help at ucar.edu
On 08/05/2013 11:39 AM, Maria Eugenia via RT wrote:
>
> Mon Aug 05 11:39:03 2013: Request 62507 was acted upon.
> Transaction: Ticket created by maria.frediani at gmail.com
> Queue: met_help
> Subject: pcp_combine output creates error in MODE
> Owner: Nobody
> Requestors: maria.frediani at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62507 >
>
>
> Hi,
>
> I'm having a problem with pcp_combine output when trying to use it
as
> MODE input.
>
> At first I thought it was due to the undefined gridcells in the
files
> I'm using as pcp_combine input. In these files (created by copygb),
> there are cells containing "NANf" values. After passing it through
> pcp_combine, the NANf receives the value "_".
>
> APCP_06 =
> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
_, _, _,
>
> As sanity check, I dumped the netcdf created by pcp_combine to an
> ascii file and replaced all "_" by "0". Then I used ncgen to create
a
> new netcdf. I used this new netcdf with zeros as input for MODE but
it
> crashed just as before:
> ERROR : process_fcst_obs_files() -> can't get data from file
> "2004-08-20_06_wrf6hAcc.nc"
>
> I know the error is coming from pcp_combine. I tested the copygb
> output directly in MODE and it works.
>
> I'm using the following command line:
>
>
> ***** Copygb *****
> frediani at yslogin4:/glade/p/work/frediani/Validation>
> $UPP_Dir/bin/copygb.exe -g"255 5 1121 881 23117 -119023 8 -105000
4763
> 4763 0 64" -x -k'4*-1 61' -i3
> /glade/p/work/frediani/WRF_Sim/USA_NE_hornet/2004-08-
20/postprd/WRFPRS_d03.006
> 2004-08-20_06-06_gridst4
>
> frediani at yslogin4:/glade/p/work/frediani/Validation>
> $UPP_Dir/bin/copygb.exe -g"255 5 1121 881 23117 -119023 8 -105000
4763
> 4763 0 64" -x -k'4*-1 61' -i3
> /glade/p/work/frediani/WRF_Sim/USA_NE_hornet/2004-08-
20/postprd/WRFPRS_d03.000
> 2004-08-20_06-00_gridst4
>
> *****pcp_combine *****
> pcp_combine -subtract 2004-08-20_06-06_gridst4 06
> 2004-08-20_06-00_gridst4 00 2004-08-20_06_wrf6hAcc.nc -v 10
>
> DEBUG 1: Reading input file: 2004-08-20_06-06_gridst4
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_Gb1".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of type "FileType_Gb1".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
> VarInfo "APCP/A6" in GRIB record 1 of GRIB file
> "2004-08-20_06-06_gridst4".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
> matching VarInfo "APCP/A6" in GRIB file "2004-08-20_06-06_gridst4".
> DEBUG 1: Reading input file: 2004-08-20_06-00_gridst4
> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created new
> Met2dDataFile object of type "FileType_Gb1".
> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
> of type "FileType_Gb1".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
> VarInfo "APCP/A0" in GRIB record 1 of GRIB file
> "2004-08-20_06-00_gridst4".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
> matching VarInfo "APCP/A0" in GRIB file "2004-08-20_06-00_gridst4".
> DEBUG 2: Performing subtraction command.
> DEBUG 1: Writing output file: 2004-08-20_06_wrf6hAcc.nc
>
> ***** Mode *****
> frediani at yslogin4:/glade/p/work/frediani/Validation> mode
> 2004-08-20_06_wrf6hAcc.nc
> GCIP-EOP-Pcp-NCEP-EMC-4KmStageIV/EOL/2004/ST4.2004082006.06h
> MODEConfig
>
> DEBUG 1: Default Config File:
>
/glade/u/apps/opt/METv4.0/4.0/intel/12.1.5/data/config/MODEConfig_default
> DEBUG 1: Match Config File: MODEConfig
> DEBUG 1: Merge Config File: MODEConfig
> DEBUG 1: Forecast File: 2004-08-20_06_wrf6hAcc.nc
> DEBUG 1: Observation File:
> GCIP-EOP-Pcp-NCEP-EMC-4KmStageIV/EOL/2004/ST4.2004082006.06h
> ERROR :
> ERROR : process_fcst_obs_files() -> can't get data from file
> "2004-08-20_06_wrf6hAcc.nc"
> ERROR :
>
> pcp_combine -subtract WRFPRS_d03.006_obsgrid 06
WRFPRS_d03.000_obsgrid
> 00 2004-08-20_06_wrf6hAcc.nc -v 10
>
> Do you have any idea what could be going wrong?
>
>
> Thanks
> Maria Eugenia B. Frediani
>
-------------------------------------------------------------------------------------
> PhD Candidate
> University of Connecticut
> School of Engineering
> 261 Glenbrook Rd
> Storrs, CT 06269
> maria.frediani at uconn.edu
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62507] pcp_combine output creates error in MODE
From: Maria Eugenia
Time: Mon Aug 05 13:25:46 2013
Hi John,
The changes in the MODEConfig file solved the problem.
I'm aware that the values in the 0-hour precipitation file are all
zeros. I'm writing a bash script to do all MET tolls steps
automatically, looping at 6h intervals for the 72h of wrf simulation.
Including the 0-hour file in the loop make the script slightly
simpler, that's why.
Thanks you very much
Maria Eugenia B. Frediani
-------------------------------------------------------------------------------------
PhD Candidate
University of Connecticut
School of Engineering
261 Glenbrook Rd
Storrs, CT 06269
maria.frediani at uconn.edu
On Mon, Aug 5, 2013 at 12:09 PM, John Halley Gotway via RT
<met_help at ucar.edu> wrote:
> Maria,
>
> Looking at the information you sent, it sounds like the problem is
that MODE isn't able to read the data from the file named "2004-08-
20_06_wrf6hAcc.nc". So the questions are... what does that file
> contain and what does the MODEConfig configuration file contain?
>
> I assume that 2004-08-20_06_wrf6hAcc.nc will contain a variable
named "APCP_06". But try running "ncdump -h 2004-08-
20_06_wrf6hAcc.nc" to make sure.
>
> And for MODE since your forecast file is NetCDF and your observation
file is GRIB, the config file should contain something like this:
>
> fcst = {
> field = {
> name = "APCP_06";
> level = "(*,*)";
> };
>
> raw_thresh = >=0.0;
> conv_radius = 5;
> conv_thresh = >=2.5;
> vld_thresh = 0.5;
> area_thresh = >=0.0;
> inten_perc_value = 100;
> inten_perc_thresh = >=0.0;
> merge_thresh = >=1.25;
> merge_flag = THRESH;
> };
> obs = {
> field = {
> name = "APCP";
> level = "A6";
> };
>
> raw_thresh = >=0.0;
> conv_radius = 5;
> conv_thresh = >=2.5;
> vld_thresh = 0.5;
> area_thresh = >=0.0;
> inten_perc_value = 100;
> inten_perc_thresh = >=0.0;
> merge_thresh = >=1.25;
> merge_flag = THRESH;
> };
>
> The "name" and "level" settings differ depending on the type for
file format being used.
>
> Please compare this to what you have in the MODE config file and
adjust as needed. If you're still stuck, you could send me the data
you're passing to MODE, and I could figure it out here. To send
> us data, you can post it on our anonymous ftp site following these
instructions:
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> Lastly, your call to pcp_combine seems a little odd. Why are you
subtracting 0-hours of precip from 6-hours of precip? The precip at
forecast hour zero should be 0's everywhere anyway. You could
> just use the 6-hour accumulation of precip you already have in the
GRIB file named "WRFPRS_d03.006_obsgrid".
>
> Thanks,
> John Halley Gotway
> met_help at ucar.edu
>
>
>
> On 08/05/2013 11:39 AM, Maria Eugenia via RT wrote:
>>
>> Mon Aug 05 11:39:03 2013: Request 62507 was acted upon.
>> Transaction: Ticket created by maria.frediani at gmail.com
>> Queue: met_help
>> Subject: pcp_combine output creates error in MODE
>> Owner: Nobody
>> Requestors: maria.frediani at gmail.com
>> Status: new
>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62507 >
>>
>>
>> Hi,
>>
>> I'm having a problem with pcp_combine output when trying to use it
as
>> MODE input.
>>
>> At first I thought it was due to the undefined gridcells in the
files
>> I'm using as pcp_combine input. In these files (created by copygb),
>> there are cells containing "NANf" values. After passing it through
>> pcp_combine, the NANf receives the value "_".
>>
>> APCP_06 =
>> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
_, _, _, _,
>>
>> As sanity check, I dumped the netcdf created by pcp_combine to an
>> ascii file and replaced all "_" by "0". Then I used ncgen to create
a
>> new netcdf. I used this new netcdf with zeros as input for MODE but
it
>> crashed just as before:
>> ERROR : process_fcst_obs_files() -> can't get data from file
>> "2004-08-20_06_wrf6hAcc.nc"
>>
>> I know the error is coming from pcp_combine. I tested the copygb
>> output directly in MODE and it works.
>>
>> I'm using the following command line:
>>
>>
>> ***** Copygb *****
>> frediani at yslogin4:/glade/p/work/frediani/Validation>
>> $UPP_Dir/bin/copygb.exe -g"255 5 1121 881 23117 -119023 8 -105000
4763
>> 4763 0 64" -x -k'4*-1 61' -i3
>> /glade/p/work/frediani/WRF_Sim/USA_NE_hornet/2004-08-
20/postprd/WRFPRS_d03.006
>> 2004-08-20_06-06_gridst4
>>
>> frediani at yslogin4:/glade/p/work/frediani/Validation>
>> $UPP_Dir/bin/copygb.exe -g"255 5 1121 881 23117 -119023 8 -105000
4763
>> 4763 0 64" -x -k'4*-1 61' -i3
>> /glade/p/work/frediani/WRF_Sim/USA_NE_hornet/2004-08-
20/postprd/WRFPRS_d03.000
>> 2004-08-20_06-00_gridst4
>>
>> *****pcp_combine *****
>> pcp_combine -subtract 2004-08-20_06-06_gridst4 06
>> 2004-08-20_06-00_gridst4 00 2004-08-20_06_wrf6hAcc.nc -v 10
>>
>> DEBUG 1: Reading input file: 2004-08-20_06-06_gridst4
>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>> Met2dDataFile object of type "FileType_Gb1".
>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>> of type "FileType_Gb1".
>> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
>> VarInfo "APCP/A6" in GRIB record 1 of GRIB file
>> "2004-08-20_06-06_gridst4".
>> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
>> matching VarInfo "APCP/A6" in GRIB file "2004-08-20_06-06_gridst4".
>> DEBUG 1: Reading input file: 2004-08-20_06-00_gridst4
>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>> Met2dDataFile object of type "FileType_Gb1".
>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>> of type "FileType_Gb1".
>> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
>> VarInfo "APCP/A0" in GRIB record 1 of GRIB file
>> "2004-08-20_06-00_gridst4".
>> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
>> matching VarInfo "APCP/A0" in GRIB file "2004-08-20_06-00_gridst4".
>> DEBUG 2: Performing subtraction command.
>> DEBUG 1: Writing output file: 2004-08-20_06_wrf6hAcc.nc
>>
>> ***** Mode *****
>> frediani at yslogin4:/glade/p/work/frediani/Validation> mode
>> 2004-08-20_06_wrf6hAcc.nc
>> GCIP-EOP-Pcp-NCEP-EMC-4KmStageIV/EOL/2004/ST4.2004082006.06h
>> MODEConfig
>>
>> DEBUG 1: Default Config File:
>>
/glade/u/apps/opt/METv4.0/4.0/intel/12.1.5/data/config/MODEConfig_default
>> DEBUG 1: Match Config File: MODEConfig
>> DEBUG 1: Merge Config File: MODEConfig
>> DEBUG 1: Forecast File: 2004-08-20_06_wrf6hAcc.nc
>> DEBUG 1: Observation File:
>> GCIP-EOP-Pcp-NCEP-EMC-4KmStageIV/EOL/2004/ST4.2004082006.06h
>> ERROR :
>> ERROR : process_fcst_obs_files() -> can't get data from file
>> "2004-08-20_06_wrf6hAcc.nc"
>> ERROR :
>>
>> pcp_combine -subtract WRFPRS_d03.006_obsgrid 06
WRFPRS_d03.000_obsgrid
>> 00 2004-08-20_06_wrf6hAcc.nc -v 10
>>
>> Do you have any idea what could be going wrong?
>>
>>
>> Thanks
>> Maria Eugenia B. Frediani
>>
-------------------------------------------------------------------------------------
>> PhD Candidate
>> University of Connecticut
>> School of Engineering
>> 261 Glenbrook Rd
>> Storrs, CT 06269
>> maria.frediani at uconn.edu
>>
>
------------------------------------------------
Subject: Re: [rt.rap.ucar.edu #62507] pcp_combine output creates error in MODE
From: John Halley Gotway
Time: Mon Aug 05 15:23:03 2013
Maria,
OK, makes sense. Just let us know if any more questions arise in your
use of MET.
Thanks,
John
On 08/05/2013 01:25 PM, Maria Eugenia via RT wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62507 >
>
> Hi John,
>
> The changes in the MODEConfig file solved the problem.
>
> I'm aware that the values in the 0-hour precipitation file are all
> zeros. I'm writing a bash script to do all MET tolls steps
> automatically, looping at 6h intervals for the 72h of wrf
simulation.
> Including the 0-hour file in the loop make the script slightly
> simpler, that's why.
>
> Thanks you very much
> Maria Eugenia B. Frediani
>
-------------------------------------------------------------------------------------
> PhD Candidate
> University of Connecticut
> School of Engineering
> 261 Glenbrook Rd
> Storrs, CT 06269
> maria.frediani at uconn.edu
>
>
> On Mon, Aug 5, 2013 at 12:09 PM, John Halley Gotway via RT
> <met_help at ucar.edu> wrote:
>> Maria,
>>
>> Looking at the information you sent, it sounds like the problem is
that MODE isn't able to read the data from the file named "2004-08-
20_06_wrf6hAcc.nc". So the questions are... what does that file
>> contain and what does the MODEConfig configuration file contain?
>>
>> I assume that 2004-08-20_06_wrf6hAcc.nc will contain a variable
named "APCP_06". But try running "ncdump -h 2004-08-
20_06_wrf6hAcc.nc" to make sure.
>>
>> And for MODE since your forecast file is NetCDF and your
observation file is GRIB, the config file should contain something
like this:
>>
>> fcst = {
>> field = {
>> name = "APCP_06";
>> level = "(*,*)";
>> };
>>
>> raw_thresh = >=0.0;
>> conv_radius = 5;
>> conv_thresh = >=2.5;
>> vld_thresh = 0.5;
>> area_thresh = >=0.0;
>> inten_perc_value = 100;
>> inten_perc_thresh = >=0.0;
>> merge_thresh = >=1.25;
>> merge_flag = THRESH;
>> };
>> obs = {
>> field = {
>> name = "APCP";
>> level = "A6";
>> };
>>
>> raw_thresh = >=0.0;
>> conv_radius = 5;
>> conv_thresh = >=2.5;
>> vld_thresh = 0.5;
>> area_thresh = >=0.0;
>> inten_perc_value = 100;
>> inten_perc_thresh = >=0.0;
>> merge_thresh = >=1.25;
>> merge_flag = THRESH;
>> };
>>
>> The "name" and "level" settings differ depending on the type for
file format being used.
>>
>> Please compare this to what you have in the MODE config file and
adjust as needed. If you're still stuck, you could send me the data
you're passing to MODE, and I could figure it out here. To send
>> us data, you can post it on our anonymous ftp site following these
instructions:
>> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>>
>> Lastly, your call to pcp_combine seems a little odd. Why are you
subtracting 0-hours of precip from 6-hours of precip? The precip at
forecast hour zero should be 0's everywhere anyway. You could
>> just use the 6-hour accumulation of precip you already have in the
GRIB file named "WRFPRS_d03.006_obsgrid".
>>
>> Thanks,
>> John Halley Gotway
>> met_help at ucar.edu
>>
>>
>>
>> On 08/05/2013 11:39 AM, Maria Eugenia via RT wrote:
>>>
>>> Mon Aug 05 11:39:03 2013: Request 62507 was acted upon.
>>> Transaction: Ticket created by maria.frediani at gmail.com
>>> Queue: met_help
>>> Subject: pcp_combine output creates error in MODE
>>> Owner: Nobody
>>> Requestors: maria.frediani at gmail.com
>>> Status: new
>>> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=62507 >
>>>
>>>
>>> Hi,
>>>
>>> I'm having a problem with pcp_combine output when trying to use it
as
>>> MODE input.
>>>
>>> At first I thought it was due to the undefined gridcells in the
files
>>> I'm using as pcp_combine input. In these files (created by
copygb),
>>> there are cells containing "NANf" values. After passing it through
>>> pcp_combine, the NANf receives the value "_".
>>>
>>> APCP_06 =
>>> _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _, _,
_, _, _, _,
>>>
>>> As sanity check, I dumped the netcdf created by pcp_combine to an
>>> ascii file and replaced all "_" by "0". Then I used ncgen to
create a
>>> new netcdf. I used this new netcdf with zeros as input for MODE
but it
>>> crashed just as before:
>>> ERROR : process_fcst_obs_files() -> can't get data from file
>>> "2004-08-20_06_wrf6hAcc.nc"
>>>
>>> I know the error is coming from pcp_combine. I tested the copygb
>>> output directly in MODE and it works.
>>>
>>> I'm using the following command line:
>>>
>>>
>>> ***** Copygb *****
>>> frediani at yslogin4:/glade/p/work/frediani/Validation>
>>> $UPP_Dir/bin/copygb.exe -g"255 5 1121 881 23117 -119023 8 -105000
4763
>>> 4763 0 64" -x -k'4*-1 61' -i3
>>> /glade/p/work/frediani/WRF_Sim/USA_NE_hornet/2004-08-
20/postprd/WRFPRS_d03.006
>>> 2004-08-20_06-06_gridst4
>>>
>>> frediani at yslogin4:/glade/p/work/frediani/Validation>
>>> $UPP_Dir/bin/copygb.exe -g"255 5 1121 881 23117 -119023 8 -105000
4763
>>> 4763 0 64" -x -k'4*-1 61' -i3
>>> /glade/p/work/frediani/WRF_Sim/USA_NE_hornet/2004-08-
20/postprd/WRFPRS_d03.000
>>> 2004-08-20_06-00_gridst4
>>>
>>> *****pcp_combine *****
>>> pcp_combine -subtract 2004-08-20_06-06_gridst4 06
>>> 2004-08-20_06-00_gridst4 00 2004-08-20_06_wrf6hAcc.nc -v 10
>>>
>>> DEBUG 1: Reading input file: 2004-08-20_06-06_gridst4
>>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>>> Met2dDataFile object of type "FileType_Gb1".
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>>> of type "FileType_Gb1".
>>> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
>>> VarInfo "APCP/A6" in GRIB record 1 of GRIB file
>>> "2004-08-20_06-06_gridst4".
>>> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
>>> matching VarInfo "APCP/A6" in GRIB file "2004-08-20_06-
06_gridst4".
>>> DEBUG 1: Reading input file: 2004-08-20_06-00_gridst4
>>> DEBUG 4: Met2dDataFileFactory::new_met_2d_data_file() -> created
new
>>> Met2dDataFile object of type "FileType_Gb1".
>>> DEBUG 4: VarInfoFactory::new_var_info() -> created new VarInfo
object
>>> of type "FileType_Gb1".
>>> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
>>> VarInfo "APCP/A0" in GRIB record 1 of GRIB file
>>> "2004-08-20_06-00_gridst4".
>>> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
>>> matching VarInfo "APCP/A0" in GRIB file "2004-08-20_06-
00_gridst4".
>>> DEBUG 2: Performing subtraction command.
>>> DEBUG 1: Writing output file: 2004-08-20_06_wrf6hAcc.nc
>>>
>>> ***** Mode *****
>>> frediani at yslogin4:/glade/p/work/frediani/Validation> mode
>>> 2004-08-20_06_wrf6hAcc.nc
>>> GCIP-EOP-Pcp-NCEP-EMC-4KmStageIV/EOL/2004/ST4.2004082006.06h
>>> MODEConfig
>>>
>>> DEBUG 1: Default Config File:
>>>
/glade/u/apps/opt/METv4.0/4.0/intel/12.1.5/data/config/MODEConfig_default
>>> DEBUG 1: Match Config File: MODEConfig
>>> DEBUG 1: Merge Config File: MODEConfig
>>> DEBUG 1: Forecast File: 2004-08-20_06_wrf6hAcc.nc
>>> DEBUG 1: Observation File:
>>> GCIP-EOP-Pcp-NCEP-EMC-4KmStageIV/EOL/2004/ST4.2004082006.06h
>>> ERROR :
>>> ERROR : process_fcst_obs_files() -> can't get data from file
>>> "2004-08-20_06_wrf6hAcc.nc"
>>> ERROR :
>>>
>>> pcp_combine -subtract WRFPRS_d03.006_obsgrid 06
WRFPRS_d03.000_obsgrid
>>> 00 2004-08-20_06_wrf6hAcc.nc -v 10
>>>
>>> Do you have any idea what could be going wrong?
>>>
>>>
>>> Thanks
>>> Maria Eugenia B. Frediani
>>>
-------------------------------------------------------------------------------------
>>> PhD Candidate
>>> University of Connecticut
>>> School of Engineering
>>> 261 Glenbrook Rd
>>> Storrs, CT 06269
>>> maria.frediani at uconn.edu
>>>
>>
------------------------------------------------
More information about the Met_help
mailing list