[Met_help] [rt.rap.ucar.edu #94991] History for MODE Analysis [NetCDF: HDF error w/StageIV + WRF Output]
John Halley Gotway via RT
met_help at ucar.edu
Thu Apr 23 14:01:37 MDT 2020
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello,
I'm trying to run ''mode" analysis for WRFv3.7/8.1 output (RAINNC+RAINC)
with respect to summed, 24hr StageIV precipitation data (over 72hrs).
Unfortunately, I encountered the following error while running
./bin/mode... According to my mode.log file, I need 3 arguments for my
forecasted dataset (I only provided 3 inputs while running mode).
[image: image.png]
Can you provide some feedback to help me resolve this problem?
Here is my procedure:
1) Get WRF total accumulation .nc file
AFWA_TOTPRECIP: "AFWA Diagnostic: Total simulation precip"
1.a) Load NCO
module load nco/4.4.0
1.b) Process data for final wrfout hour
ncks -v AFWA_TOTPRECIP wrfout_d01_2006-02-13_06:00:00 AFWA_TOTPRECIP.nc
2) Sum StageIV precipitation Analysis
2.a) Extract tar/.Z files
tar -xvf [file_name]
gunzip *.Z
2.b) Sum three, 24hrly StageIV precipitation files.
At the met-8.1.2 directory:
./bin/pcp_combine -sum 00000000_000000 24 20060213_120000 72 sum_test.nc
-pcpdir ./MODE_PROCEDURE/STAGEIV_Test/
3) Run 'mode'
At the met-8.1.2 directory:
./bin/mode ./MODE_PROCEDURE/MODE_Files_nc/AFWA_TOTPRECIP.nc
./MODE_PROCEDURE/MODE_Files_nc/sum_test.nc
./MODE_PROCEDURE/MODE_Files_nc/MODEConfig_APCP_WWE -log mode.log
I believe by adding "to_grid = FCST;" into the mode configuration file,
regridding sum_test.nc to the grid of AFWA_TOTPRECIP.nc is unnecessary.
Any help will greatly be appreciated!
I attached my configuration file, AFWA_TOTPRECIP.nc, sum_test.nc, and
mode.log if it helps!
--
*Sincerely,*
Michael S. Walters | B.S. Atmospheric Science
Graduate Research Assistant
Master's Student | Civil and Environmental Engineering
261 Glenbrook Road. Storrs, CT, 06269-3037
University of Connecticut
Phone: 774-254-6608
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: MODE Analysis [NetCDF: HDF error w/StageIV + WRF Output]
From: Julie Prestopnik
Time: Tue Apr 21 18:33:00 2020
Hi Michael.
I have assigned this ticket to John Halley Gotway, who should be able
to
assist you with this work. Please allow a few business days for a
response.
Julie
--
Julie Prestopnik
Software Engineer
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
My working day may not be your working day. Please do not feel
obliged to
reply to this email outside of your normal working hours.
------------------------------------------------
Subject: MODE Analysis [NetCDF: HDF error w/StageIV + WRF Output]
From: John Halley Gotway
Time: Wed Apr 22 23:16:17 2020
Hi Michael,
Sorry for the delay in getting back to you on this. I see you're
having
some trouble running the MODE tool. Thanks for sending your sample
data
files.
I ran mode using the files you sent was able to replicate the behavior
you
describe:
> mode AFWA_TOTPRECIP.nc sum_test.nc MODEConfig_APCP_WWE
ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane &)
const
-> needed 3 arguments for variable AFWA_TOTPRECIP, got 2
Take a look in the AFWA file at the AFWA_TOTPRECIP variable:
> ncdump -h AFWA_TOTPRECIP.nc
float AFWA_TOTPRECIP(XTIME, y, x) ;
Notice that it's indexed by 3 dimensions. In the MODE config file,
the
"level" setting must include an entry for each of those 3 dimensions.
I
updated MODEConfig_APCP_WWE as shown below:
*fcst = { field = { name = "AFWA_TOTPRECIP"; level =
"(0,*,*)"; }*
I also had to revert the colortable settings back to:
* color_table = "MET_BASE/colortables/met_default.ctable";
color_table = "MET_BASE/colortables/met_default.ctable";
color_table = "MET_BASE/colortables/mode_obj.ctable";*
And then it ran fine. See a png version of the first page of the
PostScript output from MODE.
Judging by the size of the resolved objects, I'd recommend reducing
the
convolution radius or increasing the threshold to make smaller
objects.
Hopefully this helps get you going.
Thanks,
John Halley Gotway
On Tue, Apr 21, 2020 at 6:33 PM Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991 >
>
> Hi Michael.
>
> I have assigned this ticket to John Halley Gotway, who should be
able to
> assist you with this work. Please allow a few business days for a
> response.
>
> Julie
>
> --
> Julie Prestopnik
> Software Engineer
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
> My working day may not be your working day. Please do not feel
obliged to
> reply to this email outside of your normal working hours.
>
>
------------------------------------------------
Subject: MODE Analysis [NetCDF: HDF error w/StageIV + WRF Output]
From: Michael Walters
Time: Thu Apr 23 11:03:23 2020
Thank you so very much, John! This tool is great! I have a few other
questions that I'd like your feedback for though.
My domain setup for WRF is about the size of the plotted box. Would
you
recommend masking StageIV (observation) precipitation data to the
bounds of
the WRF (forecast) domain for better results?
Also, through your experience, would you recommend changing any of the
mode
options (such as conv. radius)? A little bit of background, I ran 38
simulations (for 3 different modeling systems) for significant snow
events
based on the Regional Snowfall Index (for the Northeast U.S.) and I'm
comparing NWP qpf fields to StageIV precipitation analysis. My goal is
to
determine how well NWP performs with respect to broad areas of
precipitation (possibly linking it to banded precipitation as well)
during
snow events.
Thanks again!
<http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail>
Virus-free.
www.avg.com
<http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
On Thu, Apr 23, 2020 at 1:16 AM John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
> *Message sent from a system outside of UConn.*
>
>
> Hi Michael,
>
> Sorry for the delay in getting back to you on this. I see you're
having
> some trouble running the MODE tool. Thanks for sending your sample
data
> files.
>
> I ran mode using the files you sent was able to replicate the
behavior you
> describe:
> > mode AFWA_TOTPRECIP.nc sum_test.nc MODEConfig_APCP_WWE
> ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane &)
const
> -> needed 3 arguments for variable AFWA_TOTPRECIP, got 2
>
> Take a look in the AFWA file at the AFWA_TOTPRECIP variable:
> > ncdump -h AFWA_TOTPRECIP.nc
> float AFWA_TOTPRECIP(XTIME, y, x) ;
>
> Notice that it's indexed by 3 dimensions. In the MODE config file,
the
> "level" setting must include an entry for each of those 3
dimensions. I
> updated MODEConfig_APCP_WWE as shown below:
>
>
>
>
>
> *fcst = { field = { name = "AFWA_TOTPRECIP"; level =
> "(0,*,*)"; }*
>
> I also had to revert the colortable settings back to:
>
>
>
> * color_table = "MET_BASE/colortables/met_default.ctable";
> color_table = "MET_BASE/colortables/met_default.ctable";
> color_table = "MET_BASE/colortables/mode_obj.ctable";*
>
> And then it ran fine. See a png version of the first page of the
> PostScript output from MODE.
>
> Judging by the size of the resolved objects, I'd recommend reducing
the
> convolution radius or increasing the threshold to make smaller
objects.
>
> Hopefully this helps get you going.
>
> Thanks,
> John Halley Gotway
>
> On Tue, Apr 21, 2020 at 6:33 PM Julie Prestopnik via RT
<met_help at ucar.edu
> >
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991 >
> >
> > Hi Michael.
> >
> > I have assigned this ticket to John Halley Gotway, who should be
able to
> > assist you with this work. Please allow a few business days for a
> > response.
> >
> > Julie
> >
> > --
> > Julie Prestopnik
> > Software Engineer
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> > My working day may not be your working day. Please do not feel
obliged
> to
> > reply to this email outside of your normal working hours.
> >
> >
>
>
--
*Sincerely,*
Michael S. Walters | B.S. Atmospheric Science
Graduate Research Assistant
Master's Student | Civil and Environmental Engineering
261 Glenbrook Road. Storrs, CT, 06269-3037
University of Connecticut
Phone: 774-254-6608
<http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail>
Virus-free.
www.avg.com
<http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail>
<#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
------------------------------------------------
Subject: MODE Analysis [NetCDF: HDF error w/StageIV + WRF Output]
From: John Halley Gotway
Time: Thu Apr 23 12:33:56 2020
Michael,
Since you're already using the "regrid" configuration option to
automatically regrid the StageIV data to your WRF model domain, I
can't see
how additional masking of the StageIV data would be useful. I don't
think
it'd have any impact on the results.
As for the settings, you can modify the convolution radius and
threshold to
create objects which match the features you're trying to study. So it
really is up to you. I will say though that generally I find MODE's
results easiest to interpret when the size of the objects are
relatively
small compared to the size of the domain. In your case, with such a
small
domain in the northeast and with relatively large scale precipitation
objects, you're be left with lot's of edge effects.
If you were to run MODE on a global grid, then the MODE's object area
would
accurately represent the area of the precipitation feature. In your
case,
the precip feature will very often be truncated by the edge of the
domain.
And that makes the interpretation more difficult. When looking at
matched
fcst/obs pairs, you may see a difference in the centroid locations.
But is
that centroid offset due to a true offset in the precipitation object
or
due to what part of that object happens to fall in your domain?
Whenever
an object touches the edge of the domain, we can't really answer that
question.
Thanks,
John
On Thu, Apr 23, 2020 at 11:06 AM Michael Walters via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991 >
>
> Thank you so very much, John! This tool is great! I have a few other
> questions that I'd like your feedback for though.
>
> My domain setup for WRF is about the size of the plotted box. Would
you
> recommend masking StageIV (observation) precipitation data to the
bounds of
> the WRF (forecast) domain for better results?
>
> Also, through your experience, would you recommend changing any of
the mode
> options (such as conv. radius)? A little bit of background, I ran 38
> simulations (for 3 different modeling systems) for significant snow
events
> based on the Regional Snowfall Index (for the Northeast U.S.) and
I'm
> comparing NWP qpf fields to StageIV precipitation analysis. My goal
is to
> determine how well NWP performs with respect to broad areas of
> precipitation (possibly linking it to banded precipitation as well)
during
> snow events.
>
> Thanks again!
>
>
> <
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> >
> Virus-free.
> www.avg.com
> <
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> >
> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
> On Thu, Apr 23, 2020 at 1:16 AM John Halley Gotway via RT <
> met_help at ucar.edu>
> wrote:
>
> > *Message sent from a system outside of UConn.*
> >
> >
> > Hi Michael,
> >
> > Sorry for the delay in getting back to you on this. I see you're
having
> > some trouble running the MODE tool. Thanks for sending your
sample data
> > files.
> >
> > I ran mode using the files you sent was able to replicate the
behavior
> you
> > describe:
> > > mode AFWA_TOTPRECIP.nc sum_test.nc MODEConfig_APCP_WWE
> > ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane
&) const
> > -> needed 3 arguments for variable AFWA_TOTPRECIP, got 2
> >
> > Take a look in the AFWA file at the AFWA_TOTPRECIP variable:
> > > ncdump -h AFWA_TOTPRECIP.nc
> > float AFWA_TOTPRECIP(XTIME, y, x) ;
> >
> > Notice that it's indexed by 3 dimensions. In the MODE config
file, the
> > "level" setting must include an entry for each of those 3
dimensions. I
> > updated MODEConfig_APCP_WWE as shown below:
> >
> >
> >
> >
> >
> > *fcst = { field = { name = "AFWA_TOTPRECIP"; level =
> > "(0,*,*)"; }*
> >
> > I also had to revert the colortable settings back to:
> >
> >
> >
> > * color_table = "MET_BASE/colortables/met_default.ctable";
> > color_table = "MET_BASE/colortables/met_default.ctable";
> > color_table = "MET_BASE/colortables/mode_obj.ctable";*
> >
> > And then it ran fine. See a png version of the first page of the
> > PostScript output from MODE.
> >
> > Judging by the size of the resolved objects, I'd recommend
reducing the
> > convolution radius or increasing the threshold to make smaller
objects.
> >
> > Hopefully this helps get you going.
> >
> > Thanks,
> > John Halley Gotway
> >
> > On Tue, Apr 21, 2020 at 6:33 PM Julie Prestopnik via RT <
> met_help at ucar.edu
> > >
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991 >
> > >
> > > Hi Michael.
> > >
> > > I have assigned this ticket to John Halley Gotway, who should be
able
> to
> > > assist you with this work. Please allow a few business days for
a
> > > response.
> > >
> > > Julie
> > >
> > > --
> > > Julie Prestopnik
> > > Software Engineer
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > > My working day may not be your working day. Please do not feel
obliged
> > to
> > > reply to this email outside of your normal working hours.
> > >
> > >
> >
> >
>
> --
>
> *Sincerely,*
>
> Michael S. Walters | B.S. Atmospheric Science
> Graduate Research Assistant
> Master's Student | Civil and Environmental Engineering
> 261 Glenbrook Road. Storrs, CT, 06269-3037
> University of Connecticut
> Phone: 774-254-6608
>
> <
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> >
> Virus-free.
> www.avg.com
> <
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> >
> <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
>
------------------------------------------------
Subject: MODE Analysis [NetCDF: HDF error w/StageIV + WRF Output]
From: Michael Walters
Time: Thu Apr 23 13:45:07 2020
Thanks for your feedback and help, John.
That's a good point that I did not consider (poor centroids defined as
a
result of NWP's boundary). Given what you wrote, if this is applicable
to
what I'd like to study, I'll configure MODE specifically to what I
want to
analyze. I'll try to define any boundary related errors if they arise.
Have a great day!
<http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail>
Virus-free.
www.avg.com
<http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail>
<#m_-6519207214278552914_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
On Thu, Apr 23, 2020 at 2:34 PM John Halley Gotway via RT
<met_help at ucar.edu>
wrote:
> *Message sent from a system outside of UConn.*
>
>
> Michael,
>
> Since you're already using the "regrid" configuration option to
> automatically regrid the StageIV data to your WRF model domain, I
can't see
> how additional masking of the StageIV data would be useful. I don't
think
> it'd have any impact on the results.
>
> As for the settings, you can modify the convolution radius and
threshold to
> create objects which match the features you're trying to study. So
it
> really is up to you. I will say though that generally I find MODE's
> results easiest to interpret when the size of the objects are
relatively
> small compared to the size of the domain. In your case, with such a
small
> domain in the northeast and with relatively large scale
precipitation
> objects, you're be left with lot's of edge effects.
>
> If you were to run MODE on a global grid, then the MODE's object
area would
> accurately represent the area of the precipitation feature. In your
case,
> the precip feature will very often be truncated by the edge of the
domain.
> And that makes the interpretation more difficult. When looking at
matched
> fcst/obs pairs, you may see a difference in the centroid locations.
But is
> that centroid offset due to a true offset in the precipitation
object or
> due to what part of that object happens to fall in your domain?
Whenever
> an object touches the edge of the domain, we can't really answer
that
> question.
>
> Thanks,
> John
>
> On Thu, Apr 23, 2020 at 11:06 AM Michael Walters via RT
<met_help at ucar.edu
> >
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991 >
> >
> > Thank you so very much, John! This tool is great! I have a few
other
> > questions that I'd like your feedback for though.
> >
> > My domain setup for WRF is about the size of the plotted box.
Would you
> > recommend masking StageIV (observation) precipitation data to the
bounds
> of
> > the WRF (forecast) domain for better results?
> >
> > Also, through your experience, would you recommend changing any of
the
> mode
> > options (such as conv. radius)? A little bit of background, I ran
38
> > simulations (for 3 different modeling systems) for significant
snow
> events
> > based on the Regional Snowfall Index (for the Northeast U.S.) and
I'm
> > comparing NWP qpf fields to StageIV precipitation analysis. My
goal is to
> > determine how well NWP performs with respect to broad areas of
> > precipitation (possibly linking it to banded precipitation as
well)
> during
> > snow events.
> >
> > Thanks again!
> >
> >
> > <
> >
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> > >
> > Virus-free.
> > www.avg.com
> > <
> >
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> > >
> > <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
> >
> > On Thu, Apr 23, 2020 at 1:16 AM John Halley Gotway via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > *Message sent from a system outside of UConn.*
> > >
> > >
> > > Hi Michael,
> > >
> > > Sorry for the delay in getting back to you on this. I see
you're
> having
> > > some trouble running the MODE tool. Thanks for sending your
sample
> data
> > > files.
> > >
> > > I ran mode using the files you sent was able to replicate the
behavior
> > you
> > > describe:
> > > > mode AFWA_TOTPRECIP.nc sum_test.nc MODEConfig_APCP_WWE
> > > ERROR : NcCfFile::getData(NcVar *, const LongArray &, DataPlane
&)
> const
> > > -> needed 3 arguments for variable AFWA_TOTPRECIP, got 2
> > >
> > > Take a look in the AFWA file at the AFWA_TOTPRECIP variable:
> > > > ncdump -h AFWA_TOTPRECIP.nc
> > > float AFWA_TOTPRECIP(XTIME, y, x) ;
> > >
> > > Notice that it's indexed by 3 dimensions. In the MODE config
file, the
> > > "level" setting must include an entry for each of those 3
dimensions.
> I
> > > updated MODEConfig_APCP_WWE as shown below:
> > >
> > >
> > >
> > >
> > >
> > > *fcst = { field = { name = "AFWA_TOTPRECIP"; level
=
> > > "(0,*,*)"; }*
> > >
> > > I also had to revert the colortable settings back to:
> > >
> > >
> > >
> > > * color_table =
"MET_BASE/colortables/met_default.ctable";
> > > color_table = "MET_BASE/colortables/met_default.ctable";
> > > color_table = "MET_BASE/colortables/mode_obj.ctable";*
> > >
> > > And then it ran fine. See a png version of the first page of
the
> > > PostScript output from MODE.
> > >
> > > Judging by the size of the resolved objects, I'd recommend
reducing the
> > > convolution radius or increasing the threshold to make smaller
objects.
> > >
> > > Hopefully this helps get you going.
> > >
> > > Thanks,
> > > John Halley Gotway
> > >
> > > On Tue, Apr 21, 2020 at 6:33 PM Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > >
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991
>
> > > >
> > > > Hi Michael.
> > > >
> > > > I have assigned this ticket to John Halley Gotway, who should
be able
> > to
> > > > assist you with this work. Please allow a few business days
for a
> > > > response.
> > > >
> > > > Julie
> > > >
> > > > --
> > > > Julie Prestopnik
> > > > Software Engineer
> > > > National Center for Atmospheric Research
> > > > Research Applications Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > > > My working day may not be your working day. Please do not
feel
> obliged
> > > to
> > > > reply to this email outside of your normal working hours.
> > > >
> > > >
> > >
> > >
> >
> > --
> >
> > *Sincerely,*
> >
> > Michael S. Walters | B.S. Atmospheric Science
> > Graduate Research Assistant
> > Master's Student | Civil and Environmental Engineering
> > 261 Glenbrook Road. Storrs, CT, 06269-3037
> > University of Connecticut
> > Phone: 774-254-6608
> >
> > <
> >
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> > >
> > Virus-free.
> > www.avg.com
> > <
> >
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> > >
> > <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
> >
> >
>
>
--
*Sincerely,*
Michael S. Walters | B.S. Atmospheric Science
Graduate Research Assistant
Master's Student | Civil and Environmental Engineering
261 Glenbrook Road. Storrs, CT, 06269-3037
University of Connecticut
Phone: 774-254-6608
<http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail>
Virus-free.
www.avg.com
<http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail>
<#m_-6519207214278552914_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
------------------------------------------------
Subject: MODE Analysis [NetCDF: HDF error w/StageIV + WRF Output]
From: John Halley Gotway
Time: Thu Apr 23 14:01:35 2020
Sure thing. Hope you find MODE useful in your work.
I'll go ahead and resolve this ticket.
Thanks,
John
On Thu, Apr 23, 2020 at 1:45 PM Michael Walters via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991 >
>
> Thanks for your feedback and help, John.
>
> That's a good point that I did not consider (poor centroids defined
as a
> result of NWP's boundary). Given what you wrote, if this is
applicable to
> what I'd like to study, I'll configure MODE specifically to what I
want to
> analyze. I'll try to define any boundary related errors if they
arise.
>
> Have a great day!
>
>
>
> <
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> >
> Virus-free.
> www.avg.com
> <
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> >
> <#m_-6519207214278552914_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
> On Thu, Apr 23, 2020 at 2:34 PM John Halley Gotway via RT <
> met_help at ucar.edu>
> wrote:
>
> > *Message sent from a system outside of UConn.*
> >
> >
> > Michael,
> >
> > Since you're already using the "regrid" configuration option to
> > automatically regrid the StageIV data to your WRF model domain, I
can't
> see
> > how additional masking of the StageIV data would be useful. I
don't think
> > it'd have any impact on the results.
> >
> > As for the settings, you can modify the convolution radius and
threshold
> to
> > create objects which match the features you're trying to study.
So it
> > really is up to you. I will say though that generally I find
MODE's
> > results easiest to interpret when the size of the objects are
relatively
> > small compared to the size of the domain. In your case, with such
a
> small
> > domain in the northeast and with relatively large scale
precipitation
> > objects, you're be left with lot's of edge effects.
> >
> > If you were to run MODE on a global grid, then the MODE's object
area
> would
> > accurately represent the area of the precipitation feature. In
your
> case,
> > the precip feature will very often be truncated by the edge of the
> domain.
> > And that makes the interpretation more difficult. When looking at
> matched
> > fcst/obs pairs, you may see a difference in the centroid
locations. But
> is
> > that centroid offset due to a true offset in the precipitation
object or
> > due to what part of that object happens to fall in your domain?
Whenever
> > an object touches the edge of the domain, we can't really answer
that
> > question.
> >
> > Thanks,
> > John
> >
> > On Thu, Apr 23, 2020 at 11:06 AM Michael Walters via RT <
> met_help at ucar.edu
> > >
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991 >
> > >
> > > Thank you so very much, John! This tool is great! I have a few
other
> > > questions that I'd like your feedback for though.
> > >
> > > My domain setup for WRF is about the size of the plotted box.
Would you
> > > recommend masking StageIV (observation) precipitation data to
the
> bounds
> > of
> > > the WRF (forecast) domain for better results?
> > >
> > > Also, through your experience, would you recommend changing any
of the
> > mode
> > > options (such as conv. radius)? A little bit of background, I
ran 38
> > > simulations (for 3 different modeling systems) for significant
snow
> > events
> > > based on the Regional Snowfall Index (for the Northeast U.S.)
and I'm
> > > comparing NWP qpf fields to StageIV precipitation analysis. My
goal is
> to
> > > determine how well NWP performs with respect to broad areas of
> > > precipitation (possibly linking it to banded precipitation as
well)
> > during
> > > snow events.
> > >
> > > Thanks again!
> > >
> > >
> > > <
> > >
> >
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> > > >
> > > Virus-free.
> > > www.avg.com
> > > <
> > >
> >
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> > > >
> > > <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
> > >
> > > On Thu, Apr 23, 2020 at 1:16 AM John Halley Gotway via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > > *Message sent from a system outside of UConn.*
> > > >
> > > >
> > > > Hi Michael,
> > > >
> > > > Sorry for the delay in getting back to you on this. I see
you're
> > having
> > > > some trouble running the MODE tool. Thanks for sending your
sample
> > data
> > > > files.
> > > >
> > > > I ran mode using the files you sent was able to replicate the
> behavior
> > > you
> > > > describe:
> > > > > mode AFWA_TOTPRECIP.nc sum_test.nc MODEConfig_APCP_WWE
> > > > ERROR : NcCfFile::getData(NcVar *, const LongArray &,
DataPlane &)
> > const
> > > > -> needed 3 arguments for variable AFWA_TOTPRECIP, got 2
> > > >
> > > > Take a look in the AFWA file at the AFWA_TOTPRECIP variable:
> > > > > ncdump -h AFWA_TOTPRECIP.nc
> > > > float AFWA_TOTPRECIP(XTIME, y, x) ;
> > > >
> > > > Notice that it's indexed by 3 dimensions. In the MODE config
file,
> the
> > > > "level" setting must include an entry for each of those 3
dimensions.
> > I
> > > > updated MODEConfig_APCP_WWE as shown below:
> > > >
> > > >
> > > >
> > > >
> > > >
> > > > *fcst = { field = { name = "AFWA_TOTPRECIP";
level =
> > > > "(0,*,*)"; }*
> > > >
> > > > I also had to revert the colortable settings back to:
> > > >
> > > >
> > > >
> > > > * color_table =
"MET_BASE/colortables/met_default.ctable";
> > > > color_table = "MET_BASE/colortables/met_default.ctable";
> > > > color_table = "MET_BASE/colortables/mode_obj.ctable";*
> > > >
> > > > And then it ran fine. See a png version of the first page of
the
> > > > PostScript output from MODE.
> > > >
> > > > Judging by the size of the resolved objects, I'd recommend
reducing
> the
> > > > convolution radius or increasing the threshold to make smaller
> objects.
> > > >
> > > > Hopefully this helps get you going.
> > > >
> > > > Thanks,
> > > > John Halley Gotway
> > > >
> > > > On Tue, Apr 21, 2020 at 6:33 PM Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > >
> > > > wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=94991 >
> > > > >
> > > > > Hi Michael.
> > > > >
> > > > > I have assigned this ticket to John Halley Gotway, who
should be
> able
> > > to
> > > > > assist you with this work. Please allow a few business days
for a
> > > > > response.
> > > > >
> > > > > Julie
> > > > >
> > > > > --
> > > > > Julie Prestopnik
> > > > > Software Engineer
> > > > > National Center for Atmospheric Research
> > > > > Research Applications Laboratory
> > > > > Phone: 303.497.8399
> > > > > Email: jpresto at ucar.edu
> > > > >
> > > > > My working day may not be your working day. Please do not
feel
> > obliged
> > > > to
> > > > > reply to this email outside of your normal working hours.
> > > > >
> > > > >
> > > >
> > > >
> > >
> > > --
> > >
> > > *Sincerely,*
> > >
> > > Michael S. Walters | B.S. Atmospheric Science
> > > Graduate Research Assistant
> > > Master's Student | Civil and Environmental Engineering
> > > 261 Glenbrook Road. Storrs, CT, 06269-3037
> > > University of Connecticut
> > > Phone: 774-254-6608
> > >
> > > <
> > >
> >
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> > > >
> > > Virus-free.
> > > www.avg.com
> > > <
> > >
> >
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> > > >
> > > <#DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
> > >
> > >
> >
> >
>
> --
>
> *Sincerely,*
>
> Michael S. Walters | B.S. Atmospheric Science
> Graduate Research Assistant
> Master's Student | Civil and Environmental Engineering
> 261 Glenbrook Road. Storrs, CT, 06269-3037
> University of Connecticut
> Phone: 774-254-6608
>
> <
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> >
> Virus-free.
> www.avg.com
> <
> http://www.avg.com/email-
signature?utm_medium=email&utm_source=link&utm_campaign=sig-
email&utm_content=webmail
> >
> <#m_-6519207214278552914_DAB4FAD8-2DD7-40BB-A1B8-4E2AA1F9FDF2>
>
>
------------------------------------------------
More information about the Met_help
mailing list