[Met_help] [rt.rap.ucar.edu #94659] History for PCPCombine on cheyenne info
John Halley Gotway via RT
met_help at ucar.edu
Thu Apr 2 09:35:10 MDT 2020
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Please I need explanation on how to run pcp_combine (https://dtcenter.org/met-online-tutorial-metv8-0/grid-processing-tools/pcp-combine/run) on my data (/glade/work/gbode/WORK/8.0/MET_Tutorial/wrfMOD/daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc). Because I am only able to run the example.
Can I use the pcp_combine to
1. compute total precipitation
2. deaccumulate the precipitation from WRF output?
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: PCPCombine on cheyenne info
From: John Halley Gotway
Time: Wed Apr 01 11:06:44 2020
Hello,
I just ran across this MET-Help ticket and realize that we never sent
you a response. I apologize for letting this slip through the cracks.
Are you still having trouble running pcp_combine on cheyenne?
Please let me know what questions you have.
Thanks,
John Halley Gotway
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: John Halley Gotway
Time: Wed Apr 01 11:17:59 2020
Hello,
I just ran across this MET-Help ticket and realize that we never sent
you a response. I apologize for letting this slip through the cracks.
Are you still having trouble running pcp_combine on cheyenne?
Please let me know what questions you have.
Thanks,
John Halley Gotway
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: gbode Imoleayo
Time: Wed Apr 01 13:11:34 2020
Dear John,
Thanks for your reply.
I understand that so many things are happening at this time that may
cause you to miss emails.
Please I need explanation on how to run pcp_combine
(https://dtcenter.org/met-online-tutorial-metv8-0/grid-processing-
tools/pcp-combine/run) on my data
(/glade/work/gbode/WORK/8.0/MET_Tutorial/wrfMOD/daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc)
on Cheyenne. How can I sum or subtract specific period let's say
2001-03-01 10:30:00 to 2001-03-31 10:30:00 in the file.
Again, can I use the pcp_combine to;1. compute total precipitation
2. deaccumulate the precipitation from WRF output?
Thank you.
Best regards,Ezekiel
On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley Gotway
via RT <met_help at ucar.edu> wrote:
Hello,
I just ran across this MET-Help ticket and realize that we never sent
you a response. I apologize for letting this slip through the cracks.
Are you still having trouble running pcp_combine on cheyenne?
Please let me know what questions you have.
Thanks,
John Halley Gotway
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: gbode Imoleayo
Time: Wed Apr 01 13:41:59 2020
Dear John,
This is the error I get.
MET_Tutorial/wrfMOD> pcp_combine -sum 20010301_103000 1
20010331_103000 30 daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
-pcpdir pcp_combine_out/
DEBUG 2: Performing sum command: Init/In_Accum/Valid/Out_Accum Times =
20010301_103000/010000/20010331_103000/300000
DEBUG 2: Searching for 30 files with accumulation times of 010000 to
sum to a total accumulation time of 300000.
ERROR :
ERROR : sum_data_files() -> cannot find a file with a valid time of
20010331_103000 and accumulation time of 010000 matching the regular
expression ".*"
ERROR :
Thank you.
Best regards,Ezekiel
On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley Gotway
via RT <met_help at ucar.edu> wrote:
Hello,
I just ran across this MET-Help ticket and realize that we never sent
you a response. I apologize for letting this slip through the cracks.
Are you still having trouble running pcp_combine on cheyenne?
Please let me know what questions you have.
Thanks,
John Halley Gotway
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: John Halley Gotway
Time: Wed Apr 01 14:37:29 2020
Ezekiel,
I see you're having trouble running pcp_combine on a CF-compliant
NetCDF
file on cheyenne.
I logged onto chyenne and am working in:
*/glade/p/ral/jntp/MET/MET_Help/imoleayo_data_20200401*
Looks like your NetCDF file
(daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc)
is CF-compliant, is on a global lat/lon grid, and contains total
precip
(TP) for 306 timestemps.
Start by loading met-9.0:
*module use /glade/p/ral/jntp/MET/MET_releases/modulefiles*
*module load met/9.0*
Next, I ran ncview to get a quick look at the data:
*module load ncview*
*ncview daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc*
I'm surprised to see that your data occupies such a small tile of the
global domain. You could consider using a small grid that only covers
your
actual domain. Next, let's make sure MET can read this data using the
plot_data_plane tool:
*plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
'name="TP"; level="(100,*,*)";' -v 4*
DEBUG 4: Data plane information:
DEBUG 4: plane min: 0
DEBUG 4: plane max: 107.658
DEBUG 4: valid time: 20010609_103000
DEBUG 4: lead time: 000000
DEBUG 4: init time: 20010609_103000
DEBUG 4: accum time: 000000
This plots the 101-st time step since the time index (100) is 0-based.
A
png of the resulting image is attached. Is your data located
correctly
over western Africa? Running a debug level 4 (-v 4), I see that the
valid
time for this data is 20010609_103000. That timestamp can also be
used
when reading data from this file:
*plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
'name="TP"; level="(20010609_103000,*,*)";' -v 4*
The next question is what type of accumulation is contained in this
precip
data? I'm going to assume that each timestep contains data 1 days of
accumulated precip. You can use the "-add" option in pcp_combine, but
adding together 30 timesteps will be a very lengthy command line.
Here's
how you'd do 5 of them:
*pcp_combine -add \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010301_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010302_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010303_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010304_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010305_103000,*,*)";' \*
*precip_20010301_20010305.nc <http://precip_20010301_20010305.nc>*
Honestly, pcp_combine is probably not your best bet for summing up 1
month
of daily precip amounts. You could probably do this more efficiently
using
one of the CDO tools.
But basically, when running the "-add", "-subtract", or "-derive"
options,
you list the file name followed by a string defining how to extract
data
from that file. I can see how things would be much easier if you
could
define a range of dates like the following. But that is not currently
supported:
*pcp_combine -add \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010301_103000-20010331_103000,*,*)";' precip_march.nc
<http://precip_march.nc>*
Hope this helps get you going.
Thanks,
John
On Wed, Apr 1, 2020 at 1:42 PM gbode Imoleayo via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
>
> Dear John,
> This is the error I get.
> MET_Tutorial/wrfMOD> pcp_combine -sum 20010301_103000 1
20010331_103000 30
> daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc -pcpdir
pcp_combine_out/
> DEBUG 2: Performing sum command: Init/In_Accum/Valid/Out_Accum Times
=
> 20010301_103000/010000/20010331_103000/300000
> DEBUG 2: Searching for 30 files with accumulation times of 010000 to
sum
> to a total accumulation time of 300000.
> ERROR :
> ERROR : sum_data_files() -> cannot find a file with a valid time of
> 20010331_103000 and accumulation time of 010000 matching the regular
> expression ".*"
> ERROR :
>
>
> Thank you.
> Best regards,Ezekiel
>
>
> On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley
Gotway via
> RT <met_help at ucar.edu> wrote:
>
> Hello,
>
> I just ran across this MET-Help ticket and realize that we never
sent you
> a response. I apologize for letting this slip through the cracks.
Are you
> still having trouble running pcp_combine on cheyenne?
>
> Please let me know what questions you have.
>
> Thanks,
> John Halley Gotway
>
>
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: gbode Imoleayo
Time: Wed Apr 01 14:53:31 2020
Dear John,
Thanks for the very detailed explanation. I will go through these
steps one after the other and let you know the outcome.
Thanks for the png plot. Please how can I specify the long and lat
region for this type of png plot with plot_data_plane?
Thank you.
Best regards,Ezekiel
On Wednesday, April 1, 2020, 10:37:32 PM GMT+2, John Halley Gotway
via RT <met_help at ucar.edu> wrote:
Ezekiel,
I see you're having trouble running pcp_combine on a CF-compliant
NetCDF
file on cheyenne.
I logged onto chyenne and am working in:
*/glade/p/ral/jntp/MET/MET_Help/imoleayo_data_20200401*
Looks like your NetCDF file
(daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc)
is CF-compliant, is on a global lat/lon grid, and contains total
precip
(TP) for 306 timestemps.
Start by loading met-9.0:
*module use /glade/p/ral/jntp/MET/MET_releases/modulefiles*
*module load met/9.0*
Next, I ran ncview to get a quick look at the data:
*module load ncview*
*ncview daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc*
I'm surprised to see that your data occupies such a small tile of the
global domain. You could consider using a small grid that only covers
your
actual domain. Next, let's make sure MET can read this data using the
plot_data_plane tool:
*plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
'name="TP"; level="(100,*,*)";' -v 4*
DEBUG 4: Data plane information:
DEBUG 4: plane min: 0
DEBUG 4: plane max: 107.658
DEBUG 4: valid time: 20010609_103000
DEBUG 4: lead time: 000000
DEBUG 4: init time: 20010609_103000
DEBUG 4: accum time: 000000
This plots the 101-st time step since the time index (100) is 0-
based. A
png of the resulting image is attached. Is your data located
correctly
over western Africa? Running a debug level 4 (-v 4), I see that the
valid
time for this data is 20010609_103000. That timestamp can also be
used
when reading data from this file:
*plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
'name="TP"; level="(20010609_103000,*,*)";' -v 4*
The next question is what type of accumulation is contained in this
precip
data? I'm going to assume that each timestep contains data 1 days of
accumulated precip. You can use the "-add" option in pcp_combine, but
adding together 30 timesteps will be a very lengthy command line.
Here's
how you'd do 5 of them:
*pcp_combine -add \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010301_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010302_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010303_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010304_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010305_103000,*,*)";' \*
*precip_20010301_20010305.nc <http://precip_20010301_20010305.nc>*
Honestly, pcp_combine is probably not your best bet for summing up 1
month
of daily precip amounts. You could probably do this more efficiently
using
one of the CDO tools.
But basically, when running the "-add", "-subtract", or "-derive"
options,
you list the file name followed by a string defining how to extract
data
from that file. I can see how things would be much easier if you
could
define a range of dates like the following. But that is not currently
supported:
*pcp_combine -add \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010301_103000-20010331_103000,*,*)";' precip_march.nc
<http://precip_march.nc>*
Hope this helps get you going.
Thanks,
John
On Wed, Apr 1, 2020 at 1:42 PM gbode Imoleayo via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
>
> Dear John,
> This is the error I get.
> MET_Tutorial/wrfMOD> pcp_combine -sum 20010301_103000 1
20010331_103000 30
> daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc -pcpdir
pcp_combine_out/
> DEBUG 2: Performing sum command: Init/In_Accum/Valid/Out_Accum Times
=
> 20010301_103000/010000/20010331_103000/300000
> DEBUG 2: Searching for 30 files with accumulation times of 010000 to
sum
> to a total accumulation time of 300000.
> ERROR :
> ERROR : sum_data_files() -> cannot find a file with a valid time of
> 20010331_103000 and accumulation time of 010000 matching the regular
> expression ".*"
> ERROR :
>
>
> Thank you.
> Best regards,Ezekiel
>
>
> On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley Gotway
via
> RT <met_help at ucar.edu> wrote:
>
> Hello,
>
> I just ran across this MET-Help ticket and realize that we never
sent you
> a response. I apologize for letting this slip through the cracks.
Are you
> still having trouble running pcp_combine on cheyenne?
>
> Please let me know what questions you have.
>
> Thanks,
> John Halley Gotway
>
>
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: gbode Imoleayo
Time: Wed Apr 01 14:57:30 2020
Dear John,
Yes, my data is well situated over West Africa.
Thank you.
Best regards,Ezekiel
On Wednesday, April 1, 2020, 10:53:25 PM GMT+2, gbode Imoleayo
<gbodeimoleayo at yahoo.com> wrote:
Dear John,
Thanks for the very detailed explanation. I will go through these
steps one after the other and let you know the outcome.
Thanks for the png plot. Please how can I specify the long and lat
region for this type of png plot with plot_data_plane?
Thank you.
Best regards,Ezekiel
On Wednesday, April 1, 2020, 10:37:32 PM GMT+2, John Halley Gotway
via RT <met_help at ucar.edu> wrote:
Ezekiel,
I see you're having trouble running pcp_combine on a CF-compliant
NetCDF
file on cheyenne.
I logged onto chyenne and am working in:
*/glade/p/ral/jntp/MET/MET_Help/imoleayo_data_20200401*
Looks like your NetCDF file
(daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc)
is CF-compliant, is on a global lat/lon grid, and contains total
precip
(TP) for 306 timestemps.
Start by loading met-9.0:
*module use /glade/p/ral/jntp/MET/MET_releases/modulefiles*
*module load met/9.0*
Next, I ran ncview to get a quick look at the data:
*module load ncview*
*ncview daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc*
I'm surprised to see that your data occupies such a small tile of the
global domain. You could consider using a small grid that only covers
your
actual domain. Next, let's make sure MET can read this data using the
plot_data_plane tool:
*plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
'name="TP"; level="(100,*,*)";' -v 4*
DEBUG 4: Data plane information:
DEBUG 4: plane min: 0
DEBUG 4: plane max: 107.658
DEBUG 4: valid time: 20010609_103000
DEBUG 4: lead time: 000000
DEBUG 4: init time: 20010609_103000
DEBUG 4: accum time: 000000
This plots the 101-st time step since the time index (100) is 0-
based. A
png of the resulting image is attached. Is your data located
correctly
over western Africa? Running a debug level 4 (-v 4), I see that the
valid
time for this data is 20010609_103000. That timestamp can also be
used
when reading data from this file:
*plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
'name="TP"; level="(20010609_103000,*,*)";' -v 4*
The next question is what type of accumulation is contained in this
precip
data? I'm going to assume that each timestep contains data 1 days of
accumulated precip. You can use the "-add" option in pcp_combine, but
adding together 30 timesteps will be a very lengthy command line.
Here's
how you'd do 5 of them:
*pcp_combine -add \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010301_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010302_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010303_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010304_103000,*,*)";' \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010305_103000,*,*)";' \*
*precip_20010301_20010305.nc <http://precip_20010301_20010305.nc>*
Honestly, pcp_combine is probably not your best bet for summing up 1
month
of daily precip amounts. You could probably do this more efficiently
using
one of the CDO tools.
But basically, when running the "-add", "-subtract", or "-derive"
options,
you list the file name followed by a string defining how to extract
data
from that file. I can see how things would be much easier if you
could
define a range of dates like the following. But that is not currently
supported:
*pcp_combine -add \*
*daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
level="(20010301_103000-20010331_103000,*,*)";' precip_march.nc
<http://precip_march.nc>*
Hope this helps get you going.
Thanks,
John
On Wed, Apr 1, 2020 at 1:42 PM gbode Imoleayo via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
>
> Dear John,
> This is the error I get.
> MET_Tutorial/wrfMOD> pcp_combine -sum 20010301_103000 1
20010331_103000 30
> daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc -pcpdir
pcp_combine_out/
> DEBUG 2: Performing sum command: Init/In_Accum/Valid/Out_Accum Times
=
> 20010301_103000/010000/20010331_103000/300000
> DEBUG 2: Searching for 30 files with accumulation times of 010000 to
sum
> to a total accumulation time of 300000.
> ERROR :
> ERROR : sum_data_files() -> cannot find a file with a valid time of
> 20010331_103000 and accumulation time of 010000 matching the regular
> expression ".*"
> ERROR :
>
>
> Thank you.
> Best regards,Ezekiel
>
>
> On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley Gotway
via
> RT <met_help at ucar.edu> wrote:
>
> Hello,
>
> I just ran across this MET-Help ticket and realize that we never
sent you
> a response. I apologize for letting this slip through the cracks.
Are you
> still having trouble running pcp_combine on cheyenne?
>
> Please let me know what questions you have.
>
> Thanks,
> John Halley Gotway
>
>
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: John Halley Gotway
Time: Wed Apr 01 16:40:58 2020
Ezekial,
The plot_data_plane utility in MET just reads the projection
information
from the input NetCDF file and plots the data on that projection.
There is
no way of changing the region that's plotted by plot_data_plane. We
do not
intend it to be a fully functional plotting utility. Instead, we use
it as
a sanity check to make sure that MET is reading your data correctly
and
placing in on the correct place on the globe.
However, MET does include functionality for regridding data. This can
be
done it 2 ways:
(1) Directly in the configuration files of the tools using the
"regrid"
dictionary.
(2) Or by running the stand-alone regrid_data_plane tool.
You can find examples of this in the online tutorial.
https://dtcenter.org/met-online-tutorial-metv8-0/regridding/regrid-
data-plane
Thanks,
John
On Wed, Apr 1, 2020 at 3:06 PM gbode Imoleayo via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
>
> Dear John,
> Yes, my data is well situated over West Africa.
> Thank you.
> Best regards,Ezekiel
>
> On Wednesday, April 1, 2020, 10:53:25 PM GMT+2, gbode Imoleayo <
> gbodeimoleayo at yahoo.com> wrote:
>
> Dear John,
> Thanks for the very detailed explanation. I will go through these
steps
> one after the other and let you know the outcome.
>
> Thanks for the png plot. Please how can I specify the long and lat
region
> for this type of png plot with plot_data_plane?
> Thank you.
>
> Best regards,Ezekiel
>
>
>
> On Wednesday, April 1, 2020, 10:37:32 PM GMT+2, John Halley
Gotway via
> RT <met_help at ucar.edu> wrote:
>
> Ezekiel,
>
> I see you're having trouble running pcp_combine on a CF-compliant
NetCDF
> file on cheyenne.
>
> I logged onto chyenne and am working in:
> */glade/p/ral/jntp/MET/MET_Help/imoleayo_data_20200401*
>
> Looks like your NetCDF file
(daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc)
> is CF-compliant, is on a global lat/lon grid, and contains total
precip
> (TP) for 306 timestemps.
>
> Start by loading met-9.0:
>
> *module use /glade/p/ral/jntp/MET/MET_releases/modulefiles*
>
> *module load met/9.0*
>
> Next, I ran ncview to get a quick look at the data:
> *module load ncview*
> *ncview daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc*
>
> I'm surprised to see that your data occupies such a small tile of
the
> global domain. You could consider using a small grid that only
covers your
> actual domain. Next, let's make sure MET can read this data using
the
> plot_data_plane tool:
>
> *plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
> 'name="TP"; level="(100,*,*)";' -v 4*
> DEBUG 4: Data plane information:
> DEBUG 4: plane min: 0
> DEBUG 4: plane max: 107.658
> DEBUG 4: valid time: 20010609_103000
> DEBUG 4: lead time: 000000
> DEBUG 4: init time: 20010609_103000
> DEBUG 4: accum time: 000000
>
> This plots the 101-st time step since the time index (100) is 0-
based. A
> png of the resulting image is attached. Is your data located
correctly
> over western Africa? Running a debug level 4 (-v 4), I see that the
valid
> time for this data is 20010609_103000. That timestamp can also be
used
> when reading data from this file:
>
> *plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
> 'name="TP"; level="(20010609_103000,*,*)";' -v 4*
>
> The next question is what type of accumulation is contained in this
precip
> data? I'm going to assume that each timestep contains data 1 days
of
> accumulated precip. You can use the "-add" option in pcp_combine,
but
> adding together 30 timesteps will be a very lengthy command line.
Here's
> how you'd do 5 of them:
>
> *pcp_combine -add \*
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010301_103000,*,*)";' \*
>
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010302_103000,*,*)";' \*
>
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010303_103000,*,*)";' \*
>
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010304_103000,*,*)";' \*
>
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010305_103000,*,*)";' \*
> *precip_20010301_20010305.nc <http://precip_20010301_20010305.nc>*
>
> Honestly, pcp_combine is probably not your best bet for summing up 1
month
> of daily precip amounts. You could probably do this more
efficiently using
> one of the CDO tools.
>
> But basically, when running the "-add", "-subtract", or "-derive"
options,
> you list the file name followed by a string defining how to extract
data
> from that file. I can see how things would be much easier if you
could
> define a range of dates like the following. But that is not
currently
> supported:
> *pcp_combine -add \*
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010301_103000-20010331_103000,*,*)";' precip_march.nc
> <http://precip_march.nc>*
>
> Hope this helps get you going.
>
> Thanks,
> John
>
> On Wed, Apr 1, 2020 at 1:42 PM gbode Imoleayo via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
> >
> > Dear John,
> > This is the error I get.
> > MET_Tutorial/wrfMOD> pcp_combine -sum 20010301_103000 1
20010331_103000
> 30
> > daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc -pcpdir
pcp_combine_out/
> > DEBUG 2: Performing sum command: Init/In_Accum/Valid/Out_Accum
Times =
> > 20010301_103000/010000/20010331_103000/300000
> > DEBUG 2: Searching for 30 files with accumulation times of 010000
to sum
> > to a total accumulation time of 300000.
> > ERROR :
> > ERROR : sum_data_files() -> cannot find a file with a valid time
of
> > 20010331_103000 and accumulation time of 010000 matching the
regular
> > expression ".*"
> > ERROR :
> >
> >
> > Thank you.
> > Best regards,Ezekiel
> >
> >
> > On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley
Gotway via
> > RT <met_help at ucar.edu> wrote:
> >
> > Hello,
> >
> > I just ran across this MET-Help ticket and realize that we never
sent you
> > a response. I apologize for letting this slip through the cracks.
Are you
> > still having trouble running pcp_combine on cheyenne?
> >
> > Please let me know what questions you have.
> >
> > Thanks,
> > John Halley Gotway
> >
> >
>
>
>
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: gbode Imoleayo
Time: Wed Apr 01 23:50:37 2020
Dear John,
Thanks for the explanation.
Is it possible to specify the region (e.g. West Africa) when running
"grid_stat"?
Thank you.
Best regards,Ezekiel
On Thursday, April 2, 2020, 12:41:00 AM GMT+2, John Halley Gotway
via RT <met_help at ucar.edu> wrote:
Ezekial,
The plot_data_plane utility in MET just reads the projection
information
from the input NetCDF file and plots the data on that projection.
There is
no way of changing the region that's plotted by plot_data_plane. We
do not
intend it to be a fully functional plotting utility. Instead, we use
it as
a sanity check to make sure that MET is reading your data correctly
and
placing in on the correct place on the globe.
However, MET does include functionality for regridding data. This can
be
done it 2 ways:
(1) Directly in the configuration files of the tools using the
"regrid"
dictionary.
(2) Or by running the stand-alone regrid_data_plane tool.
You can find examples of this in the online tutorial.
https://dtcenter.org/met-online-tutorial-metv8-0/regridding/regrid-
data-plane
Thanks,
John
On Wed, Apr 1, 2020 at 3:06 PM gbode Imoleayo via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
>
> Dear John,
> Yes, my data is well situated over West Africa.
> Thank you.
> Best regards,Ezekiel
>
> On Wednesday, April 1, 2020, 10:53:25 PM GMT+2, gbode Imoleayo <
> gbodeimoleayo at yahoo.com> wrote:
>
> Dear John,
> Thanks for the very detailed explanation. I will go through these
steps
> one after the other and let you know the outcome.
>
> Thanks for the png plot. Please how can I specify the long and lat
region
> for this type of png plot with plot_data_plane?
> Thank you.
>
> Best regards,Ezekiel
>
>
>
> On Wednesday, April 1, 2020, 10:37:32 PM GMT+2, John Halley
Gotway via
> RT <met_help at ucar.edu> wrote:
>
> Ezekiel,
>
> I see you're having trouble running pcp_combine on a CF-compliant
NetCDF
> file on cheyenne.
>
> I logged onto chyenne and am working in:
> */glade/p/ral/jntp/MET/MET_Help/imoleayo_data_20200401*
>
> Looks like your NetCDF file
(daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc)
> is CF-compliant, is on a global lat/lon grid, and contains total
precip
> (TP) for 306 timestemps.
>
> Start by loading met-9.0:
>
> *module use /glade/p/ral/jntp/MET/MET_releases/modulefiles*
>
> *module load met/9.0*
>
> Next, I ran ncview to get a quick look at the data:
> *module load ncview*
> *ncview daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc*
>
> I'm surprised to see that your data occupies such a small tile of
the
> global domain. You could consider using a small grid that only
covers your
> actual domain. Next, let's make sure MET can read this data using
the
> plot_data_plane tool:
>
> *plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
> 'name="TP"; level="(100,*,*)";' -v 4*
> DEBUG 4: Data plane information:
> DEBUG 4: plane min: 0
> DEBUG 4: plane max: 107.658
> DEBUG 4: valid time: 20010609_103000
> DEBUG 4: lead time: 000000
> DEBUG 4: init time: 20010609_103000
> DEBUG 4: accum time: 000000
>
> This plots the 101-st time step since the time index (100) is 0-
based. A
> png of the resulting image is attached. Is your data located
correctly
> over western Africa? Running a debug level 4 (-v 4), I see that the
valid
> time for this data is 20010609_103000. That timestamp can also be
used
> when reading data from this file:
>
> *plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
> 'name="TP"; level="(20010609_103000,*,*)";' -v 4*
>
> The next question is what type of accumulation is contained in this
precip
> data? I'm going to assume that each timestep contains data 1 days
of
> accumulated precip. You can use the "-add" option in pcp_combine,
but
> adding together 30 timesteps will be a very lengthy command line.
Here's
> how you'd do 5 of them:
>
> *pcp_combine -add \*
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010301_103000,*,*)";' \*
>
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010302_103000,*,*)";' \*
>
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010303_103000,*,*)";' \*
>
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010304_103000,*,*)";' \*
>
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010305_103000,*,*)";' \*
> *precip_20010301_20010305.nc <http://precip_20010301_20010305.nc>*
>
> Honestly, pcp_combine is probably not your best bet for summing up 1
month
> of daily precip amounts. You could probably do this more
efficiently using
> one of the CDO tools.
>
> But basically, when running the "-add", "-subtract", or "-derive"
options,
> you list the file name followed by a string defining how to extract
data
> from that file. I can see how things would be much easier if you
could
> define a range of dates like the following. But that is not
currently
> supported:
> *pcp_combine -add \*
> *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> level="(20010301_103000-20010331_103000,*,*)";' precip_march.nc
> <http://precip_march.nc>*
>
> Hope this helps get you going.
>
> Thanks,
> John
>
> On Wed, Apr 1, 2020 at 1:42 PM gbode Imoleayo via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
> >
> > Dear John,
> > This is the error I get.
> > MET_Tutorial/wrfMOD> pcp_combine -sum 20010301_103000 1
20010331_103000
> 30
> > daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc -pcpdir
pcp_combine_out/
> > DEBUG 2: Performing sum command: Init/In_Accum/Valid/Out_Accum
Times =
> > 20010301_103000/010000/20010331_103000/300000
> > DEBUG 2: Searching for 30 files with accumulation times of 010000
to sum
> > to a total accumulation time of 300000.
> > ERROR :
> > ERROR : sum_data_files() -> cannot find a file with a valid time
of
> > 20010331_103000 and accumulation time of 010000 matching the
regular
> > expression ".*"
> > ERROR :
> >
> >
> > Thank you.
> > Best regards,Ezekiel
> >
> >
> > On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley
Gotway via
> > RT <met_help at ucar.edu> wrote:
> >
> > Hello,
> >
> > I just ran across this MET-Help ticket and realize that we never
sent you
> > a response. I apologize for letting this slip through the cracks.
Are you
> > still having trouble running pcp_combine on cheyenne?
> >
> > Please let me know what questions you have.
> >
> > Thanks,
> > John Halley Gotway
> >
> >
>
>
>
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: John Halley Gotway
Time: Thu Apr 02 07:15:24 2020
Sure, we call defining spatial sub-regions “masking”. And you
configure
Grid-Stat to compute statistics for however many masking regions you’d
like. You just list them in the “mask” dictionary of the Grid-Stat
configuration file.
In order to generate those masks, I’d recommend running the Gen-Vx-
Mask
tool. It supports multiple methods for defining them. Take a look at
the
online tutorial for examples:
https://dtcenter.org/met-online-tutorial-metv8-0/grid-processing-
tools/gen-vx-mask
Thanks
John
On Wed, Apr 1, 2020 at 11:51 PM gbode Imoleayo via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
>
> Dear John,
> Thanks for the explanation.
> Is it possible to specify the region (e.g. West Africa) when running
> "grid_stat"?
> Thank you.
> Best regards,Ezekiel
>
> On Thursday, April 2, 2020, 12:41:00 AM GMT+2, John Halley
Gotway via
> RT <met_help at ucar.edu> wrote:
>
> Ezekial,
>
> The plot_data_plane utility in MET just reads the projection
information
> from the input NetCDF file and plots the data on that projection.
There is
> no way of changing the region that's plotted by plot_data_plane. We
do not
> intend it to be a fully functional plotting utility. Instead, we
use it as
> a sanity check to make sure that MET is reading your data correctly
and
> placing in on the correct place on the globe.
>
> However, MET does include functionality for regridding data. This
can be
> done it 2 ways:
>
> (1) Directly in the configuration files of the tools using the
"regrid"
> dictionary.
> (2) Or by running the stand-alone regrid_data_plane tool.
>
> You can find examples of this in the online tutorial.
>
> https://dtcenter.org/met-online-tutorial-metv8-0/regridding/regrid-
data-plane
>
> Thanks,
> John
>
> On Wed, Apr 1, 2020 at 3:06 PM gbode Imoleayo via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
> >
> > Dear John,
> > Yes, my data is well situated over West Africa.
> > Thank you.
> > Best regards,Ezekiel
> >
> > On Wednesday, April 1, 2020, 10:53:25 PM GMT+2, gbode Imoleayo
<
> > gbodeimoleayo at yahoo.com> wrote:
> >
> > Dear John,
> > Thanks for the very detailed explanation. I will go through these
steps
> > one after the other and let you know the outcome.
> >
> > Thanks for the png plot. Please how can I specify the long and lat
region
> > for this type of png plot with plot_data_plane?
> > Thank you.
> >
> > Best regards,Ezekiel
> >
> >
> >
> > On Wednesday, April 1, 2020, 10:37:32 PM GMT+2, John Halley
Gotway via
> > RT <met_help at ucar.edu> wrote:
> >
> > Ezekiel,
> >
> > I see you're having trouble running pcp_combine on a CF-compliant
NetCDF
> > file on cheyenne.
> >
> > I logged onto chyenne and am working in:
> > */glade/p/ral/jntp/MET/MET_Help/imoleayo_data_20200401*
> >
> > Looks like your NetCDF file
> (daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc)
> > is CF-compliant, is on a global lat/lon grid, and contains total
precip
> > (TP) for 306 timestemps.
> >
> > Start by loading met-9.0:
> >
> > *module use /glade/p/ral/jntp/MET/MET_releases/modulefiles*
> >
> > *module load met/9.0*
> >
> > Next, I ran ncview to get a quick look at the data:
> > *module load ncview*
> > *ncview daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc*
> >
> > I'm surprised to see that your data occupies such a small tile of
the
> > global domain. You could consider using a small grid that only
covers
> your
> > actual domain. Next, let's make sure MET can read this data using
the
> > plot_data_plane tool:
> >
> > *plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
> > 'name="TP"; level="(100,*,*)";' -v 4*
> > DEBUG 4: Data plane information:
> > DEBUG 4: plane min: 0
> > DEBUG 4: plane max: 107.658
> > DEBUG 4: valid time: 20010609_103000
> > DEBUG 4: lead time: 000000
> > DEBUG 4: init time: 20010609_103000
> > DEBUG 4: accum time: 000000
> >
> > This plots the 101-st time step since the time index (100) is 0-
based. A
> > png of the resulting image is attached. Is your data located
correctly
> > over western Africa? Running a debug level 4 (-v 4), I see that
the
> valid
> > time for this data is 20010609_103000. That timestamp can also be
used
> > when reading data from this file:
> >
> > *plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
> > 'name="TP"; level="(20010609_103000,*,*)";' -v 4*
> >
> > The next question is what type of accumulation is contained in
this
> precip
> > data? I'm going to assume that each timestep contains data 1 days
of
> > accumulated precip. You can use the "-add" option in pcp_combine,
but
> > adding together 30 timesteps will be a very lengthy command line.
Here's
> > how you'd do 5 of them:
> >
> > *pcp_combine -add \*
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010301_103000,*,*)";' \*
> >
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010302_103000,*,*)";' \*
> >
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010303_103000,*,*)";' \*
> >
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010304_103000,*,*)";' \*
> >
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010305_103000,*,*)";' \*
> > *precip_20010301_20010305.nc <http://precip_20010301_20010305.nc>*
> >
> > Honestly, pcp_combine is probably not your best bet for summing up
1
> month
> > of daily precip amounts. You could probably do this more
efficiently
> using
> > one of the CDO tools.
> >
> > But basically, when running the "-add", "-subtract", or "-derive"
> options,
> > you list the file name followed by a string defining how to
extract data
> > from that file. I can see how things would be much easier if you
could
> > define a range of dates like the following. But that is not
currently
> > supported:
> > *pcp_combine -add \*
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010301_103000-20010331_103000,*,*)";' precip_march.nc
> > <http://precip_march.nc>*
> >
> > Hope this helps get you going.
> >
> > Thanks,
> > John
> >
> > On Wed, Apr 1, 2020 at 1:42 PM gbode Imoleayo via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
> > >
> > > Dear John,
> > > This is the error I get.
> > > MET_Tutorial/wrfMOD> pcp_combine -sum 20010301_103000 1
20010331_103000
> > 30
> > > daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc -pcpdir
pcp_combine_out/
> > > DEBUG 2: Performing sum command: Init/In_Accum/Valid/Out_Accum
Times =
> > > 20010301_103000/010000/20010331_103000/300000
> > > DEBUG 2: Searching for 30 files with accumulation times of
010000 to
> sum
> > > to a total accumulation time of 300000.
> > > ERROR :
> > > ERROR : sum_data_files() -> cannot find a file with a valid
time of
> > > 20010331_103000 and accumulation time of 010000 matching the
regular
> > > expression ".*"
> > > ERROR :
> > >
> > >
> > > Thank you.
> > > Best regards,Ezekiel
> > >
> > >
> > > On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley
Gotway
> via
> > > RT <met_help at ucar.edu> wrote:
> > >
> > > Hello,
> > >
> > > I just ran across this MET-Help ticket and realize that we never
sent
> you
> > > a response. I apologize for letting this slip through the
cracks. Are
> you
> > > still having trouble running pcp_combine on cheyenne?
> > >
> > > Please let me know what questions you have.
> > >
> > > Thanks,
> > > John Halley Gotway
> > >
> > >
> >
> >
> >
>
>
>
------------------------------------------------
Subject: PCPCombine on cheyenne info
From: gbode Imoleayo
Time: Thu Apr 02 07:18:36 2020
Dear John,
Thanks for your help and the suggestion.
Best regards,Ezekiel
On Thursday, April 2, 2020, 3:15:26 PM GMT+2, John Halley Gotway
via RT <met_help at ucar.edu> wrote:
Sure, we call defining spatial sub-regions “masking”. And you
configure
Grid-Stat to compute statistics for however many masking regions you’d
like. You just list them in the “mask” dictionary of the Grid-Stat
configuration file.
In order to generate those masks, I’d recommend running the Gen-Vx-
Mask
tool. It supports multiple methods for defining them. Take a look at
the
online tutorial for examples:
https://dtcenter.org/met-online-tutorial-metv8-0/grid-processing-
tools/gen-vx-mask
Thanks
John
On Wed, Apr 1, 2020 at 11:51 PM gbode Imoleayo via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
>
> Dear John,
> Thanks for the explanation.
> Is it possible to specify the region (e.g. West Africa) when running
> "grid_stat"?
> Thank you.
> Best regards,Ezekiel
>
> On Thursday, April 2, 2020, 12:41:00 AM GMT+2, John Halley Gotway
via
> RT <met_help at ucar.edu> wrote:
>
> Ezekial,
>
> The plot_data_plane utility in MET just reads the projection
information
> from the input NetCDF file and plots the data on that projection.
There is
> no way of changing the region that's plotted by plot_data_plane. We
do not
> intend it to be a fully functional plotting utility. Instead, we
use it as
> a sanity check to make sure that MET is reading your data correctly
and
> placing in on the correct place on the globe.
>
> However, MET does include functionality for regridding data. This
can be
> done it 2 ways:
>
> (1) Directly in the configuration files of the tools using the
"regrid"
> dictionary.
> (2) Or by running the stand-alone regrid_data_plane tool.
>
> You can find examples of this in the online tutorial.
>
> https://dtcenter.org/met-online-tutorial-metv8-0/regridding/regrid-
data-plane
>
> Thanks,
> John
>
> On Wed, Apr 1, 2020 at 3:06 PM gbode Imoleayo via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
> >
> > Dear John,
> > Yes, my data is well situated over West Africa.
> > Thank you.
> > Best regards,Ezekiel
> >
> > On Wednesday, April 1, 2020, 10:53:25 PM GMT+2, gbode Imoleayo
<
> > gbodeimoleayo at yahoo.com> wrote:
> >
> > Dear John,
> > Thanks for the very detailed explanation. I will go through these
steps
> > one after the other and let you know the outcome.
> >
> > Thanks for the png plot. Please how can I specify the long and lat
region
> > for this type of png plot with plot_data_plane?
> > Thank you.
> >
> > Best regards,Ezekiel
> >
> >
> >
> > On Wednesday, April 1, 2020, 10:37:32 PM GMT+2, John Halley
Gotway via
> > RT <met_help at ucar.edu> wrote:
> >
> > Ezekiel,
> >
> > I see you're having trouble running pcp_combine on a CF-compliant
NetCDF
> > file on cheyenne.
> >
> > I logged onto chyenne and am working in:
> > */glade/p/ral/jntp/MET/MET_Help/imoleayo_data_20200401*
> >
> > Looks like your NetCDF file
> (daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc)
> > is CF-compliant, is on a global lat/lon grid, and contains total
precip
> > (TP) for 306 timestemps.
> >
> > Start by loading met-9.0:
> >
> > *module use /glade/p/ral/jntp/MET/MET_releases/modulefiles*
> >
> > *module load met/9.0*
> >
> > Next, I ran ncview to get a quick look at the data:
> > *module load ncview*
> > *ncview daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc*
> >
> > I'm surprised to see that your data occupies such a small tile of
the
> > global domain. You could consider using a small grid that only
covers
> your
> > actual domain. Next, let's make sure MET can read this data using
the
> > plot_data_plane tool:
> >
> > *plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
> > 'name="TP"; level="(100,*,*)";' -v 4*
> > DEBUG 4: Data plane information:
> > DEBUG 4: plane min: 0
> > DEBUG 4: plane max: 107.658
> > DEBUG 4: valid time: 20010609_103000
> > DEBUG 4: lead time: 000000
> > DEBUG 4: init time: 20010609_103000
> > DEBUG 4: accum time: 000000
> >
> > This plots the 101-st time step since the time index (100) is 0-
based. A
> > png of the resulting image is attached. Is your data located
correctly
> > over western Africa? Running a debug level 4 (-v 4), I see that
the
> valid
> > time for this data is 20010609_103000. That timestamp can also be
used
> > when reading data from this file:
> >
> > *plot_data_plane daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc
plot_TP.ps
> > 'name="TP"; level="(20010609_103000,*,*)";' -v 4*
> >
> > The next question is what type of accumulation is contained in
this
> precip
> > data? I'm going to assume that each timestep contains data 1 days
of
> > accumulated precip. You can use the "-add" option in pcp_combine,
but
> > adding together 30 timesteps will be a very lengthy command line.
Here's
> > how you'd do 5 of them:
> >
> > *pcp_combine -add \*
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010301_103000,*,*)";' \*
> >
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010302_103000,*,*)";' \*
> >
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010303_103000,*,*)";' \*
> >
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010304_103000,*,*)";' \*
> >
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010305_103000,*,*)";' \*
> > *precip_20010301_20010305.nc <http://precip_20010301_20010305.nc>*
> >
> > Honestly, pcp_combine is probably not your best bet for summing up
1
> month
> > of daily precip amounts. You could probably do this more
efficiently
> using
> > one of the CDO tools.
> >
> > But basically, when running the "-add", "-subtract", or "-derive"
> options,
> > you list the file name followed by a string defining how to
extract data
> > from that file. I can see how things would be much easier if you
could
> > define a range of dates like the following. But that is not
currently
> > supported:
> > *pcp_combine -add \*
> > *daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc 'name="TP";
> > level="(20010301_103000-20010331_103000,*,*)";' precip_march.nc
> > <http://precip_march.nc>*
> >
> > Hope this helps get you going.
> >
> > Thanks,
> > John
> >
> > On Wed, Apr 1, 2020 at 1:42 PM gbode Imoleayo via RT
<met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94659 >
> > >
> > > Dear John,
> > > This is the error I get.
> > > MET_Tutorial/wrfMOD> pcp_combine -sum 20010301_103000 1
20010331_103000
> > 30
> > > daily_TotalPrecip_mar_dec_2001_GD_MYJ_BMJ.nc -pcpdir
pcp_combine_out/
> > > DEBUG 2: Performing sum command: Init/In_Accum/Valid/Out_Accum
Times =
> > > 20010301_103000/010000/20010331_103000/300000
> > > DEBUG 2: Searching for 30 files with accumulation times of
010000 to
> sum
> > > to a total accumulation time of 300000.
> > > ERROR :
> > > ERROR : sum_data_files() -> cannot find a file with a valid
time of
> > > 20010331_103000 and accumulation time of 010000 matching the
regular
> > > expression ".*"
> > > ERROR :
> > >
> > >
> > > Thank you.
> > > Best regards,Ezekiel
> > >
> > >
> > > On Wednesday, April 1, 2020, 7:26:01 PM GMT+2, John Halley
Gotway
> via
> > > RT <met_help at ucar.edu> wrote:
> > >
> > > Hello,
> > >
> > > I just ran across this MET-Help ticket and realize that we never
sent
> you
> > > a response. I apologize for letting this slip through the
cracks. Are
> you
> > > still having trouble running pcp_combine on cheyenne?
> > >
> > > Please let me know what questions you have.
> > >
> > > Thanks,
> > > John Halley Gotway
> > >
> > >
> >
> >
> >
>
>
>
------------------------------------------------
More information about the Met_help
mailing list