[Met_help] [rt.rap.ucar.edu #74467] History for point_stat for WRF output question
Julie Prestopnik via RT
met_help at ucar.edu
Wed Feb 10 09:36:12 MST 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Hello there,
I got point_stat run through for the tutorial case, and would like to try
on the different domain like Nepal, where you have kindly helped creating
the mask file. I edit the config file as attached. The error shows up as
follow:
~/met-5.1_bugfix$ ./bin/point_stat \
>
/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-20_18:00:00
\
> tutorial/out/ascii2nc/nepal.nc \
> tutorial/config/PointStatConfig_Nepal \
> -outdir tutorial/out/point_stat \
> -v 2
DEBUG 1: Default Config File:
/home/vsurapipith/met-5.1_bugfix/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
ERROR :
ERROR : Trouble reading forecast file
"/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-20_18:00:00"
ERROR :
There seems to be some error about reading the WRF output file. Could you
give me some clue what kind of error that might be, such that I may
correct, please. Is there something to set in the config file? I am not
sure what could be wrong with the WRF output, which has worked okay with
NCL/other tools.
Best wishes,
Vanisa
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: point_stat for WRF output question
From: Julie Prestopnik
Time: Tue Dec 22 11:10:06 2015
Hi Vanisa.
It is difficult to say without having access to the file you are
using.
Could you please upload your input forecast and observation file to
our ftp
site using the directions on this page:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
In an effort to try to help before having access to the files, I will
say
that we have seen this problem before when the first four characters
in the
file are not "GRIB". That is called the "magic cookie" to indicate
the
file type.
You could use the command "od -c wrfout_d02_2013-06-20_18:00:00 |
more" to
determine what the first four characters are. You should see
something
like this:
0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0 \0 001 P
´
If the first four characters are not "GRIB", you could get around this
by
changing the file's suffix to be ".grb" or by using the "file_type =
GRIB1;" option in the config file (see data/config/README for more
information).
I will be out of the office starting today at 2:45pm (US Mountain
time)
through January 4th. I hope I can help you with this problem before I
am
out of the office.
Julie
On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via RT
<met_help at ucar.edu
> wrote:
>
> Tue Dec 22 03:18:04 2015: Request 74467 was acted upon.
> Transaction: Ticket created by vanisas at gmail.com
> Queue: met_help
> Subject: point_stat for WRF output question
> Owner: Nobody
> Requestors: vanisas at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>
>
> Hello there,
>
> I got point_stat run through for the tutorial case, and would like
to try
> on the different domain like Nepal, where you have kindly helped
creating
> the mask file. I edit the config file as attached. The error shows
up as
> follow:
>
> ~/met-5.1_bugfix$ ./bin/point_stat \
> >
>
> /home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00
> \
> > tutorial/out/ascii2nc/nepal.nc \
> > tutorial/config/PointStatConfig_Nepal \
> > -outdir tutorial/out/point_stat \
> > -v 2
> DEBUG 1: Default Config File:
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> ERROR :
> ERROR : Trouble reading forecast file
>
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> ERROR :
>
> There seems to be some error about reading the WRF output file.
Could you
> give me some clue what kind of error that might be, such that I may
> correct, please. Is there something to set in the config file? I am
not
> sure what could be wrong with the WRF output, which has worked okay
with
> NCL/other tools.
>
> Best wishes,
> Vanisa
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: RE: [rt.rap.ucar.edu #74467] point_stat for WRF output question
From: Vanisa Surapipith
Time: Wed Dec 23 12:52:36 2015
Hi Julie,
Thank you very much for your kind guidance. I am sorry that I just
read this e-mail, since there was no internet connection along my
travel way, whole day, towards home in Thailand. I really appreciate
your effort to help me before the holidays begin; I'm only now at home
in Thailand with internet, but my Nepal workstation with WRF and MET
seems to down now, after midnight of December 23rd in Kathmandu. I can
ask the IT admin to check what happen, but he would only be able to
resolve this in the next 8 hours (Christmas is actually not a holiday
in Nepal, but I'm taking end-of-year leave rather).
Never mind; let's celebrate the Season & give the best to people at
home; I'll find out as you recommend in this while, and report soon,
but I'll be okay to find you response later. Wish there a warm
Christmas Eve & Happy coming New Year.
Best wishes,
Vanisa
-----Original Message-----
From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
Sent: 22 December 2015 23:55
To: vanisas at gmail.com
Cc: Vanisa Surapipith
Subject: Re: [rt.rap.ucar.edu #74467] point_stat for WRF output
question
Hi Vanisa.
It is difficult to say without having access to the file you are
using.
Could you please upload your input forecast and observation file to
our ftp site using the directions on this page:
http://www.dtcenter.org/met/users/support/met_help.php#ftp
In an effort to try to help before having access to the files, I will
say that we have seen this problem before when the first four
characters in the file are not "GRIB". That is called the "magic
cookie" to indicate the file type.
You could use the command "od -c wrfout_d02_2013-06-20_18:00:00 |
more" to determine what the first four characters are. You should see
something like this:
0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0 \0 001 P
´
If the first four characters are not "GRIB", you could get around this
by changing the file's suffix to be ".grb" or by using the "file_type
= GRIB1;" option in the config file (see data/config/README for more
information).
I will be out of the office starting today at 2:45pm (US Mountain
time) through January 4th. I hope I can help you with this problem
before I am out of the office.
Julie
On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via RT
<met_help at ucar.edu
> wrote:
>
> Tue Dec 22 03:18:04 2015: Request 74467 was acted upon.
> Transaction: Ticket created by vanisas at gmail.com
> Queue: met_help
> Subject: point_stat for WRF output question
> Owner: Nobody
> Requestors: vanisas at gmail.com
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> >
>
>
> Hello there,
>
> I got point_stat run through for the tutorial case, and would like
to
> try on the different domain like Nepal, where you have kindly helped
> creating the mask file. I edit the config file as attached. The
error
> shows up as
> follow:
>
> ~/met-5.1_bugfix$ ./bin/point_stat \
> >
>
> /home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-
> 20_18:00:00
> \
> > tutorial/out/ascii2nc/nepal.nc \
> > tutorial/config/PointStatConfig_Nepal \ -outdir
> > tutorial/out/point_stat \ -v 2
> DEBUG 1: Default Config File:
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> ult DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> ERROR :
> ERROR : Trouble reading forecast file
>
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> ERROR :
>
> There seems to be some error about reading the WRF output file.
Could
> you give me some clue what kind of error that might be, such that I
> may correct, please. Is there something to set in the config file? I
> am not sure what could be wrong with the WRF output, which has
worked
> okay with NCL/other tools.
>
> Best wishes,
> Vanisa
>
>
--
Julie Prestopnik
National Center for Atmospheric Research Research Applications
Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Sun Dec 27 23:37:17 2015
Dear Julie & all there,
I just manage to try as Julie recommended right before the Christmas
holidays, i.e. to put the input file I use on the ftp site. I am not
sure
if it works, or not. What shows up say:
220 brew.rap.ucar.edu FTP server ready.
Name (ftp.rap.ucar.edu:vsurapipith): anonymous
331 Guest login ok, send your complete e-mail address as password.
Password:
230 Guest login ok, access restrictions apply.
Remote system type is UNIX.
Using binary mode to transfer files.
ftp> cd incoming/irap/met_help
250 CWD command successful.
ftp> mkdir surapipith_data
257 "/incoming/irap/met_help/surapipith_data" new directory created.
ftp> cd surapipith_data
250 CWD command successful.
ftp> put wrfout_d02_2013-06-20_18:00:00
local: wrfout_d02_2013-06-20_18:00:00 remote: wrfout_d02_2013-06-
20_18:00:00
No control connection for command: Success
Passive mode refused.
ftp> bye
Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00 | more", there
were
many lines running up, and I could not find the first 4 line if it has
GRIB
word or not. I will find the way to try again. Yet I hope you may see
the
WRF output that I upload to the ftp site of ucar. If not, kindly let
me
know do it right, please.
Thanks very much.
Best wishes,
Vanisa
On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via RT
<met_help at ucar.edu
> wrote:
> Hi Julie,
>
> Thank you very much for your kind guidance. I am sorry that I just
read
> this e-mail, since there was no internet connection along my travel
way,
> whole day, towards home in Thailand. I really appreciate your effort
to
> help me before the holidays begin; I'm only now at home in Thailand
with
> internet, but my Nepal workstation with WRF and MET seems to down
now,
> after midnight of December 23rd in Kathmandu. I can ask the IT admin
to
> check what happen, but he would only be able to resolve this in the
next 8
> hours (Christmas is actually not a holiday in Nepal, but I'm taking
> end-of-year leave rather).
>
> Never mind; let's celebrate the Season & give the best to people at
home;
> I'll find out as you recommend in this while, and report soon, but
I'll be
> okay to find you response later. Wish there a warm Christmas Eve &
Happy
> coming New Year.
>
> Best wishes,
> Vanisa
>
> -----Original Message-----
> From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
> Sent: 22 December 2015 23:55
> To: vanisas at gmail.com
> Cc: Vanisa Surapipith
> Subject: Re: [rt.rap.ucar.edu #74467] point_stat for WRF output
question
>
> Hi Vanisa.
>
> It is difficult to say without having access to the file you are
using.
> Could you please upload your input forecast and observation file to
our
> ftp site using the directions on this page:
>
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> In an effort to try to help before having access to the files, I
will say
> that we have seen this problem before when the first four characters
in the
> file are not "GRIB". That is called the "magic cookie" to indicate
the
> file type.
>
> You could use the command "od -c wrfout_d02_2013-06-20_18:00:00 |
more"
> to determine what the first four characters are. You should see
something
> like this:
>
> 0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0 \0 001 P
´
>
> If the first four characters are not "GRIB", you could get around
this by
> changing the file's suffix to be ".grb" or by using the "file_type =
> GRIB1;" option in the config file (see data/config/README for more
> information).
>
> I will be out of the office starting today at 2:45pm (US Mountain
time)
> through January 4th. I hope I can help you with this problem before
I am
> out of the office.
>
> Julie
>
> On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > Tue Dec 22 03:18:04 2015: Request 74467 was acted upon.
> > Transaction: Ticket created by vanisas at gmail.com
> > Queue: met_help
> > Subject: point_stat for WRF output question
> > Owner: Nobody
> > Requestors: vanisas at gmail.com
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > >
> >
> >
> > Hello there,
> >
> > I got point_stat run through for the tutorial case, and would like
to
> > try on the different domain like Nepal, where you have kindly
helped
> > creating the mask file. I edit the config file as attached. The
error
> > shows up as
> > follow:
> >
> > ~/met-5.1_bugfix$ ./bin/point_stat \
> > >
> >
> >
/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-
> > 20_18:00:00
> > \
> > > tutorial/out/ascii2nc/nepal.nc \
> > > tutorial/config/PointStatConfig_Nepal \ -outdir
> > > tutorial/out/point_stat \ -v 2
> > DEBUG 1: Default Config File:
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> > ult DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > ERROR :
> > ERROR : Trouble reading forecast file
> >
> >
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> > ERROR :
> >
> > There seems to be some error about reading the WRF output file.
Could
> > you give me some clue what kind of error that might be, such that
I
> > may correct, please. Is there something to set in the config file?
I
> > am not sure what could be wrong with the WRF output, which has
worked
> > okay with NCL/other tools.
> >
> > Best wishes,
> > Vanisa
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research Research Applications
Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>
>
>
------------------------------------------------
Subject: point_stat for WRF output question
From: Julie Prestopnik
Time: Mon Jan 04 15:42:31 2016
Hi Vanisa. The file you uploaded may have aged off of our ftp site.
I do
not see it there currently.
I neglected to ask you if your WRF output file was post-processed?
MET is
unable to read the output of WRF directly. Instead, you must first
post-process it using either the Unified PostProcessor (UPP) or the
pinterp
tool (for ARW only). We strongly recommend using UPP since MET's
support
is most robust for the GRIB1 format that UPP creates.
If you have already post-processed your file, please trying uploading
to
our ftp site again. If you have not, please post-process the file and
give
it another try. Thanks!
Happy New Year!
Julie
On Sun, Dec 27, 2015 at 11:37 PM, vanisas at gmail.com via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>
> Dear Julie & all there,
>
> I just manage to try as Julie recommended right before the Christmas
> holidays, i.e. to put the input file I use on the ftp site. I am not
sure
> if it works, or not. What shows up say:
>
> 220 brew.rap.ucar.edu FTP server ready.
> Name (ftp.rap.ucar.edu:vsurapipith): anonymous
> 331 Guest login ok, send your complete e-mail address as password.
> Password:
> 230 Guest login ok, access restrictions apply.
> Remote system type is UNIX.
> Using binary mode to transfer files.
> ftp> cd incoming/irap/met_help
> 250 CWD command successful.
> ftp> mkdir surapipith_data
> 257 "/incoming/irap/met_help/surapipith_data" new directory created.
> ftp> cd surapipith_data
> 250 CWD command successful.
> ftp> put wrfout_d02_2013-06-20_18:00:00
> local: wrfout_d02_2013-06-20_18:00:00 remote:
> wrfout_d02_2013-06-20_18:00:00
> No control connection for command: Success
> Passive mode refused.
> ftp> bye
>
> Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00 | more",
there were
> many lines running up, and I could not find the first 4 line if it
has GRIB
> word or not. I will find the way to try again. Yet I hope you may
see the
> WRF output that I upload to the ftp site of ucar. If not, kindly let
me
> know do it right, please.
>
> Thanks very much.
> Best wishes,
> Vanisa
>
> On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via RT <
> met_help at ucar.edu
> > wrote:
>
> > Hi Julie,
> >
> > Thank you very much for your kind guidance. I am sorry that I just
read
> > this e-mail, since there was no internet connection along my
travel way,
> > whole day, towards home in Thailand. I really appreciate your
effort to
> > help me before the holidays begin; I'm only now at home in
Thailand with
> > internet, but my Nepal workstation with WRF and MET seems to down
now,
> > after midnight of December 23rd in Kathmandu. I can ask the IT
admin to
> > check what happen, but he would only be able to resolve this in
the next
> 8
> > hours (Christmas is actually not a holiday in Nepal, but I'm
taking
> > end-of-year leave rather).
> >
> > Never mind; let's celebrate the Season & give the best to people
at home;
> > I'll find out as you recommend in this while, and report soon, but
I'll
> be
> > okay to find you response later. Wish there a warm Christmas Eve &
Happy
> > coming New Year.
> >
> > Best wishes,
> > Vanisa
> >
> > -----Original Message-----
> > From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
> > Sent: 22 December 2015 23:55
> > To: vanisas at gmail.com
> > Cc: Vanisa Surapipith
> > Subject: Re: [rt.rap.ucar.edu #74467] point_stat for WRF output
question
> >
> > Hi Vanisa.
> >
> > It is difficult to say without having access to the file you are
using.
> > Could you please upload your input forecast and observation file
to our
> > ftp site using the directions on this page:
> >
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > In an effort to try to help before having access to the files, I
will say
> > that we have seen this problem before when the first four
characters in
> the
> > file are not "GRIB". That is called the "magic cookie" to
indicate the
> > file type.
> >
> > You could use the command "od -c wrfout_d02_2013-06-20_18:00:00 |
more"
> > to determine what the first four characters are. You should see
> something
> > like this:
> >
> > 0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0 \0 001
P ´
> >
> > If the first four characters are not "GRIB", you could get around
this by
> > changing the file's suffix to be ".grb" or by using the "file_type
=
> > GRIB1;" option in the config file (see data/config/README for more
> > information).
> >
> > I will be out of the office starting today at 2:45pm (US Mountain
time)
> > through January 4th. I hope I can help you with this problem
before I am
> > out of the office.
> >
> > Julie
> >
> > On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > >
> > > Tue Dec 22 03:18:04 2015: Request 74467 was acted upon.
> > > Transaction: Ticket created by vanisas at gmail.com
> > > Queue: met_help
> > > Subject: point_stat for WRF output question
> > > Owner: Nobody
> > > Requestors: vanisas at gmail.com
> > > Status: new
> > > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > > >
> > >
> > >
> > > Hello there,
> > >
> > > I got point_stat run through for the tutorial case, and would
like to
> > > try on the different domain like Nepal, where you have kindly
helped
> > > creating the mask file. I edit the config file as attached. The
error
> > > shows up as
> > > follow:
> > >
> > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > >
> > >
> > >
/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-
> > > 20_18:00:00
> > > \
> > > > tutorial/out/ascii2nc/nepal.nc \
> > > > tutorial/config/PointStatConfig_Nepal \ -outdir
> > > > tutorial/out/point_stat \ -v 2
> > > DEBUG 1: Default Config File:
> > > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> > > ult DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > ERROR :
> > > ERROR : Trouble reading forecast file
> > >
> > >
> >
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> > > ERROR :
> > >
> > > There seems to be some error about reading the WRF output file.
Could
> > > you give me some clue what kind of error that might be, such
that I
> > > may correct, please. Is there something to set in the config
file? I
> > > am not sure what could be wrong with the WRF output, which has
worked
> > > okay with NCL/other tools.
> > >
> > > Best wishes,
> > > Vanisa
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research Research Applications
Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
> >
> >
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Jan 05 04:13:42 2016
Hi Julie,
Happy New Year...:) Thank you very much for your reply upon back to
office.
Indeed, I did not use UPP for processing the WRF (ARW) output.
Therefore,
today I look into how to install UPP on my workstation for this
purpose.
However, there is a problem that unipost.exe was not generated after
compilation, while other files came through. I was wondered if it has
to do
with liking to JASPER. I have not been able to spot why, and now
attach the
compile_upp.log to you here, in case you might be able to see what is
wrong
in there.
Thanks again; update more soon.
Best wishes,
Vanisa
On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Vanisa. The file you uploaded may have aged off of our ftp site.
I do
> not see it there currently.
>
> I neglected to ask you if your WRF output file was post-processed?
MET is
> unable to read the output of WRF directly. Instead, you must first
> post-process it using either the Unified PostProcessor (UPP) or the
pinterp
> tool (for ARW only). We strongly recommend using UPP since MET's
support
> is most robust for the GRIB1 format that UPP creates.
>
> If you have already post-processed your file, please trying
uploading to
> our ftp site again. If you have not, please post-process the file
and give
> it another try. Thanks!
>
> Happy New Year!
>
> Julie
>
> On Sun, Dec 27, 2015 at 11:37 PM, vanisas at gmail.com via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >
> > Dear Julie & all there,
> >
> > I just manage to try as Julie recommended right before the
Christmas
> > holidays, i.e. to put the input file I use on the ftp site. I am
not sure
> > if it works, or not. What shows up say:
> >
> > 220 brew.rap.ucar.edu FTP server ready.
> > Name (ftp.rap.ucar.edu:vsurapipith): anonymous
> > 331 Guest login ok, send your complete e-mail address as password.
> > Password:
> > 230 Guest login ok, access restrictions apply.
> > Remote system type is UNIX.
> > Using binary mode to transfer files.
> > ftp> cd incoming/irap/met_help
> > 250 CWD command successful.
> > ftp> mkdir surapipith_data
> > 257 "/incoming/irap/met_help/surapipith_data" new directory
created.
> > ftp> cd surapipith_data
> > 250 CWD command successful.
> > ftp> put wrfout_d02_2013-06-20_18:00:00
> > local: wrfout_d02_2013-06-20_18:00:00 remote:
> > wrfout_d02_2013-06-20_18:00:00
> > No control connection for command: Success
> > Passive mode refused.
> > ftp> bye
> >
> > Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00 | more",
there
> were
> > many lines running up, and I could not find the first 4 line if it
has
> GRIB
> > word or not. I will find the way to try again. Yet I hope you may
see the
> > WRF output that I upload to the ftp site of ucar. If not, kindly
let me
> > know do it right, please.
> >
> > Thanks very much.
> > Best wishes,
> > Vanisa
> >
> > On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Hi Julie,
> > >
> > > Thank you very much for your kind guidance. I am sorry that I
just read
> > > this e-mail, since there was no internet connection along my
travel
> way,
> > > whole day, towards home in Thailand. I really appreciate your
effort to
> > > help me before the holidays begin; I'm only now at home in
Thailand
> with
> > > internet, but my Nepal workstation with WRF and MET seems to
down now,
> > > after midnight of December 23rd in Kathmandu. I can ask the IT
admin to
> > > check what happen, but he would only be able to resolve this in
the
> next
> > 8
> > > hours (Christmas is actually not a holiday in Nepal, but I'm
taking
> > > end-of-year leave rather).
> > >
> > > Never mind; let's celebrate the Season & give the best to people
at
> home;
> > > I'll find out as you recommend in this while, and report soon,
but I'll
> > be
> > > okay to find you response later. Wish there a warm Christmas Eve
&
> Happy
> > > coming New Year.
> > >
> > > Best wishes,
> > > Vanisa
> > >
> > > -----Original Message-----
> > > From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
> > > Sent: 22 December 2015 23:55
> > > To: vanisas at gmail.com
> > > Cc: Vanisa Surapipith
> > > Subject: Re: [rt.rap.ucar.edu #74467] point_stat for WRF output
> question
> > >
> > > Hi Vanisa.
> > >
> > > It is difficult to say without having access to the file you are
using.
> > > Could you please upload your input forecast and observation file
to our
> > > ftp site using the directions on this page:
> > >
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > In an effort to try to help before having access to the files, I
will
> say
> > > that we have seen this problem before when the first four
characters in
> > the
> > > file are not "GRIB". That is called the "magic cookie" to
indicate the
> > > file type.
> > >
> > > You could use the command "od -c wrfout_d02_2013-06-20_18:00:00
|
> more"
> > > to determine what the first four characters are. You should see
> > something
> > > like this:
> > >
> > > 0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0 \0 001
P ´
> > >
> > > If the first four characters are not "GRIB", you could get
around this
> by
> > > changing the file's suffix to be ".grb" or by using the
"file_type =
> > > GRIB1;" option in the config file (see data/config/README for
more
> > > information).
> > >
> > > I will be out of the office starting today at 2:45pm (US
Mountain time)
> > > through January 4th. I hope I can help you with this problem
before I
> am
> > > out of the office.
> > >
> > > Julie
> > >
> > > On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > >
> > > > Tue Dec 22 03:18:04 2015: Request 74467 was acted upon.
> > > > Transaction: Ticket created by vanisas at gmail.com
> > > > Queue: met_help
> > > > Subject: point_stat for WRF output question
> > > > Owner: Nobody
> > > > Requestors: vanisas at gmail.com
> > > > Status: new
> > > > Ticket <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > > > >
> > > >
> > > >
> > > > Hello there,
> > > >
> > > > I got point_stat run through for the tutorial case, and would
like to
> > > > try on the different domain like Nepal, where you have kindly
helped
> > > > creating the mask file. I edit the config file as attached.
The error
> > > > shows up as
> > > > follow:
> > > >
> > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > >
> > > >
> > > >
> /home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-
> > > > 20_18:00:00
> > > > \
> > > > > tutorial/out/ascii2nc/nepal.nc \
> > > > > tutorial/config/PointStatConfig_Nepal \ -outdir
> > > > > tutorial/out/point_stat \ -v 2
> > > > DEBUG 1: Default Config File:
> > > >
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> > > > ult DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > > ERROR :
> > > > ERROR : Trouble reading forecast file
> > > >
> > > >
> > >
> >
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> > > > ERROR :
> > > >
> > > > There seems to be some error about reading the WRF output
file. Could
> > > > you give me some clue what kind of error that might be, such
that I
> > > > may correct, please. Is there something to set in the config
file? I
> > > > am not sure what could be wrong with the WRF output, which has
worked
> > > > okay with NCL/other tools.
> > > >
> > > > Best wishes,
> > > > Vanisa
> > > >
> > > >
> > >
> > >
> > > --
> > > Julie Prestopnik
> > > National Center for Atmospheric Research Research Applications
> Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> > >
> > >
> > >
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>
--
"Our ignorance is not so vast as our failure to use what we know."
M. King Hubbert
October 5th, 1903 -- October 11th, 1989
------------------------------------------------
Subject: point_stat for WRF output question
From: Julie Prestopnik
Time: Tue Jan 05 10:21:29 2016
Hi Vanisa.
While I don't have any experience installing UPP myself, I did take a
look
at your log file where I saw:
enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No such file or
directory compilation terminated.
So, it looks like UPP is having trouble finding jasper.h.
Support for UPP is provided via wrfhelp at ucar.edu. If you're unable
to
resolve the jasper problem, I would suggest writing to them regarding
this
issue.
Thanks,
Julie
On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>
> Hi Julie,
>
> Happy New Year...:) Thank you very much for your reply upon back to
office.
> Indeed, I did not use UPP for processing the WRF (ARW) output.
Therefore,
> today I look into how to install UPP on my workstation for this
purpose.
> However, there is a problem that unipost.exe was not generated after
> compilation, while other files came through. I was wondered if it
has to do
> with liking to JASPER. I have not been able to spot why, and now
attach the
> compile_upp.log to you here, in case you might be able to see what
is wrong
> in there.
>
> Thanks again; update more soon.
>
> Best wishes,
> Vanisa
>
> On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT
<met_help at ucar.edu
> >
> wrote:
>
> > Hi Vanisa. The file you uploaded may have aged off of our ftp
site. I
> do
> > not see it there currently.
> >
> > I neglected to ask you if your WRF output file was post-processed?
MET
> is
> > unable to read the output of WRF directly. Instead, you must
first
> > post-process it using either the Unified PostProcessor (UPP) or
the
> pinterp
> > tool (for ARW only). We strongly recommend using UPP since MET's
support
> > is most robust for the GRIB1 format that UPP creates.
> >
> > If you have already post-processed your file, please trying
uploading to
> > our ftp site again. If you have not, please post-process the file
and
> give
> > it another try. Thanks!
> >
> > Happy New Year!
> >
> > Julie
> >
> > On Sun, Dec 27, 2015 at 11:37 PM, vanisas at gmail.com via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > >
> > > Dear Julie & all there,
> > >
> > > I just manage to try as Julie recommended right before the
Christmas
> > > holidays, i.e. to put the input file I use on the ftp site. I am
not
> sure
> > > if it works, or not. What shows up say:
> > >
> > > 220 brew.rap.ucar.edu FTP server ready.
> > > Name (ftp.rap.ucar.edu:vsurapipith): anonymous
> > > 331 Guest login ok, send your complete e-mail address as
password.
> > > Password:
> > > 230 Guest login ok, access restrictions apply.
> > > Remote system type is UNIX.
> > > Using binary mode to transfer files.
> > > ftp> cd incoming/irap/met_help
> > > 250 CWD command successful.
> > > ftp> mkdir surapipith_data
> > > 257 "/incoming/irap/met_help/surapipith_data" new directory
created.
> > > ftp> cd surapipith_data
> > > 250 CWD command successful.
> > > ftp> put wrfout_d02_2013-06-20_18:00:00
> > > local: wrfout_d02_2013-06-20_18:00:00 remote:
> > > wrfout_d02_2013-06-20_18:00:00
> > > No control connection for command: Success
> > > Passive mode refused.
> > > ftp> bye
> > >
> > > Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00 | more",
there
> > were
> > > many lines running up, and I could not find the first 4 line if
it has
> > GRIB
> > > word or not. I will find the way to try again. Yet I hope you
may see
> the
> > > WRF output that I upload to the ftp site of ucar. If not, kindly
let me
> > > know do it right, please.
> > >
> > > Thanks very much.
> > > Best wishes,
> > > Vanisa
> > >
> > > On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hi Julie,
> > > >
> > > > Thank you very much for your kind guidance. I am sorry that I
just
> read
> > > > this e-mail, since there was no internet connection along my
travel
> > way,
> > > > whole day, towards home in Thailand. I really appreciate your
effort
> to
> > > > help me before the holidays begin; I'm only now at home in
Thailand
> > with
> > > > internet, but my Nepal workstation with WRF and MET seems to
down
> now,
> > > > after midnight of December 23rd in Kathmandu. I can ask the IT
admin
> to
> > > > check what happen, but he would only be able to resolve this
in the
> > next
> > > 8
> > > > hours (Christmas is actually not a holiday in Nepal, but I'm
taking
> > > > end-of-year leave rather).
> > > >
> > > > Never mind; let's celebrate the Season & give the best to
people at
> > home;
> > > > I'll find out as you recommend in this while, and report soon,
but
> I'll
> > > be
> > > > okay to find you response later. Wish there a warm Christmas
Eve &
> > Happy
> > > > coming New Year.
> > > >
> > > > Best wishes,
> > > > Vanisa
> > > >
> > > > -----Original Message-----
> > > > From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
> > > > Sent: 22 December 2015 23:55
> > > > To: vanisas at gmail.com
> > > > Cc: Vanisa Surapipith
> > > > Subject: Re: [rt.rap.ucar.edu #74467] point_stat for WRF
output
> > question
> > > >
> > > > Hi Vanisa.
> > > >
> > > > It is difficult to say without having access to the file you
are
> using.
> > > > Could you please upload your input forecast and observation
file to
> our
> > > > ftp site using the directions on this page:
> > > >
> > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > In an effort to try to help before having access to the files,
I will
> > say
> > > > that we have seen this problem before when the first four
characters
> in
> > > the
> > > > file are not "GRIB". That is called the "magic cookie" to
indicate
> the
> > > > file type.
> > > >
> > > > You could use the command "od -c wrfout_d02_2013-06-
20_18:00:00 |
> > more"
> > > > to determine what the first four characters are. You should
see
> > > something
> > > > like this:
> > > >
> > > > 0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0 \0
001 P
> ´
> > > >
> > > > If the first four characters are not "GRIB", you could get
around
> this
> > by
> > > > changing the file's suffix to be ".grb" or by using the
"file_type =
> > > > GRIB1;" option in the config file (see data/config/README for
more
> > > > information).
> > > >
> > > > I will be out of the office starting today at 2:45pm (US
Mountain
> time)
> > > > through January 4th. I hope I can help you with this problem
before
> I
> > am
> > > > out of the office.
> > > >
> > > > Julie
> > > >
> > > > On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > >
> > > > > Tue Dec 22 03:18:04 2015: Request 74467 was acted upon.
> > > > > Transaction: Ticket created by vanisas at gmail.com
> > > > > Queue: met_help
> > > > > Subject: point_stat for WRF output question
> > > > > Owner: Nobody
> > > > > Requestors: vanisas at gmail.com
> > > > > Status: new
> > > > > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > > > > >
> > > > >
> > > > >
> > > > > Hello there,
> > > > >
> > > > > I got point_stat run through for the tutorial case, and
would like
> to
> > > > > try on the different domain like Nepal, where you have
kindly
> helped
> > > > > creating the mask file. I edit the config file as attached.
The
> error
> > > > > shows up as
> > > > > follow:
> > > > >
> > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > >
> > > > >
> > > > >
> >
/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-
> > > > > 20_18:00:00
> > > > > \
> > > > > > tutorial/out/ascii2nc/nepal.nc \
> > > > > > tutorial/config/PointStatConfig_Nepal \ -outdir
> > > > > > tutorial/out/point_stat \ -v 2
> > > > > DEBUG 1: Default Config File:
> > > > >
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> > > > > ult DEBUG 1: User Config File:
> tutorial/config/PointStatConfig_Nepal
> > > > > ERROR :
> > > > > ERROR : Trouble reading forecast file
> > > > >
> > > > >
> > > >
> > >
> >
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> > > > > ERROR :
> > > > >
> > > > > There seems to be some error about reading the WRF output
file.
> Could
> > > > > you give me some clue what kind of error that might be, such
that I
> > > > > may correct, please. Is there something to set in the config
file?
> I
> > > > > am not sure what could be wrong with the WRF output, which
has
> worked
> > > > > okay with NCL/other tools.
> > > > >
> > > > > Best wishes,
> > > > > Vanisa
> > > > >
> > > > >
> > > >
> > > >
> > > > --
> > > > Julie Prestopnik
> > > > National Center for Atmospheric Research Research Applications
> > Laboratory
> > > > Phone: 303.497.8399
> > > > Email: jpresto at ucar.edu
> > > >
> > > >
> > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>
> --
> "Our ignorance is not so vast as our failure to use what we know."
> M. King Hubbert
> October 5th, 1903 -- October 11th, 1989
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Jan 26 04:01:54 2016
Hi Julie,
I'm back....after a couple of weeks spending on stumbling over the
installation of UPP until I got it through, now I have the GRIB1 files
systematically converted from the WRF output generated here.
I try the GRIB1 files with the PointStatConfig_Nepal that I sent
earlier to
do Point_Stat with our these WRF output, but currently the error shows
up
as:
~/met-5.1_bugfix$ ./bin/point_stat
/home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
tutorial/out/ascii2nc/nepal.nc tutorial/config/PointStatConfig_Nepal
-outdir tutorial/out/point_stat -v 2
DEBUG 1: Default Config File:
/home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
WARNING:
WARNING: parse_conf_interp() -> For neareast neighbor interpolation
method,
resetting method from "UW_MEAN" to "NEAREST" since the interpolation
width
is 1.
WARNING:
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744072070604630
DEBUG 1: Forecast File:
/home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SPFH/P500.
ERROR :
ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is not
defined.
ERROR :
It looks to me like the mask file - nepal.nc, that I deliberately took
from
what you generated earlier - was not seen/defined well by MET here. I
have
actually just put the WRF output files in GRIB1 format (checked with
as you
recommend) on the rap.ucar.edu ftp site (surapipith_data)...Would you
kindly let me know what should be done to correct this, please?
Best wishes,
Vanisa
On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Vanisa.
>
> While I don't have any experience installing UPP myself, I did take
a look
> at your log file where I saw:
>
> enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No such file or
> directory compilation terminated.
>
> So, it looks like UPP is having trouble finding jasper.h.
>
> Support for UPP is provided via wrfhelp at ucar.edu. If you're
unable to
> resolve the jasper problem, I would suggest writing to them
regarding this
> issue.
>
> Thanks,
> Julie
>
>
> On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via RT <
> met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >
> > Hi Julie,
> >
> > Happy New Year...:) Thank you very much for your reply upon back
to
> office.
> > Indeed, I did not use UPP for processing the WRF (ARW) output.
Therefore,
> > today I look into how to install UPP on my workstation for this
purpose.
> > However, there is a problem that unipost.exe was not generated
after
> > compilation, while other files came through. I was wondered if it
has to
> do
> > with liking to JASPER. I have not been able to spot why, and now
attach
> the
> > compile_upp.log to you here, in case you might be able to see what
is
> wrong
> > in there.
> >
> > Thanks again; update more soon.
> >
> > Best wishes,
> > Vanisa
> >
> > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT <
> met_help at ucar.edu
> > >
> > wrote:
> >
> > > Hi Vanisa. The file you uploaded may have aged off of our ftp
site. I
> > do
> > > not see it there currently.
> > >
> > > I neglected to ask you if your WRF output file was post-
processed? MET
> > is
> > > unable to read the output of WRF directly. Instead, you must
first
> > > post-process it using either the Unified PostProcessor (UPP) or
the
> > pinterp
> > > tool (for ARW only). We strongly recommend using UPP since
MET's
> support
> > > is most robust for the GRIB1 format that UPP creates.
> > >
> > > If you have already post-processed your file, please trying
uploading
> to
> > > our ftp site again. If you have not, please post-process the
file and
> > give
> > > it another try. Thanks!
> > >
> > > Happy New Year!
> > >
> > > Julie
> > >
> > > On Sun, Dec 27, 2015 at 11:37 PM, vanisas at gmail.com via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
>
> > > >
> > > > Dear Julie & all there,
> > > >
> > > > I just manage to try as Julie recommended right before the
Christmas
> > > > holidays, i.e. to put the input file I use on the ftp site. I
am not
> > sure
> > > > if it works, or not. What shows up say:
> > > >
> > > > 220 brew.rap.ucar.edu FTP server ready.
> > > > Name (ftp.rap.ucar.edu:vsurapipith): anonymous
> > > > 331 Guest login ok, send your complete e-mail address as
password.
> > > > Password:
> > > > 230 Guest login ok, access restrictions apply.
> > > > Remote system type is UNIX.
> > > > Using binary mode to transfer files.
> > > > ftp> cd incoming/irap/met_help
> > > > 250 CWD command successful.
> > > > ftp> mkdir surapipith_data
> > > > 257 "/incoming/irap/met_help/surapipith_data" new directory
created.
> > > > ftp> cd surapipith_data
> > > > 250 CWD command successful.
> > > > ftp> put wrfout_d02_2013-06-20_18:00:00
> > > > local: wrfout_d02_2013-06-20_18:00:00 remote:
> > > > wrfout_d02_2013-06-20_18:00:00
> > > > No control connection for command: Success
> > > > Passive mode refused.
> > > > ftp> bye
> > > >
> > > > Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00 |
more", there
> > > were
> > > > many lines running up, and I could not find the first 4 line
if it
> has
> > > GRIB
> > > > word or not. I will find the way to try again. Yet I hope you
may see
> > the
> > > > WRF output that I upload to the ftp site of ucar. If not,
kindly let
> me
> > > > know do it right, please.
> > > >
> > > > Thanks very much.
> > > > Best wishes,
> > > > Vanisa
> > > >
> > > > On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via RT <
> > > > met_help at ucar.edu
> > > > > wrote:
> > > >
> > > > > Hi Julie,
> > > > >
> > > > > Thank you very much for your kind guidance. I am sorry that
I just
> > read
> > > > > this e-mail, since there was no internet connection along my
travel
> > > way,
> > > > > whole day, towards home in Thailand. I really appreciate
your
> effort
> > to
> > > > > help me before the holidays begin; I'm only now at home in
Thailand
> > > with
> > > > > internet, but my Nepal workstation with WRF and MET seems to
down
> > now,
> > > > > after midnight of December 23rd in Kathmandu. I can ask the
IT
> admin
> > to
> > > > > check what happen, but he would only be able to resolve this
in the
> > > next
> > > > 8
> > > > > hours (Christmas is actually not a holiday in Nepal, but I'm
taking
> > > > > end-of-year leave rather).
> > > > >
> > > > > Never mind; let's celebrate the Season & give the best to
people at
> > > home;
> > > > > I'll find out as you recommend in this while, and report
soon, but
> > I'll
> > > > be
> > > > > okay to find you response later. Wish there a warm Christmas
Eve &
> > > Happy
> > > > > coming New Year.
> > > > >
> > > > > Best wishes,
> > > > > Vanisa
> > > > >
> > > > > -----Original Message-----
> > > > > From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
> > > > > Sent: 22 December 2015 23:55
> > > > > To: vanisas at gmail.com
> > > > > Cc: Vanisa Surapipith
> > > > > Subject: Re: [rt.rap.ucar.edu #74467] point_stat for WRF
output
> > > question
> > > > >
> > > > > Hi Vanisa.
> > > > >
> > > > > It is difficult to say without having access to the file you
are
> > using.
> > > > > Could you please upload your input forecast and observation
file to
> > our
> > > > > ftp site using the directions on this page:
> > > > >
> > > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > >
> > > > > In an effort to try to help before having access to the
files, I
> will
> > > say
> > > > > that we have seen this problem before when the first four
> characters
> > in
> > > > the
> > > > > file are not "GRIB". That is called the "magic cookie" to
indicate
> > the
> > > > > file type.
> > > > >
> > > > > You could use the command "od -c wrfout_d02_2013-06-
20_18:00:00 |
> > > more"
> > > > > to determine what the first four characters are. You should
see
> > > > something
> > > > > like this:
> > > > >
> > > > > 0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0 \0
001 P
> > ´
> > > > >
> > > > > If the first four characters are not "GRIB", you could get
around
> > this
> > > by
> > > > > changing the file's suffix to be ".grb" or by using the
"file_type
> =
> > > > > GRIB1;" option in the config file (see data/config/README
for more
> > > > > information).
> > > > >
> > > > > I will be out of the office starting today at 2:45pm (US
Mountain
> > time)
> > > > > through January 4th. I hope I can help you with this
problem
> before
> > I
> > > am
> > > > > out of the office.
> > > > >
> > > > > Julie
> > > > >
> > > > > On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via RT <
> > > > > met_help at ucar.edu
> > > > > > wrote:
> > > > >
> > > > > >
> > > > > > Tue Dec 22 03:18:04 2015: Request 74467 was acted upon.
> > > > > > Transaction: Ticket created by vanisas at gmail.com
> > > > > > Queue: met_help
> > > > > > Subject: point_stat for WRF output question
> > > > > > Owner: Nobody
> > > > > > Requestors: vanisas at gmail.com
> > > > > > Status: new
> > > > > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > > > > > >
> > > > > >
> > > > > >
> > > > > > Hello there,
> > > > > >
> > > > > > I got point_stat run through for the tutorial case, and
would
> like
> > to
> > > > > > try on the different domain like Nepal, where you have
kindly
> > helped
> > > > > > creating the mask file. I edit the config file as
attached. The
> > error
> > > > > > shows up as
> > > > > > follow:
> > > > > >
> > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > > >
> > > > > >
> > > > > >
> > >
/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-
> > > > > > 20_18:00:00
> > > > > > \
> > > > > > > tutorial/out/ascii2nc/nepal.nc \
> > > > > > > tutorial/config/PointStatConfig_Nepal \ -outdir
> > > > > > > tutorial/out/point_stat \ -v 2
> > > > > > DEBUG 1: Default Config File:
> > > > > >
> > > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> > > > > > ult DEBUG 1: User Config File:
> > tutorial/config/PointStatConfig_Nepal
> > > > > > ERROR :
> > > > > > ERROR : Trouble reading forecast file
> > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> > > > > > ERROR :
> > > > > >
> > > > > > There seems to be some error about reading the WRF output
file.
> > Could
> > > > > > you give me some clue what kind of error that might be,
such
> that I
> > > > > > may correct, please. Is there something to set in the
config
> file?
> > I
> > > > > > am not sure what could be wrong with the WRF output, which
has
> > worked
> > > > > > okay with NCL/other tools.
> > > > > >
> > > > > > Best wishes,
> > > > > > Vanisa
> > > > > >
> > > > > >
> > > > >
> > > > >
> > > > > --
> > > > > Julie Prestopnik
> > > > > National Center for Atmospheric Research Research
Applications
> > > Laboratory
> > > > > Phone: 303.497.8399
> > > > > Email: jpresto at ucar.edu
> > > > >
>
------------------------------------------------
Subject: point_stat for WRF output question
From: John Halley Gotway
Time: Tue Jan 26 09:59:36 2016
Hello Vanisa,
My name is John, and I work with Julie who is out of the office today.
I
see that you're having some trouble getting Point-Stat to run:
ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is
not
defined.
This should be pretty easy to resolve. You've worked with Julie to
generate a mask file for Nepal, and that mask is contained in some
NetCDF
file (the output of the gen_vx_mask tool). For example, let's say
it's in
"/home/vsurapipith/nepal_mask.nc". You'll need to put the path to
that
file in the "mask" section of the Point-Stat configuration file:
mask = {
grid = [];
poly = [ "/home/vsurapipith/nepal_mask.nc" ]; // or whatever the
actual path is to that data file
sid = [];
}
And then try running Point-Stat again. For a more complete
description of
the masking options available, please see pages 54-55 of the MET
User's
Guide:
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
Please let me know how it goes.
Thanks,
John
On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>
> Hi Julie,
>
> I'm back....after a couple of weeks spending on stumbling over the
> installation of UPP until I got it through, now I have the GRIB1
files
> systematically converted from the WRF output generated here.
>
> I try the GRIB1 files with the PointStatConfig_Nepal that I sent
earlier to
> do Point_Stat with our these WRF output, but currently the error
shows up
> as:
>
> ~/met-5.1_bugfix$ ./bin/point_stat
> /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> tutorial/out/ascii2nc/nepal.nc tutorial/config/PointStatConfig_Nepal
> -outdir tutorial/out/point_stat -v 2
> DEBUG 1: Default Config File:
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> WARNING:
> WARNING: parse_conf_interp() -> For neareast neighbor interpolation
method,
> resetting method from "UW_MEAN" to "NEAREST" since the interpolation
width
> is 1.
> WARNING:
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744072070604630
> DEBUG 1: Forecast File:
> /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> DEBUG 2:
> DEBUG 2:
>
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for SPFH/P500.
> ERROR :
> ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is not
> defined.
> ERROR :
>
> It looks to me like the mask file - nepal.nc, that I deliberately
took
> from
> what you generated earlier - was not seen/defined well by MET here.
I have
> actually just put the WRF output files in GRIB1 format (checked with
as you
> recommend) on the rap.ucar.edu ftp site (surapipith_data)...Would
you
> kindly let me know what should be done to correct this, please?
>
> Best wishes,
> Vanisa
>
> On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Vanisa.
> >
> > While I don't have any experience installing UPP myself, I did
take a
> look
> > at your log file where I saw:
> >
> > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No such file or
> > directory compilation terminated.
> >
> > So, it looks like UPP is having trouble finding jasper.h.
> >
> > Support for UPP is provided via wrfhelp at ucar.edu. If you're
unable
> to
> > resolve the jasper problem, I would suggest writing to them
regarding
> this
> > issue.
> >
> > Thanks,
> > Julie
> >
> >
> > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > >
> > > Hi Julie,
> > >
> > > Happy New Year...:) Thank you very much for your reply upon back
to
> > office.
> > > Indeed, I did not use UPP for processing the WRF (ARW) output.
> Therefore,
> > > today I look into how to install UPP on my workstation for this
> purpose.
> > > However, there is a problem that unipost.exe was not generated
after
> > > compilation, while other files came through. I was wondered if
it has
> to
> > do
> > > with liking to JASPER. I have not been able to spot why, and now
attach
> > the
> > > compile_upp.log to you here, in case you might be able to see
what is
> > wrong
> > > in there.
> > >
> > > Thanks again; update more soon.
> > >
> > > Best wishes,
> > > Vanisa
> > >
> > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > >
> > > wrote:
> > >
> > > > Hi Vanisa. The file you uploaded may have aged off of our ftp
> site. I
> > > do
> > > > not see it there currently.
> > > >
> > > > I neglected to ask you if your WRF output file was post-
processed?
> MET
> > > is
> > > > unable to read the output of WRF directly. Instead, you must
first
> > > > post-process it using either the Unified PostProcessor (UPP)
or the
> > > pinterp
> > > > tool (for ARW only). We strongly recommend using UPP since
MET's
> > support
> > > > is most robust for the GRIB1 format that UPP creates.
> > > >
> > > > If you have already post-processed your file, please trying
uploading
> > to
> > > > our ftp site again. If you have not, please post-process the
file
> and
> > > give
> > > > it another try. Thanks!
> > > >
> > > > Happy New Year!
> > > >
> > > > Julie
> > > >
> > > > On Sun, Dec 27, 2015 at 11:37 PM, vanisas at gmail.com via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > > > >
> > > > > Dear Julie & all there,
> > > > >
> > > > > I just manage to try as Julie recommended right before the
> Christmas
> > > > > holidays, i.e. to put the input file I use on the ftp site.
I am
> not
> > > sure
> > > > > if it works, or not. What shows up say:
> > > > >
> > > > > 220 brew.rap.ucar.edu FTP server ready.
> > > > > Name (ftp.rap.ucar.edu:vsurapipith): anonymous
> > > > > 331 Guest login ok, send your complete e-mail address as
password.
> > > > > Password:
> > > > > 230 Guest login ok, access restrictions apply.
> > > > > Remote system type is UNIX.
> > > > > Using binary mode to transfer files.
> > > > > ftp> cd incoming/irap/met_help
> > > > > 250 CWD command successful.
> > > > > ftp> mkdir surapipith_data
> > > > > 257 "/incoming/irap/met_help/surapipith_data" new directory
> created.
> > > > > ftp> cd surapipith_data
> > > > > 250 CWD command successful.
> > > > > ftp> put wrfout_d02_2013-06-20_18:00:00
> > > > > local: wrfout_d02_2013-06-20_18:00:00 remote:
> > > > > wrfout_d02_2013-06-20_18:00:00
> > > > > No control connection for command: Success
> > > > > Passive mode refused.
> > > > > ftp> bye
> > > > >
> > > > > Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00 |
more",
> there
> > > > were
> > > > > many lines running up, and I could not find the first 4 line
if it
> > has
> > > > GRIB
> > > > > word or not. I will find the way to try again. Yet I hope
you may
> see
> > > the
> > > > > WRF output that I upload to the ftp site of ucar. If not,
kindly
> let
> > me
> > > > > know do it right, please.
> > > > >
> > > > > Thanks very much.
> > > > > Best wishes,
> > > > > Vanisa
> > > > >
> > > > > On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via RT <
> > > > > met_help at ucar.edu
> > > > > > wrote:
> > > > >
> > > > > > Hi Julie,
> > > > > >
> > > > > > Thank you very much for your kind guidance. I am sorry
that I
> just
> > > read
> > > > > > this e-mail, since there was no internet connection along
my
> travel
> > > > way,
> > > > > > whole day, towards home in Thailand. I really appreciate
your
> > effort
> > > to
> > > > > > help me before the holidays begin; I'm only now at home in
> Thailand
> > > > with
> > > > > > internet, but my Nepal workstation with WRF and MET seems
to down
> > > now,
> > > > > > after midnight of December 23rd in Kathmandu. I can ask
the IT
> > admin
> > > to
> > > > > > check what happen, but he would only be able to resolve
this in
> the
> > > > next
> > > > > 8
> > > > > > hours (Christmas is actually not a holiday in Nepal, but
I'm
> taking
> > > > > > end-of-year leave rather).
> > > > > >
> > > > > > Never mind; let's celebrate the Season & give the best to
people
> at
> > > > home;
> > > > > > I'll find out as you recommend in this while, and report
soon,
> but
> > > I'll
> > > > > be
> > > > > > okay to find you response later. Wish there a warm
Christmas Eve
> &
> > > > Happy
> > > > > > coming New Year.
> > > > > >
> > > > > > Best wishes,
> > > > > > Vanisa
> > > > > >
> > > > > > -----Original Message-----
> > > > > > From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
> > > > > > Sent: 22 December 2015 23:55
> > > > > > To: vanisas at gmail.com
> > > > > > Cc: Vanisa Surapipith
> > > > > > Subject: Re: [rt.rap.ucar.edu #74467] point_stat for WRF
output
> > > > question
> > > > > >
> > > > > > Hi Vanisa.
> > > > > >
> > > > > > It is difficult to say without having access to the file
you are
> > > using.
> > > > > > Could you please upload your input forecast and
observation file
> to
> > > our
> > > > > > ftp site using the directions on this page:
> > > > > >
> > > > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > > >
> > > > > > In an effort to try to help before having access to the
files, I
> > will
> > > > say
> > > > > > that we have seen this problem before when the first four
> > characters
> > > in
> > > > > the
> > > > > > file are not "GRIB". That is called the "magic cookie" to
> indicate
> > > the
> > > > > > file type.
> > > > > >
> > > > > > You could use the command "od -c wrfout_d02_2013-06-
20_18:00:00
> |
> > > > more"
> > > > > > to determine what the first four characters are. You
should see
> > > > > something
> > > > > > like this:
> > > > > >
> > > > > > 0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0
\0 001
> P
> > > ´
> > > > > >
> > > > > > If the first four characters are not "GRIB", you could get
around
> > > this
> > > > by
> > > > > > changing the file's suffix to be ".grb" or by using the
> "file_type
> > =
> > > > > > GRIB1;" option in the config file (see data/config/README
for
> more
> > > > > > information).
> > > > > >
> > > > > > I will be out of the office starting today at 2:45pm (US
Mountain
> > > time)
> > > > > > through January 4th. I hope I can help you with this
problem
> > before
> > > I
> > > > am
> > > > > > out of the office.
> > > > > >
> > > > > > Julie
> > > > > >
> > > > > > On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via RT
<
> > > > > > met_help at ucar.edu
> > > > > > > wrote:
> > > > > >
> > > > > > >
> > > > > > > Tue Dec 22 03:18:04 2015: Request 74467 was acted upon.
> > > > > > > Transaction: Ticket created by vanisas at gmail.com
> > > > > > > Queue: met_help
> > > > > > > Subject: point_stat for WRF output question
> > > > > > > Owner: Nobody
> > > > > > > Requestors: vanisas at gmail.com
> > > > > > > Status: new
> > > > > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > Hello there,
> > > > > > >
> > > > > > > I got point_stat run through for the tutorial case, and
would
> > like
> > > to
> > > > > > > try on the different domain like Nepal, where you have
kindly
> > > helped
> > > > > > > creating the mask file. I edit the config file as
attached. The
> > > error
> > > > > > > shows up as
> > > > > > > follow:
> > > > > > >
> > > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > > > >
> > > > > > >
> > > > > > >
> > > >
> /home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-
> > > > > > > 20_18:00:00
> > > > > > > \
> > > > > > > > tutorial/out/ascii2nc/nepal.nc \
> > > > > > > > tutorial/config/PointStatConfig_Nepal \ -outdir
> > > > > > > > tutorial/out/point_stat \ -v 2
> > > > > > > DEBUG 1: Default Config File:
> > > > > > >
> > > >
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> > > > > > > ult DEBUG 1: User Config File:
> > > tutorial/config/PointStatConfig_Nepal
> > > > > > > ERROR :
> > > > > > > ERROR : Trouble reading forecast file
> > > > > > >
> > > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> > > > > > > ERROR :
> > > > > > >
> > > > > > > There seems to be some error about reading the WRF
output file.
> > > Could
> > > > > > > you give me some clue what kind of error that might be,
such
> > that I
> > > > > > > may correct, please. Is there something to set in the
config
> > file?
> > > I
> > > > > > > am not sure what could be wrong with the WRF output,
which has
> > > worked
> > > > > > > okay with NCL/other tools.
> > > > > > >
> > > > > > > Best wishes,
> > > > > > > Vanisa
> > > > > > >
> > > > > > >
> > > > > >
> > > > > >
> > > > > > --
> > > > > > Julie Prestopnik
> > > > > > National Center for Atmospheric Research Research
Applications
> > > > Laboratory
> > > > > > Phone: 303.497.8399
> > > > > > Email: jpresto at ucar.edu
> > > > > >
> >
>
>
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Jan 26 22:03:41 2016
Hi John,
thank you for your reply; I recall your help earlier, and I am only
gradually getting to know MET. I am grateful to the help as to what
Julie
was responding patiently as well. I notice the (presumably) mask files
that
Julie helped generating for Nepal, and I try them. However, the
similar
error still shows up, and I am trying to understand why:
~/met-5.1_bugfix$ ./bin/point_stat \
> /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> tutorial/out/ascii2nc/nepal.nc tutorial/config/PointStatConfig_Nepal
\
> -outdir tutorial/out/point_stat \
> -v 2
DEBUG 1: Default Config File:
/home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
WARNING:
WARNING: parse_conf_interp() -> For neareast neighbor interpolation
method,
resetting method from "UW_MEAN" to "NEAREST" since the interpolation
width
is 1.
WARNING:
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=1682548822
DEBUG 1: Forecast File:
/home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SPFH/P500.
ERROR :
ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is not
defined.
ERROR :
Kindly see in the surapipith_data folder in rap.ucar.edu site the
latest
config file, as well as the suite of Nepal mask files that Julie made
and I
use. Hope some progress would be before the end of the day, or once
some
hints are next sent from there...
Best wishes,
Vanisa
On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway via RT <
met_help at ucar.edu> wrote:
> Hello Vanisa,
>
> My name is John, and I work with Julie who is out of the office
today. I
> see that you're having some trouble getting Point-Stat to run:
>
> ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is
not
> defined.
>
> This should be pretty easy to resolve. You've worked with Julie to
> generate a mask file for Nepal, and that mask is contained in some
NetCDF
> file (the output of the gen_vx_mask tool). For example, let's say
it's in
> "/home/vsurapipith/nepal_mask.nc". You'll need to put the path to
that
> file in the "mask" section of the Point-Stat configuration file:
>
> mask = {
> grid = [];
> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; // or whatever
the
> actual path is to that data file
> sid = [];
> }
>
> And then try running Point-Stat again. For a more complete
description of
> the masking options available, please see pages 54-55 of the MET
User's
> Guide:
>
>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
>
> Please let me know how it goes.
>
> Thanks,
> John
>
>
>
> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via RT <
> met_help at ucar.edu
> > wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >
> > Hi Julie,
> >
> > I'm back....after a couple of weeks spending on stumbling over the
> > installation of UPP until I got it through, now I have the GRIB1
files
> > systematically converted from the WRF output generated here.
> >
> > I try the GRIB1 files with the PointStatConfig_Nepal that I sent
earlier
> to
> > do Point_Stat with our these WRF output, but currently the error
shows up
> > as:
> >
> > ~/met-5.1_bugfix$ ./bin/point_stat
> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > tutorial/out/ascii2nc/nepal.nc
tutorial/config/PointStatConfig_Nepal
> > -outdir tutorial/out/point_stat -v 2
> > DEBUG 1: Default Config File:
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> > WARNING:
> > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> method,
> > resetting method from "UW_MEAN" to "NEAREST" since the
interpolation
> width
> > is 1.
> > WARNING:
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=18446744072070604630
> > DEBUG 1: Forecast File:
> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > DEBUG 2:
> > DEBUG 2:
> >
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for SPFH/P500.
> > ERROR :
> > ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is
not
> > defined.
> > ERROR :
> >
> > It looks to me like the mask file - nepal.nc, that I deliberately
took
> > from
> > what you generated earlier - was not seen/defined well by MET
here. I
> have
> > actually just put the WRF output files in GRIB1 format (checked
with as
> you
> > recommend) on the rap.ucar.edu ftp site (surapipith_data)...Would
you
> > kindly let me know what should be done to correct this, please?
> >
> > Best wishes,
> > Vanisa
> >
> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > > Hi Vanisa.
> > >
> > > While I don't have any experience installing UPP myself, I did
take a
> > look
> > > at your log file where I saw:
> > >
> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No such file
or
> > > directory compilation terminated.
> > >
> > > So, it looks like UPP is having trouble finding jasper.h.
> > >
> > > Support for UPP is provided via wrfhelp at ucar.edu. If you're
unable
> > to
> > > resolve the jasper problem, I would suggest writing to them
regarding
> > this
> > > issue.
> > >
> > > Thanks,
> > > Julie
> > >
> > >
> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via RT <
> > > met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
>
> > > >
> > > > Hi Julie,
> > > >
> > > > Happy New Year...:) Thank you very much for your reply upon
back to
> > > office.
> > > > Indeed, I did not use UPP for processing the WRF (ARW) output.
> > Therefore,
> > > > today I look into how to install UPP on my workstation for
this
> > purpose.
> > > > However, there is a problem that unipost.exe was not generated
after
> > > > compilation, while other files came through. I was wondered if
it has
> > to
> > > do
> > > > with liking to JASPER. I have not been able to spot why, and
now
> attach
> > > the
> > > > compile_upp.log to you here, in case you might be able to see
what is
> > > wrong
> > > > in there.
> > > >
> > > > Thanks again; update more soon.
> > > >
> > > > Best wishes,
> > > > Vanisa
> > > >
> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > >
> > > > wrote:
> > > >
> > > > > Hi Vanisa. The file you uploaded may have aged off of our
ftp
> > site. I
> > > > do
> > > > > not see it there currently.
> > > > >
> > > > > I neglected to ask you if your WRF output file was post-
processed?
> > MET
> > > > is
> > > > > unable to read the output of WRF directly. Instead, you
must first
> > > > > post-process it using either the Unified PostProcessor (UPP)
or the
> > > > pinterp
> > > > > tool (for ARW only). We strongly recommend using UPP since
MET's
> > > support
> > > > > is most robust for the GRIB1 format that UPP creates.
> > > > >
> > > > > If you have already post-processed your file, please trying
> uploading
> > > to
> > > > > our ftp site again. If you have not, please post-process
the file
> > and
> > > > give
> > > > > it another try. Thanks!
> > > > >
> > > > > Happy New Year!
> > > > >
> > > > > Julie
> > > > >
> > > > > On Sun, Dec 27, 2015 at 11:37 PM, vanisas at gmail.com via RT <
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > > >
> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > > > > >
> > > > > > Dear Julie & all there,
> > > > > >
> > > > > > I just manage to try as Julie recommended right before the
> > Christmas
> > > > > > holidays, i.e. to put the input file I use on the ftp
site. I am
> > not
> > > > sure
> > > > > > if it works, or not. What shows up say:
> > > > > >
> > > > > > 220 brew.rap.ucar.edu FTP server ready.
> > > > > > Name (ftp.rap.ucar.edu:vsurapipith): anonymous
> > > > > > 331 Guest login ok, send your complete e-mail address as
> password.
> > > > > > Password:
> > > > > > 230 Guest login ok, access restrictions apply.
> > > > > > Remote system type is UNIX.
> > > > > > Using binary mode to transfer files.
> > > > > > ftp> cd incoming/irap/met_help
> > > > > > 250 CWD command successful.
> > > > > > ftp> mkdir surapipith_data
> > > > > > 257 "/incoming/irap/met_help/surapipith_data" new
directory
> > created.
> > > > > > ftp> cd surapipith_data
> > > > > > 250 CWD command successful.
> > > > > > ftp> put wrfout_d02_2013-06-20_18:00:00
> > > > > > local: wrfout_d02_2013-06-20_18:00:00 remote:
> > > > > > wrfout_d02_2013-06-20_18:00:00
> > > > > > No control connection for command: Success
> > > > > > Passive mode refused.
> > > > > > ftp> bye
> > > > > >
> > > > > > Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00 |
more",
> > there
> > > > > were
> > > > > > many lines running up, and I could not find the first 4
line if
> it
> > > has
> > > > > GRIB
> > > > > > word or not. I will find the way to try again. Yet I hope
you may
> > see
> > > > the
> > > > > > WRF output that I upload to the ftp site of ucar. If not,
kindly
> > let
> > > me
> > > > > > know do it right, please.
> > > > > >
> > > > > > Thanks very much.
> > > > > > Best wishes,
> > > > > > Vanisa
> > > > > >
> > > > > > On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via RT
<
> > > > > > met_help at ucar.edu
> > > > > > > wrote:
> > > > > >
> > > > > > > Hi Julie,
> > > > > > >
> > > > > > > Thank you very much for your kind guidance. I am sorry
that I
> > just
> > > > read
> > > > > > > this e-mail, since there was no internet connection
along my
> > travel
> > > > > way,
> > > > > > > whole day, towards home in Thailand. I really appreciate
your
> > > effort
> > > > to
> > > > > > > help me before the holidays begin; I'm only now at home
in
> > Thailand
> > > > > with
> > > > > > > internet, but my Nepal workstation with WRF and MET
seems to
> down
> > > > now,
> > > > > > > after midnight of December 23rd in Kathmandu. I can ask
the IT
> > > admin
> > > > to
> > > > > > > check what happen, but he would only be able to resolve
this in
> > the
> > > > > next
> > > > > > 8
> > > > > > > hours (Christmas is actually not a holiday in Nepal, but
I'm
> > taking
> > > > > > > end-of-year leave rather).
> > > > > > >
> > > > > > > Never mind; let's celebrate the Season & give the best
to
> people
> > at
> > > > > home;
> > > > > > > I'll find out as you recommend in this while, and report
soon,
> > but
> > > > I'll
> > > > > > be
> > > > > > > okay to find you response later. Wish there a warm
Christmas
> Eve
> > &
> > > > > Happy
> > > > > > > coming New Year.
> > > > > > >
> > > > > > > Best wishes,
> > > > > > > Vanisa
> > > > > > >
> > > > > > > -----Original Message-----
> > > > > > > From: Julie Prestopnik via RT [mailto:met_help at ucar.edu]
> > > > > > > Sent: 22 December 2015 23:55
> > > > > > > To: vanisas at gmail.com
> > > > > > > Cc: Vanisa Surapipith
> > > > > > > Subject: Re: [rt.rap.ucar.edu #74467] point_stat for WRF
> output
> > > > > question
> > > > > > >
> > > > > > > Hi Vanisa.
> > > > > > >
> > > > > > > It is difficult to say without having access to the file
you
> are
> > > > using.
> > > > > > > Could you please upload your input forecast and
observation
> file
> > to
> > > > our
> > > > > > > ftp site using the directions on this page:
> > > > > > >
> > > > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > > > > >
> > > > > > > In an effort to try to help before having access to the
files,
> I
> > > will
> > > > > say
> > > > > > > that we have seen this problem before when the first
four
> > > characters
> > > > in
> > > > > > the
> > > > > > > file are not "GRIB". That is called the "magic cookie"
to
> > indicate
> > > > the
> > > > > > > file type.
> > > > > > >
> > > > > > > You could use the command "od -c
> wrfout_d02_2013-06-20_18:00:00
> > |
> > > > > more"
> > > > > > > to determine what the first four characters are. You
should
> see
> > > > > > something
> > > > > > > like this:
> > > > > > >
> > > > > > > 0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0
\0 001
> > P
> > > > ´
> > > > > > >
> > > > > > > If the first four characters are not "GRIB", you could
get
> around
> > > > this
> > > > > by
> > > > > > > changing the file's suffix to be ".grb" or by using the
> > "file_type
> > > =
> > > > > > > GRIB1;" option in the config file (see
data/config/README for
> > more
> > > > > > > information).
> > > > > > >
> > > > > > > I will be out of the office starting today at 2:45pm (US
> Mountain
> > > > time)
> > > > > > > through January 4th. I hope I can help you with this
problem
> > > before
> > > > I
> > > > > am
> > > > > > > out of the office.
> > > > > > >
> > > > > > > Julie
> > > > > > >
> > > > > > > On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via
RT <
> > > > > > > met_help at ucar.edu
> > > > > > > > wrote:
> > > > > > >
> > > > > > > >
> > > > > > > > Tue Dec 22 03:18:04 2015: Request 74467 was acted
upon.
> > > > > > > > Transaction: Ticket created by vanisas at gmail.com
> > > > > > > > Queue: met_help
> > > > > > > > Subject: point_stat for WRF output question
> > > > > > > > Owner: Nobody
> > > > > > > > Requestors: vanisas at gmail.com
> > > > > > > > Status: new
> > > > > > > > Ticket <URL:
> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > > > > > Hello there,
> > > > > > > >
> > > > > > > > I got point_stat run through for the tutorial case,
and would
> > > like
> > > > to
> > > > > > > > try on the different domain like Nepal, where you have
kindly
> > > > helped
> > > > > > > > creating the mask file. I edit the config file as
attached.
> The
> > > > error
> > > > > > > > shows up as
> > > > > > > > follow:
> > > > > > > >
> > > > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > > > > >
> > > > > > > >
> > > > > > > >
> > > > >
> >
/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-
> > > > > > > > 20_18:00:00
> > > > > > > > \
> > > > > > > > > tutorial/out/ascii2nc/nepal.nc \
> > > > > > > > > tutorial/config/PointStatConfig_Nepal \ -outdir
> > > > > > > > > tutorial/out/point_stat \ -v 2
> > > > > > > > DEBUG 1: Default Config File:
> > > > > > > >
> > > > >
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> > > > > > > > ult DEBUG 1: User Config File:
> > > > tutorial/config/PointStatConfig_Nepal
> > > > > > > > ERROR :
> > > > > > > > ERROR : Trouble reading forecast file
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > >
> > > > >
> > > >
> > >
> >
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> > > > > > > > ERROR :
> > > > > > > >
> > > > > > > > There seems to be some error about reading the WRF
output
> file.
> > > > Could
> > > > > > > > you give me some clue what kind of error that might
be, such
> > > that I
> > > > > > > > may correct, please. Is there something to set in the
config
> > > file?
> > > > I
> > > > > > > > am not sure what could be wrong with the WRF output,
which
> has
> > > > worked
> > > > > > > > okay with NCL/other tools.
> > > > > > > >
> > > > > > > > Best wishes,
> > > > > > > > Vanisa
> > > > > > > >
> > > > > > > >
> > > > > > >
> > > > > > >
> > > > > > > --
> > > > > > > Julie Prestopnik
> > > > > > > National Center for Atmospheric Research Research
Applications
> > > > > Laboratory
> > > > > > > Phone: 303.497.8399
> > > > > > > Email: jpresto at ucar.edu
> > > > > > >
> > >
> >
> >
>
>
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Jan 26 23:04:21 2016
Hi John,
apologies, I just notice that I did not remove the word "nepal" in the
bracket ([]) behind grid = within mask section in the PointConfig
file. I
have then removed that, and found a different error (also rename the
ascii
file containing the ADPUPA observation from our stations in Nepal to
be
called "nepal_ascii.nc" to avoid confusion) :
~/met-5.1_bugfix$ ./bin/point_stat \
> /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> tutorial/out/ascii2nc/nepal_ascii.nc \
> tutorial/config/PointStatConfig_Nepal \
> -outdir tutorial/out/point_stat \
> -v 2
DEBUG 1: Default Config File:
/home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
WARNING:
WARNING: parse_conf_interp() -> For neareast neighbor interpolation
method,
resetting method from "UW_MEAN" to "NEAREST" since the interpolation
width
is 1.
WARNING:
GSL_RNG_TYPE=mt19937
GSL_RNG_SEED=18446744072507402326
DEBUG 1: Forecast File:
/home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal_ascii.nc
DEBUG 2:
DEBUG 2:
--------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for SPFH/P500.
ERROR :
ERROR : parse_poly_mask() -> The masking and verification grids do
not
match:
ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll: -90.000 lon_ll:
-0.000
delta_lat: 0.500 delta_lon: 0.500 !=
ERROR : Projection: Mercator Nx: 270 Ny: 180 Lat_LL_radians: 0.4178
Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423 Lon_UR_radians: -1.5680
Mx:
1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
ERROR :
It looks like I need to set for some correct projection somewhere.
Could
you give a bit of guideline (in the while, I'll try to read the
Met_Users_Guide_v5.1.pdf in more detail; may send update again before
end
of the day here)...
Best wishes,
Vanisa
On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith
<vanisas at gmail.com>
wrote:
> Hi John,
>
> thank you for your reply; I recall your help earlier, and I am only
> gradually getting to know MET. I am grateful to the help as to what
Julie
> was responding patiently as well. I notice the (presumably) mask
files that
> Julie helped generating for Nepal, and I try them. However, the
similar
> error still shows up, and I am trying to understand why:
>
> ~/met-5.1_bugfix$ ./bin/point_stat \
> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > tutorial/out/ascii2nc/nepal.nc
tutorial/config/PointStatConfig_Nepal \
> > -outdir tutorial/out/point_stat \
> > -v 2
> DEBUG 1: Default Config File:
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> WARNING:
> WARNING: parse_conf_interp() -> For neareast neighbor interpolation
> method, resetting method from "UW_MEAN" to "NEAREST" since the
> interpolation width is 1.
> WARNING:
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=1682548822
> DEBUG 1: Forecast File:
> /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for SPFH/P500.
> ERROR :
> ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is not
> defined.
> ERROR :
>
> Kindly see in the surapipith_data folder in rap.ucar.edu site the
latest
> config file, as well as the suite of Nepal mask files that Julie
made and I
> use. Hope some progress would be before the end of the day, or once
some
> hints are next sent from there...
>
> Best wishes,
> Vanisa
>
>
> On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway via RT <
> met_help at ucar.edu> wrote:
>
>> Hello Vanisa,
>>
>> My name is John, and I work with Julie who is out of the office
today. I
>> see that you're having some trouble getting Point-Stat to run:
>>
>> ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is
not
>> defined.
>>
>> This should be pretty easy to resolve. You've worked with Julie to
>> generate a mask file for Nepal, and that mask is contained in some
NetCDF
>> file (the output of the gen_vx_mask tool). For example, let's say
it's in
>> "/home/vsurapipith/nepal_mask.nc". You'll need to put the path to
that
>> file in the "mask" section of the Point-Stat configuration file:
>>
>> mask = {
>> grid = [];
>> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; // or whatever
the
>> actual path is to that data file
>> sid = [];
>> }
>>
>> And then try running Point-Stat again. For a more complete
description of
>> the masking options available, please see pages 54-55 of the MET
User's
>> Guide:
>>
>>
>>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
>>
>> Please let me know how it goes.
>>
>> Thanks,
>> John
>>
>>
>>
>> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via RT <
>> met_help at ucar.edu
>> > wrote:
>>
>> >
>> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>> >
>> > Hi Julie,
>> >
>> > I'm back....after a couple of weeks spending on stumbling over
the
>> > installation of UPP until I got it through, now I have the GRIB1
files
>> > systematically converted from the WRF output generated here.
>> >
>> > I try the GRIB1 files with the PointStatConfig_Nepal that I sent
>> earlier to
>> > do Point_Stat with our these WRF output, but currently the error
shows
>> up
>> > as:
>> >
>> > ~/met-5.1_bugfix$ ./bin/point_stat
>> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
>> > tutorial/out/ascii2nc/nepal.nc
tutorial/config/PointStatConfig_Nepal
>> > -outdir tutorial/out/point_stat -v 2
>> > DEBUG 1: Default Config File:
>> >
>> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
>> > DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
>> > WARNING:
>> > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
>> method,
>> > resetting method from "UW_MEAN" to "NEAREST" since the
interpolation
>> width
>> > is 1.
>> > WARNING:
>> > GSL_RNG_TYPE=mt19937
>> > GSL_RNG_SEED=18446744072070604630
>> > DEBUG 1: Forecast File:
>> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
>> > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
>> > DEBUG 2:
>> > DEBUG 2:
>> >
>> >
>>
--------------------------------------------------------------------------------
>> > DEBUG 2:
>> > DEBUG 2: Reading data for SPFH/P500.
>> > ERROR :
>> > ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is
not
>> > defined.
>> > ERROR :
>> >
>> > It looks to me like the mask file - nepal.nc, that I deliberately
took
>> > from
>> > what you generated earlier - was not seen/defined well by MET
here. I
>> have
>> > actually just put the WRF output files in GRIB1 format (checked
with as
>> you
>> > recommend) on the rap.ucar.edu ftp site (surapipith_data)...Would
you
>> > kindly let me know what should be done to correct this, please?
>> >
>> > Best wishes,
>> > Vanisa
>> >
>> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via RT <
>> > met_help at ucar.edu>
>> > wrote:
>> >
>> > > Hi Vanisa.
>> > >
>> > > While I don't have any experience installing UPP myself, I did
take a
>> > look
>> > > at your log file where I saw:
>> > >
>> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No such file
or
>> > > directory compilation terminated.
>> > >
>> > > So, it looks like UPP is having trouble finding jasper.h.
>> > >
>> > > Support for UPP is provided via wrfhelp at ucar.edu. If you're
>> unable
>> > to
>> > > resolve the jasper problem, I would suggest writing to them
regarding
>> > this
>> > > issue.
>> > >
>> > > Thanks,
>> > > Julie
>> > >
>> > >
>> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via RT <
>> > > met_help at ucar.edu>
>> > > wrote:
>> > >
>> > > >
>> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
>
>> > > >
>> > > > Hi Julie,
>> > > >
>> > > > Happy New Year...:) Thank you very much for your reply upon
back to
>> > > office.
>> > > > Indeed, I did not use UPP for processing the WRF (ARW)
output.
>> > Therefore,
>> > > > today I look into how to install UPP on my workstation for
this
>> > purpose.
>> > > > However, there is a problem that unipost.exe was not
generated after
>> > > > compilation, while other files came through. I was wondered
if it
>> has
>> > to
>> > > do
>> > > > with liking to JASPER. I have not been able to spot why, and
now
>> attach
>> > > the
>> > > > compile_upp.log to you here, in case you might be able to see
what
>> is
>> > > wrong
>> > > > in there.
>> > > >
>> > > > Thanks again; update more soon.
>> > > >
>> > > > Best wishes,
>> > > > Vanisa
>> > > >
>> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT <
>> > > met_help at ucar.edu
>> > > > >
>> > > > wrote:
>> > > >
>> > > > > Hi Vanisa. The file you uploaded may have aged off of our
ftp
>> > site. I
>> > > > do
>> > > > > not see it there currently.
>> > > > >
>> > > > > I neglected to ask you if your WRF output file was post-
processed?
>> > MET
>> > > > is
>> > > > > unable to read the output of WRF directly. Instead, you
must
>> first
>> > > > > post-process it using either the Unified PostProcessor
(UPP) or
>> the
>> > > > pinterp
>> > > > > tool (for ARW only). We strongly recommend using UPP since
MET's
>> > > support
>> > > > > is most robust for the GRIB1 format that UPP creates.
>> > > > >
>> > > > > If you have already post-processed your file, please trying
>> uploading
>> > > to
>> > > > > our ftp site again. If you have not, please post-process
the file
>> > and
>> > > > give
>> > > > > it another try. Thanks!
>> > > > >
>> > > > > Happy New Year!
>> > > > >
>> > > > > Julie
>> > > > >
>> > > > > On Sun, Dec 27, 2015 at 11:37 PM, vanisas at gmail.com via RT
<
>> > > > > met_help at ucar.edu> wrote:
>> > > > >
>> > > > > >
>> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>> > > > > >
>> > > > > > Dear Julie & all there,
>> > > > > >
>> > > > > > I just manage to try as Julie recommended right before
the
>> > Christmas
>> > > > > > holidays, i.e. to put the input file I use on the ftp
site. I am
>> > not
>> > > > sure
>> > > > > > if it works, or not. What shows up say:
>> > > > > >
>> > > > > > 220 brew.rap.ucar.edu FTP server ready.
>> > > > > > Name (ftp.rap.ucar.edu:vsurapipith): anonymous
>> > > > > > 331 Guest login ok, send your complete e-mail address as
>> password.
>> > > > > > Password:
>> > > > > > 230 Guest login ok, access restrictions apply.
>> > > > > > Remote system type is UNIX.
>> > > > > > Using binary mode to transfer files.
>> > > > > > ftp> cd incoming/irap/met_help
>> > > > > > 250 CWD command successful.
>> > > > > > ftp> mkdir surapipith_data
>> > > > > > 257 "/incoming/irap/met_help/surapipith_data" new
directory
>> > created.
>> > > > > > ftp> cd surapipith_data
>> > > > > > 250 CWD command successful.
>> > > > > > ftp> put wrfout_d02_2013-06-20_18:00:00
>> > > > > > local: wrfout_d02_2013-06-20_18:00:00 remote:
>> > > > > > wrfout_d02_2013-06-20_18:00:00
>> > > > > > No control connection for command: Success
>> > > > > > Passive mode refused.
>> > > > > > ftp> bye
>> > > > > >
>> > > > > > Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00 |
more",
>> > there
>> > > > > were
>> > > > > > many lines running up, and I could not find the first 4
line if
>> it
>> > > has
>> > > > > GRIB
>> > > > > > word or not. I will find the way to try again. Yet I hope
you
>> may
>> > see
>> > > > the
>> > > > > > WRF output that I upload to the ftp site of ucar. If not,
kindly
>> > let
>> > > me
>> > > > > > know do it right, please.
>> > > > > >
>> > > > > > Thanks very much.
>> > > > > > Best wishes,
>> > > > > > Vanisa
>> > > > > >
>> > > > > > On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via RT
<
>> > > > > > met_help at ucar.edu
>> > > > > > > wrote:
>> > > > > >
>> > > > > > > Hi Julie,
>> > > > > > >
>> > > > > > > Thank you very much for your kind guidance. I am sorry
that I
>> > just
>> > > > read
>> > > > > > > this e-mail, since there was no internet connection
along my
>> > travel
>> > > > > way,
>> > > > > > > whole day, towards home in Thailand. I really
appreciate your
>> > > effort
>> > > > to
>> > > > > > > help me before the holidays begin; I'm only now at home
in
>> > Thailand
>> > > > > with
>> > > > > > > internet, but my Nepal workstation with WRF and MET
seems to
>> down
>> > > > now,
>> > > > > > > after midnight of December 23rd in Kathmandu. I can ask
the IT
>> > > admin
>> > > > to
>> > > > > > > check what happen, but he would only be able to resolve
this
>> in
>> > the
>> > > > > next
>> > > > > > 8
>> > > > > > > hours (Christmas is actually not a holiday in Nepal,
but I'm
>> > taking
>> > > > > > > end-of-year leave rather).
>> > > > > > >
>> > > > > > > Never mind; let's celebrate the Season & give the best
to
>> people
>> > at
>> > > > > home;
>> > > > > > > I'll find out as you recommend in this while, and
report soon,
>> > but
>> > > > I'll
>> > > > > > be
>> > > > > > > okay to find you response later. Wish there a warm
Christmas
>> Eve
>> > &
>> > > > > Happy
>> > > > > > > coming New Year.
>> > > > > > >
>> > > > > > > Best wishes,
>> > > > > > > Vanisa
>> > > > > > >
>> > > > > > > -----Original Message-----
>> > > > > > > From: Julie Prestopnik via RT
[mailto:met_help at ucar.edu]
>> > > > > > > Sent: 22 December 2015 23:55
>> > > > > > > To: vanisas at gmail.com
>> > > > > > > Cc: Vanisa Surapipith
>> > > > > > > Subject: Re: [rt.rap.ucar.edu #74467] point_stat for
WRF
>> output
>> > > > > question
>> > > > > > >
>> > > > > > > Hi Vanisa.
>> > > > > > >
>> > > > > > > It is difficult to say without having access to the
file you
>> are
>> > > > using.
>> > > > > > > Could you please upload your input forecast and
observation
>> file
>> > to
>> > > > our
>> > > > > > > ftp site using the directions on this page:
>> > > > > > >
>> > > > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
>> > > > > > >
>> > > > > > > In an effort to try to help before having access to the
>> files, I
>> > > will
>> > > > > say
>> > > > > > > that we have seen this problem before when the first
four
>> > > characters
>> > > > in
>> > > > > > the
>> > > > > > > file are not "GRIB". That is called the "magic cookie"
to
>> > indicate
>> > > > the
>> > > > > > > file type.
>> > > > > > >
>> > > > > > > You could use the command "od -c
>> wrfout_d02_2013-06-20_18:00:00
>> > |
>> > > > > more"
>> > > > > > > to determine what the first four characters are. You
should
>> see
>> > > > > > something
>> > > > > > > like this:
>> > > > > > >
>> > > > > > > 0000000 G R I B \0 \0 \0 002 \0 \0 \0 \0
\0
>> 001
>> > P
>> > > > ´
>> > > > > > >
>> > > > > > > If the first four characters are not "GRIB", you could
get
>> around
>> > > > this
>> > > > > by
>> > > > > > > changing the file's suffix to be ".grb" or by using the
>> > "file_type
>> > > =
>> > > > > > > GRIB1;" option in the config file (see
data/config/README for
>> > more
>> > > > > > > information).
>> > > > > > >
>> > > > > > > I will be out of the office starting today at 2:45pm
(US
>> Mountain
>> > > > time)
>> > > > > > > through January 4th. I hope I can help you with this
problem
>> > > before
>> > > > I
>> > > > > am
>> > > > > > > out of the office.
>> > > > > > >
>> > > > > > > Julie
>> > > > > > >
>> > > > > > > On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com via
RT <
>> > > > > > > met_help at ucar.edu
>> > > > > > > > wrote:
>> > > > > > >
>> > > > > > > >
>> > > > > > > > Tue Dec 22 03:18:04 2015: Request 74467 was acted
upon.
>> > > > > > > > Transaction: Ticket created by vanisas at gmail.com
>> > > > > > > > Queue: met_help
>> > > > > > > > Subject: point_stat for WRF output question
>> > > > > > > > Owner: Nobody
>> > > > > > > > Requestors: vanisas at gmail.com
>> > > > > > > > Status: new
>> > > > > > > > Ticket <URL:
>> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
>> > > > > > > > >
>> > > > > > > >
>> > > > > > > >
>> > > > > > > > Hello there,
>> > > > > > > >
>> > > > > > > > I got point_stat run through for the tutorial case,
and
>> would
>> > > like
>> > > > to
>> > > > > > > > try on the different domain like Nepal, where you
have
>> kindly
>> > > > helped
>> > > > > > > > creating the mask file. I edit the config file as
attached.
>> The
>> > > > error
>> > > > > > > > shows up as
>> > > > > > > > follow:
>> > > > > > > >
>> > > > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
>> > > > > > > > >
>> > > > > > > >
>> > > > > > > >
>> > > > >
>> >
/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-
>> > > > > > > > 20_18:00:00
>> > > > > > > > \
>> > > > > > > > > tutorial/out/ascii2nc/nepal.nc \
>> > > > > > > > > tutorial/config/PointStatConfig_Nepal \ -outdir
>> > > > > > > > > tutorial/out/point_stat \ -v 2
>> > > > > > > > DEBUG 1: Default Config File:
>> > > > > > > >
>> > > > >
>> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
>> > > > > > > > ult DEBUG 1: User Config File:
>> > > > tutorial/config/PointStatConfig_Nepal
>> > > > > > > > ERROR :
>> > > > > > > > ERROR : Trouble reading forecast file
>> > > > > > > >
>> > > > > > > >
>> > > > > > >
>> > > > > >
>> > > > >
>> > > >
>> > >
>> >
>>
"/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
>> > > > > > > > ERROR :
>> > > > > > > >
>> > > > > > > > There seems to be some error about reading the WRF
output
>> file.
>> > > > Could
>> > > > > > > > you give me some clue what kind of error that might
be, such
>> > > that I
>> > > > > > > > may correct, please. Is there something to set in the
config
>> > > file?
>> > > > I
>> > > > > > > > am not sure what could be wrong with the WRF output,
which
>> has
>> > > > worked
>> > > > > > > > okay with NCL/other tools.
>> > > > > > > >
>> > > > > > > > Best wishes,
>> > > > > > > > Vanisa
>> > > > > > > >
>> > > > > > > >
>> > > > > > >
>> > > > > > >
>> > > > > > > --
>> > > > > > > Julie Prestopnik
>> > > > > > > National Center for Atmospheric Research Research
Applications
>> > > > > Laboratory
>> > > > > > > Phone: 303.497.8399
>> > > > > > > Email: jpresto at ucar.edu
>> > > > > > >
>> > >
>> >
>> >
>>
>>
>
>
>
>
--
"Our ignorance is not so vast as our failure to use what we know."
M. King Hubbert
October 5th, 1903 -- October 11th, 1989
------------------------------------------------
Subject: point_stat for WRF output question
From: Julie Prestopnik
Time: Wed Jan 27 11:15:39 2016
Hi Vanisa. This is Julie again.
I see that you are encountering a problem due to the masking and
verification grids not matching.
One projection is a Lat/Lon grid and the other is a Mercator.
The mask file I created for you using gen_vx_mask used a Lat/Lon grid
file
(gfs_2012040900_F003.grib). If you would like to get your netcdf mask
file
on a Mercator grid, you could run gen_vx_mask using your file that
contains
the grid that you would like to verify, which I believe is on a
Mercator
projection (WRFPRS_d02.00?) and the nepal poly file which will then
produce
a netcdf file on that same projection.
So, for example, something like:
gen_vx_mask WRFPRS_d02.00 nepal_detail.poly nepal_detail.nc
Then, you can rerun point_stat using the new netcdf mask file on the
Mercator grid and that should fix the problem.
However, if you're not going to be running thousands of cases of Point
Stat, you don't necessarily need to run gen_vx_mask, but rather you
could
use the poly file directly with Point Stat in the configuration file
in the
mask section that John described earlier. Using the poly file
directly can
be slow, so if you're running point_stat on a large number of files,
this
is likely not the best way to use point_stat.
I hope this helps. Please let me know if you have any questions.
Thanks,
Julie
On Tue, Jan 26, 2016 at 11:04 PM, vanisas at gmail.com via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>
> Hi John,
>
> apologies, I just notice that I did not remove the word "nepal" in
the
> bracket ([]) behind grid = within mask section in the PointConfig
file. I
> have then removed that, and found a different error (also rename the
ascii
> file containing the ADPUPA observation from our stations in Nepal to
be
> called "nepal_ascii.nc" to avoid confusion) :
>
> ~/met-5.1_bugfix$ ./bin/point_stat \
> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > tutorial/out/ascii2nc/nepal_ascii.nc \
> > tutorial/config/PointStatConfig_Nepal \
> > -outdir tutorial/out/point_stat \
> > -v 2
> DEBUG 1: Default Config File:
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> WARNING:
> WARNING: parse_conf_interp() -> For neareast neighbor interpolation
method,
> resetting method from "UW_MEAN" to "NEAREST" since the interpolation
width
> is 1.
> WARNING:
> GSL_RNG_TYPE=mt19937
> GSL_RNG_SEED=18446744072507402326
> DEBUG 1: Forecast File:
> /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal_ascii.nc
> DEBUG 2:
> DEBUG 2:
>
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for SPFH/P500.
> ERROR :
> ERROR : parse_poly_mask() -> The masking and verification grids do
not
> match:
> ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll: -90.000 lon_ll:
-0.000
> delta_lat: 0.500 delta_lon: 0.500 !=
> ERROR : Projection: Mercator Nx: 270 Ny: 180 Lat_LL_radians: 0.4178
> Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423 Lon_UR_radians:
-1.5680 Mx:
> 1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
> ERROR :
>
> It looks like I need to set for some correct projection somewhere.
Could
> you give a bit of guideline (in the while, I'll try to read the
> Met_Users_Guide_v5.1.pdf in more detail; may send update again
before end
> of the day here)...
>
> Best wishes,
> Vanisa
>
> On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith
<vanisas at gmail.com>
> wrote:
>
> > Hi John,
> >
> > thank you for your reply; I recall your help earlier, and I am
only
> > gradually getting to know MET. I am grateful to the help as to
what Julie
> > was responding patiently as well. I notice the (presumably) mask
files
> that
> > Julie helped generating for Nepal, and I try them. However, the
similar
> > error still shows up, and I am trying to understand why:
> >
> > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > > tutorial/out/ascii2nc/nepal.nc
tutorial/config/PointStatConfig_Nepal \
> > > -outdir tutorial/out/point_stat \
> > > -v 2
> > DEBUG 1: Default Config File:
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> > WARNING:
> > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> > method, resetting method from "UW_MEAN" to "NEAREST" since the
> > interpolation width is 1.
> > WARNING:
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=1682548822
> > DEBUG 1: Forecast File:
> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > DEBUG 2:
> > DEBUG 2:
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for SPFH/P500.
> > ERROR :
> > ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is
not
> > defined.
> > ERROR :
> >
> > Kindly see in the surapipith_data folder in rap.ucar.edu site the
latest
> > config file, as well as the suite of Nepal mask files that Julie
made
> and I
> > use. Hope some progress would be before the end of the day, or
once some
> > hints are next sent from there...
> >
> > Best wishes,
> > Vanisa
> >
> >
> > On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway via RT <
> > met_help at ucar.edu> wrote:
> >
> >> Hello Vanisa,
> >>
> >> My name is John, and I work with Julie who is out of the office
today.
> I
> >> see that you're having some trouble getting Point-Stat to run:
> >>
> >> ERROR : parse_grid_mask() -> the mask_grid requested "nepal"
is not
> >> defined.
> >>
> >> This should be pretty easy to resolve. You've worked with Julie
to
> >> generate a mask file for Nepal, and that mask is contained in
some
> NetCDF
> >> file (the output of the gen_vx_mask tool). For example, let's
say it's
> in
> >> "/home/vsurapipith/nepal_mask.nc". You'll need to put the path
to that
> >> file in the "mask" section of the Point-Stat configuration file:
> >>
> >> mask = {
> >> grid = [];
> >> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; // or
whatever the
> >> actual path is to that data file
> >> sid = [];
> >> }
> >>
> >> And then try running Point-Stat again. For a more complete
description
> of
> >> the masking options available, please see pages 54-55 of the MET
User's
> >> Guide:
> >>
> >>
> >>
>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
> >>
> >> Please let me know how it goes.
> >>
> >> Thanks,
> >> John
> >>
> >>
> >>
> >> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via RT <
> >> met_help at ucar.edu
> >> > wrote:
> >>
> >> >
> >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >> >
> >> > Hi Julie,
> >> >
> >> > I'm back....after a couple of weeks spending on stumbling over
the
> >> > installation of UPP until I got it through, now I have the
GRIB1 files
> >> > systematically converted from the WRF output generated here.
> >> >
> >> > I try the GRIB1 files with the PointStatConfig_Nepal that I
sent
> >> earlier to
> >> > do Point_Stat with our these WRF output, but currently the
error shows
> >> up
> >> > as:
> >> >
> >> > ~/met-5.1_bugfix$ ./bin/point_stat
> >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> >> > tutorial/out/ascii2nc/nepal.nc
tutorial/config/PointStatConfig_Nepal
> >> > -outdir tutorial/out/point_stat -v 2
> >> > DEBUG 1: Default Config File:
> >> >
> >>
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> >> > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> >> > WARNING:
> >> > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> >> method,
> >> > resetting method from "UW_MEAN" to "NEAREST" since the
interpolation
> >> width
> >> > is 1.
> >> > WARNING:
> >> > GSL_RNG_TYPE=mt19937
> >> > GSL_RNG_SEED=18446744072070604630
> >> > DEBUG 1: Forecast File:
> >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> >> > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> >> > DEBUG 2:
> >> > DEBUG 2:
> >> >
> >> >
> >>
>
--------------------------------------------------------------------------------
> >> > DEBUG 2:
> >> > DEBUG 2: Reading data for SPFH/P500.
> >> > ERROR :
> >> > ERROR : parse_grid_mask() -> the mask_grid requested "nepal"
is not
> >> > defined.
> >> > ERROR :
> >> >
> >> > It looks to me like the mask file - nepal.nc, that I
deliberately
> took
> >> > from
> >> > what you generated earlier - was not seen/defined well by MET
here. I
> >> have
> >> > actually just put the WRF output files in GRIB1 format (checked
with
> as
> >> you
> >> > recommend) on the rap.ucar.edu ftp site
(surapipith_data)...Would you
> >> > kindly let me know what should be done to correct this, please?
> >> >
> >> > Best wishes,
> >> > Vanisa
> >> >
> >> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via RT <
> >> > met_help at ucar.edu>
> >> > wrote:
> >> >
> >> > > Hi Vanisa.
> >> > >
> >> > > While I don't have any experience installing UPP myself, I
did take
> a
> >> > look
> >> > > at your log file where I saw:
> >> > >
> >> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No such
file or
> >> > > directory compilation terminated.
> >> > >
> >> > > So, it looks like UPP is having trouble finding jasper.h.
> >> > >
> >> > > Support for UPP is provided via wrfhelp at ucar.edu. If
you're
> >> unable
> >> > to
> >> > > resolve the jasper problem, I would suggest writing to them
> regarding
> >> > this
> >> > > issue.
> >> > >
> >> > > Thanks,
> >> > > Julie
> >> > >
> >> > >
> >> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via RT <
> >> > > met_help at ucar.edu>
> >> > > wrote:
> >> > >
> >> > > >
> >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >> > > >
> >> > > > Hi Julie,
> >> > > >
> >> > > > Happy New Year...:) Thank you very much for your reply upon
back
> to
> >> > > office.
> >> > > > Indeed, I did not use UPP for processing the WRF (ARW)
output.
> >> > Therefore,
> >> > > > today I look into how to install UPP on my workstation for
this
> >> > purpose.
> >> > > > However, there is a problem that unipost.exe was not
generated
> after
> >> > > > compilation, while other files came through. I was wondered
if it
> >> has
> >> > to
> >> > > do
> >> > > > with liking to JASPER. I have not been able to spot why,
and now
> >> attach
> >> > > the
> >> > > > compile_upp.log to you here, in case you might be able to
see what
> >> is
> >> > > wrong
> >> > > > in there.
> >> > > >
> >> > > > Thanks again; update more soon.
> >> > > >
> >> > > > Best wishes,
> >> > > > Vanisa
> >> > > >
> >> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT <
> >> > > met_help at ucar.edu
> >> > > > >
> >> > > > wrote:
> >> > > >
> >> > > > > Hi Vanisa. The file you uploaded may have aged off of
our ftp
> >> > site. I
> >> > > > do
> >> > > > > not see it there currently.
> >> > > > >
> >> > > > > I neglected to ask you if your WRF output file was
> post-processed?
> >> > MET
> >> > > > is
> >> > > > > unable to read the output of WRF directly. Instead, you
must
> >> first
> >> > > > > post-process it using either the Unified PostProcessor
(UPP) or
> >> the
> >> > > > pinterp
> >> > > > > tool (for ARW only). We strongly recommend using UPP
since
> MET's
> >> > > support
> >> > > > > is most robust for the GRIB1 format that UPP creates.
> >> > > > >
> >> > > > > If you have already post-processed your file, please
trying
> >> uploading
> >> > > to
> >> > > > > our ftp site again. If you have not, please post-process
the
> file
> >> > and
> >> > > > give
> >> > > > > it another try. Thanks!
> >> > > > >
> >> > > > > Happy New Year!
> >> > > > >
> >> > > > > Julie
> >> > > > >
> >> > > > > On Sun, Dec 27, 2015 at 11:37 PM, vanisas at gmail.com via
RT <
> >> > > > > met_help at ucar.edu> wrote:
> >> > > > >
> >> > > > > >
> >> > > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> >
> >> > > > > >
> >> > > > > > Dear Julie & all there,
> >> > > > > >
> >> > > > > > I just manage to try as Julie recommended right before
the
> >> > Christmas
> >> > > > > > holidays, i.e. to put the input file I use on the ftp
site. I
> am
> >> > not
> >> > > > sure
> >> > > > > > if it works, or not. What shows up say:
> >> > > > > >
> >> > > > > > 220 brew.rap.ucar.edu FTP server ready.
> >> > > > > > Name (ftp.rap.ucar.edu:vsurapipith): anonymous
> >> > > > > > 331 Guest login ok, send your complete e-mail address
as
> >> password.
> >> > > > > > Password:
> >> > > > > > 230 Guest login ok, access restrictions apply.
> >> > > > > > Remote system type is UNIX.
> >> > > > > > Using binary mode to transfer files.
> >> > > > > > ftp> cd incoming/irap/met_help
> >> > > > > > 250 CWD command successful.
> >> > > > > > ftp> mkdir surapipith_data
> >> > > > > > 257 "/incoming/irap/met_help/surapipith_data" new
directory
> >> > created.
> >> > > > > > ftp> cd surapipith_data
> >> > > > > > 250 CWD command successful.
> >> > > > > > ftp> put wrfout_d02_2013-06-20_18:00:00
> >> > > > > > local: wrfout_d02_2013-06-20_18:00:00 remote:
> >> > > > > > wrfout_d02_2013-06-20_18:00:00
> >> > > > > > No control connection for command: Success
> >> > > > > > Passive mode refused.
> >> > > > > > ftp> bye
> >> > > > > >
> >> > > > > > Also while I try "oc -d wrfout_d02_2013-06-20_18:00:00
|
> more",
> >> > there
> >> > > > > were
> >> > > > > > many lines running up, and I could not find the first 4
line
> if
> >> it
> >> > > has
> >> > > > > GRIB
> >> > > > > > word or not. I will find the way to try again. Yet I
hope you
> >> may
> >> > see
> >> > > > the
> >> > > > > > WRF output that I upload to the ftp site of ucar. If
not,
> kindly
> >> > let
> >> > > me
> >> > > > > > know do it right, please.
> >> > > > > >
> >> > > > > > Thanks very much.
> >> > > > > > Best wishes,
> >> > > > > > Vanisa
> >> > > > > >
> >> > > > > > On Thu, Dec 24, 2015 at 1:37 AM, Vanisa Surapipith via
RT <
> >> > > > > > met_help at ucar.edu
> >> > > > > > > wrote:
> >> > > > > >
> >> > > > > > > Hi Julie,
> >> > > > > > >
> >> > > > > > > Thank you very much for your kind guidance. I am
sorry that
> I
> >> > just
> >> > > > read
> >> > > > > > > this e-mail, since there was no internet connection
along my
> >> > travel
> >> > > > > way,
> >> > > > > > > whole day, towards home in Thailand. I really
appreciate
> your
> >> > > effort
> >> > > > to
> >> > > > > > > help me before the holidays begin; I'm only now at
home in
> >> > Thailand
> >> > > > > with
> >> > > > > > > internet, but my Nepal workstation with WRF and MET
seems to
> >> down
> >> > > > now,
> >> > > > > > > after midnight of December 23rd in Kathmandu. I can
ask the
> IT
> >> > > admin
> >> > > > to
> >> > > > > > > check what happen, but he would only be able to
resolve this
> >> in
> >> > the
> >> > > > > next
> >> > > > > > 8
> >> > > > > > > hours (Christmas is actually not a holiday in Nepal,
but I'm
> >> > taking
> >> > > > > > > end-of-year leave rather).
> >> > > > > > >
> >> > > > > > > Never mind; let's celebrate the Season & give the
best to
> >> people
> >> > at
> >> > > > > home;
> >> > > > > > > I'll find out as you recommend in this while, and
report
> soon,
> >> > but
> >> > > > I'll
> >> > > > > > be
> >> > > > > > > okay to find you response later. Wish there a warm
Christmas
> >> Eve
> >> > &
> >> > > > > Happy
> >> > > > > > > coming New Year.
> >> > > > > > >
> >> > > > > > > Best wishes,
> >> > > > > > > Vanisa
> >> > > > > > >
> >> > > > > > > -----Original Message-----
> >> > > > > > > From: Julie Prestopnik via RT
[mailto:met_help at ucar.edu]
> >> > > > > > > Sent: 22 December 2015 23:55
> >> > > > > > > To: vanisas at gmail.com
> >> > > > > > > Cc: Vanisa Surapipith
> >> > > > > > > Subject: Re: [rt.rap.ucar.edu #74467] point_stat for
WRF
> >> output
> >> > > > > question
> >> > > > > > >
> >> > > > > > > Hi Vanisa.
> >> > > > > > >
> >> > > > > > > It is difficult to say without having access to the
file you
> >> are
> >> > > > using.
> >> > > > > > > Could you please upload your input forecast and
observation
> >> file
> >> > to
> >> > > > our
> >> > > > > > > ftp site using the directions on this page:
> >> > > > > > >
> >> > > > > > >
http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >> > > > > > >
> >> > > > > > > In an effort to try to help before having access to
the
> >> files, I
> >> > > will
> >> > > > > say
> >> > > > > > > that we have seen this problem before when the first
four
> >> > > characters
> >> > > > in
> >> > > > > > the
> >> > > > > > > file are not "GRIB". That is called the "magic
cookie" to
> >> > indicate
> >> > > > the
> >> > > > > > > file type.
> >> > > > > > >
> >> > > > > > > You could use the command "od -c
> >> wrfout_d02_2013-06-20_18:00:00
> >> > |
> >> > > > > more"
> >> > > > > > > to determine what the first four characters are. You
should
> >> see
> >> > > > > > something
> >> > > > > > > like this:
> >> > > > > > >
> >> > > > > > > 0000000 G R I B \0 \0 \0 002 \0 \0 \0
\0 \0
> >> 001
> >> > P
> >> > > > ´
> >> > > > > > >
> >> > > > > > > If the first four characters are not "GRIB", you
could get
> >> around
> >> > > > this
> >> > > > > by
> >> > > > > > > changing the file's suffix to be ".grb" or by using
the
> >> > "file_type
> >> > > =
> >> > > > > > > GRIB1;" option in the config file (see
data/config/README
> for
> >> > more
> >> > > > > > > information).
> >> > > > > > >
> >> > > > > > > I will be out of the office starting today at 2:45pm
(US
> >> Mountain
> >> > > > time)
> >> > > > > > > through January 4th. I hope I can help you with this
> problem
> >> > > before
> >> > > > I
> >> > > > > am
> >> > > > > > > out of the office.
> >> > > > > > >
> >> > > > > > > Julie
> >> > > > > > >
> >> > > > > > > On Tue, Dec 22, 2015 at 3:18 AM, vanisas at gmail.com
via RT <
> >> > > > > > > met_help at ucar.edu
> >> > > > > > > > wrote:
> >> > > > > > >
> >> > > > > > > >
> >> > > > > > > > Tue Dec 22 03:18:04 2015: Request 74467 was acted
upon.
> >> > > > > > > > Transaction: Ticket created by vanisas at gmail.com
> >> > > > > > > > Queue: met_help
> >> > > > > > > > Subject: point_stat for WRF output question
> >> > > > > > > > Owner: Nobody
> >> > > > > > > > Requestors: vanisas at gmail.com
> >> > > > > > > > Status: new
> >> > > > > > > > Ticket <URL:
> >> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> >> > > > > > > > >
> >> > > > > > > >
> >> > > > > > > >
> >> > > > > > > > Hello there,
> >> > > > > > > >
> >> > > > > > > > I got point_stat run through for the tutorial case,
and
> >> would
> >> > > like
> >> > > > to
> >> > > > > > > > try on the different domain like Nepal, where you
have
> >> kindly
> >> > > > helped
> >> > > > > > > > creating the mask file. I edit the config file as
> attached.
> >> The
> >> > > > error
> >> > > > > > > > shows up as
> >> > > > > > > > follow:
> >> > > > > > > >
> >> > > > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> >> > > > > > > > >
> >> > > > > > > >
> >> > > > > > > >
> >> > > > >
> >> >
/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-06-
> >> > > > > > > > 20_18:00:00
> >> > > > > > > > \
> >> > > > > > > > > tutorial/out/ascii2nc/nepal.nc \
> >> > > > > > > > > tutorial/config/PointStatConfig_Nepal \ -outdir
> >> > > > > > > > > tutorial/out/point_stat \ -v 2
> >> > > > > > > > DEBUG 1: Default Config File:
> >> > > > > > > >
> >> > > > >
> >> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_defa
> >> > > > > > > > ult DEBUG 1: User Config File:
> >> > > > tutorial/config/PointStatConfig_Nepal
> >> > > > > > > > ERROR :
> >> > > > > > > > ERROR : Trouble reading forecast file
> >> > > > > > > >
> >> > > > > > > >
> >> > > > > > >
> >> > > > > >
> >> > > > >
> >> > > >
> >> > >
> >> >
> >>
> "/home/vsurapipith/AMET/model_data/MET/wrfoutSusKAT/wrfout_d02_2013-
06-20_18:00:00"
> >> > > > > > > > ERROR :
> >> > > > > > > >
> >> > > > > > > > There seems to be some error about reading the WRF
output
> >> file.
> >> > > > Could
> >> > > > > > > > you give me some clue what kind of error that might
be,
> such
> >> > > that I
> >> > > > > > > > may correct, please. Is there something to set in
the
> config
> >> > > file?
> >> > > > I
> >> > > > > > > > am not sure what could be wrong with the WRF
output, which
> >> has
> >> > > > worked
> >> > > > > > > > okay with NCL/other tools.
> >> > > > > > > >
> >> > > > > > > > Best wishes,
> >> > > > > > > > Vanisa
> >> > > > > > > >
> >> > > > > > > >
> >> > > > > > >
> >> > > > > > >
> >> > > > > > > --
> >> > > > > > > Julie Prestopnik
> >> > > > > > > National Center for Atmospheric Research Research
> Applications
> >> > > > > Laboratory
> >> > > > > > > Phone: 303.497.8399
> >> > > > > > > Email: jpresto at ucar.edu
> >> > > > > > >
> >> > >
> >> >
> >> >
> >>
> >>
> >
> >
> >
> >
>
>
> --
> "Our ignorance is not so vast as our failure to use what we know."
> M. King Hubbert
> October 5th, 1903 -- October 11th, 1989
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Feb 02 02:32:01 2016
Hi Julie,
thank you for your reply last week; apologies for my delay in
reporting.
What you suggest works. I can now generate the output files using
Point_Stat. However, I do not yet understand how to interpret or make
use
of the output...I have been trying to read what described in the
MET_User_Guide. Frankly speaking, I still do not get it. Could you
kindly
see in my folder in the rap.ucar.edu ftp site in
incoming/irap/met_help
under surapipith_data. I attached here some of the files. It looks
like
only the headers appear...or what I should do to generate the
statistical
evaluation output. What surprise me was that the example files I got
from
the US example (follow the online tutorial) look similar. Could you
guide
me where to read to comprehend MET?
One other question I have is that probably this is because I use only
one
hour of WRF output for this try. Do I need to generate a GRIB1 file
that
combine all the required hours...otherwise, can MET read in more than
one
GRIB1 files at a time for Point_Stat analysis?
Best wishes,
Vanisa
On Thu, Jan 28, 2016 at 12:00 AM, Julie Prestopnik via RT
<met_help at ucar.edu
> wrote:
> Hi Vanisa. This is Julie again.
>
> I see that you are encountering a problem due to the masking and
> verification grids not matching.
>
> One projection is a Lat/Lon grid and the other is a Mercator.
>
> The mask file I created for you using gen_vx_mask used a Lat/Lon
grid file
> (gfs_2012040900_F003.grib). If you would like to get your netcdf
mask
> file
> on a Mercator grid, you could run gen_vx_mask using your file that
contains
> the grid that you would like to verify, which I believe is on a
Mercator
> projection (WRFPRS_d02.00?) and the nepal poly file which will then
produce
> a netcdf file on that same projection.
>
> So, for example, something like:
>
> gen_vx_mask WRFPRS_d02.00 nepal_detail.poly nepal_detail.nc
>
> Then, you can rerun point_stat using the new netcdf mask file on
the
> Mercator grid and that should fix the problem.
>
> However, if you're not going to be running thousands of cases of
Point
> Stat, you don't necessarily need to run gen_vx_mask, but rather you
could
> use the poly file directly with Point Stat in the configuration file
in the
> mask section that John described earlier. Using the poly file
directly can
> be slow, so if you're running point_stat on a large number of files,
this
> is likely not the best way to use point_stat.
>
> I hope this helps. Please let me know if you have any questions.
>
> Thanks,
> Julie
>
> On Tue, Jan 26, 2016 at 11:04 PM, vanisas at gmail.com via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >
> > Hi John,
> >
> > apologies, I just notice that I did not remove the word "nepal" in
the
> > bracket ([]) behind grid = within mask section in the PointConfig
file. I
> > have then removed that, and found a different error (also rename
the
> ascii
> > file containing the ADPUPA observation from our stations in Nepal
to be
> > called "nepal_ascii.nc" to avoid confusion) :
> >
> > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > > tutorial/out/ascii2nc/nepal_ascii.nc \
> > > tutorial/config/PointStatConfig_Nepal \
> > > -outdir tutorial/out/point_stat \
> > > -v 2
> > DEBUG 1: Default Config File:
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> > WARNING:
> > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> method,
> > resetting method from "UW_MEAN" to "NEAREST" since the
interpolation
> width
> > is 1.
> > WARNING:
> > GSL_RNG_TYPE=mt19937
> > GSL_RNG_SEED=18446744072507402326
> > DEBUG 1: Forecast File:
> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal_ascii.nc
> > DEBUG 2:
> > DEBUG 2:
> >
> >
>
--------------------------------------------------------------------------------
> > DEBUG 2:
> > DEBUG 2: Reading data for SPFH/P500.
> > ERROR :
> > ERROR : parse_poly_mask() -> The masking and verification grids
do not
> > match:
> > ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll: -90.000
lon_ll:
> -0.000
> > delta_lat: 0.500 delta_lon: 0.500 !=
> > ERROR : Projection: Mercator Nx: 270 Ny: 180 Lat_LL_radians:
0.4178
> > Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423 Lon_UR_radians:
-1.5680
> Mx:
> > 1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
> > ERROR :
> >
> > It looks like I need to set for some correct projection somewhere.
Could
> > you give a bit of guideline (in the while, I'll try to read the
> > Met_Users_Guide_v5.1.pdf in more detail; may send update again
before end
> > of the day here)...
> >
> > Best wishes,
> > Vanisa
> >
> > On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith
<vanisas at gmail.com>
> > wrote:
> >
> > > Hi John,
> > >
> > > thank you for your reply; I recall your help earlier, and I am
only
> > > gradually getting to know MET. I am grateful to the help as to
what
> Julie
> > > was responding patiently as well. I notice the (presumably) mask
files
> > that
> > > Julie helped generating for Nepal, and I try them. However, the
similar
> > > error still shows up, and I am trying to understand why:
> > >
> > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > > > tutorial/out/ascii2nc/nepal.nc
> tutorial/config/PointStatConfig_Nepal \
> > > > -outdir tutorial/out/point_stat \
> > > > -v 2
> > > DEBUG 1: Default Config File:
> > >
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> > > WARNING:
> > > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> > > method, resetting method from "UW_MEAN" to "NEAREST" since the
> > > interpolation width is 1.
> > > WARNING:
> > > GSL_RNG_TYPE=mt19937
> > > GSL_RNG_SEED=1682548822
> > > DEBUG 1: Forecast File:
> > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > > DEBUG 2:
> > > DEBUG 2:
> > >
> >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for SPFH/P500.
> > > ERROR :
> > > ERROR : parse_grid_mask() -> the mask_grid requested "nepal" is
not
> > > defined.
> > > ERROR :
> > >
> > > Kindly see in the surapipith_data folder in rap.ucar.edu site
the
> latest
> > > config file, as well as the suite of Nepal mask files that Julie
made
> > and I
> > > use. Hope some progress would be before the end of the day, or
once
> some
> > > hints are next sent from there...
> > >
> > > Best wishes,
> > > Vanisa
> > >
> > >
> > > On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >> Hello Vanisa,
> > >>
> > >> My name is John, and I work with Julie who is out of the office
today.
> > I
> > >> see that you're having some trouble getting Point-Stat to run:
> > >>
> > >> ERROR : parse_grid_mask() -> the mask_grid requested
"nepal" is
> not
> > >> defined.
> > >>
> > >> This should be pretty easy to resolve. You've worked with
Julie to
> > >> generate a mask file for Nepal, and that mask is contained in
some
> > NetCDF
> > >> file (the output of the gen_vx_mask tool). For example, let's
say
> it's
> > in
> > >> "/home/vsurapipith/nepal_mask.nc". You'll need to put the path
to
> that
> > >> file in the "mask" section of the Point-Stat configuration
file:
> > >>
> > >> mask = {
> > >> grid = [];
> > >> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; // or
whatever
> the
> > >> actual path is to that data file
> > >> sid = [];
> > >> }
> > >>
> > >> And then try running Point-Stat again. For a more complete
> description
> > of
> > >> the masking options available, please see pages 54-55 of the
MET
> User's
> > >> Guide:
> > >>
> > >>
> > >>
> >
>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
> > >>
> > >> Please let me know how it goes.
> > >>
> > >> Thanks,
> > >> John
> > >>
> > >>
> > >>
> > >> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via RT <
> > >> met_help at ucar.edu
> > >> > wrote:
> > >>
> > >> >
> > >> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
>
> > >> >
> > >> > Hi Julie,
> > >> >
> > >> > I'm back....after a couple of weeks spending on stumbling
over the
> > >> > installation of UPP until I got it through, now I have the
GRIB1
> files
> > >> > systematically converted from the WRF output generated here.
> > >> >
> > >> > I try the GRIB1 files with the PointStatConfig_Nepal that I
sent
> > >> earlier to
> > >> > do Point_Stat with our these WRF output, but currently the
error
> shows
> > >> up
> > >> > as:
> > >> >
> > >> > ~/met-5.1_bugfix$ ./bin/point_stat
> > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > >> > tutorial/out/ascii2nc/nepal.nc
> tutorial/config/PointStatConfig_Nepal
> > >> > -outdir tutorial/out/point_stat -v 2
> > >> > DEBUG 1: Default Config File:
> > >> >
> > >>
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > >> > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > >> > WARNING:
> > >> > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> > >> method,
> > >> > resetting method from "UW_MEAN" to "NEAREST" since the
interpolation
> > >> width
> > >> > is 1.
> > >> > WARNING:
> > >> > GSL_RNG_TYPE=mt19937
> > >> > GSL_RNG_SEED=18446744072070604630
> > >> > DEBUG 1: Forecast File:
> > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > >> > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > >> > DEBUG 2:
> > >> > DEBUG 2:
> > >> >
> > >> >
> > >>
> >
>
--------------------------------------------------------------------------------
> > >> > DEBUG 2:
> > >> > DEBUG 2: Reading data for SPFH/P500.
> > >> > ERROR :
> > >> > ERROR : parse_grid_mask() -> the mask_grid requested "nepal"
is not
> > >> > defined.
> > >> > ERROR :
> > >> >
> > >> > It looks to me like the mask file - nepal.nc, that I
deliberately
> > took
> > >> > from
> > >> > what you generated earlier - was not seen/defined well by MET
here.
> I
> > >> have
> > >> > actually just put the WRF output files in GRIB1 format
(checked with
> > as
> > >> you
> > >> > recommend) on the rap.ucar.edu ftp site
(surapipith_data)...Would
> you
> > >> > kindly let me know what should be done to correct this,
please?
> > >> >
> > >> > Best wishes,
> > >> > Vanisa
> > >> >
> > >> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via RT <
> > >> > met_help at ucar.edu>
> > >> > wrote:
> > >> >
> > >> > > Hi Vanisa.
> > >> > >
> > >> > > While I don't have any experience installing UPP myself, I
did
> take
> > a
> > >> > look
> > >> > > at your log file where I saw:
> > >> > >
> > >> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No such
file or
> > >> > > directory compilation terminated.
> > >> > >
> > >> > > So, it looks like UPP is having trouble finding jasper.h.
> > >> > >
> > >> > > Support for UPP is provided via wrfhelp at ucar.edu. If
you're
> > >> unable
> > >> > to
> > >> > > resolve the jasper problem, I would suggest writing to them
> > regarding
> > >> > this
> > >> > > issue.
> > >> > >
> > >> > > Thanks,
> > >> > > Julie
> > >> > >
> > >> > >
> > >> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via RT <
> > >> > > met_help at ucar.edu>
> > >> > > wrote:
> > >> > >
> > >> > > >
> > >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > >> > > >
> > >> > > > Hi Julie,
> > >> > > >
> > >> > > > Happy New Year...:) Thank you very much for your reply
upon back
> > to
> > >> > > office.
> > >> > > > Indeed, I did not use UPP for processing the WRF (ARW)
output.
> > >> > Therefore,
> > >> > > > today I look into how to install UPP on my workstation
for this
> > >> > purpose.
> > >> > > > However, there is a problem that unipost.exe was not
generated
> > after
> > >> > > > compilation, while other files came through. I was
wondered if
> it
> > >> has
> > >> > to
> > >> > > do
> > >> > > > with liking to JASPER. I have not been able to spot why,
and now
> > >> attach
> > >> > > the
> > >> > > > compile_upp.log to you here, in case you might be able to
see
> what
> > >> is
> > >> > > wrong
> > >> > > > in there.
> > >> > > >
> > >> > > > Thanks again; update more soon.
> > >> > > >
> > >> > > > Best wishes,
> > >> > > > Vanisa
> > >> > > >
> > >> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT <
> > >> > > met_help at ucar.edu
> > >> > > > >
> > >> > > > wrote:
> > >> > > >
> > >> > > > > Hi Vanisa. The file you uploaded may have aged off of
our ftp
> > >> > site. I
> > >> > > > do
> > >> > > > > not see it there currently.
> > >> > > > >
> > >> > > > > I neglected to ask you if your WRF output file was
> > post-processed?
> > >> > MET
> > >> > > > is
> > >> > > > > unable to read the output of WRF directly. Instead,
you must
> > >> first
> > >> > > > > post-process it using either the Unified PostProcessor
(UPP)
> or
> > >> the
> > >> > > > pinterp
> > >> > > > > tool (for ARW only). We strongly recommend using UPP
since
> > MET's
> > >> > > support
> > >> > > > > is most robust for the GRIB1 format that UPP creates.
> > >> > > > >
> > >> > > > > If you have already post-processed your file, please
trying
> > >> uploading
> > >> > > to
> > >> > > > > our ftp site again. If you have not, please post-
process the
> > file
> > >> > and
> > >> > > > give
> > >> > > > > it another try. Thanks!
> > >> > > > >
> > >> > > > > Happy New Year!
> > >> > > > >
> > >> > > > > Julie
> > >> > > > >
>
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Feb 02 02:32:01 2016
VERSION MODEL FCST_LEAD FCST_VALID_BEG FCST_VALID_END OBS_LEAD
OBS_VALID_BEG OBS_VALID_END FCST_VAR FCST_LEV OBS_VAR OBS_LEV OBTYPE
VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH OBS_THRESH COV_THRESH
ALPHA LINE_TYPE TOTAL FBAR FBAR_NCL FBAR_NCU FBAR_BCL FBAR_BCU FSTDEV
FSTDEV_NCL FSTDEV_NCU FSTDEV_BCL FSTDEV_BCU OBAR OBAR_NCL OBAR_NCU
OBAR_BCL OBAR_BCU OSTDEV OSTDEV_NCL OSTDEV_NCU OSTDEV_BCL OSTDEV_BCU
PR_CORR PR_CORR_NCL PR_CORR_NCU PR_CORR_BCL PR_CORR_BCU SP_CORR
KT_CORR RANKS FRANK_TIES ORANK_TIES ME ME_NCL ME_NCU ME_BCL ME_BCU
ESTDEV ESTDEV_NCL ESTDEV_NCU ESTDEV_BCL ESTDEV_BCU MBIAS MBIAS_BCL
MBIAS_BCU MAE MAE_BCL MAE_BCU MSE MSE_BCL MSE_BCU BCMSE BCMSE_BCL
BCMSE_BCU RMSE RMSE_BCL RMSE_BCU E10 E10_BCL E10_BCU E25 E25_BCL
E25_BCU E50 E50_BCL E50_BCU E75 E75_BCL E75_BCU E90 E90_BCL E90_BCU
EIQR EIQR_BCL EIQR_BCU MAD MAD_BCL MAD_BCU ANOM_CORR ANOM_CORR_NCL
ANOM_CORR_NCU ANOM_CORR_BCL ANOM_CORR_BCU ME2 ME2_BCL ME2_BCU MSESS
MSESS_BCL MSESS_BCU
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Feb 02 02:32:01 2016
VERSION MODEL FCST_LEAD FCST_VALID_BEG FCST_VALID_END OBS_LEAD
OBS_VALID_BEG OBS_VALID_END FCST_VAR FCST_LEV OBS_VAR OBS_LEV OBTYPE
VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH OBS_THRESH COV_THRESH
ALPHA LINE_TYPE TOTAL FY_OY FY_ON FN_OY FN_ON
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Feb 02 02:32:01 2016
VERSION MODEL FCST_LEAD FCST_VALID_BEG FCST_VALID_END OBS_LEAD
OBS_VALID_BEG OBS_VALID_END FCST_VAR FCST_LEV OBS_VAR OBS_LEV OBTYPE
VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH OBS_THRESH COV_THRESH
ALPHA LINE_TYPE TOTAL FBAR FBAR_NCL FBAR_NCU FBAR_BCL FBAR_BCU FSTDEV
FSTDEV_NCL FSTDEV_NCU FSTDEV_BCL FSTDEV_BCU OBAR OBAR_NCL OBAR_NCU
OBAR_BCL OBAR_BCU OSTDEV OSTDEV_NCL OSTDEV_NCU OSTDEV_BCL OSTDEV_BCU
PR_CORR PR_CORR_NCL PR_CORR_NCU PR_CORR_BCL PR_CORR_BCU SP_CORR
KT_CORR RANKS FRANK_TIES ORANK_TIES ME ME_NCL ME_NCU ME_BCL ME_BCU
ESTDEV ESTDEV_NCL ESTDEV_NCU ESTDEV_BCL ESTDEV_BCU MBIAS MBIAS_BCL
MBIAS_BCU MAE MAE_BCL MAE_BCU MSE MSE_BCL MSE_BCU BCMSE BCMSE_BCL
BCMSE_BCU RMSE RMSE_BCL RMSE_BCU E10 E10_BCL E10_BCU E25 E25_BCL
E25_BCU E50 E50_BCL E50_BCU E75 E75_BCL E75_BCU E90 E90_BCL E90_BCU
EIQR EIQR_BCL EIQR_BCU MAD MAD_BCL MAD_BCU ANOM_CORR ANOM_CORR_NCL
ANOM_CORR_NCU ANOM_CORR_BCL ANOM_CORR_BCU ME2 ME2_BCL ME2_BCU MSESS
MSESS_BCL MSESS_BCU
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Tue Feb 02 02:32:01 2016
VERSION MODEL FCST_LEAD FCST_VALID_BEG FCST_VALID_END OBS_LEAD
OBS_VALID_BEG OBS_VALID_END FCST_VAR FCST_LEV OBS_VAR OBS_LEV OBTYPE
VX_MASK INTERP_MTHD INTERP_PNTS FCST_THRESH OBS_THRESH COV_THRESH
ALPHA LINE_TYPE TOTAL FY_OY FY_ON FN_OY FN_ON
------------------------------------------------
Subject: point_stat for WRF output question
From: Julie Prestopnik
Time: Wed Feb 03 10:56:06 2016
Hi Vanisa.
Could you please upload your input forecast file, input observation
file,
masking file, and your configuration file? Also, could you please
send me
the command line argument for point_stat that you are running?
I can take a closer look and let you know what I find out.
Thanks,
Julie
On Tue, Feb 2, 2016 at 2:32 AM, vanisas at gmail.com via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>
> Hi Julie,
>
> thank you for your reply last week; apologies for my delay in
reporting.
> What you suggest works. I can now generate the output files using
> Point_Stat. However, I do not yet understand how to interpret or
make use
> of the output...I have been trying to read what described in the
> MET_User_Guide. Frankly speaking, I still do not get it. Could you
kindly
> see in my folder in the rap.ucar.edu ftp site in
incoming/irap/met_help
> under surapipith_data. I attached here some of the files. It looks
like
> only the headers appear...or what I should do to generate the
statistical
> evaluation output. What surprise me was that the example files I got
from
> the US example (follow the online tutorial) look similar. Could you
guide
> me where to read to comprehend MET?
>
> One other question I have is that probably this is because I use
only one
> hour of WRF output for this try. Do I need to generate a GRIB1 file
that
> combine all the required hours...otherwise, can MET read in more
than one
> GRIB1 files at a time for Point_Stat analysis?
>
> Best wishes,
> Vanisa
>
> On Thu, Jan 28, 2016 at 12:00 AM, Julie Prestopnik via RT <
> met_help at ucar.edu
> > wrote:
>
> > Hi Vanisa. This is Julie again.
> >
> > I see that you are encountering a problem due to the masking and
> > verification grids not matching.
> >
> > One projection is a Lat/Lon grid and the other is a Mercator.
> >
> > The mask file I created for you using gen_vx_mask used a Lat/Lon
grid
> file
> > (gfs_2012040900_F003.grib). If you would like to get your netcdf
mask
> > file
> > on a Mercator grid, you could run gen_vx_mask using your file that
> contains
> > the grid that you would like to verify, which I believe is on a
Mercator
> > projection (WRFPRS_d02.00?) and the nepal poly file which will
then
> produce
> > a netcdf file on that same projection.
> >
> > So, for example, something like:
> >
> > gen_vx_mask WRFPRS_d02.00 nepal_detail.poly nepal_detail.nc
> >
> > Then, you can rerun point_stat using the new netcdf mask file on
the
> > Mercator grid and that should fix the problem.
> >
> > However, if you're not going to be running thousands of cases of
Point
> > Stat, you don't necessarily need to run gen_vx_mask, but rather
you could
> > use the poly file directly with Point Stat in the configuration
file in
> the
> > mask section that John described earlier. Using the poly file
directly
> can
> > be slow, so if you're running point_stat on a large number of
files, this
> > is likely not the best way to use point_stat.
> >
> > I hope this helps. Please let me know if you have any questions.
> >
> > Thanks,
> > Julie
> >
> > On Tue, Jan 26, 2016 at 11:04 PM, vanisas at gmail.com via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > >
> > > Hi John,
> > >
> > > apologies, I just notice that I did not remove the word "nepal"
in the
> > > bracket ([]) behind grid = within mask section in the
PointConfig
> file. I
> > > have then removed that, and found a different error (also rename
the
> > ascii
> > > file containing the ADPUPA observation from our stations in
Nepal to be
> > > called "nepal_ascii.nc" to avoid confusion) :
> > >
> > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > > > tutorial/out/ascii2nc/nepal_ascii.nc \
> > > > tutorial/config/PointStatConfig_Nepal \
> > > > -outdir tutorial/out/point_stat \
> > > > -v 2
> > > DEBUG 1: Default Config File:
> > >
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > DEBUG 1: User Config File: tutorial/config/PointStatConfig_Nepal
> > > WARNING:
> > > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> > method,
> > > resetting method from "UW_MEAN" to "NEAREST" since the
interpolation
> > width
> > > is 1.
> > > WARNING:
> > > GSL_RNG_TYPE=mt19937
> > > GSL_RNG_SEED=18446744072507402326
> > > DEBUG 1: Forecast File:
> > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal_ascii.nc
> > > DEBUG 2:
> > > DEBUG 2:
> > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for SPFH/P500.
> > > ERROR :
> > > ERROR : parse_poly_mask() -> The masking and verification grids
do not
> > > match:
> > > ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll: -90.000
lon_ll:
> > -0.000
> > > delta_lat: 0.500 delta_lon: 0.500 !=
> > > ERROR : Projection: Mercator Nx: 270 Ny: 180 Lat_LL_radians:
0.4178
> > > Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423 Lon_UR_radians:
-1.5680
> > Mx:
> > > 1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
> > > ERROR :
> > >
> > > It looks like I need to set for some correct projection
somewhere.
> Could
> > > you give a bit of guideline (in the while, I'll try to read the
> > > Met_Users_Guide_v5.1.pdf in more detail; may send update again
before
> end
> > > of the day here)...
> > >
> > > Best wishes,
> > > Vanisa
> > >
> > > On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith
<vanisas at gmail.com
> >
> > > wrote:
> > >
> > > > Hi John,
> > > >
> > > > thank you for your reply; I recall your help earlier, and I am
only
> > > > gradually getting to know MET. I am grateful to the help as to
what
> > Julie
> > > > was responding patiently as well. I notice the (presumably)
mask
> files
> > > that
> > > > Julie helped generating for Nepal, and I try them. However,
the
> similar
> > > > error still shows up, and I am trying to understand why:
> > > >
> > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > > > > tutorial/out/ascii2nc/nepal.nc
> > tutorial/config/PointStatConfig_Nepal \
> > > > > -outdir tutorial/out/point_stat \
> > > > > -v 2
> > > > DEBUG 1: Default Config File:
> > > >
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > > WARNING:
> > > > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> > > > method, resetting method from "UW_MEAN" to "NEAREST" since the
> > > > interpolation width is 1.
> > > > WARNING:
> > > > GSL_RNG_TYPE=mt19937
> > > > GSL_RNG_SEED=1682548822
> > > > DEBUG 1: Forecast File:
> > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > DEBUG 2:
> > > > DEBUG 2: Reading data for SPFH/P500.
> > > > ERROR :
> > > > ERROR : parse_grid_mask() -> the mask_grid requested "nepal"
is not
> > > > defined.
> > > > ERROR :
> > > >
> > > > Kindly see in the surapipith_data folder in rap.ucar.edu site
the
> > latest
> > > > config file, as well as the suite of Nepal mask files that
Julie made
> > > and I
> > > > use. Hope some progress would be before the end of the day, or
once
> > some
> > > > hints are next sent from there...
> > > >
> > > > Best wishes,
> > > > Vanisa
> > > >
> > > >
> > > > On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > >> Hello Vanisa,
> > > >>
> > > >> My name is John, and I work with Julie who is out of the
office
> today.
> > > I
> > > >> see that you're having some trouble getting Point-Stat to
run:
> > > >>
> > > >> ERROR : parse_grid_mask() -> the mask_grid requested
"nepal" is
> > not
> > > >> defined.
> > > >>
> > > >> This should be pretty easy to resolve. You've worked with
Julie to
> > > >> generate a mask file for Nepal, and that mask is contained in
some
> > > NetCDF
> > > >> file (the output of the gen_vx_mask tool). For example,
let's say
> > it's
> > > in
> > > >> "/home/vsurapipith/nepal_mask.nc". You'll need to put the
path to
> > that
> > > >> file in the "mask" section of the Point-Stat configuration
file:
> > > >>
> > > >> mask = {
> > > >> grid = [];
> > > >> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; // or
whatever
> > the
> > > >> actual path is to that data file
> > > >> sid = [];
> > > >> }
> > > >>
> > > >> And then try running Point-Stat again. For a more complete
> > description
> > > of
> > > >> the masking options available, please see pages 54-55 of the
MET
> > User's
> > > >> Guide:
> > > >>
> > > >>
> > > >>
> > >
> >
>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
> > > >>
> > > >> Please let me know how it goes.
> > > >>
> > > >> Thanks,
> > > >> John
> > > >>
> > > >>
> > > >>
> > > >> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via RT <
> > > >> met_help at ucar.edu
> > > >> > wrote:
> > > >>
> > > >> >
> > > >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > > >> >
> > > >> > Hi Julie,
> > > >> >
> > > >> > I'm back....after a couple of weeks spending on stumbling
over the
> > > >> > installation of UPP until I got it through, now I have the
GRIB1
> > files
> > > >> > systematically converted from the WRF output generated
here.
> > > >> >
> > > >> > I try the GRIB1 files with the PointStatConfig_Nepal that I
sent
> > > >> earlier to
> > > >> > do Point_Stat with our these WRF output, but currently the
error
> > shows
> > > >> up
> > > >> > as:
> > > >> >
> > > >> > ~/met-5.1_bugfix$ ./bin/point_stat
> > > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > >> > tutorial/out/ascii2nc/nepal.nc
> > tutorial/config/PointStatConfig_Nepal
> > > >> > -outdir tutorial/out/point_stat -v 2
> > > >> > DEBUG 1: Default Config File:
> > > >> >
> > > >>
> > >
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > >> > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > >> > WARNING:
> > > >> > WARNING: parse_conf_interp() -> For neareast neighbor
> interpolation
> > > >> method,
> > > >> > resetting method from "UW_MEAN" to "NEAREST" since the
> interpolation
> > > >> width
> > > >> > is 1.
> > > >> > WARNING:
> > > >> > GSL_RNG_TYPE=mt19937
> > > >> > GSL_RNG_SEED=18446744072070604630
> > > >> > DEBUG 1: Forecast File:
> > > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > >> > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > > >> > DEBUG 2:
> > > >> > DEBUG 2:
> > > >> >
> > > >> >
> > > >>
> > >
> >
>
--------------------------------------------------------------------------------
> > > >> > DEBUG 2:
> > > >> > DEBUG 2: Reading data for SPFH/P500.
> > > >> > ERROR :
> > > >> > ERROR : parse_grid_mask() -> the mask_grid requested
"nepal" is
> not
> > > >> > defined.
> > > >> > ERROR :
> > > >> >
> > > >> > It looks to me like the mask file - nepal.nc, that I
deliberately
> > > took
> > > >> > from
> > > >> > what you generated earlier - was not seen/defined well by
MET
> here.
> > I
> > > >> have
> > > >> > actually just put the WRF output files in GRIB1 format
(checked
> with
> > > as
> > > >> you
> > > >> > recommend) on the rap.ucar.edu ftp site
(surapipith_data)...Would
> > you
> > > >> > kindly let me know what should be done to correct this,
please?
> > > >> >
> > > >> > Best wishes,
> > > >> > Vanisa
> > > >> >
> > > >> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via RT <
> > > >> > met_help at ucar.edu>
> > > >> > wrote:
> > > >> >
> > > >> > > Hi Vanisa.
> > > >> > >
> > > >> > > While I don't have any experience installing UPP myself,
I did
> > take
> > > a
> > > >> > look
> > > >> > > at your log file where I saw:
> > > >> > >
> > > >> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No
such file
> or
> > > >> > > directory compilation terminated.
> > > >> > >
> > > >> > > So, it looks like UPP is having trouble finding jasper.h.
> > > >> > >
> > > >> > > Support for UPP is provided via wrfhelp at ucar.edu. If
you're
> > > >> unable
> > > >> > to
> > > >> > > resolve the jasper problem, I would suggest writing to
them
> > > regarding
> > > >> > this
> > > >> > > issue.
> > > >> > >
> > > >> > > Thanks,
> > > >> > > Julie
> > > >> > >
> > > >> > >
> > > >> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via RT
<
> > > >> > > met_help at ucar.edu>
> > > >> > > wrote:
> > > >> > >
> > > >> > > >
> > > >> > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> >
> > > >> > > >
> > > >> > > > Hi Julie,
> > > >> > > >
> > > >> > > > Happy New Year...:) Thank you very much for your reply
upon
> back
> > > to
> > > >> > > office.
> > > >> > > > Indeed, I did not use UPP for processing the WRF (ARW)
output.
> > > >> > Therefore,
> > > >> > > > today I look into how to install UPP on my workstation
for
> this
> > > >> > purpose.
> > > >> > > > However, there is a problem that unipost.exe was not
generated
> > > after
> > > >> > > > compilation, while other files came through. I was
wondered if
> > it
> > > >> has
> > > >> > to
> > > >> > > do
> > > >> > > > with liking to JASPER. I have not been able to spot
why, and
> now
> > > >> attach
> > > >> > > the
> > > >> > > > compile_upp.log to you here, in case you might be able
to see
> > what
> > > >> is
> > > >> > > wrong
> > > >> > > > in there.
> > > >> > > >
> > > >> > > > Thanks again; update more soon.
> > > >> > > >
> > > >> > > > Best wishes,
> > > >> > > > Vanisa
> > > >> > > >
> > > >> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via RT
<
> > > >> > > met_help at ucar.edu
> > > >> > > > >
> > > >> > > > wrote:
> > > >> > > >
> > > >> > > > > Hi Vanisa. The file you uploaded may have aged off
of our
> ftp
> > > >> > site. I
> > > >> > > > do
> > > >> > > > > not see it there currently.
> > > >> > > > >
> > > >> > > > > I neglected to ask you if your WRF output file was
> > > post-processed?
> > > >> > MET
> > > >> > > > is
> > > >> > > > > unable to read the output of WRF directly. Instead,
you
> must
> > > >> first
> > > >> > > > > post-process it using either the Unified
PostProcessor (UPP)
> > or
> > > >> the
> > > >> > > > pinterp
> > > >> > > > > tool (for ARW only). We strongly recommend using UPP
since
> > > MET's
> > > >> > > support
> > > >> > > > > is most robust for the GRIB1 format that UPP creates.
> > > >> > > > >
> > > >> > > > > If you have already post-processed your file, please
trying
> > > >> uploading
> > > >> > > to
> > > >> > > > > our ftp site again. If you have not, please post-
process
> the
> > > file
> > > >> > and
> > > >> > > > give
> > > >> > > > > it another try. Thanks!
> > > >> > > > >
> > > >> > > > > Happy New Year!
> > > >> > > > >
> > > >> > > > > Julie
> > > >> > > > >
> >
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Thu Feb 04 05:47:37 2016
Hi Julie,
I have just uploaded the 4 files I used for running point_stat that
generated the output files that I put there earlier. Hope these are
comprehensible. the command lines that I used are as follows:
~/met-5.1_bugfix$ ./bin/point_stat \
> /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> tutorial/out/ascii2nc/nepal_ascii.nc \
> tutorial/config/PointStatConfig_Nepal \
> -outdir tutorial/out/point_stat \
> -v 2
Thanks alot for your kind help. Wish to see this through soon.
Best wishes,
Vanisa
On Wed, Feb 3, 2016 at 11:41 PM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Vanisa.
>
> Could you please upload your input forecast file, input observation
file,
> masking file, and your configuration file? Also, could you please
send me
> the command line argument for point_stat that you are running?
>
> I can take a closer look and let you know what I find out.
>
> Thanks,
> Julie
>
> On Tue, Feb 2, 2016 at 2:32 AM, vanisas at gmail.com via RT <
> met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >
> > Hi Julie,
> >
> > thank you for your reply last week; apologies for my delay in
reporting.
> > What you suggest works. I can now generate the output files using
> > Point_Stat. However, I do not yet understand how to interpret or
make use
> > of the output...I have been trying to read what described in the
> > MET_User_Guide. Frankly speaking, I still do not get it. Could you
kindly
> > see in my folder in the rap.ucar.edu ftp site in
incoming/irap/met_help
> > under surapipith_data. I attached here some of the files. It looks
like
> > only the headers appear...or what I should do to generate the
statistical
> > evaluation output. What surprise me was that the example files I
got from
> > the US example (follow the online tutorial) look similar. Could
you guide
> > me where to read to comprehend MET?
> >
> > One other question I have is that probably this is because I use
only one
> > hour of WRF output for this try. Do I need to generate a GRIB1
file that
> > combine all the required hours...otherwise, can MET read in more
than one
> > GRIB1 files at a time for Point_Stat analysis?
> >
> > Best wishes,
> > Vanisa
> >
> > On Thu, Jan 28, 2016 at 12:00 AM, Julie Prestopnik via RT <
> > met_help at ucar.edu
> > > wrote:
> >
> > > Hi Vanisa. This is Julie again.
> > >
> > > I see that you are encountering a problem due to the masking and
> > > verification grids not matching.
> > >
> > > One projection is a Lat/Lon grid and the other is a Mercator.
> > >
> > > The mask file I created for you using gen_vx_mask used a Lat/Lon
grid
> > file
> > > (gfs_2012040900_F003.grib). If you would like to get your
netcdf mask
> > > file
> > > on a Mercator grid, you could run gen_vx_mask using your file
that
> > contains
> > > the grid that you would like to verify, which I believe is on a
> Mercator
> > > projection (WRFPRS_d02.00?) and the nepal poly file which will
then
> > produce
> > > a netcdf file on that same projection.
> > >
> > > So, for example, something like:
> > >
> > > gen_vx_mask WRFPRS_d02.00 nepal_detail.poly nepal_detail.nc
> > >
> > > Then, you can rerun point_stat using the new netcdf mask file
on the
> > > Mercator grid and that should fix the problem.
> > >
> > > However, if you're not going to be running thousands of cases of
Point
> > > Stat, you don't necessarily need to run gen_vx_mask, but rather
you
> could
> > > use the poly file directly with Point Stat in the configuration
file in
> > the
> > > mask section that John described earlier. Using the poly file
directly
> > can
> > > be slow, so if you're running point_stat on a large number of
files,
> this
> > > is likely not the best way to use point_stat.
> > >
> > > I hope this helps. Please let me know if you have any
questions.
> > >
> > > Thanks,
> > > Julie
> > >
> > > On Tue, Jan 26, 2016 at 11:04 PM, vanisas at gmail.com via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
>
> > > >
> > > > Hi John,
> > > >
> > > > apologies, I just notice that I did not remove the word
"nepal" in
> the
> > > > bracket ([]) behind grid = within mask section in the
PointConfig
> > file. I
> > > > have then removed that, and found a different error (also
rename the
> > > ascii
> > > > file containing the ADPUPA observation from our stations in
Nepal to
> be
> > > > called "nepal_ascii.nc" to avoid confusion) :
> > > >
> > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > > > > tutorial/out/ascii2nc/nepal_ascii.nc \
> > > > > tutorial/config/PointStatConfig_Nepal \
> > > > > -outdir tutorial/out/point_stat \
> > > > > -v 2
> > > > DEBUG 1: Default Config File:
> > > >
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > > WARNING:
> > > > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> > > method,
> > > > resetting method from "UW_MEAN" to "NEAREST" since the
interpolation
> > > width
> > > > is 1.
> > > > WARNING:
> > > > GSL_RNG_TYPE=mt19937
> > > > GSL_RNG_SEED=18446744072507402326
> > > > DEBUG 1: Forecast File:
> > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal_ascii.nc
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > DEBUG 2:
> > > > DEBUG 2: Reading data for SPFH/P500.
> > > > ERROR :
> > > > ERROR : parse_poly_mask() -> The masking and verification
grids do
> not
> > > > match:
> > > > ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll: -90.000
lon_ll:
> > > -0.000
> > > > delta_lat: 0.500 delta_lon: 0.500 !=
> > > > ERROR : Projection: Mercator Nx: 270 Ny: 180 Lat_LL_radians:
0.4178
> > > > Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423 Lon_UR_radians:
> -1.5680
> > > Mx:
> > > > 1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
> > > > ERROR :
> > > >
> > > > It looks like I need to set for some correct projection
somewhere.
> > Could
> > > > you give a bit of guideline (in the while, I'll try to read
the
> > > > Met_Users_Guide_v5.1.pdf in more detail; may send update again
before
> > end
> > > > of the day here)...
> > > >
> > > > Best wishes,
> > > > Vanisa
> > > >
> > > > On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith <
> vanisas at gmail.com
> > >
> > > > wrote:
> > > >
> > > > > Hi John,
> > > > >
> > > > > thank you for your reply; I recall your help earlier, and I
am only
> > > > > gradually getting to know MET. I am grateful to the help as
to what
> > > Julie
> > > > > was responding patiently as well. I notice the (presumably)
mask
> > files
> > > > that
> > > > > Julie helped generating for Nepal, and I try them. However,
the
> > similar
> > > > > error still shows up, and I am trying to understand why:
> > > > >
> > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > > > > > tutorial/out/ascii2nc/nepal.nc
> > > tutorial/config/PointStatConfig_Nepal \
> > > > > > -outdir tutorial/out/point_stat \
> > > > > > -v 2
> > > > > DEBUG 1: Default Config File:
> > > > >
> > >
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > > > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > > > WARNING:
> > > > > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> > > > > method, resetting method from "UW_MEAN" to "NEAREST" since
the
> > > > > interpolation width is 1.
> > > > > WARNING:
> > > > > GSL_RNG_TYPE=mt19937
> > > > > GSL_RNG_SEED=1682548822
> > > > > DEBUG 1: Forecast File:
> > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > DEBUG 2:
> > > > > DEBUG 2: Reading data for SPFH/P500.
> > > > > ERROR :
> > > > > ERROR : parse_grid_mask() -> the mask_grid requested
"nepal" is
> not
> > > > > defined.
> > > > > ERROR :
> > > > >
> > > > > Kindly see in the surapipith_data folder in rap.ucar.edu
site the
> > > latest
> > > > > config file, as well as the suite of Nepal mask files that
Julie
> made
> > > > and I
> > > > > use. Hope some progress would be before the end of the day,
or once
> > > some
> > > > > hints are next sent from there...
> > > > >
> > > > > Best wishes,
> > > > > Vanisa
> > > > >
> > > > >
> > > > > On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway via RT
<
> > > > > met_help at ucar.edu> wrote:
> > > > >
> > > > >> Hello Vanisa,
> > > > >>
> > > > >> My name is John, and I work with Julie who is out of the
office
> > today.
> > > > I
> > > > >> see that you're having some trouble getting Point-Stat to
run:
> > > > >>
> > > > >> ERROR : parse_grid_mask() -> the mask_grid requested
"nepal"
> is
> > > not
> > > > >> defined.
> > > > >>
> > > > >> This should be pretty easy to resolve. You've worked with
Julie
> to
> > > > >> generate a mask file for Nepal, and that mask is contained
in some
> > > > NetCDF
> > > > >> file (the output of the gen_vx_mask tool). For example,
let's say
> > > it's
> > > > in
> > > > >> "/home/vsurapipith/nepal_mask.nc". You'll need to put the
path
> to
> > > that
> > > > >> file in the "mask" section of the Point-Stat configuration
file:
> > > > >>
> > > > >> mask = {
> > > > >> grid = [];
> > > > >> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; // or
> whatever
> > > the
> > > > >> actual path is to that data file
> > > > >> sid = [];
> > > > >> }
> > > > >>
> > > > >> And then try running Point-Stat again. For a more complete
> > > description
> > > > of
> > > > >> the masking options available, please see pages 54-55 of
the MET
> > > User's
> > > > >> Guide:
> > > > >>
> > > > >>
> > > > >>
> > > >
> > >
> >
>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
> > > > >>
> > > > >> Please let me know how it goes.
> > > > >>
> > > > >> Thanks,
> > > > >> John
> > > > >>
> > > > >>
> > > > >>
> > > > >> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via RT <
> > > > >> met_help at ucar.edu
> > > > >> > wrote:
> > > > >>
> > > > >> >
> > > > >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > > > >> >
> > > > >> > Hi Julie,
> > > > >> >
> > > > >> > I'm back....after a couple of weeks spending on stumbling
over
> the
> > > > >> > installation of UPP until I got it through, now I have
the GRIB1
> > > files
> > > > >> > systematically converted from the WRF output generated
here.
> > > > >> >
> > > > >> > I try the GRIB1 files with the PointStatConfig_Nepal that
I sent
> > > > >> earlier to
> > > > >> > do Point_Stat with our these WRF output, but currently
the error
> > > shows
> > > > >> up
> > > > >> > as:
> > > > >> >
> > > > >> > ~/met-5.1_bugfix$ ./bin/point_stat
> > > > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > >> > tutorial/out/ascii2nc/nepal.nc
> > > tutorial/config/PointStatConfig_Nepal
> > > > >> > -outdir tutorial/out/point_stat -v 2
> > > > >> > DEBUG 1: Default Config File:
> > > > >> >
> > > > >>
> > > >
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > > >> > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > > >> > WARNING:
> > > > >> > WARNING: parse_conf_interp() -> For neareast neighbor
> > interpolation
> > > > >> method,
> > > > >> > resetting method from "UW_MEAN" to "NEAREST" since the
> > interpolation
> > > > >> width
> > > > >> > is 1.
> > > > >> > WARNING:
> > > > >> > GSL_RNG_TYPE=mt19937
> > > > >> > GSL_RNG_SEED=18446744072070604630
> > > > >> > DEBUG 1: Forecast File:
> > > > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > >> > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > > > >> > DEBUG 2:
> > > > >> > DEBUG 2:
> > > > >> >
> > > > >> >
> > > > >>
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > >> > DEBUG 2:
> > > > >> > DEBUG 2: Reading data for SPFH/P500.
> > > > >> > ERROR :
> > > > >> > ERROR : parse_grid_mask() -> the mask_grid requested
"nepal" is
> > not
> > > > >> > defined.
> > > > >> > ERROR :
> > > > >> >
> > > > >> > It looks to me like the mask file - nepal.nc, that I
> deliberately
> > > > took
> > > > >> > from
> > > > >> > what you generated earlier - was not seen/defined well by
MET
> > here.
> > > I
> > > > >> have
> > > > >> > actually just put the WRF output files in GRIB1 format
(checked
> > with
> > > > as
> > > > >> you
> > > > >> > recommend) on the rap.ucar.edu ftp site
> (surapipith_data)...Would
> > > you
> > > > >> > kindly let me know what should be done to correct this,
please?
> > > > >> >
> > > > >> > Best wishes,
> > > > >> > Vanisa
> > > > >> >
> > > > >> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via RT
<
> > > > >> > met_help at ucar.edu>
> > > > >> > wrote:
> > > > >> >
> > > > >> > > Hi Vanisa.
> > > > >> > >
> > > > >> > > While I don't have any experience installing UPP
myself, I did
> > > take
> > > > a
> > > > >> > look
> > > > >> > > at your log file where I saw:
> > > > >> > >
> > > > >> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No
such
> file
> > or
> > > > >> > > directory compilation terminated.
> > > > >> > >
> > > > >> > > So, it looks like UPP is having trouble finding
jasper.h.
> > > > >> > >
> > > > >> > > Support for UPP is provided via wrfhelp at ucar.edu.
If
> you're
> > > > >> unable
> > > > >> > to
> > > > >> > > resolve the jasper problem, I would suggest writing to
them
> > > > regarding
> > > > >> > this
> > > > >> > > issue.
> > > > >> > >
> > > > >> > > Thanks,
> > > > >> > > Julie
> > > > >> > >
> > > > >> > >
> > > > >> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via
RT <
> > > > >> > > met_help at ucar.edu>
> > > > >> > > wrote:
> > > > >> > >
> > > > >> > > >
> > > > >> > > > <URL:
> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > >
> > > > >> > > >
> > > > >> > > > Hi Julie,
> > > > >> > > >
> > > > >> > > > Happy New Year...:) Thank you very much for your
reply upon
> > back
> > > > to
> > > > >> > > office.
> > > > >> > > > Indeed, I did not use UPP for processing the WRF
(ARW)
> output.
> > > > >> > Therefore,
> > > > >> > > > today I look into how to install UPP on my
workstation for
> > this
> > > > >> > purpose.
> > > > >> > > > However, there is a problem that unipost.exe was not
> generated
> > > > after
> > > > >> > > > compilation, while other files came through. I was
wondered
> if
> > > it
> > > > >> has
> > > > >> > to
> > > > >> > > do
> > > > >> > > > with liking to JASPER. I have not been able to spot
why, and
> > now
> > > > >> attach
> > > > >> > > the
> > > > >> > > > compile_upp.log to you here, in case you might be
able to
> see
> > > what
> > > > >> is
> > > > >> > > wrong
> > > > >> > > > in there.
> > > > >> > > >
> > > > >> > > > Thanks again; update more soon.
> > > > >> > > >
> > > > >> > > > Best wishes,
> > > > >> > > > Vanisa
> > > > >> > > >
> > > > >> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik via
RT <
> > > > >> > > met_help at ucar.edu
> > > > >> > > > >
> > > > >> > > > wrote:
> > > > >> > > >
> > > > >> > > > > Hi Vanisa. The file you uploaded may have aged off
of our
> > ftp
> > > > >> > site. I
> > > > >> > > > do
> > > > >> > > > > not see it there currently.
> > > > >> > > > >
> > > > >> > > > > I neglected to ask you if your WRF output file was
> > > > post-processed?
> > > > >> > MET
> > > > >> > > > is
> > > > >> > > > > unable to read the output of WRF directly.
Instead, you
> > must
> > > > >> first
> > > > >> > > > > post-process it using either the Unified
PostProcessor
> (UPP)
> > > or
> > > > >> the
> > > > >> > > > pinterp
> > > > >> > > > > tool (for ARW only). We strongly recommend using
UPP
> since
> > > > MET's
> > > > >> > > support
> > > > >> > > > > is most robust for the GRIB1 format that UPP
creates.
> > > > >> > > > >
> > > > >> > > > > If you have already post-processed your file,
please
> trying
> > > > >> uploading
> > > > >> > > to
> > > > >> > > > > our ftp site again. If you have not, please post-
process
> > the
> > > > file
> > > > >> > and
> > > > >> > > > give
> > > > >> > > > > it another try. Thanks!
> > > > >> > > > >
> > > > >> > > > > Happy New Year!
> > > > >> > > > >
> > > > >> > > > > Julie
> > > > >> > > > >
> > >
> >
> >
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>
------------------------------------------------
Subject: point_stat for WRF output question
From: Julie Prestopnik
Time: Fri Feb 05 10:32:00 2016
Hi Vanisa. I apologize for the delayed response. I was out of the
office
yesterday.
Unfortunately, I need some help from John to figure out what exactly
is
going on here, but he is out of the office today. We will revisit
this on
Monday and get back to you as soon as we can.
Thank you,
Julie
On Thu, Feb 4, 2016 at 5:47 AM, vanisas at gmail.com via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>
> Hi Julie,
>
> I have just uploaded the 4 files I used for running point_stat that
> generated the output files that I put there earlier. Hope these are
> comprehensible. the command lines that I used are as follows:
>
> ~/met-5.1_bugfix$ ./bin/point_stat \
>
> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
>
> > tutorial/out/ascii2nc/nepal_ascii.nc \
>
> > tutorial/config/PointStatConfig_Nepal \
>
> > -outdir tutorial/out/point_stat \
>
> > -v 2
>
>
> Thanks alot for your kind help. Wish to see this through soon.
>
>
> Best wishes,
> Vanisa
>
> On Wed, Feb 3, 2016 at 11:41 PM, Julie Prestopnik via RT <
> met_help at ucar.edu>
> wrote:
>
> > Hi Vanisa.
> >
> > Could you please upload your input forecast file, input
observation file,
> > masking file, and your configuration file? Also, could you please
send
> me
> > the command line argument for point_stat that you are running?
> >
> > I can take a closer look and let you know what I find out.
> >
> > Thanks,
> > Julie
> >
> > On Tue, Feb 2, 2016 at 2:32 AM, vanisas at gmail.com via RT <
> > met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > >
> > > Hi Julie,
> > >
> > > thank you for your reply last week; apologies for my delay in
> reporting.
> > > What you suggest works. I can now generate the output files
using
> > > Point_Stat. However, I do not yet understand how to interpret or
make
> use
> > > of the output...I have been trying to read what described in the
> > > MET_User_Guide. Frankly speaking, I still do not get it. Could
you
> kindly
> > > see in my folder in the rap.ucar.edu ftp site in
> incoming/irap/met_help
> > > under surapipith_data. I attached here some of the files. It
looks like
> > > only the headers appear...or what I should do to generate the
> statistical
> > > evaluation output. What surprise me was that the example files I
got
> from
> > > the US example (follow the online tutorial) look similar. Could
you
> guide
> > > me where to read to comprehend MET?
> > >
> > > One other question I have is that probably this is because I use
only
> one
> > > hour of WRF output for this try. Do I need to generate a GRIB1
file
> that
> > > combine all the required hours...otherwise, can MET read in more
than
> one
> > > GRIB1 files at a time for Point_Stat analysis?
> > >
> > > Best wishes,
> > > Vanisa
> > >
> > > On Thu, Jan 28, 2016 at 12:00 AM, Julie Prestopnik via RT <
> > > met_help at ucar.edu
> > > > wrote:
> > >
> > > > Hi Vanisa. This is Julie again.
> > > >
> > > > I see that you are encountering a problem due to the masking
and
> > > > verification grids not matching.
> > > >
> > > > One projection is a Lat/Lon grid and the other is a Mercator.
> > > >
> > > > The mask file I created for you using gen_vx_mask used a
Lat/Lon grid
> > > file
> > > > (gfs_2012040900_F003.grib). If you would like to get your
netcdf
> mask
> > > > file
> > > > on a Mercator grid, you could run gen_vx_mask using your file
that
> > > contains
> > > > the grid that you would like to verify, which I believe is on
a
> > Mercator
> > > > projection (WRFPRS_d02.00?) and the nepal poly file which will
then
> > > produce
> > > > a netcdf file on that same projection.
> > > >
> > > > So, for example, something like:
> > > >
> > > > gen_vx_mask WRFPRS_d02.00 nepal_detail.poly nepal_detail.nc
> > > >
> > > > Then, you can rerun point_stat using the new netcdf mask file
on the
> > > > Mercator grid and that should fix the problem.
> > > >
> > > > However, if you're not going to be running thousands of cases
of
> Point
> > > > Stat, you don't necessarily need to run gen_vx_mask, but
rather you
> > could
> > > > use the poly file directly with Point Stat in the
configuration file
> in
> > > the
> > > > mask section that John described earlier. Using the poly file
> directly
> > > can
> > > > be slow, so if you're running point_stat on a large number of
files,
> > this
> > > > is likely not the best way to use point_stat.
> > > >
> > > > I hope this helps. Please let me know if you have any
questions.
> > > >
> > > > Thanks,
> > > > Julie
> > > >
> > > > On Tue, Jan 26, 2016 at 11:04 PM, vanisas at gmail.com via RT <
> > > > met_help at ucar.edu> wrote:
> > > >
> > > > >
> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > > > >
> > > > > Hi John,
> > > > >
> > > > > apologies, I just notice that I did not remove the word
"nepal" in
> > the
> > > > > bracket ([]) behind grid = within mask section in the
PointConfig
> > > file. I
> > > > > have then removed that, and found a different error (also
rename
> the
> > > > ascii
> > > > > file containing the ADPUPA observation from our stations in
Nepal
> to
> > be
> > > > > called "nepal_ascii.nc" to avoid confusion) :
> > > > >
> > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > > > > > tutorial/out/ascii2nc/nepal_ascii.nc \
> > > > > > tutorial/config/PointStatConfig_Nepal \
> > > > > > -outdir tutorial/out/point_stat \
> > > > > > -v 2
> > > > > DEBUG 1: Default Config File:
> > > > >
> > >
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > > > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > > > WARNING:
> > > > > WARNING: parse_conf_interp() -> For neareast neighbor
interpolation
> > > > method,
> > > > > resetting method from "UW_MEAN" to "NEAREST" since the
> interpolation
> > > > width
> > > > > is 1.
> > > > > WARNING:
> > > > > GSL_RNG_TYPE=mt19937
> > > > > GSL_RNG_SEED=18446744072507402326
> > > > > DEBUG 1: Forecast File:
> > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal_ascii.nc
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > DEBUG 2:
> > > > > DEBUG 2: Reading data for SPFH/P500.
> > > > > ERROR :
> > > > > ERROR : parse_poly_mask() -> The masking and verification
grids do
> > not
> > > > > match:
> > > > > ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll: -90.000
> lon_ll:
> > > > -0.000
> > > > > delta_lat: 0.500 delta_lon: 0.500 !=
> > > > > ERROR : Projection: Mercator Nx: 270 Ny: 180
Lat_LL_radians:
> 0.4178
> > > > > Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423
Lon_UR_radians:
> > -1.5680
> > > > Mx:
> > > > > 1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
> > > > > ERROR :
> > > > >
> > > > > It looks like I need to set for some correct projection
somewhere.
> > > Could
> > > > > you give a bit of guideline (in the while, I'll try to read
the
> > > > > Met_Users_Guide_v5.1.pdf in more detail; may send update
again
> before
> > > end
> > > > > of the day here)...
> > > > >
> > > > > Best wishes,
> > > > > Vanisa
> > > > >
> > > > > On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith <
> > vanisas at gmail.com
> > > >
> > > > > wrote:
> > > > >
> > > > > > Hi John,
> > > > > >
> > > > > > thank you for your reply; I recall your help earlier, and
I am
> only
> > > > > > gradually getting to know MET. I am grateful to the help
as to
> what
> > > > Julie
> > > > > > was responding patiently as well. I notice the
(presumably) mask
> > > files
> > > > > that
> > > > > > Julie helped generating for Nepal, and I try them.
However, the
> > > similar
> > > > > > error still shows up, and I am trying to understand why:
> > > > > >
> > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
\
> > > > > > > tutorial/out/ascii2nc/nepal.nc
> > > > tutorial/config/PointStatConfig_Nepal \
> > > > > > > -outdir tutorial/out/point_stat \
> > > > > > > -v 2
> > > > > > DEBUG 1: Default Config File:
> > > > > >
> > > >
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > > > > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> > > > > > WARNING:
> > > > > > WARNING: parse_conf_interp() -> For neareast neighbor
> interpolation
> > > > > > method, resetting method from "UW_MEAN" to "NEAREST" since
the
> > > > > > interpolation width is 1.
> > > > > > WARNING:
> > > > > > GSL_RNG_TYPE=mt19937
> > > > > > GSL_RNG_SEED=1682548822
> > > > > > DEBUG 1: Forecast File:
> > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > > > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
> > > > > > DEBUG 2:
> > > > > > DEBUG 2:
> > > > > >
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > > DEBUG 2:
> > > > > > DEBUG 2: Reading data for SPFH/P500.
> > > > > > ERROR :
> > > > > > ERROR : parse_grid_mask() -> the mask_grid requested
"nepal" is
> > not
> > > > > > defined.
> > > > > > ERROR :
> > > > > >
> > > > > > Kindly see in the surapipith_data folder in rap.ucar.edu
site
> the
> > > > latest
> > > > > > config file, as well as the suite of Nepal mask files that
Julie
> > made
> > > > > and I
> > > > > > use. Hope some progress would be before the end of the
day, or
> once
> > > > some
> > > > > > hints are next sent from there...
> > > > > >
> > > > > > Best wishes,
> > > > > > Vanisa
> > > > > >
> > > > > >
> > > > > > On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway via
RT <
> > > > > > met_help at ucar.edu> wrote:
> > > > > >
> > > > > >> Hello Vanisa,
> > > > > >>
> > > > > >> My name is John, and I work with Julie who is out of the
office
> > > today.
> > > > > I
> > > > > >> see that you're having some trouble getting Point-Stat to
run:
> > > > > >>
> > > > > >> ERROR : parse_grid_mask() -> the mask_grid requested
"nepal"
> > is
> > > > not
> > > > > >> defined.
> > > > > >>
> > > > > >> This should be pretty easy to resolve. You've worked
with Julie
> > to
> > > > > >> generate a mask file for Nepal, and that mask is
contained in
> some
> > > > > NetCDF
> > > > > >> file (the output of the gen_vx_mask tool). For example,
let's
> say
> > > > it's
> > > > > in
> > > > > >> "/home/vsurapipith/nepal_mask.nc". You'll need to put
the path
> > to
> > > > that
> > > > > >> file in the "mask" section of the Point-Stat
configuration file:
> > > > > >>
> > > > > >> mask = {
> > > > > >> grid = [];
> > > > > >> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; // or
> > whatever
> > > > the
> > > > > >> actual path is to that data file
> > > > > >> sid = [];
> > > > > >> }
> > > > > >>
> > > > > >> And then try running Point-Stat again. For a more
complete
> > > > description
> > > > > of
> > > > > >> the masking options available, please see pages 54-55 of
the MET
> > > > User's
> > > > > >> Guide:
> > > > > >>
> > > > > >>
> > > > > >>
> > > > >
> > > >
> > >
> >
>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
> > > > > >>
> > > > > >> Please let me know how it goes.
> > > > > >>
> > > > > >> Thanks,
> > > > > >> John
> > > > > >>
> > > > > >>
> > > > > >>
> > > > > >> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via RT
<
> > > > > >> met_help at ucar.edu
> > > > > >> > wrote:
> > > > > >>
> > > > > >> >
> > > > > >> > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> >
> > > > > >> >
> > > > > >> > Hi Julie,
> > > > > >> >
> > > > > >> > I'm back....after a couple of weeks spending on
stumbling over
> > the
> > > > > >> > installation of UPP until I got it through, now I have
the
> GRIB1
> > > > files
> > > > > >> > systematically converted from the WRF output generated
here.
> > > > > >> >
> > > > > >> > I try the GRIB1 files with the PointStatConfig_Nepal
that I
> sent
> > > > > >> earlier to
> > > > > >> > do Point_Stat with our these WRF output, but currently
the
> error
> > > > shows
> > > > > >> up
> > > > > >> > as:
> > > > > >> >
> > > > > >> > ~/met-5.1_bugfix$ ./bin/point_stat
> > > > > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > > >> > tutorial/out/ascii2nc/nepal.nc
> > > > tutorial/config/PointStatConfig_Nepal
> > > > > >> > -outdir tutorial/out/point_stat -v 2
> > > > > >> > DEBUG 1: Default Config File:
> > > > > >> >
> > > > > >>
> > > > >
> > >
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > > > > >> > DEBUG 1: User Config File:
> tutorial/config/PointStatConfig_Nepal
> > > > > >> > WARNING:
> > > > > >> > WARNING: parse_conf_interp() -> For neareast neighbor
> > > interpolation
> > > > > >> method,
> > > > > >> > resetting method from "UW_MEAN" to "NEAREST" since the
> > > interpolation
> > > > > >> width
> > > > > >> > is 1.
> > > > > >> > WARNING:
> > > > > >> > GSL_RNG_TYPE=mt19937
> > > > > >> > GSL_RNG_SEED=18446744072070604630
> > > > > >> > DEBUG 1: Forecast File:
> > > > > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > > > > >> > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal.nc
> > > > > >> > DEBUG 2:
> > > > > >> > DEBUG 2:
> > > > > >> >
> > > > > >> >
> > > > > >>
> > > > >
> > > >
> > >
> >
>
--------------------------------------------------------------------------------
> > > > > >> > DEBUG 2:
> > > > > >> > DEBUG 2: Reading data for SPFH/P500.
> > > > > >> > ERROR :
> > > > > >> > ERROR : parse_grid_mask() -> the mask_grid requested
"nepal"
> is
> > > not
> > > > > >> > defined.
> > > > > >> > ERROR :
> > > > > >> >
> > > > > >> > It looks to me like the mask file - nepal.nc, that I
> > deliberately
> > > > > took
> > > > > >> > from
> > > > > >> > what you generated earlier - was not seen/defined well
by MET
> > > here.
> > > > I
> > > > > >> have
> > > > > >> > actually just put the WRF output files in GRIB1 format
> (checked
> > > with
> > > > > as
> > > > > >> you
> > > > > >> > recommend) on the rap.ucar.edu ftp site
> > (surapipith_data)...Would
> > > > you
> > > > > >> > kindly let me know what should be done to correct this,
> please?
> > > > > >> >
> > > > > >> > Best wishes,
> > > > > >> > Vanisa
> > > > > >> >
> > > > > >> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via
RT <
> > > > > >> > met_help at ucar.edu>
> > > > > >> > wrote:
> > > > > >> >
> > > > > >> > > Hi Vanisa.
> > > > > >> > >
> > > > > >> > > While I don't have any experience installing UPP
myself, I
> did
> > > > take
> > > > > a
> > > > > >> > look
> > > > > >> > > at your log file where I saw:
> > > > > >> > >
> > > > > >> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h: No
such
> > file
> > > or
> > > > > >> > > directory compilation terminated.
> > > > > >> > >
> > > > > >> > > So, it looks like UPP is having trouble finding
jasper.h.
> > > > > >> > >
> > > > > >> > > Support for UPP is provided via wrfhelp at ucar.edu.
If
> > you're
> > > > > >> unable
> > > > > >> > to
> > > > > >> > > resolve the jasper problem, I would suggest writing
to them
> > > > > regarding
> > > > > >> > this
> > > > > >> > > issue.
> > > > > >> > >
> > > > > >> > > Thanks,
> > > > > >> > > Julie
> > > > > >> > >
> > > > > >> > >
> > > > > >> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com via
RT <
> > > > > >> > > met_help at ucar.edu>
> > > > > >> > > wrote:
> > > > > >> > >
> > > > > >> > > >
> > > > > >> > > > <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > > >
> > > > > >> > > >
> > > > > >> > > > Hi Julie,
> > > > > >> > > >
> > > > > >> > > > Happy New Year...:) Thank you very much for your
reply
> upon
> > > back
> > > > > to
> > > > > >> > > office.
> > > > > >> > > > Indeed, I did not use UPP for processing the WRF
(ARW)
> > output.
> > > > > >> > Therefore,
> > > > > >> > > > today I look into how to install UPP on my
workstation for
> > > this
> > > > > >> > purpose.
> > > > > >> > > > However, there is a problem that unipost.exe was
not
> > generated
> > > > > after
> > > > > >> > > > compilation, while other files came through. I was
> wondered
> > if
> > > > it
> > > > > >> has
> > > > > >> > to
> > > > > >> > > do
> > > > > >> > > > with liking to JASPER. I have not been able to spot
why,
> and
> > > now
> > > > > >> attach
> > > > > >> > > the
> > > > > >> > > > compile_upp.log to you here, in case you might be
able to
> > see
> > > > what
> > > > > >> is
> > > > > >> > > wrong
> > > > > >> > > > in there.
> > > > > >> > > >
> > > > > >> > > > Thanks again; update more soon.
> > > > > >> > > >
> > > > > >> > > > Best wishes,
> > > > > >> > > > Vanisa
> > > > > >> > > >
> > > > > >> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik
via RT <
> > > > > >> > > met_help at ucar.edu
> > > > > >> > > > >
> > > > > >> > > > wrote:
> > > > > >> > > >
> > > > > >> > > > > Hi Vanisa. The file you uploaded may have aged
off of
> our
> > > ftp
> > > > > >> > site. I
> > > > > >> > > > do
> > > > > >> > > > > not see it there currently.
> > > > > >> > > > >
> > > > > >> > > > > I neglected to ask you if your WRF output file
was
> > > > > post-processed?
> > > > > >> > MET
> > > > > >> > > > is
> > > > > >> > > > > unable to read the output of WRF directly.
Instead, you
> > > must
> > > > > >> first
> > > > > >> > > > > post-process it using either the Unified
PostProcessor
> > (UPP)
> > > > or
> > > > > >> the
> > > > > >> > > > pinterp
> > > > > >> > > > > tool (for ARW only). We strongly recommend using
UPP
> > since
> > > > > MET's
> > > > > >> > > support
> > > > > >> > > > > is most robust for the GRIB1 format that UPP
creates.
> > > > > >> > > > >
> > > > > >> > > > > If you have already post-processed your file,
please
> > trying
> > > > > >> uploading
> > > > > >> > > to
> > > > > >> > > > > our ftp site again. If you have not, please
> post-process
> > > the
> > > > > file
> > > > > >> > and
> > > > > >> > > > give
> > > > > >> > > > > it another try. Thanks!
> > > > > >> > > > >
> > > > > >> > > > > Happy New Year!
> > > > > >> > > > >
> > > > > >> > > > > Julie
> > > > > >> > > > >
> > > >
> > >
> > >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: point_stat for WRF output question
From: Julie Prestopnik
Time: Mon Feb 08 15:58:35 2016
Hi Vanisa. I was able to speak with John today. He helped me better
understand why you are not getting the desired output.
Looking in the nepal_ascii.nc, you can see in the "obs_arr" that the
second
column (grib code) contains only values of 11. That grib code is for
2-m
TMP. This explains why point_stat was showing us that for all other
fields
there were rejected observations due to the grib code, since the
observation file only contains observations for the 2-m TMP:
DEBUG 3: Observations processed = 33816
DEBUG 3: Rejected: SID exclusion = 0
DEBUG 3: Rejected: GRIB code = 33816
With regard to verifying the 2-m TMP, I saw that the observations were
rejected due to the valid time:
DEBUG 3: Observations processed = 33816
DEBUG 3: Rejected: SID exclusion = 0
DEBUG 3: Rejected: GRIB code = 0
DEBUG 3: Rejected: valid time = 33816
I first ran "wgrib WRFPRS_d02.00 | grep TMP" to get the record number
for
the 2-m TMP:
...
296:22645020:d=13010118:TMP:kpds5=11:kpds6=105:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:2
m above gnd:anl:NAve=0
...
Then I ran "wgrib -d 296 -V WRFPRS_d02.00" to get the date/time
information
for the 2-m TMP and got:
rec 296:22645020:date 2013010118 TMP kpds5=11 kpds6=105 kpds7=2
levels=(0,2) grid=255 2 m above gnd anl:
TMP=Temp. [K]
timerange 0 P1 0 P2 0 TimeU 1 nx 270 ny 180 GDS grid 1 num_in_ave 0
missing 0
center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
Mercator: lat 23.937000 to 31.070000 by 5.000000 km nxny 48600
long 77.740000 to 89.838000 by 5.000000 km, (270 x 180) scan
64
mode 136 Latin 0.000000 bdsgrid 1
min/max data 238.877 291.292 num bits 16 BDS_Ref 238877 DecScale
3
BinScale 0
You can see the date/time is 2013010118. Looking for that time in the
observation file, nepal_ascii.nc produced no results. However, there
were
results for 2013010100, 2013010106, and 2013010112. So, as we can
see,
just as the debug statement showed, the observations were rejected
because
the valid time didn't match up. Now you wouldn't necessarily want to
verify the 18 hour with the 12 hour, 6 hour, or 0 hour, but to show
that
you will get MET output, you could extend the obs_window "beg" and/or
"end"
values in this case and you would get some matched pair output.
Also, to answer your question "Do I need to generate a GRIB1 file that
combine all the required hours...otherwise, can MET read in more than
one
GRIB1 files at a time for Point_Stat analysis?" Most users use a grib
file
for each time and run point stat once for each separate output time.
I hope this helps.
Please let us know if you have additional questions.
Thanks,
Julie
On Fri, Feb 5, 2016 at 10:31 AM, Julie Prestopnik <jpresto at ucar.edu>
wrote:
> Hi Vanisa. I apologize for the delayed response. I was out of the
office
> yesterday.
>
> Unfortunately, I need some help from John to figure out what exactly
is
> going on here, but he is out of the office today. We will revisit
this on
> Monday and get back to you as soon as we can.
>
> Thank you,
> Julie
>
> On Thu, Feb 4, 2016 at 5:47 AM, vanisas at gmail.com via RT <
> met_help at ucar.edu> wrote:
>
>>
>> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>>
>> Hi Julie,
>>
>> I have just uploaded the 4 files I used for running point_stat that
>> generated the output files that I put there earlier. Hope these are
>> comprehensible. the command lines that I used are as follows:
>>
>> ~/met-5.1_bugfix$ ./bin/point_stat \
>>
>> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
>>
>> > tutorial/out/ascii2nc/nepal_ascii.nc \
>>
>> > tutorial/config/PointStatConfig_Nepal \
>>
>> > -outdir tutorial/out/point_stat \
>>
>> > -v 2
>>
>>
>> Thanks alot for your kind help. Wish to see this through soon.
>>
>>
>> Best wishes,
>> Vanisa
>>
>> On Wed, Feb 3, 2016 at 11:41 PM, Julie Prestopnik via RT <
>> met_help at ucar.edu>
>> wrote:
>>
>> > Hi Vanisa.
>> >
>> > Could you please upload your input forecast file, input
observation
>> file,
>> > masking file, and your configuration file? Also, could you
please send
>> me
>> > the command line argument for point_stat that you are running?
>> >
>> > I can take a closer look and let you know what I find out.
>> >
>> > Thanks,
>> > Julie
>> >
>> > On Tue, Feb 2, 2016 at 2:32 AM, vanisas at gmail.com via RT <
>> > met_help at ucar.edu>
>> > wrote:
>> >
>> > >
>> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>> > >
>> > > Hi Julie,
>> > >
>> > > thank you for your reply last week; apologies for my delay in
>> reporting.
>> > > What you suggest works. I can now generate the output files
using
>> > > Point_Stat. However, I do not yet understand how to interpret
or make
>> use
>> > > of the output...I have been trying to read what described in
the
>> > > MET_User_Guide. Frankly speaking, I still do not get it. Could
you
>> kindly
>> > > see in my folder in the rap.ucar.edu ftp site in
>> incoming/irap/met_help
>> > > under surapipith_data. I attached here some of the files. It
looks
>> like
>> > > only the headers appear...or what I should do to generate the
>> statistical
>> > > evaluation output. What surprise me was that the example files
I got
>> from
>> > > the US example (follow the online tutorial) look similar. Could
you
>> guide
>> > > me where to read to comprehend MET?
>> > >
>> > > One other question I have is that probably this is because I
use only
>> one
>> > > hour of WRF output for this try. Do I need to generate a GRIB1
file
>> that
>> > > combine all the required hours...otherwise, can MET read in
more than
>> one
>> > > GRIB1 files at a time for Point_Stat analysis?
>> > >
>> > > Best wishes,
>> > > Vanisa
>> > >
>> > > On Thu, Jan 28, 2016 at 12:00 AM, Julie Prestopnik via RT <
>> > > met_help at ucar.edu
>> > > > wrote:
>> > >
>> > > > Hi Vanisa. This is Julie again.
>> > > >
>> > > > I see that you are encountering a problem due to the masking
and
>> > > > verification grids not matching.
>> > > >
>> > > > One projection is a Lat/Lon grid and the other is a Mercator.
>> > > >
>> > > > The mask file I created for you using gen_vx_mask used a
Lat/Lon
>> grid
>> > > file
>> > > > (gfs_2012040900_F003.grib). If you would like to get your
netcdf
>> mask
>> > > > file
>> > > > on a Mercator grid, you could run gen_vx_mask using your file
that
>> > > contains
>> > > > the grid that you would like to verify, which I believe is on
a
>> > Mercator
>> > > > projection (WRFPRS_d02.00?) and the nepal poly file which
will then
>> > > produce
>> > > > a netcdf file on that same projection.
>> > > >
>> > > > So, for example, something like:
>> > > >
>> > > > gen_vx_mask WRFPRS_d02.00 nepal_detail.poly nepal_detail.nc
>> > > >
>> > > > Then, you can rerun point_stat using the new netcdf mask
file on
>> the
>> > > > Mercator grid and that should fix the problem.
>> > > >
>> > > > However, if you're not going to be running thousands of cases
of
>> Point
>> > > > Stat, you don't necessarily need to run gen_vx_mask, but
rather you
>> > could
>> > > > use the poly file directly with Point Stat in the
configuration
>> file in
>> > > the
>> > > > mask section that John described earlier. Using the poly
file
>> directly
>> > > can
>> > > > be slow, so if you're running point_stat on a large number of
files,
>> > this
>> > > > is likely not the best way to use point_stat.
>> > > >
>> > > > I hope this helps. Please let me know if you have any
questions.
>> > > >
>> > > > Thanks,
>> > > > Julie
>> > > >
>> > > > On Tue, Jan 26, 2016 at 11:04 PM, vanisas at gmail.com via RT <
>> > > > met_help at ucar.edu> wrote:
>> > > >
>> > > > >
>> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>> > > > >
>> > > > > Hi John,
>> > > > >
>> > > > > apologies, I just notice that I did not remove the word
"nepal" in
>> > the
>> > > > > bracket ([]) behind grid = within mask section in the
PointConfig
>> > > file. I
>> > > > > have then removed that, and found a different error (also
rename
>> the
>> > > > ascii
>> > > > > file containing the ADPUPA observation from our stations in
Nepal
>> to
>> > be
>> > > > > called "nepal_ascii.nc" to avoid confusion) :
>> > > > >
>> > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
>> > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
>> > > > > > tutorial/out/ascii2nc/nepal_ascii.nc \
>> > > > > > tutorial/config/PointStatConfig_Nepal \
>> > > > > > -outdir tutorial/out/point_stat \
>> > > > > > -v 2
>> > > > > DEBUG 1: Default Config File:
>> > > > >
>> > >
>> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
>> > > > > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
>> > > > > WARNING:
>> > > > > WARNING: parse_conf_interp() -> For neareast neighbor
>> interpolation
>> > > > method,
>> > > > > resetting method from "UW_MEAN" to "NEAREST" since the
>> interpolation
>> > > > width
>> > > > > is 1.
>> > > > > WARNING:
>> > > > > GSL_RNG_TYPE=mt19937
>> > > > > GSL_RNG_SEED=18446744072507402326
>> > > > > DEBUG 1: Forecast File:
>> > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
>> > > > > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal_ascii.nc
>> > > > > DEBUG 2:
>> > > > > DEBUG 2:
>> > > > >
>> > > > >
>> > > >
>> > >
>> >
>>
--------------------------------------------------------------------------------
>> > > > > DEBUG 2:
>> > > > > DEBUG 2: Reading data for SPFH/P500.
>> > > > > ERROR :
>> > > > > ERROR : parse_poly_mask() -> The masking and verification
grids
>> do
>> > not
>> > > > > match:
>> > > > > ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll:
-90.000
>> lon_ll:
>> > > > -0.000
>> > > > > delta_lat: 0.500 delta_lon: 0.500 !=
>> > > > > ERROR : Projection: Mercator Nx: 270 Ny: 180
Lat_LL_radians:
>> 0.4178
>> > > > > Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423
Lon_UR_radians:
>> > -1.5680
>> > > > Mx:
>> > > > > 1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
>> > > > > ERROR :
>> > > > >
>> > > > > It looks like I need to set for some correct projection
somewhere.
>> > > Could
>> > > > > you give a bit of guideline (in the while, I'll try to read
the
>> > > > > Met_Users_Guide_v5.1.pdf in more detail; may send update
again
>> before
>> > > end
>> > > > > of the day here)...
>> > > > >
>> > > > > Best wishes,
>> > > > > Vanisa
>> > > > >
>> > > > > On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith <
>> > vanisas at gmail.com
>> > > >
>> > > > > wrote:
>> > > > >
>> > > > > > Hi John,
>> > > > > >
>> > > > > > thank you for your reply; I recall your help earlier, and
I am
>> only
>> > > > > > gradually getting to know MET. I am grateful to the help
as to
>> what
>> > > > Julie
>> > > > > > was responding patiently as well. I notice the
(presumably) mask
>> > > files
>> > > > > that
>> > > > > > Julie helped generating for Nepal, and I try them.
However, the
>> > > similar
>> > > > > > error still shows up, and I am trying to understand why:
>> > > > > >
>> > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
>> > > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
\
>> > > > > > > tutorial/out/ascii2nc/nepal.nc
>> > > > tutorial/config/PointStatConfig_Nepal \
>> > > > > > > -outdir tutorial/out/point_stat \
>> > > > > > > -v 2
>> > > > > > DEBUG 1: Default Config File:
>> > > > > >
>> > > >
>> >
>> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
>> > > > > > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
>> > > > > > WARNING:
>> > > > > > WARNING: parse_conf_interp() -> For neareast neighbor
>> interpolation
>> > > > > > method, resetting method from "UW_MEAN" to "NEAREST"
since the
>> > > > > > interpolation width is 1.
>> > > > > > WARNING:
>> > > > > > GSL_RNG_TYPE=mt19937
>> > > > > > GSL_RNG_SEED=1682548822
>> > > > > > DEBUG 1: Forecast File:
>> > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
>> > > > > > DEBUG 1: Observation File: tutorial/out/ascii2nc/nepal.nc
>> > > > > > DEBUG 2:
>> > > > > > DEBUG 2:
>> > > > > >
>> > > > >
>> > > >
>> > >
>> >
>>
--------------------------------------------------------------------------------
>> > > > > > DEBUG 2:
>> > > > > > DEBUG 2: Reading data for SPFH/P500.
>> > > > > > ERROR :
>> > > > > > ERROR : parse_grid_mask() -> the mask_grid requested
"nepal" is
>> > not
>> > > > > > defined.
>> > > > > > ERROR :
>> > > > > >
>> > > > > > Kindly see in the surapipith_data folder in rap.ucar.edu
site
>> the
>> > > > latest
>> > > > > > config file, as well as the suite of Nepal mask files
that Julie
>> > made
>> > > > > and I
>> > > > > > use. Hope some progress would be before the end of the
day, or
>> once
>> > > > some
>> > > > > > hints are next sent from there...
>> > > > > >
>> > > > > > Best wishes,
>> > > > > > Vanisa
>> > > > > >
>> > > > > >
>> > > > > > On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway via
RT <
>> > > > > > met_help at ucar.edu> wrote:
>> > > > > >
>> > > > > >> Hello Vanisa,
>> > > > > >>
>> > > > > >> My name is John, and I work with Julie who is out of the
office
>> > > today.
>> > > > > I
>> > > > > >> see that you're having some trouble getting Point-Stat
to run:
>> > > > > >>
>> > > > > >> ERROR : parse_grid_mask() -> the mask_grid requested
>> "nepal"
>> > is
>> > > > not
>> > > > > >> defined.
>> > > > > >>
>> > > > > >> This should be pretty easy to resolve. You've worked
with
>> Julie
>> > to
>> > > > > >> generate a mask file for Nepal, and that mask is
contained in
>> some
>> > > > > NetCDF
>> > > > > >> file (the output of the gen_vx_mask tool). For example,
let's
>> say
>> > > > it's
>> > > > > in
>> > > > > >> "/home/vsurapipith/nepal_mask.nc". You'll need to put
the
>> path
>> > to
>> > > > that
>> > > > > >> file in the "mask" section of the Point-Stat
configuration
>> file:
>> > > > > >>
>> > > > > >> mask = {
>> > > > > >> grid = [];
>> > > > > >> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; //
or
>> > whatever
>> > > > the
>> > > > > >> actual path is to that data file
>> > > > > >> sid = [];
>> > > > > >> }
>> > > > > >>
>> > > > > >> And then try running Point-Stat again. For a more
complete
>> > > > description
>> > > > > of
>> > > > > >> the masking options available, please see pages 54-55 of
the
>> MET
>> > > > User's
>> > > > > >> Guide:
>> > > > > >>
>> > > > > >>
>> > > > > >>
>> > > > >
>> > > >
>> > >
>> >
>>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
>> > > > > >>
>> > > > > >> Please let me know how it goes.
>> > > > > >>
>> > > > > >> Thanks,
>> > > > > >> John
>> > > > > >>
>> > > > > >>
>> > > > > >>
>> > > > > >> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via
RT <
>> > > > > >> met_help at ucar.edu
>> > > > > >> > wrote:
>> > > > > >>
>> > > > > >> >
>> > > > > >> > <URL:
>> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>> > > > > >> >
>> > > > > >> > Hi Julie,
>> > > > > >> >
>> > > > > >> > I'm back....after a couple of weeks spending on
stumbling
>> over
>> > the
>> > > > > >> > installation of UPP until I got it through, now I have
the
>> GRIB1
>> > > > files
>> > > > > >> > systematically converted from the WRF output generated
here.
>> > > > > >> >
>> > > > > >> > I try the GRIB1 files with the PointStatConfig_Nepal
that I
>> sent
>> > > > > >> earlier to
>> > > > > >> > do Point_Stat with our these WRF output, but currently
the
>> error
>> > > > shows
>> > > > > >> up
>> > > > > >> > as:
>> > > > > >> >
>> > > > > >> > ~/met-5.1_bugfix$ ./bin/point_stat
>> > > > > >> > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00
>> > > > > >> > tutorial/out/ascii2nc/nepal.nc
>> > > > tutorial/config/PointStatConfig_Nepal
>> > > > > >> > -outdir tutorial/out/point_stat -v 2
>> > > > > >> > DEBUG 1: Default Config File:
>> > > > > >> >
>> > > > > >>
>> > > > >
>> > >
>> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
>> > > > > >> > DEBUG 1: User Config File:
>> tutorial/config/PointStatConfig_Nepal
>> > > > > >> > WARNING:
>> > > > > >> > WARNING: parse_conf_interp() -> For neareast neighbor
>> > > interpolation
>> > > > > >> method,
>> > > > > >> > resetting method from "UW_MEAN" to "NEAREST" since the
>> > > interpolation
>> > > > > >> width
>> > > > > >> > is 1.
>> > > > > >> > WARNING:
>> > > > > >> > GSL_RNG_TYPE=mt19937
>> > > > > >> > GSL_RNG_SEED=18446744072070604630
>> > > > > >> > DEBUG 1: Forecast File:
>> > > > > >> > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00
>> > > > > >> > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal.nc
>> > > > > >> > DEBUG 2:
>> > > > > >> > DEBUG 2:
>> > > > > >> >
>> > > > > >> >
>> > > > > >>
>> > > > >
>> > > >
>> > >
>> >
>>
--------------------------------------------------------------------------------
>> > > > > >> > DEBUG 2:
>> > > > > >> > DEBUG 2: Reading data for SPFH/P500.
>> > > > > >> > ERROR :
>> > > > > >> > ERROR : parse_grid_mask() -> the mask_grid requested
>> "nepal" is
>> > > not
>> > > > > >> > defined.
>> > > > > >> > ERROR :
>> > > > > >> >
>> > > > > >> > It looks to me like the mask file - nepal.nc, that I
>> > deliberately
>> > > > > took
>> > > > > >> > from
>> > > > > >> > what you generated earlier - was not seen/defined well
by MET
>> > > here.
>> > > > I
>> > > > > >> have
>> > > > > >> > actually just put the WRF output files in GRIB1 format
>> (checked
>> > > with
>> > > > > as
>> > > > > >> you
>> > > > > >> > recommend) on the rap.ucar.edu ftp site
>> > (surapipith_data)...Would
>> > > > you
>> > > > > >> > kindly let me know what should be done to correct
this,
>> please?
>> > > > > >> >
>> > > > > >> > Best wishes,
>> > > > > >> > Vanisa
>> > > > > >> >
>> > > > > >> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik via
RT <
>> > > > > >> > met_help at ucar.edu>
>> > > > > >> > wrote:
>> > > > > >> >
>> > > > > >> > > Hi Vanisa.
>> > > > > >> > >
>> > > > > >> > > While I don't have any experience installing UPP
myself, I
>> did
>> > > > take
>> > > > > a
>> > > > > >> > look
>> > > > > >> > > at your log file where I saw:
>> > > > > >> > >
>> > > > > >> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h:
No such
>> > file
>> > > or
>> > > > > >> > > directory compilation terminated.
>> > > > > >> > >
>> > > > > >> > > So, it looks like UPP is having trouble finding
jasper.h.
>> > > > > >> > >
>> > > > > >> > > Support for UPP is provided via wrfhelp at ucar.edu.
If
>> > you're
>> > > > > >> unable
>> > > > > >> > to
>> > > > > >> > > resolve the jasper problem, I would suggest writing
to them
>> > > > > regarding
>> > > > > >> > this
>> > > > > >> > > issue.
>> > > > > >> > >
>> > > > > >> > > Thanks,
>> > > > > >> > > Julie
>> > > > > >> > >
>> > > > > >> > >
>> > > > > >> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com
via RT <
>> > > > > >> > > met_help at ucar.edu>
>> > > > > >> > > wrote:
>> > > > > >> > >
>> > > > > >> > > >
>> > > > > >> > > > <URL:
>> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
>> > > >
>> > > > > >> > > >
>> > > > > >> > > > Hi Julie,
>> > > > > >> > > >
>> > > > > >> > > > Happy New Year...:) Thank you very much for your
reply
>> upon
>> > > back
>> > > > > to
>> > > > > >> > > office.
>> > > > > >> > > > Indeed, I did not use UPP for processing the WRF
(ARW)
>> > output.
>> > > > > >> > Therefore,
>> > > > > >> > > > today I look into how to install UPP on my
workstation
>> for
>> > > this
>> > > > > >> > purpose.
>> > > > > >> > > > However, there is a problem that unipost.exe was
not
>> > generated
>> > > > > after
>> > > > > >> > > > compilation, while other files came through. I was
>> wondered
>> > if
>> > > > it
>> > > > > >> has
>> > > > > >> > to
>> > > > > >> > > do
>> > > > > >> > > > with liking to JASPER. I have not been able to
spot why,
>> and
>> > > now
>> > > > > >> attach
>> > > > > >> > > the
>> > > > > >> > > > compile_upp.log to you here, in case you might be
able to
>> > see
>> > > > what
>> > > > > >> is
>> > > > > >> > > wrong
>> > > > > >> > > > in there.
>> > > > > >> > > >
>> > > > > >> > > > Thanks again; update more soon.
>> > > > > >> > > >
>> > > > > >> > > > Best wishes,
>> > > > > >> > > > Vanisa
>> > > > > >> > > >
>> > > > > >> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik
via RT <
>> > > > > >> > > met_help at ucar.edu
>> > > > > >> > > > >
>> > > > > >> > > > wrote:
>> > > > > >> > > >
>> > > > > >> > > > > Hi Vanisa. The file you uploaded may have aged
off of
>> our
>> > > ftp
>> > > > > >> > site. I
>> > > > > >> > > > do
>> > > > > >> > > > > not see it there currently.
>> > > > > >> > > > >
>> > > > > >> > > > > I neglected to ask you if your WRF output file
was
>> > > > > post-processed?
>> > > > > >> > MET
>> > > > > >> > > > is
>> > > > > >> > > > > unable to read the output of WRF directly.
Instead,
>> you
>> > > must
>> > > > > >> first
>> > > > > >> > > > > post-process it using either the Unified
PostProcessor
>> > (UPP)
>> > > > or
>> > > > > >> the
>> > > > > >> > > > pinterp
>> > > > > >> > > > > tool (for ARW only). We strongly recommend
using UPP
>> > since
>> > > > > MET's
>> > > > > >> > > support
>> > > > > >> > > > > is most robust for the GRIB1 format that UPP
creates.
>> > > > > >> > > > >
>> > > > > >> > > > > If you have already post-processed your file,
please
>> > trying
>> > > > > >> uploading
>> > > > > >> > > to
>> > > > > >> > > > > our ftp site again. If you have not, please
>> post-process
>> > > the
>> > > > > file
>> > > > > >> > and
>> > > > > >> > > > give
>> > > > > >> > > > > it another try. Thanks!
>> > > > > >> > > > >
>> > > > > >> > > > > Happy New Year!
>> > > > > >> > > > >
>> > > > > >> > > > > Julie
>> > > > > >> > > > >
>> > > >
>> > >
>> > >
>> >
>> >
>> > --
>> > Julie Prestopnik
>> > National Center for Atmospheric Research
>> > Research Applications Laboratory
>> > Phone: 303.497.8399
>> > Email: jpresto at ucar.edu
>> >
>> >
>>
>>
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
Subject: point_stat for WRF output question
From: vanisas at gmail.com
Time: Wed Feb 10 04:01:54 2016
Hi Julie and John,
thank you so much for your help identifying the problem and giving
detailed
explanation. I got this solved now with adjusting the obs_window. I am
working into understanding the analysis using MET, and may send next
questions on another topic soon.
Best wishes,
Vanisa
On Tue, Feb 9, 2016 at 4:43 AM, Julie Prestopnik via RT
<met_help at ucar.edu>
wrote:
> Hi Vanisa. I was able to speak with John today. He helped me better
> understand why you are not getting the desired output.
>
> Looking in the nepal_ascii.nc, you can see in the "obs_arr" that the
> second
> column (grib code) contains only values of 11. That grib code is
for 2-m
> TMP. This explains why point_stat was showing us that for all other
fields
> there were rejected observations due to the grib code, since the
> observation file only contains observations for the 2-m TMP:
>
> DEBUG 3: Observations processed = 33816
> DEBUG 3: Rejected: SID exclusion = 0
> DEBUG 3: Rejected: GRIB code = 33816
>
> With regard to verifying the 2-m TMP, I saw that the observations
were
> rejected due to the valid time:
>
> DEBUG 3: Observations processed = 33816
> DEBUG 3: Rejected: SID exclusion = 0
> DEBUG 3: Rejected: GRIB code = 0
> DEBUG 3: Rejected: valid time = 33816
>
> I first ran "wgrib WRFPRS_d02.00 | grep TMP" to get the record
number for
> the 2-m TMP:
> ...
>
>
296:22645020:d=13010118:TMP:kpds5=11:kpds6=105:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:2
> m above gnd:anl:NAve=0
> ...
>
> Then I ran "wgrib -d 296 -V WRFPRS_d02.00" to get the date/time
information
> for the 2-m TMP and got:
>
> rec 296:22645020:date 2013010118 TMP kpds5=11 kpds6=105 kpds7=2
> levels=(0,2) grid=255 2 m above gnd anl:
> TMP=Temp. [K]
> timerange 0 P1 0 P2 0 TimeU 1 nx 270 ny 180 GDS grid 1 num_in_ave
0
> missing 0
> center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
> Mercator: lat 23.937000 to 31.070000 by 5.000000 km nxny 48600
> long 77.740000 to 89.838000 by 5.000000 km, (270 x 180)
scan 64
> mode 136 Latin 0.000000 bdsgrid 1
> min/max data 238.877 291.292 num bits 16 BDS_Ref 238877
DecScale 3
> BinScale 0
>
>
> You can see the date/time is 2013010118. Looking for that time in
the
> observation file, nepal_ascii.nc produced no results. However,
there were
> results for 2013010100, 2013010106, and 2013010112. So, as we can
see,
> just as the debug statement showed, the observations were rejected
because
> the valid time didn't match up. Now you wouldn't necessarily want
to
> verify the 18 hour with the 12 hour, 6 hour, or 0 hour, but to show
that
> you will get MET output, you could extend the obs_window "beg"
and/or "end"
> values in this case and you would get some matched pair output.
>
> Also, to answer your question "Do I need to generate a GRIB1 file
that
> combine all the required hours...otherwise, can MET read in more
than one
> GRIB1 files at a time for Point_Stat analysis?" Most users use a
grib file
> for each time and run point stat once for each separate output time.
>
> I hope this helps.
>
> Please let us know if you have additional questions.
>
> Thanks,
> Julie
>
>
> On Fri, Feb 5, 2016 at 10:31 AM, Julie Prestopnik <jpresto at ucar.edu>
> wrote:
>
> > Hi Vanisa. I apologize for the delayed response. I was out of
the
> office
> > yesterday.
> >
> > Unfortunately, I need some help from John to figure out what
exactly is
> > going on here, but he is out of the office today. We will revisit
this
> on
> > Monday and get back to you as soon as we can.
> >
> > Thank you,
> > Julie
> >
> > On Thu, Feb 4, 2016 at 5:47 AM, vanisas at gmail.com via RT <
> > met_help at ucar.edu> wrote:
> >
> >>
> >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >>
> >> Hi Julie,
> >>
> >> I have just uploaded the 4 files I used for running point_stat
that
> >> generated the output files that I put there earlier. Hope these
are
> >> comprehensible. the command lines that I used are as follows:
> >>
> >> ~/met-5.1_bugfix$ ./bin/point_stat \
> >>
> >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> >>
> >> > tutorial/out/ascii2nc/nepal_ascii.nc \
> >>
> >> > tutorial/config/PointStatConfig_Nepal \
> >>
> >> > -outdir tutorial/out/point_stat \
> >>
> >> > -v 2
> >>
> >>
> >> Thanks alot for your kind help. Wish to see this through soon.
> >>
> >>
> >> Best wishes,
> >> Vanisa
> >>
> >> On Wed, Feb 3, 2016 at 11:41 PM, Julie Prestopnik via RT <
> >> met_help at ucar.edu>
> >> wrote:
> >>
> >> > Hi Vanisa.
> >> >
> >> > Could you please upload your input forecast file, input
observation
> >> file,
> >> > masking file, and your configuration file? Also, could you
please
> send
> >> me
> >> > the command line argument for point_stat that you are running?
> >> >
> >> > I can take a closer look and let you know what I find out.
> >> >
> >> > Thanks,
> >> > Julie
> >> >
> >> > On Tue, Feb 2, 2016 at 2:32 AM, vanisas at gmail.com via RT <
> >> > met_help at ucar.edu>
> >> > wrote:
> >> >
> >> > >
> >> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
>
> >> > >
> >> > > Hi Julie,
> >> > >
> >> > > thank you for your reply last week; apologies for my delay in
> >> reporting.
> >> > > What you suggest works. I can now generate the output files
using
> >> > > Point_Stat. However, I do not yet understand how to interpret
or
> make
> >> use
> >> > > of the output...I have been trying to read what described in
the
> >> > > MET_User_Guide. Frankly speaking, I still do not get it.
Could you
> >> kindly
> >> > > see in my folder in the rap.ucar.edu ftp site in
> >> incoming/irap/met_help
> >> > > under surapipith_data. I attached here some of the files. It
looks
> >> like
> >> > > only the headers appear...or what I should do to generate the
> >> statistical
> >> > > evaluation output. What surprise me was that the example
files I got
> >> from
> >> > > the US example (follow the online tutorial) look similar.
Could you
> >> guide
> >> > > me where to read to comprehend MET?
> >> > >
> >> > > One other question I have is that probably this is because I
use
> only
> >> one
> >> > > hour of WRF output for this try. Do I need to generate a
GRIB1 file
> >> that
> >> > > combine all the required hours...otherwise, can MET read in
more
> than
> >> one
> >> > > GRIB1 files at a time for Point_Stat analysis?
> >> > >
> >> > > Best wishes,
> >> > > Vanisa
> >> > >
> >> > > On Thu, Jan 28, 2016 at 12:00 AM, Julie Prestopnik via RT <
> >> > > met_help at ucar.edu
> >> > > > wrote:
> >> > >
> >> > > > Hi Vanisa. This is Julie again.
> >> > > >
> >> > > > I see that you are encountering a problem due to the
masking and
> >> > > > verification grids not matching.
> >> > > >
> >> > > > One projection is a Lat/Lon grid and the other is a
Mercator.
> >> > > >
> >> > > > The mask file I created for you using gen_vx_mask used a
Lat/Lon
> >> grid
> >> > > file
> >> > > > (gfs_2012040900_F003.grib). If you would like to get your
netcdf
> >> mask
> >> > > > file
> >> > > > on a Mercator grid, you could run gen_vx_mask using your
file that
> >> > > contains
> >> > > > the grid that you would like to verify, which I believe is
on a
> >> > Mercator
> >> > > > projection (WRFPRS_d02.00?) and the nepal poly file which
will
> then
> >> > > produce
> >> > > > a netcdf file on that same projection.
> >> > > >
> >> > > > So, for example, something like:
> >> > > >
> >> > > > gen_vx_mask WRFPRS_d02.00 nepal_detail.poly nepal_detail.nc
> >> > > >
> >> > > > Then, you can rerun point_stat using the new netcdf mask
file on
> >> the
> >> > > > Mercator grid and that should fix the problem.
> >> > > >
> >> > > > However, if you're not going to be running thousands of
cases of
> >> Point
> >> > > > Stat, you don't necessarily need to run gen_vx_mask, but
rather
> you
> >> > could
> >> > > > use the poly file directly with Point Stat in the
configuration
> >> file in
> >> > > the
> >> > > > mask section that John described earlier. Using the poly
file
> >> directly
> >> > > can
> >> > > > be slow, so if you're running point_stat on a large number
of
> files,
> >> > this
> >> > > > is likely not the best way to use point_stat.
> >> > > >
> >> > > > I hope this helps. Please let me know if you have any
questions.
> >> > > >
> >> > > > Thanks,
> >> > > > Julie
> >> > > >
> >> > > > On Tue, Jan 26, 2016 at 11:04 PM, vanisas at gmail.com via RT
<
> >> > > > met_help at ucar.edu> wrote:
> >> > > >
> >> > > > >
> >> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >> > > > >
> >> > > > > Hi John,
> >> > > > >
> >> > > > > apologies, I just notice that I did not remove the word
"nepal"
> in
> >> > the
> >> > > > > bracket ([]) behind grid = within mask section in the
> PointConfig
> >> > > file. I
> >> > > > > have then removed that, and found a different error (also
rename
> >> the
> >> > > > ascii
> >> > > > > file containing the ADPUPA observation from our stations
in
> Nepal
> >> to
> >> > be
> >> > > > > called "nepal_ascii.nc" to avoid confusion) :
> >> > > > >
> >> > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> >> > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
\
> >> > > > > > tutorial/out/ascii2nc/nepal_ascii.nc \
> >> > > > > > tutorial/config/PointStatConfig_Nepal \
> >> > > > > > -outdir tutorial/out/point_stat \
> >> > > > > > -v 2
> >> > > > > DEBUG 1: Default Config File:
> >> > > > >
> >> > >
> >>
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> >> > > > > DEBUG 1: User Config File:
tutorial/config/PointStatConfig_Nepal
> >> > > > > WARNING:
> >> > > > > WARNING: parse_conf_interp() -> For neareast neighbor
> >> interpolation
> >> > > > method,
> >> > > > > resetting method from "UW_MEAN" to "NEAREST" since the
> >> interpolation
> >> > > > width
> >> > > > > is 1.
> >> > > > > WARNING:
> >> > > > > GSL_RNG_TYPE=mt19937
> >> > > > > GSL_RNG_SEED=18446744072507402326
> >> > > > > DEBUG 1: Forecast File:
> >> > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> >> > > > > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal_ascii.nc
> >> > > > > DEBUG 2:
> >> > > > > DEBUG 2:
> >> > > > >
> >> > > > >
> >> > > >
> >> > >
> >> >
> >>
>
--------------------------------------------------------------------------------
> >> > > > > DEBUG 2:
> >> > > > > DEBUG 2: Reading data for SPFH/P500.
> >> > > > > ERROR :
> >> > > > > ERROR : parse_poly_mask() -> The masking and
verification grids
> >> do
> >> > not
> >> > > > > match:
> >> > > > > ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll:
-90.000
> >> lon_ll:
> >> > > > -0.000
> >> > > > > delta_lat: 0.500 delta_lon: 0.500 !=
> >> > > > > ERROR : Projection: Mercator Nx: 270 Ny: 180
Lat_LL_radians:
> >> 0.4178
> >> > > > > Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423
Lon_UR_radians:
> >> > -1.5680
> >> > > > Mx:
> >> > > > > 1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
> >> > > > > ERROR :
> >> > > > >
> >> > > > > It looks like I need to set for some correct projection
> somewhere.
> >> > > Could
> >> > > > > you give a bit of guideline (in the while, I'll try to
read the
> >> > > > > Met_Users_Guide_v5.1.pdf in more detail; may send update
again
> >> before
> >> > > end
> >> > > > > of the day here)...
> >> > > > >
> >> > > > > Best wishes,
> >> > > > > Vanisa
> >> > > > >
> >> > > > > On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith <
> >> > vanisas at gmail.com
> >> > > >
> >> > > > > wrote:
> >> > > > >
> >> > > > > > Hi John,
> >> > > > > >
> >> > > > > > thank you for your reply; I recall your help earlier,
and I am
> >> only
> >> > > > > > gradually getting to know MET. I am grateful to the
help as to
> >> what
> >> > > > Julie
> >> > > > > > was responding patiently as well. I notice the
(presumably)
> mask
> >> > > files
> >> > > > > that
> >> > > > > > Julie helped generating for Nepal, and I try them.
However,
> the
> >> > > similar
> >> > > > > > error still shows up, and I am trying to understand
why:
> >> > > > > >
> >> > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> >> > > > > > > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> >> > > > > > > tutorial/out/ascii2nc/nepal.nc
> >> > > > tutorial/config/PointStatConfig_Nepal \
> >> > > > > > > -outdir tutorial/out/point_stat \
> >> > > > > > > -v 2
> >> > > > > > DEBUG 1: Default Config File:
> >> > > > > >
> >> > > >
> >> >
> >>
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> >> > > > > > DEBUG 1: User Config File:
> tutorial/config/PointStatConfig_Nepal
> >> > > > > > WARNING:
> >> > > > > > WARNING: parse_conf_interp() -> For neareast neighbor
> >> interpolation
> >> > > > > > method, resetting method from "UW_MEAN" to "NEAREST"
since the
> >> > > > > > interpolation width is 1.
> >> > > > > > WARNING:
> >> > > > > > GSL_RNG_TYPE=mt19937
> >> > > > > > GSL_RNG_SEED=1682548822
> >> > > > > > DEBUG 1: Forecast File:
> >> > > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> >> > > > > > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal.nc
> >> > > > > > DEBUG 2:
> >> > > > > > DEBUG 2:
> >> > > > > >
> >> > > > >
> >> > > >
> >> > >
> >> >
> >>
>
--------------------------------------------------------------------------------
> >> > > > > > DEBUG 2:
> >> > > > > > DEBUG 2: Reading data for SPFH/P500.
> >> > > > > > ERROR :
> >> > > > > > ERROR : parse_grid_mask() -> the mask_grid requested
"nepal"
> is
> >> > not
> >> > > > > > defined.
> >> > > > > > ERROR :
> >> > > > > >
> >> > > > > > Kindly see in the surapipith_data folder in
rap.ucar.edu site
> >> the
> >> > > > latest
> >> > > > > > config file, as well as the suite of Nepal mask files
that
> Julie
> >> > made
> >> > > > > and I
> >> > > > > > use. Hope some progress would be before the end of the
day, or
> >> once
> >> > > > some
> >> > > > > > hints are next sent from there...
> >> > > > > >
> >> > > > > > Best wishes,
> >> > > > > > Vanisa
> >> > > > > >
> >> > > > > >
> >> > > > > > On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway
via RT <
> >> > > > > > met_help at ucar.edu> wrote:
> >> > > > > >
> >> > > > > >> Hello Vanisa,
> >> > > > > >>
> >> > > > > >> My name is John, and I work with Julie who is out of
the
> office
> >> > > today.
> >> > > > > I
> >> > > > > >> see that you're having some trouble getting Point-Stat
to
> run:
> >> > > > > >>
> >> > > > > >> ERROR : parse_grid_mask() -> the mask_grid
requested
> >> "nepal"
> >> > is
> >> > > > not
> >> > > > > >> defined.
> >> > > > > >>
> >> > > > > >> This should be pretty easy to resolve. You've worked
with
> >> Julie
> >> > to
> >> > > > > >> generate a mask file for Nepal, and that mask is
contained in
> >> some
> >> > > > > NetCDF
> >> > > > > >> file (the output of the gen_vx_mask tool). For
example,
> let's
> >> say
> >> > > > it's
> >> > > > > in
> >> > > > > >> "/home/vsurapipith/nepal_mask.nc". You'll need to put
the
> >> path
> >> > to
> >> > > > that
> >> > > > > >> file in the "mask" section of the Point-Stat
configuration
> >> file:
> >> > > > > >>
> >> > > > > >> mask = {
> >> > > > > >> grid = [];
> >> > > > > >> poly = [ "/home/vsurapipith/nepal_mask.nc" ]; //
or
> >> > whatever
> >> > > > the
> >> > > > > >> actual path is to that data file
> >> > > > > >> sid = [];
> >> > > > > >> }
> >> > > > > >>
> >> > > > > >> And then try running Point-Stat again. For a more
complete
> >> > > > description
> >> > > > > of
> >> > > > > >> the masking options available, please see pages 54-55
of the
> >> MET
> >> > > > User's
> >> > > > > >> Guide:
> >> > > > > >>
> >> > > > > >>
> >> > > > > >>
> >> > > > >
> >> > > >
> >> > >
> >> >
> >>
>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
> >> > > > > >>
> >> > > > > >> Please let me know how it goes.
> >> > > > > >>
> >> > > > > >> Thanks,
> >> > > > > >> John
> >> > > > > >>
> >> > > > > >>
> >> > > > > >>
> >> > > > > >> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com via
RT <
> >> > > > > >> met_help at ucar.edu
> >> > > > > >> > wrote:
> >> > > > > >>
> >> > > > > >> >
> >> > > > > >> > <URL:
> >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> >> > > > > >> >
> >> > > > > >> > Hi Julie,
> >> > > > > >> >
> >> > > > > >> > I'm back....after a couple of weeks spending on
stumbling
> >> over
> >> > the
> >> > > > > >> > installation of UPP until I got it through, now I
have the
> >> GRIB1
> >> > > > files
> >> > > > > >> > systematically converted from the WRF output
generated
> here.
> >> > > > > >> >
> >> > > > > >> > I try the GRIB1 files with the PointStatConfig_Nepal
that I
> >> sent
> >> > > > > >> earlier to
> >> > > > > >> > do Point_Stat with our these WRF output, but
currently the
> >> error
> >> > > > shows
> >> > > > > >> up
> >> > > > > >> > as:
> >> > > > > >> >
> >> > > > > >> > ~/met-5.1_bugfix$ ./bin/point_stat
> >> > > > > >> > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00
> >> > > > > >> > tutorial/out/ascii2nc/nepal.nc
> >> > > > tutorial/config/PointStatConfig_Nepal
> >> > > > > >> > -outdir tutorial/out/point_stat -v 2
> >> > > > > >> > DEBUG 1: Default Config File:
> >> > > > > >> >
> >> > > > > >>
> >> > > > >
> >> > >
> >>
> /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> >> > > > > >> > DEBUG 1: User Config File:
> >> tutorial/config/PointStatConfig_Nepal
> >> > > > > >> > WARNING:
> >> > > > > >> > WARNING: parse_conf_interp() -> For neareast
neighbor
> >> > > interpolation
> >> > > > > >> method,
> >> > > > > >> > resetting method from "UW_MEAN" to "NEAREST" since
the
> >> > > interpolation
> >> > > > > >> width
> >> > > > > >> > is 1.
> >> > > > > >> > WARNING:
> >> > > > > >> > GSL_RNG_TYPE=mt19937
> >> > > > > >> > GSL_RNG_SEED=18446744072070604630
> >> > > > > >> > DEBUG 1: Forecast File:
> >> > > > > >> > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00
> >> > > > > >> > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal.nc
> >> > > > > >> > DEBUG 2:
> >> > > > > >> > DEBUG 2:
> >> > > > > >> >
> >> > > > > >> >
> >> > > > > >>
> >> > > > >
> >> > > >
> >> > >
> >> >
> >>
>
--------------------------------------------------------------------------------
> >> > > > > >> > DEBUG 2:
> >> > > > > >> > DEBUG 2: Reading data for SPFH/P500.
> >> > > > > >> > ERROR :
> >> > > > > >> > ERROR : parse_grid_mask() -> the mask_grid
requested
> >> "nepal" is
> >> > > not
> >> > > > > >> > defined.
> >> > > > > >> > ERROR :
> >> > > > > >> >
> >> > > > > >> > It looks to me like the mask file - nepal.nc, that I
> >> > deliberately
> >> > > > > took
> >> > > > > >> > from
> >> > > > > >> > what you generated earlier - was not seen/defined
well by
> MET
> >> > > here.
> >> > > > I
> >> > > > > >> have
> >> > > > > >> > actually just put the WRF output files in GRIB1
format
> >> (checked
> >> > > with
> >> > > > > as
> >> > > > > >> you
> >> > > > > >> > recommend) on the rap.ucar.edu ftp site
> >> > (surapipith_data)...Would
> >> > > > you
> >> > > > > >> > kindly let me know what should be done to correct
this,
> >> please?
> >> > > > > >> >
> >> > > > > >> > Best wishes,
> >> > > > > >> > Vanisa
> >> > > > > >> >
> >> > > > > >> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik
via RT <
> >> > > > > >> > met_help at ucar.edu>
> >> > > > > >> > wrote:
> >> > > > > >> >
> >> > > > > >> > > Hi Vanisa.
> >> > > > > >> > >
> >> > > > > >> > > While I don't have any experience installing UPP
myself,
> I
> >> did
> >> > > > take
> >> > > > > a
> >> > > > > >> > look
> >> > > > > >> > > at your log file where I saw:
> >> > > > > >> > >
> >> > > > > >> > > enc_jpeg2000.c:4:27: fatal error: jasper/jasper.h:
No
> such
> >> > file
> >> > > or
> >> > > > > >> > > directory compilation terminated.
> >> > > > > >> > >
> >> > > > > >> > > So, it looks like UPP is having trouble finding
jasper.h.
> >> > > > > >> > >
> >> > > > > >> > > Support for UPP is provided via wrfhelp at
ucar.edu. If
> >> > you're
> >> > > > > >> unable
> >> > > > > >> > to
> >> > > > > >> > > resolve the jasper problem, I would suggest
writing to
> them
> >> > > > > regarding
> >> > > > > >> > this
> >> > > > > >> > > issue.
> >> > > > > >> > >
> >> > > > > >> > > Thanks,
> >> > > > > >> > > Julie
> >> > > > > >> > >
> >> > > > > >> > >
> >> > > > > >> > > On Tue, Jan 5, 2016 at 4:13 AM, vanisas at gmail.com
via
> RT <
> >> > > > > >> > > met_help at ucar.edu>
> >> > > > > >> > > wrote:
> >> > > > > >> > >
> >> > > > > >> > > >
> >> > > > > >> > > > <URL:
> >> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> >> > > >
> >> > > > > >> > > >
> >> > > > > >> > > > Hi Julie,
> >> > > > > >> > > >
> >> > > > > >> > > > Happy New Year...:) Thank you very much for your
reply
> >> upon
> >> > > back
> >> > > > > to
> >> > > > > >> > > office.
> >> > > > > >> > > > Indeed, I did not use UPP for processing the WRF
(ARW)
> >> > output.
> >> > > > > >> > Therefore,
> >> > > > > >> > > > today I look into how to install UPP on my
workstation
> >> for
> >> > > this
> >> > > > > >> > purpose.
> >> > > > > >> > > > However, there is a problem that unipost.exe was
not
> >> > generated
> >> > > > > after
> >> > > > > >> > > > compilation, while other files came through. I
was
> >> wondered
> >> > if
> >> > > > it
> >> > > > > >> has
> >> > > > > >> > to
> >> > > > > >> > > do
> >> > > > > >> > > > with liking to JASPER. I have not been able to
spot
> why,
> >> and
> >> > > now
> >> > > > > >> attach
> >> > > > > >> > > the
> >> > > > > >> > > > compile_upp.log to you here, in case you might
be able
> to
> >> > see
> >> > > > what
> >> > > > > >> is
> >> > > > > >> > > wrong
> >> > > > > >> > > > in there.
> >> > > > > >> > > >
> >> > > > > >> > > > Thanks again; update more soon.
> >> > > > > >> > > >
> >> > > > > >> > > > Best wishes,
> >> > > > > >> > > > Vanisa
> >> > > > > >> > > >
> >> > > > > >> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie Prestopnik
via
> RT <
> >> > > > > >> > > met_help at ucar.edu
> >> > > > > >> > > > >
> >> > > > > >> > > > wrote:
> >> > > > > >> > > >
> >> > > > > >> > > > > Hi Vanisa. The file you uploaded may have
aged off
> of
> >> our
> >> > > ftp
> >> > > > > >> > site. I
> >> > > > > >> > > > do
> >> > > > > >> > > > > not see it there currently.
> >> > > > > >> > > > >
> >> > > > > >> > > > > I neglected to ask you if your WRF output file
was
> >> > > > > post-processed?
> >> > > > > >> > MET
> >> > > > > >> > > > is
> >> > > > > >> > > > > unable to read the output of WRF directly.
Instead,
> >> you
> >> > > must
> >> > > > > >> first
> >> > > > > >> > > > > post-process it using either the Unified
> PostProcessor
> >> > (UPP)
> >> > > > or
> >> > > > > >> the
> >> > > > > >> > > > pinterp
> >> > > > > >> > > > > tool (for ARW only). We strongly recommend
using UPP
> >> > since
> >> > > > > MET's
> >> > > > > >> > > support
> >> > > > > >> > > > > is most robust for the GRIB1 format that UPP
creates.
> >> > > > > >> > > > >
> >> > > > > >> > > > > If you have already post-processed your file,
please
> >> > trying
> >> > > > > >> uploading
> >> > > > > >> > > to
> >> > > > > >> > > > > our ftp site again. If you have not, please
> >> post-process
> >> > > the
> >> > > > > file
> >> > > > > >> > and
> >> > > > > >> > > > give
> >> > > > > >> > > > > it another try. Thanks!
> >> > > > > >> > > > >
> >> > > > > >> > > > > Happy New Year!
> >> > > > > >> > > > >
> >> > > > > >> > > > > Julie
> >> > > > > >> > > > >
> >> > > >
> >> > >
> >> > >
> >> >
> >> >
> >> > --
> >> > Julie Prestopnik
> >> > National Center for Atmospheric Research
> >> > Research Applications Laboratory
> >> > Phone: 303.497.8399
> >> > Email: jpresto at ucar.edu
> >> >
> >> >
> >>
> >>
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
>
>
>
> --
> Julie Prestopnik
> National Center for Atmospheric Research
> Research Applications Laboratory
> Phone: 303.497.8399
> Email: jpresto at ucar.edu
>
>
------------------------------------------------
Subject: point_stat for WRF output question
From: Julie Prestopnik
Time: Wed Feb 10 09:36:11 2016
Hi Vanisa. I'm glad that we were able to help. I'll go ahead and
resolve
this ticket, but please start a new one with any other questions you
may
have.
Thanks,
Julie
On Wed, Feb 10, 2016 at 4:01 AM, vanisas at gmail.com via RT
<met_help at ucar.edu
> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
>
> Hi Julie and John,
>
> thank you so much for your help identifying the problem and giving
detailed
> explanation. I got this solved now with adjusting the obs_window. I
am
> working into understanding the analysis using MET, and may send next
> questions on another topic soon.
>
> Best wishes,
> Vanisa
>
> On Tue, Feb 9, 2016 at 4:43 AM, Julie Prestopnik via RT
<met_help at ucar.edu
> >
> wrote:
>
> > Hi Vanisa. I was able to speak with John today. He helped me
better
> > understand why you are not getting the desired output.
> >
> > Looking in the nepal_ascii.nc, you can see in the "obs_arr" that
the
> > second
> > column (grib code) contains only values of 11. That grib code is
for 2-m
> > TMP. This explains why point_stat was showing us that for all
other
> fields
> > there were rejected observations due to the grib code, since the
> > observation file only contains observations for the 2-m TMP:
> >
> > DEBUG 3: Observations processed = 33816
> > DEBUG 3: Rejected: SID exclusion = 0
> > DEBUG 3: Rejected: GRIB code = 33816
> >
> > With regard to verifying the 2-m TMP, I saw that the observations
were
> > rejected due to the valid time:
> >
> > DEBUG 3: Observations processed = 33816
> > DEBUG 3: Rejected: SID exclusion = 0
> > DEBUG 3: Rejected: GRIB code = 0
> > DEBUG 3: Rejected: valid time = 33816
> >
> > I first ran "wgrib WRFPRS_d02.00 | grep TMP" to get the record
number for
> > the 2-m TMP:
> > ...
> >
> >
>
296:22645020:d=13010118:TMP:kpds5=11:kpds6=105:kpds7=2:TR=0:P1=0:P2=0:TimeU=1:2
> > m above gnd:anl:NAve=0
> > ...
> >
> > Then I ran "wgrib -d 296 -V WRFPRS_d02.00" to get the date/time
> information
> > for the 2-m TMP and got:
> >
> > rec 296:22645020:date 2013010118 TMP kpds5=11 kpds6=105 kpds7=2
> > levels=(0,2) grid=255 2 m above gnd anl:
> > TMP=Temp. [K]
> > timerange 0 P1 0 P2 0 TimeU 1 nx 270 ny 180 GDS grid 1
num_in_ave 0
> > missing 0
> > center 7 subcenter 0 process 125 Table 2 scan: WE:SN winds(grid)
> > Mercator: lat 23.937000 to 31.070000 by 5.000000 km nxny 48600
> > long 77.740000 to 89.838000 by 5.000000 km, (270 x 180)
scan 64
> > mode 136 Latin 0.000000 bdsgrid 1
> > min/max data 238.877 291.292 num bits 16 BDS_Ref 238877
DecScale 3
> > BinScale 0
> >
> >
> > You can see the date/time is 2013010118. Looking for that time in
the
> > observation file, nepal_ascii.nc produced no results. However,
there
> were
> > results for 2013010100, 2013010106, and 2013010112. So, as we can
see,
> > just as the debug statement showed, the observations were rejected
> because
> > the valid time didn't match up. Now you wouldn't necessarily want
to
> > verify the 18 hour with the 12 hour, 6 hour, or 0 hour, but to
show that
> > you will get MET output, you could extend the obs_window "beg"
and/or
> "end"
> > values in this case and you would get some matched pair output.
> >
> > Also, to answer your question "Do I need to generate a GRIB1 file
that
> > combine all the required hours...otherwise, can MET read in more
than one
> > GRIB1 files at a time for Point_Stat analysis?" Most users use a
grib
> file
> > for each time and run point stat once for each separate output
time.
> >
> > I hope this helps.
> >
> > Please let us know if you have additional questions.
> >
> > Thanks,
> > Julie
> >
> >
> > On Fri, Feb 5, 2016 at 10:31 AM, Julie Prestopnik
<jpresto at ucar.edu>
> > wrote:
> >
> > > Hi Vanisa. I apologize for the delayed response. I was out of
the
> > office
> > > yesterday.
> > >
> > > Unfortunately, I need some help from John to figure out what
exactly is
> > > going on here, but he is out of the office today. We will
revisit this
> > on
> > > Monday and get back to you as soon as we can.
> > >
> > > Thank you,
> > > Julie
> > >
> > > On Thu, Feb 4, 2016 at 5:47 AM, vanisas at gmail.com via RT <
> > > met_help at ucar.edu> wrote:
> > >
> > >>
> > >> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > >>
> > >> Hi Julie,
> > >>
> > >> I have just uploaded the 4 files I used for running point_stat
that
> > >> generated the output files that I put there earlier. Hope these
are
> > >> comprehensible. the command lines that I used are as follows:
> > >>
> > >> ~/met-5.1_bugfix$ ./bin/point_stat \
> > >>
> > >> > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > >>
> > >> > tutorial/out/ascii2nc/nepal_ascii.nc \
> > >>
> > >> > tutorial/config/PointStatConfig_Nepal \
> > >>
> > >> > -outdir tutorial/out/point_stat \
> > >>
> > >> > -v 2
> > >>
> > >>
> > >> Thanks alot for your kind help. Wish to see this through soon.
> > >>
> > >>
> > >> Best wishes,
> > >> Vanisa
> > >>
> > >> On Wed, Feb 3, 2016 at 11:41 PM, Julie Prestopnik via RT <
> > >> met_help at ucar.edu>
> > >> wrote:
> > >>
> > >> > Hi Vanisa.
> > >> >
> > >> > Could you please upload your input forecast file, input
observation
> > >> file,
> > >> > masking file, and your configuration file? Also, could you
please
> > send
> > >> me
> > >> > the command line argument for point_stat that you are
running?
> > >> >
> > >> > I can take a closer look and let you know what I find out.
> > >> >
> > >> > Thanks,
> > >> > Julie
> > >> >
> > >> > On Tue, Feb 2, 2016 at 2:32 AM, vanisas at gmail.com via RT <
> > >> > met_help at ucar.edu>
> > >> > wrote:
> > >> >
> > >> > >
> > >> > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > >> > >
> > >> > > Hi Julie,
> > >> > >
> > >> > > thank you for your reply last week; apologies for my delay
in
> > >> reporting.
> > >> > > What you suggest works. I can now generate the output files
using
> > >> > > Point_Stat. However, I do not yet understand how to
interpret or
> > make
> > >> use
> > >> > > of the output...I have been trying to read what described
in the
> > >> > > MET_User_Guide. Frankly speaking, I still do not get it.
Could you
> > >> kindly
> > >> > > see in my folder in the rap.ucar.edu ftp site in
> > >> incoming/irap/met_help
> > >> > > under surapipith_data. I attached here some of the files.
It looks
> > >> like
> > >> > > only the headers appear...or what I should do to generate
the
> > >> statistical
> > >> > > evaluation output. What surprise me was that the example
files I
> got
> > >> from
> > >> > > the US example (follow the online tutorial) look similar.
Could
> you
> > >> guide
> > >> > > me where to read to comprehend MET?
> > >> > >
> > >> > > One other question I have is that probably this is because
I use
> > only
> > >> one
> > >> > > hour of WRF output for this try. Do I need to generate a
GRIB1
> file
> > >> that
> > >> > > combine all the required hours...otherwise, can MET read in
more
> > than
> > >> one
> > >> > > GRIB1 files at a time for Point_Stat analysis?
> > >> > >
> > >> > > Best wishes,
> > >> > > Vanisa
> > >> > >
> > >> > > On Thu, Jan 28, 2016 at 12:00 AM, Julie Prestopnik via RT <
> > >> > > met_help at ucar.edu
> > >> > > > wrote:
> > >> > >
> > >> > > > Hi Vanisa. This is Julie again.
> > >> > > >
> > >> > > > I see that you are encountering a problem due to the
masking and
> > >> > > > verification grids not matching.
> > >> > > >
> > >> > > > One projection is a Lat/Lon grid and the other is a
Mercator.
> > >> > > >
> > >> > > > The mask file I created for you using gen_vx_mask used a
Lat/Lon
> > >> grid
> > >> > > file
> > >> > > > (gfs_2012040900_F003.grib). If you would like to get
your
> netcdf
> > >> mask
> > >> > > > file
> > >> > > > on a Mercator grid, you could run gen_vx_mask using your
file
> that
> > >> > > contains
> > >> > > > the grid that you would like to verify, which I believe
is on a
> > >> > Mercator
> > >> > > > projection (WRFPRS_d02.00?) and the nepal poly file which
will
> > then
> > >> > > produce
> > >> > > > a netcdf file on that same projection.
> > >> > > >
> > >> > > > So, for example, something like:
> > >> > > >
> > >> > > > gen_vx_mask WRFPRS_d02.00 nepal_detail.poly
nepal_detail.nc
> > >> > > >
> > >> > > > Then, you can rerun point_stat using the new netcdf mask
file
> on
> > >> the
> > >> > > > Mercator grid and that should fix the problem.
> > >> > > >
> > >> > > > However, if you're not going to be running thousands of
cases of
> > >> Point
> > >> > > > Stat, you don't necessarily need to run gen_vx_mask, but
rather
> > you
> > >> > could
> > >> > > > use the poly file directly with Point Stat in the
configuration
> > >> file in
> > >> > > the
> > >> > > > mask section that John described earlier. Using the poly
file
> > >> directly
> > >> > > can
> > >> > > > be slow, so if you're running point_stat on a large
number of
> > files,
> > >> > this
> > >> > > > is likely not the best way to use point_stat.
> > >> > > >
> > >> > > > I hope this helps. Please let me know if you have any
> questions.
> > >> > > >
> > >> > > > Thanks,
> > >> > > > Julie
> > >> > > >
> > >> > > > On Tue, Jan 26, 2016 at 11:04 PM, vanisas at gmail.com via
RT <
> > >> > > > met_help at ucar.edu> wrote:
> > >> > > >
> > >> > > > >
> > >> > > > > <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> >
> > >> > > > >
> > >> > > > > Hi John,
> > >> > > > >
> > >> > > > > apologies, I just notice that I did not remove the word
> "nepal"
> > in
> > >> > the
> > >> > > > > bracket ([]) behind grid = within mask section in the
> > PointConfig
> > >> > > file. I
> > >> > > > > have then removed that, and found a different error
(also
> rename
> > >> the
> > >> > > > ascii
> > >> > > > > file containing the ADPUPA observation from our
stations in
> > Nepal
> > >> to
> > >> > be
> > >> > > > > called "nepal_ascii.nc" to avoid confusion) :
> > >> > > > >
> > >> > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > >> > > > > > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > >> > > > > > tutorial/out/ascii2nc/nepal_ascii.nc \
> > >> > > > > > tutorial/config/PointStatConfig_Nepal \
> > >> > > > > > -outdir tutorial/out/point_stat \
> > >> > > > > > -v 2
> > >> > > > > DEBUG 1: Default Config File:
> > >> > > > >
> > >> > >
> > >>
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > >> > > > > DEBUG 1: User Config File:
> tutorial/config/PointStatConfig_Nepal
> > >> > > > > WARNING:
> > >> > > > > WARNING: parse_conf_interp() -> For neareast neighbor
> > >> interpolation
> > >> > > > method,
> > >> > > > > resetting method from "UW_MEAN" to "NEAREST" since the
> > >> interpolation
> > >> > > > width
> > >> > > > > is 1.
> > >> > > > > WARNING:
> > >> > > > > GSL_RNG_TYPE=mt19937
> > >> > > > > GSL_RNG_SEED=18446744072507402326
> > >> > > > > DEBUG 1: Forecast File:
> > >> > > > > /home/vsurapipith/met-5.1_bugfix/data/wrf/WRFPRS_d02.00
> > >> > > > > DEBUG 1: Observation File: tutorial/out/ascii2nc/
> nepal_ascii.nc
> > >> > > > > DEBUG 2:
> > >> > > > > DEBUG 2:
> > >> > > > >
> > >> > > > >
> > >> > > >
> > >> > >
> > >> >
> > >>
> >
>
--------------------------------------------------------------------------------
> > >> > > > > DEBUG 2:
> > >> > > > > DEBUG 2: Reading data for SPFH/P500.
> > >> > > > > ERROR :
> > >> > > > > ERROR : parse_poly_mask() -> The masking and
verification
> grids
> > >> do
> > >> > not
> > >> > > > > match:
> > >> > > > > ERROR : Projection: Lat/Lon Nx: 720 Ny: 361 lat_ll:
-90.000
> > >> lon_ll:
> > >> > > > -0.000
> > >> > > > > delta_lat: 0.500 delta_lon: 0.500 !=
> > >> > > > > ERROR : Projection: Mercator Nx: 270 Ny: 180
Lat_LL_radians:
> > >> 0.4178
> > >> > > > > Lon_LL_radians: -1.3568 Lat_UR_radians: 0.5423
Lon_UR_radians:
> > >> > -1.5680
> > >> > > > Mx:
> > >> > > > > 1273.9763 My: 1274.0471 Bx: -1728.5551 By: -548.4663
> > >> > > > > ERROR :
> > >> > > > >
> > >> > > > > It looks like I need to set for some correct projection
> > somewhere.
> > >> > > Could
> > >> > > > > you give a bit of guideline (in the while, I'll try to
read
> the
> > >> > > > > Met_Users_Guide_v5.1.pdf in more detail; may send
update again
> > >> before
> > >> > > end
> > >> > > > > of the day here)...
> > >> > > > >
> > >> > > > > Best wishes,
> > >> > > > > Vanisa
> > >> > > > >
> > >> > > > > On Wed, Jan 27, 2016 at 10:48 AM, Vanisa Surapipith <
> > >> > vanisas at gmail.com
> > >> > > >
> > >> > > > > wrote:
> > >> > > > >
> > >> > > > > > Hi John,
> > >> > > > > >
> > >> > > > > > thank you for your reply; I recall your help earlier,
and I
> am
> > >> only
> > >> > > > > > gradually getting to know MET. I am grateful to the
help as
> to
> > >> what
> > >> > > > Julie
> > >> > > > > > was responding patiently as well. I notice the
(presumably)
> > mask
> > >> > > files
> > >> > > > > that
> > >> > > > > > Julie helped generating for Nepal, and I try them.
However,
> > the
> > >> > > similar
> > >> > > > > > error still shows up, and I am trying to understand
why:
> > >> > > > > >
> > >> > > > > > ~/met-5.1_bugfix$ ./bin/point_stat \
> > >> > > > > > > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00 \
> > >> > > > > > > tutorial/out/ascii2nc/nepal.nc
> > >> > > > tutorial/config/PointStatConfig_Nepal \
> > >> > > > > > > -outdir tutorial/out/point_stat \
> > >> > > > > > > -v 2
> > >> > > > > > DEBUG 1: Default Config File:
> > >> > > > > >
> > >> > > >
> > >> >
> > >>
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > >> > > > > > DEBUG 1: User Config File:
> > tutorial/config/PointStatConfig_Nepal
> > >> > > > > > WARNING:
> > >> > > > > > WARNING: parse_conf_interp() -> For neareast neighbor
> > >> interpolation
> > >> > > > > > method, resetting method from "UW_MEAN" to "NEAREST"
since
> the
> > >> > > > > > interpolation width is 1.
> > >> > > > > > WARNING:
> > >> > > > > > GSL_RNG_TYPE=mt19937
> > >> > > > > > GSL_RNG_SEED=1682548822
> > >> > > > > > DEBUG 1: Forecast File:
> > >> > > > > > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00
> > >> > > > > > DEBUG 1: Observation File:
tutorial/out/ascii2nc/nepal.nc
> > >> > > > > > DEBUG 2:
> > >> > > > > > DEBUG 2:
> > >> > > > > >
> > >> > > > >
> > >> > > >
> > >> > >
> > >> >
> > >>
> >
>
--------------------------------------------------------------------------------
> > >> > > > > > DEBUG 2:
> > >> > > > > > DEBUG 2: Reading data for SPFH/P500.
> > >> > > > > > ERROR :
> > >> > > > > > ERROR : parse_grid_mask() -> the mask_grid requested
> "nepal"
> > is
> > >> > not
> > >> > > > > > defined.
> > >> > > > > > ERROR :
> > >> > > > > >
> > >> > > > > > Kindly see in the surapipith_data folder in
rap.ucar.edu
> site
> > >> the
> > >> > > > latest
> > >> > > > > > config file, as well as the suite of Nepal mask files
that
> > Julie
> > >> > made
> > >> > > > > and I
> > >> > > > > > use. Hope some progress would be before the end of
the day,
> or
> > >> once
> > >> > > > some
> > >> > > > > > hints are next sent from there...
> > >> > > > > >
> > >> > > > > > Best wishes,
> > >> > > > > > Vanisa
> > >> > > > > >
> > >> > > > > >
> > >> > > > > > On Tue, Jan 26, 2016 at 10:44 PM, John Halley Gotway
via RT
> <
> > >> > > > > > met_help at ucar.edu> wrote:
> > >> > > > > >
> > >> > > > > >> Hello Vanisa,
> > >> > > > > >>
> > >> > > > > >> My name is John, and I work with Julie who is out of
the
> > office
> > >> > > today.
> > >> > > > > I
> > >> > > > > >> see that you're having some trouble getting Point-
Stat to
> > run:
> > >> > > > > >>
> > >> > > > > >> ERROR : parse_grid_mask() -> the mask_grid
requested
> > >> "nepal"
> > >> > is
> > >> > > > not
> > >> > > > > >> defined.
> > >> > > > > >>
> > >> > > > > >> This should be pretty easy to resolve. You've
worked with
> > >> Julie
> > >> > to
> > >> > > > > >> generate a mask file for Nepal, and that mask is
contained
> in
> > >> some
> > >> > > > > NetCDF
> > >> > > > > >> file (the output of the gen_vx_mask tool). For
example,
> > let's
> > >> say
> > >> > > > it's
> > >> > > > > in
> > >> > > > > >> "/home/vsurapipith/nepal_mask.nc". You'll need to
put the
> > >> path
> > >> > to
> > >> > > > that
> > >> > > > > >> file in the "mask" section of the Point-Stat
configuration
> > >> file:
> > >> > > > > >>
> > >> > > > > >> mask = {
> > >> > > > > >> grid = [];
> > >> > > > > >> poly = [ "/home/vsurapipith/nepal_mask.nc" ];
// or
> > >> > whatever
> > >> > > > the
> > >> > > > > >> actual path is to that data file
> > >> > > > > >> sid = [];
> > >> > > > > >> }
> > >> > > > > >>
> > >> > > > > >> And then try running Point-Stat again. For a more
complete
> > >> > > > description
> > >> > > > > of
> > >> > > > > >> the masking options available, please see pages 54-
55 of
> the
> > >> MET
> > >> > > > User's
> > >> > > > > >> Guide:
> > >> > > > > >>
> > >> > > > > >>
> > >> > > > > >>
> > >> > > > >
> > >> > > >
> > >> > >
> > >> >
> > >>
> >
>
http://www.dtcenter.org/met/users/docs/users_guide/MET_Users_Guide_v5.1.pdf
> > >> > > > > >>
> > >> > > > > >> Please let me know how it goes.
> > >> > > > > >>
> > >> > > > > >> Thanks,
> > >> > > > > >> John
> > >> > > > > >>
> > >> > > > > >>
> > >> > > > > >>
> > >> > > > > >> On Tue, Jan 26, 2016 at 4:01 AM, vanisas at gmail.com
via RT
> <
> > >> > > > > >> met_help at ucar.edu
> > >> > > > > >> > wrote:
> > >> > > > > >>
> > >> > > > > >> >
> > >> > > > > >> > <URL:
> > >> https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467 >
> > >> > > > > >> >
> > >> > > > > >> > Hi Julie,
> > >> > > > > >> >
> > >> > > > > >> > I'm back....after a couple of weeks spending on
stumbling
> > >> over
> > >> > the
> > >> > > > > >> > installation of UPP until I got it through, now I
have
> the
> > >> GRIB1
> > >> > > > files
> > >> > > > > >> > systematically converted from the WRF output
generated
> > here.
> > >> > > > > >> >
> > >> > > > > >> > I try the GRIB1 files with the
PointStatConfig_Nepal
> that I
> > >> sent
> > >> > > > > >> earlier to
> > >> > > > > >> > do Point_Stat with our these WRF output, but
currently
> the
> > >> error
> > >> > > > shows
> > >> > > > > >> up
> > >> > > > > >> > as:
> > >> > > > > >> >
> > >> > > > > >> > ~/met-5.1_bugfix$ ./bin/point_stat
> > >> > > > > >> > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00
> > >> > > > > >> > tutorial/out/ascii2nc/nepal.nc
> > >> > > > tutorial/config/PointStatConfig_Nepal
> > >> > > > > >> > -outdir tutorial/out/point_stat -v 2
> > >> > > > > >> > DEBUG 1: Default Config File:
> > >> > > > > >> >
> > >> > > > > >>
> > >> > > > >
> > >> > >
> > >>
> > /home/vsurapipith/met-
5.1_bugfix/share/met/config/PointStatConfig_default
> > >> > > > > >> > DEBUG 1: User Config File:
> > >> tutorial/config/PointStatConfig_Nepal
> > >> > > > > >> > WARNING:
> > >> > > > > >> > WARNING: parse_conf_interp() -> For neareast
neighbor
> > >> > > interpolation
> > >> > > > > >> method,
> > >> > > > > >> > resetting method from "UW_MEAN" to "NEAREST" since
the
> > >> > > interpolation
> > >> > > > > >> width
> > >> > > > > >> > is 1.
> > >> > > > > >> > WARNING:
> > >> > > > > >> > GSL_RNG_TYPE=mt19937
> > >> > > > > >> > GSL_RNG_SEED=18446744072070604630
> > >> > > > > >> > DEBUG 1: Forecast File:
> > >> > > > > >> > /home/vsurapipith/met-
5.1_bugfix/data/wrf/WRFPRS_d02.00
> > >> > > > > >> > DEBUG 1: Observation File: tutorial/out/ascii2nc/
> nepal.nc
> > >> > > > > >> > DEBUG 2:
> > >> > > > > >> > DEBUG 2:
> > >> > > > > >> >
> > >> > > > > >> >
> > >> > > > > >>
> > >> > > > >
> > >> > > >
> > >> > >
> > >> >
> > >>
> >
>
--------------------------------------------------------------------------------
> > >> > > > > >> > DEBUG 2:
> > >> > > > > >> > DEBUG 2: Reading data for SPFH/P500.
> > >> > > > > >> > ERROR :
> > >> > > > > >> > ERROR : parse_grid_mask() -> the mask_grid
requested
> > >> "nepal" is
> > >> > > not
> > >> > > > > >> > defined.
> > >> > > > > >> > ERROR :
> > >> > > > > >> >
> > >> > > > > >> > It looks to me like the mask file - nepal.nc, that
I
> > >> > deliberately
> > >> > > > > took
> > >> > > > > >> > from
> > >> > > > > >> > what you generated earlier - was not seen/defined
well by
> > MET
> > >> > > here.
> > >> > > > I
> > >> > > > > >> have
> > >> > > > > >> > actually just put the WRF output files in GRIB1
format
> > >> (checked
> > >> > > with
> > >> > > > > as
> > >> > > > > >> you
> > >> > > > > >> > recommend) on the rap.ucar.edu ftp site
> > >> > (surapipith_data)...Would
> > >> > > > you
> > >> > > > > >> > kindly let me know what should be done to correct
this,
> > >> please?
> > >> > > > > >> >
> > >> > > > > >> > Best wishes,
> > >> > > > > >> > Vanisa
> > >> > > > > >> >
> > >> > > > > >> > On Tue, Jan 5, 2016 at 11:06 PM, Julie Prestopnik
via RT
> <
> > >> > > > > >> > met_help at ucar.edu>
> > >> > > > > >> > wrote:
> > >> > > > > >> >
> > >> > > > > >> > > Hi Vanisa.
> > >> > > > > >> > >
> > >> > > > > >> > > While I don't have any experience installing UPP
> myself,
> > I
> > >> did
> > >> > > > take
> > >> > > > > a
> > >> > > > > >> > look
> > >> > > > > >> > > at your log file where I saw:
> > >> > > > > >> > >
> > >> > > > > >> > > enc_jpeg2000.c:4:27: fatal error:
jasper/jasper.h: No
> > such
> > >> > file
> > >> > > or
> > >> > > > > >> > > directory compilation terminated.
> > >> > > > > >> > >
> > >> > > > > >> > > So, it looks like UPP is having trouble finding
> jasper.h.
> > >> > > > > >> > >
> > >> > > > > >> > > Support for UPP is provided via wrfhelp at
ucar.edu.
> If
> > >> > you're
> > >> > > > > >> unable
> > >> > > > > >> > to
> > >> > > > > >> > > resolve the jasper problem, I would suggest
writing to
> > them
> > >> > > > > regarding
> > >> > > > > >> > this
> > >> > > > > >> > > issue.
> > >> > > > > >> > >
> > >> > > > > >> > > Thanks,
> > >> > > > > >> > > Julie
> > >> > > > > >> > >
> > >> > > > > >> > >
> > >> > > > > >> > > On Tue, Jan 5, 2016 at 4:13 AM,
vanisas at gmail.com via
> > RT <
> > >> > > > > >> > > met_help at ucar.edu>
> > >> > > > > >> > > wrote:
> > >> > > > > >> > >
> > >> > > > > >> > > >
> > >> > > > > >> > > > <URL:
> > >> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=74467
> > >> > > >
> > >> > > > > >> > > >
> > >> > > > > >> > > > Hi Julie,
> > >> > > > > >> > > >
> > >> > > > > >> > > > Happy New Year...:) Thank you very much for
your
> reply
> > >> upon
> > >> > > back
> > >> > > > > to
> > >> > > > > >> > > office.
> > >> > > > > >> > > > Indeed, I did not use UPP for processing the
WRF
> (ARW)
> > >> > output.
> > >> > > > > >> > Therefore,
> > >> > > > > >> > > > today I look into how to install UPP on my
> workstation
> > >> for
> > >> > > this
> > >> > > > > >> > purpose.
> > >> > > > > >> > > > However, there is a problem that unipost.exe
was not
> > >> > generated
> > >> > > > > after
> > >> > > > > >> > > > compilation, while other files came through. I
was
> > >> wondered
> > >> > if
> > >> > > > it
> > >> > > > > >> has
> > >> > > > > >> > to
> > >> > > > > >> > > do
> > >> > > > > >> > > > with liking to JASPER. I have not been able to
spot
> > why,
> > >> and
> > >> > > now
> > >> > > > > >> attach
> > >> > > > > >> > > the
> > >> > > > > >> > > > compile_upp.log to you here, in case you might
be
> able
> > to
> > >> > see
> > >> > > > what
> > >> > > > > >> is
> > >> > > > > >> > > wrong
> > >> > > > > >> > > > in there.
> > >> > > > > >> > > >
> > >> > > > > >> > > > Thanks again; update more soon.
> > >> > > > > >> > > >
> > >> > > > > >> > > > Best wishes,
> > >> > > > > >> > > > Vanisa
> > >> > > > > >> > > >
> > >> > > > > >> > > > On Tue, Jan 5, 2016 at 4:27 AM, Julie
Prestopnik via
> > RT <
> > >> > > > > >> > > met_help at ucar.edu
> > >> > > > > >> > > > >
> > >> > > > > >> > > > wrote:
> > >> > > > > >> > > >
> > >> > > > > >> > > > > Hi Vanisa. The file you uploaded may have
aged off
> > of
> > >> our
> > >> > > ftp
> > >> > > > > >> > site. I
> > >> > > > > >> > > > do
> > >> > > > > >> > > > > not see it there currently.
> > >> > > > > >> > > > >
> > >> > > > > >> > > > > I neglected to ask you if your WRF output
file was
> > >> > > > > post-processed?
> > >> > > > > >> > MET
> > >> > > > > >> > > > is
> > >> > > > > >> > > > > unable to read the output of WRF directly.
> Instead,
> > >> you
> > >> > > must
> > >> > > > > >> first
> > >> > > > > >> > > > > post-process it using either the Unified
> > PostProcessor
> > >> > (UPP)
> > >> > > > or
> > >> > > > > >> the
> > >> > > > > >> > > > pinterp
> > >> > > > > >> > > > > tool (for ARW only). We strongly recommend
using
> UPP
> > >> > since
> > >> > > > > MET's
> > >> > > > > >> > > support
> > >> > > > > >> > > > > is most robust for the GRIB1 format that UPP
> creates.
> > >> > > > > >> > > > >
> > >> > > > > >> > > > > If you have already post-processed your
file,
> please
> > >> > trying
> > >> > > > > >> uploading
> > >> > > > > >> > > to
> > >> > > > > >> > > > > our ftp site again. If you have not, please
> > >> post-process
> > >> > > the
> > >> > > > > file
> > >> > > > > >> > and
> > >> > > > > >> > > > give
> > >> > > > > >> > > > > it another try. Thanks!
> > >> > > > > >> > > > >
> > >> > > > > >> > > > > Happy New Year!
> > >> > > > > >> > > > >
> > >> > > > > >> > > > > Julie
> > >> > > > > >> > > > >
> > >> > > >
> > >> > >
> > >> > >
> > >> >
> > >> >
> > >> > --
> > >> > Julie Prestopnik
> > >> > National Center for Atmospheric Research
> > >> > Research Applications Laboratory
> > >> > Phone: 303.497.8399
> > >> > Email: jpresto at ucar.edu
> > >> >
> > >> >
> > >>
> > >>
> > >
> > >
> > > --
> > > Julie Prestopnik
> > > National Center for Atmospheric Research
> > > Research Applications Laboratory
> > > Phone: 303.497.8399
> > > Email: jpresto at ucar.edu
> > >
> >
> >
> >
> > --
> > Julie Prestopnik
> > National Center for Atmospheric Research
> > Research Applications Laboratory
> > Phone: 303.497.8399
> > Email: jpresto at ucar.edu
> >
> >
>
>
--
Julie Prestopnik
National Center for Atmospheric Research
Research Applications Laboratory
Phone: 303.497.8399
Email: jpresto at ucar.edu
------------------------------------------------
More information about the Met_help
mailing list