[Met_help] [rt.rap.ucar.edu #76349] History for pointstat problem
John Halley Gotway via RT
met_help at ucar.edu
Thu May 19 13:53:25 MDT 2016
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
Dear All,
I'm facing some troubles with pointstat routine, in detail I couldn't obtain matching pairs using my data.
Using verbose debugging I see that Grid definition is not correct because for example
as lower left (ll) latitude is indicated a value that in effects is my upper level value , as you can see in the following screenshots that are taken from pointstat log the first
and a call to cdo libraries the latter.
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
DEBUG 2:
DEBUG 2: Reading data for TMP/Z2.
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match for VarInfo "TMP/Z2" in GRIB record 181 of GRIB file "/lustre/arpa/gallaii/scratch/wrf4oper//wr
f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB records matching VarInfo "TMP/Z2" in GRIB file "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
DEBUG 3: Use the matching forecast and observation grids.
DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90 lat_ll: 46.920 lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025
DEBUG 2: For TMP/Z2 found 1 forecast levels and 0 climatology levels.
DEBUG 2:
DEBUG 2: --------------------------------------------------------------------------------
.....
Whereas interrogating the same GRIB file with cdo I obtain
Grid coordinates :
1 : lonlat > size : dim = 8100 nx = 90 ny = 90
lon : first = 11.888 last = 14.149 inc = 0.0254044944 degrees_east
lat : first = 45.322 last = 46.92 inc = 0.0179550562 degrees_north
I'm missing something obvious? GRIB files come from UPP (remapping the grid from lambert conformal to regular lon lat)
It's more evident looking at the image produced with plot_data_plane : gray/blue colors should be the Alps in Italy
Whereas political boundaries show part of Denmark , Sweden ...
Any advice?
Thanks
Irene
+------------------------------------+
WE ALL are environmentally responsible. By not printing this email you've helped save paper, ink, and millions of trees.
+------------------------------------+
Irene Gallai
Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
SOC Direzione Tecnica Scientifica
CRMA - Centro Regionale di Modellistica Ambientale
Via Cairoli 14
I-33057 Palmanova (UD) ITALY
tel pers.: +39 0432 1918044
tel centr.: +39 0432 1918120
e-mail : irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
+------------------------------------+
AVVISO DI RISERVATEZZA
"Informazioni riservate possono essere contenute nel messaggio o nei suoi allegati. Se non siete i destinatari indicati nel messaggio, o responsabili per la sua consegna alla persona, o se avete ricevuto il messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed i suoi allegati. Grazie."
CONFIDENTIALITY NOTICE
"Confidential information may be contained in this message or in its attachments. If you are not the addressee indicated in this message, or responsible for message delivering to that person, or if you have received this message in error, you may not transcribe, copy or deliver this message to anyone. In that case, you should delete/destroy this message and its attachments. Thank you .
AVVISO DI RISERVATEZZA Informazioni riservate possono essere contenute nel messaggio o nei suoi allegati. Se non siete i destinatari indicati nel messaggio, o responsabili per la sua consegna alla persona, o se avete ricevuto il messaggio per errore, siete pregati di non trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo a cancellare il messaggio ed i suoi allegati. Grazie.
CONFIDENTIALITY NOTICE Confidential information may be contained in this message or in its attachments. If you are not the addressee indicated in this message, or responsible for message delivering to that person, or if you have received this message in error, you may not transcribe, copy or deliver this message to anyone. In that case, you should delete this message and its attachments. Thank you.
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: pointstat problem
From: John Halley Gotway
Time: Thu May 12 14:40:19 2016
Hello Irene,
I see you're having trouble with the grid definition in a GRIB file.
Could
you please post a sample GRIB file to our anonymous ftp site following
these instructions?
http://www.dtcenter.org/met/users/support/met_help.php#ftp
I'd like to look at the grid definition in the file, plot is using
plot_data_plane (as you've done), and also plot it using Unidata's IDV
tool. Also, do you have a plot showing the "intended" domain?
Hopefully that info will help us track down the discrepancy you're
seeing.
Lastly, please tell me the version of MET you're using.
Thanks,
John
On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> Transaction: Ticket created by irene.gallai at arpa.fvg.it
> Queue: met_help
> Subject: pointstat problem
> Owner: Nobody
> Requestors: irene.gallai at arpa.fvg.it
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
>
>
>
> Dear All,
>
> I'm facing some troubles with pointstat routine, in detail I
couldn't
> obtain matching pairs using my data.
> Using verbose debugging I see that Grid definition is not correct
because
> for example
> as lower left (ll) latitude is indicated a value that in effects is
my
> upper level value , as you can see in the following screenshots that
are
> taken from pointstat log the first
> and a call to cdo libraries the latter.
>
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> DEBUG 2:
> DEBUG 2: Reading data for TMP/Z2.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
> VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
> matching VarInfo "TMP/Z2" in GRIB file
> "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> DEBUG 3: Use the matching forecast and observation grids.
> DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90 lat_ll:
46.920
> lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025
> DEBUG 2: For TMP/Z2 found 1 forecast levels and 0 climatology
levels.
> DEBUG 2:
> DEBUG 2:
>
--------------------------------------------------------------------------------
> .....
>
> Whereas interrogating the same GRIB file with cdo I obtain
> Grid coordinates :
> 1 : lonlat > size : dim = 8100 nx = 90 ny = 90
> lon : first = 11.888 last = 14.149
inc =
> 0.0254044944 degrees_east
> lat : first = 45.322 last = 46.92
inc =
> 0.0179550562 degrees_north
>
>
> I'm missing something obvious? GRIB files come from UPP (remapping
the
> grid from lambert conformal to regular lon lat)
>
> It's more evident looking at the image produced with plot_data_plane
:
> gray/blue colors should be the Alps in Italy
> Whereas political boundaries show part of Denmark , Sweden ...
> Any advice?
>
> Thanks
> Irene
>
>
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail :
irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
suoi
> allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili
> per la sua consegna alla persona, o se avete ricevuto il messaggio
per
> errore, non dovete trascriverlo, copiarlo o inviarlo a nessuno. In
tal
> caso, dovete cancellare/distruggere il messaggio ed i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
or
> responsible for message delivering to that person, or if you have
received
> this message in error, you may not transcribe, copy or deliver this
message
> to anyone. In that case, you should delete/destroy this message and
its
> attachments. Thank you .
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute nel
> messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel
> messaggio, o responsabili per la sua consegna alla persona, o se
avete
> ricevuto il messaggio per errore, siete pregati di non trascriverlo,
> copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo a cancellare
il
> messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
this
> message or in its attachments. If you are not the addressee
indicated in
> this message, or responsible for message delivering to that person,
or if
> you have received this message in error, you may not transcribe,
copy or
> deliver this message to anyone. In that case, you should delete this
> message and its attachments. Thank you.
>
>
------------------------------------------------
Subject: R: [rt.rap.ucar.edu #76349] pointstat problem
From: Gallai Irene
Time: Fri May 13 03:52:05 2016
Hello,
As requested I posted 3 GRIB files (one for each of our domains) and 2
png images.
The one named config_domain show the 3 domains.
d02.png instead is been obtained plotting with plot_data_plane the
GRIB file for the second domain obtained with UPP starting from the
same wrf output files but leaving the coordinate in lambert conformal
and not transforming it in regular lat lon: in this case domain
localization is correct
Concerning the version of MET I'm using is the 5.1 with bug-fix
available at the date of 7 january 2016.
Thanks
Irene
+------------------------------------+
WE ALL are environmentally responsible. By not printing this email
you've helped save paper, ink, and millions of trees.
+------------------------------------+
Irene Gallai
Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
SOC Direzione Tecnica Scientifica
CRMA - Centro Regionale di Modellistica Ambientale
Via Cairoli 14
I-33057 Palmanova (UD) ITALY
tel pers.: +39 0432 1918044
tel centr.: +39 0432 1918120
e-mail : irene.gallai at arpa.fvg.it
+------------------------------------+
AVVISO DI RISERVATEZZA
"Informazioni riservate possono essere contenute nel messaggio o nei
suoi allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili per la sua consegna alla persona, o se avete ricevuto il
messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a
nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed i
suoi allegati. Grazie."
CONFIDENTIALITY NOTICE
"Confidential information may be contained in this message or in its
attachments. If you are not the addressee indicated in this message,
or responsible for message delivering to that person, or if you have
received this message in error, you may not transcribe, copy or
deliver this message to anyone. In that case, you should
delete/destroy this message and its attachments. Thank you .
-----Messaggio originale-----
Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Inviato: giovedì 12 maggio 2016 22:40
A: Gallai Irene
Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
Hello Irene,
I see you're having trouble with the grid definition in a GRIB file.
Could you please post a sample GRIB file to our anonymous ftp site
following these instructions?
http://www.dtcenter.org/met/users/support/met_help.php#ftp
I'd like to look at the grid definition in the file, plot is using
plot_data_plane (as you've done), and also plot it using Unidata's IDV
tool. Also, do you have a plot showing the "intended" domain?
Hopefully that info will help us track down the discrepancy you're
seeing.
Lastly, please tell me the version of MET you're using.
Thanks,
John
On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> Transaction: Ticket created by irene.gallai at arpa.fvg.it
> Queue: met_help
> Subject: pointstat problem
> Owner: Nobody
> Requestors: irene.gallai at arpa.fvg.it
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
> >
>
>
>
> Dear All,
>
> I'm facing some troubles with pointstat routine, in detail I
couldn't
> obtain matching pairs using my data.
> Using verbose debugging I see that Grid definition is not correct
> because for example as lower left (ll) latitude is indicated a
value
> that in effects is my upper level value , as you can see in the
> following screenshots that are taken from pointstat log the first
and
> a call to cdo libraries the latter.
>
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> DEBUG 2:
> DEBUG 2: Reading data for TMP/Z2.
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
> VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
> matching VarInfo "TMP/Z2" in GRIB file
> "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> DEBUG 3: Use the matching forecast and observation grids.
> DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90 lat_ll:
> 46.920
> lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025 DEBUG 2: For
TMP/Z2
> found 1 forecast levels and 0 climatology levels.
> DEBUG 2:
> DEBUG 2:
>
----------------------------------------------------------------------
> ----------
> .....
>
> Whereas interrogating the same GRIB file with cdo I obtain
> Grid coordinates :
> 1 : lonlat > size : dim = 8100 nx = 90 ny = 90
> lon : first = 11.888 last = 14.149
inc =
> 0.0254044944 degrees_east
> lat : first = 45.322 last = 46.92
inc =
> 0.0179550562 degrees_north
>
>
> I'm missing something obvious? GRIB files come from UPP (remapping
the
> grid from lambert conformal to regular lon lat)
>
> It's more evident looking at the image produced with plot_data_plane
:
> gray/blue colors should be the Alps in Italy Whereas political
> boundaries show part of Denmark , Sweden ...
> Any advice?
>
> Thanks
> Irene
>
>
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail :
irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
> suoi allegati. Se non siete i destinatari indicati nel messaggio, o
> responsabili per la sua consegna alla persona, o se avete ricevuto
il
> messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a
> nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed
i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
> or responsible for message delivering to that person, or if you have
> received this message in error, you may not transcribe, copy or
> deliver this message to anyone. In that case, you should
> delete/destroy this message and its attachments. Thank you .
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute
> nel messaggio o nei suoi allegati. Se non siete i destinatari
indicati
> nel messaggio, o responsabili per la sua consegna alla persona, o se
> avete ricevuto il messaggio per errore, siete pregati di non
> trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi
invitiamo
> a cancellare il messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
> this message or in its attachments. If you are not the addressee
> indicated in this message, or responsible for message delivering to
> that person, or if you have received this message in error, you may
> not transcribe, copy or deliver this message to anyone. In that
case,
> you should delete this message and its attachments. Thank you.
>
>
AVVISO DI RISERVATEZZA Informazioni riservate possono essere contenute
nel messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel messaggio, o responsabili per la sua consegna alla persona, o se
avete ricevuto il messaggio per errore, siete pregati di non
trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo
a cancellare il messaggio ed i suoi allegati. Grazie.
CONFIDENTIALITY NOTICE Confidential information may be contained in
this message or in its attachments. If you are not the addressee
indicated in this message, or responsible for message delivering to
that person, or if you have received this message in error, you may
not transcribe, copy or deliver this message to anyone. In that case,
you should delete this message and its attachments. Thank you.
------------------------------------------------
Subject: pointstat problem
From: John Halley Gotway
Time: Fri May 13 17:18:23 2016
Irene,
Thank you for sending me those sample data files. I took a look and
am
able to reproduce the problem you're seeing. I ran your data file
both
through MET and IDV and see that IDV plots it correctly while MET does
not.
I strongly suspect that MET has a bug in the processing of the GRIB
grid
description section. I'll work on narrowing down the problem, fixing
it,
and posting a bugfix next week.
I will email when a bugfix for this issue has been posted to the MET
website.
Thanks for letting us know about this problem!
John
On Fri, May 13, 2016 at 3:52 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
>
> Hello,
> As requested I posted 3 GRIB files (one for each of our domains) and
2 png
> images.
> The one named config_domain show the 3 domains.
> d02.png instead is been obtained plotting with plot_data_plane the
GRIB
> file for the second domain obtained with UPP starting from the same
wrf
> output files but leaving the coordinate in lambert conformal and not
> transforming it in regular lat lon: in this case domain localization
is
> correct
>
> Concerning the version of MET I'm using is the 5.1 with bug-fix
available
> at the date of 7 january 2016.
>
> Thanks
> Irene
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail : irene.gallai at arpa.fvg.it
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
suoi
> allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili
> per la sua consegna alla persona, o se avete ricevuto il messaggio
per
> errore, non dovete trascriverlo, copiarlo o inviarlo a nessuno. In
tal
> caso, dovete cancellare/distruggere il messaggio ed i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
or
> responsible for message delivering to that person, or if you have
received
> this message in error, you may not transcribe, copy or deliver this
message
> to anyone. In that case, you should delete/destroy this message and
its
> attachments. Thank you .
>
>
> -----Messaggio originale-----
> Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Inviato: giovedì 12 maggio 2016 22:40
> A: Gallai Irene
> Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
>
> Hello Irene,
>
> I see you're having trouble with the grid definition in a GRIB file.
> Could you please post a sample GRIB file to our anonymous ftp site
> following these instructions?
>
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> I'd like to look at the grid definition in the file, plot is using
> plot_data_plane (as you've done), and also plot it using Unidata's
IDV
> tool. Also, do you have a plot showing the "intended" domain?
>
> Hopefully that info will help us track down the discrepancy you're
seeing.
>
> Lastly, please tell me the version of MET you're using.
>
> Thanks,
> John
>
>
> On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> > Transaction: Ticket created by irene.gallai at arpa.fvg.it
> > Queue: met_help
> > Subject: pointstat problem
> > Owner: Nobody
> > Requestors: irene.gallai at arpa.fvg.it
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
> > >
> >
> >
> >
> > Dear All,
> >
> > I'm facing some troubles with pointstat routine, in detail I
couldn't
> > obtain matching pairs using my data.
> > Using verbose debugging I see that Grid definition is not correct
> > because for example as lower left (ll) latitude is indicated a
value
> > that in effects is my upper level value , as you can see in the
> > following screenshots that are taken from pointstat log the first
and
> > a call to cdo libraries the latter.
> >
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > DEBUG 2:
> > DEBUG 2: Reading data for TMP/Z2.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
for
> > VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> > "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> > f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
records
> > matching VarInfo "TMP/Z2" in GRIB file
> > "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> > 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > DEBUG 3: Use the matching forecast and observation grids.
> > DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90
lat_ll:
> > 46.920
> > lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025 DEBUG 2: For
TMP/Z2
> > found 1 forecast levels and 0 climatology levels.
> > DEBUG 2:
> > DEBUG 2:
> >
----------------------------------------------------------------------
> > ----------
> > .....
> >
> > Whereas interrogating the same GRIB file with cdo I obtain
> > Grid coordinates :
> > 1 : lonlat > size : dim = 8100 nx = 90 ny = 90
> > lon : first = 11.888 last = 14.149
inc =
> > 0.0254044944 degrees_east
> > lat : first = 45.322 last = 46.92
inc =
> > 0.0179550562 degrees_north
> >
> >
> > I'm missing something obvious? GRIB files come from UPP (remapping
the
> > grid from lambert conformal to regular lon lat)
> >
> > It's more evident looking at the image produced with
plot_data_plane :
> > gray/blue colors should be the Alps in Italy Whereas political
> > boundaries show part of Denmark , Sweden ...
> > Any advice?
> >
> > Thanks
> > Irene
> >
> >
> > +------------------------------------+
> > WE ALL are environmentally responsible. By not printing this
email
> > you've helped save paper, ink, and millions of trees.
> > +------------------------------------+
> >
> > Irene Gallai
> >
> > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > SOC Direzione Tecnica Scientifica
> > CRMA - Centro Regionale di Modellistica Ambientale
> >
> > Via Cairoli 14
> > I-33057 Palmanova (UD) ITALY
> > tel pers.: +39 0432 1918044
> > tel centr.: +39 0432 1918120
> > e-mail :
irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> > +------------------------------------+
> >
> >
> > AVVISO DI RISERVATEZZA
> > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > suoi allegati. Se non siete i destinatari indicati nel messaggio,
o
> > responsabili per la sua consegna alla persona, o se avete ricevuto
il
> > messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo
a
> > nessuno. In tal caso, dovete cancellare/distruggere il messaggio
ed i
> suoi allegati.
> > Grazie."
> >
> > CONFIDENTIALITY NOTICE
> > "Confidential information may be contained in this message or in
its
> > attachments. If you are not the addressee indicated in this
message,
> > or responsible for message delivering to that person, or if you
have
> > received this message in error, you may not transcribe, copy or
> > deliver this message to anyone. In that case, you should
> > delete/destroy this message and its attachments. Thank you .
> >
> >
> > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute
> > nel messaggio o nei suoi allegati. Se non siete i destinatari
indicati
> > nel messaggio, o responsabili per la sua consegna alla persona, o
se
> > avete ricevuto il messaggio per errore, siete pregati di non
> > trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi
invitiamo
> > a cancellare il messaggio ed i suoi allegati. Grazie.
> > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > this message or in its attachments. If you are not the addressee
> > indicated in this message, or responsible for message delivering
to
> > that person, or if you have received this message in error, you
may
> > not transcribe, copy or deliver this message to anyone. In that
case,
> > you should delete this message and its attachments. Thank you.
> >
> >
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute nel
> messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel
> messaggio, o responsabili per la sua consegna alla persona, o se
avete
> ricevuto il messaggio per errore, siete pregati di non trascriverlo,
> copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo a cancellare
il
> messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
this
> message or in its attachments. If you are not the addressee
indicated in
> this message, or responsible for message delivering to that person,
or if
> you have received this message in error, you may not transcribe,
copy or
> deliver this message to anyone. In that case, you should delete this
> message and its attachments. Thank you.
>
>
>
------------------------------------------------
Subject: R: [rt.rap.ucar.edu #76349] pointstat problem
From: Gallai Irene
Time: Mon May 16 01:40:35 2016
John,
many thanks to you,
I'm happy to contribute to MET development, even if only discovering
bugs...
I will be waiting for the fix.
Thank you again
Irene
+------------------------------------+
WE ALL are environmentally responsible. By not printing this email
you've helped save paper, ink, and millions of trees.
+------------------------------------+
Irene Gallai
Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
SOC Direzione Tecnica Scientifica
CRMA - Centro Regionale di Modellistica Ambientale
Via Cairoli 14
I-33057 Palmanova (UD) ITALY
tel pers.: +39 0432 1918044
tel centr.: +39 0432 1918120
e-mail : irene.gallai at arpa.fvg.it
+------------------------------------+
AVVISO DI RISERVATEZZA
"Informazioni riservate possono essere contenute nel messaggio o nei
suoi allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili per la sua consegna alla persona, o se avete ricevuto il
messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a
nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed i
suoi allegati. Grazie."
CONFIDENTIALITY NOTICE
"Confidential information may be contained in this message or in its
attachments. If you are not the addressee indicated in this message,
or responsible for message delivering to that person, or if you have
received this message in error, you may not transcribe, copy or
deliver this message to anyone. In that case, you should
delete/destroy this message and its attachments. Thank you .
-----Messaggio originale-----
Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Inviato: sabato 14 maggio 2016 01:18
A: Gallai Irene
Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
Irene,
Thank you for sending me those sample data files. I took a look and
am able to reproduce the problem you're seeing. I ran your data file
both through MET and IDV and see that IDV plots it correctly while MET
does not.
I strongly suspect that MET has a bug in the processing of the GRIB
grid description section. I'll work on narrowing down the problem,
fixing it, and posting a bugfix next week.
I will email when a bugfix for this issue has been posted to the MET
website.
Thanks for letting us know about this problem!
John
On Fri, May 13, 2016 at 3:52 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
>
> Hello,
> As requested I posted 3 GRIB files (one for each of our domains) and
2
> png images.
> The one named config_domain show the 3 domains.
> d02.png instead is been obtained plotting with plot_data_plane the
> GRIB file for the second domain obtained with UPP starting from the
> same wrf output files but leaving the coordinate in lambert
conformal
> and not transforming it in regular lat lon: in this case domain
> localization is correct
>
> Concerning the version of MET I'm using is the 5.1 with bug-fix
> available at the date of 7 january 2016.
>
> Thanks
> Irene
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail : irene.gallai at arpa.fvg.it
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
> suoi allegati. Se non siete i destinatari indicati nel messaggio, o
> responsabili per la sua consegna alla persona, o se avete ricevuto
il
> messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a
> nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed
i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
> or responsible for message delivering to that person, or if you have
> received this message in error, you may not transcribe, copy or
> deliver this message to anyone. In that case, you should
> delete/destroy this message and its attachments. Thank you .
>
>
> -----Messaggio originale-----
> Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Inviato: giovedì 12 maggio 2016 22:40
> A: Gallai Irene
> Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
>
> Hello Irene,
>
> I see you're having trouble with the grid definition in a GRIB file.
> Could you please post a sample GRIB file to our anonymous ftp site
> following these instructions?
>
> http://www.dtcenter.org/met/users/support/met_help.php#ftp
>
> I'd like to look at the grid definition in the file, plot is using
> plot_data_plane (as you've done), and also plot it using Unidata's
IDV
> tool. Also, do you have a plot showing the "intended" domain?
>
> Hopefully that info will help us track down the discrepancy you're
seeing.
>
> Lastly, please tell me the version of MET you're using.
>
> Thanks,
> John
>
>
> On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
> <met_help at ucar.edu>
> wrote:
>
> >
> > Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> > Transaction: Ticket created by irene.gallai at arpa.fvg.it
> > Queue: met_help
> > Subject: pointstat problem
> > Owner: Nobody
> > Requestors: irene.gallai at arpa.fvg.it
> > Status: new
> > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
> > >
> >
> >
> >
> > Dear All,
> >
> > I'm facing some troubles with pointstat routine, in detail I
> > couldn't obtain matching pairs using my data.
> > Using verbose debugging I see that Grid definition is not correct
> > because for example as lower left (ll) latitude is indicated a
> > value that in effects is my upper level value , as you can see in
> > the following screenshots that are taken from pointstat log the
> > first and a call to cdo libraries the latter.
> >
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------
> > --
> > ----------
> > DEBUG 2:
> > DEBUG 2: Reading data for TMP/Z2.
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range match
> > for VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> > "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> > f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
> > records matching VarInfo "TMP/Z2" in GRIB file
> > "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> > 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > DEBUG 3: Use the matching forecast and observation grids.
> > DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90
lat_ll:
> > 46.920
> > lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025 DEBUG 2: For
> > TMP/Z2 found 1 forecast levels and 0 climatology levels.
> > DEBUG 2:
> > DEBUG 2:
> >
--------------------------------------------------------------------
> > --
> > ----------
> > .....
> >
> > Whereas interrogating the same GRIB file with cdo I obtain
> > Grid coordinates :
> > 1 : lonlat > size : dim = 8100 nx = 90 ny = 90
> > lon : first = 11.888 last = 14.149
inc =
> > 0.0254044944 degrees_east
> > lat : first = 45.322 last = 46.92
inc =
> > 0.0179550562 degrees_north
> >
> >
> > I'm missing something obvious? GRIB files come from UPP (remapping
> > the grid from lambert conformal to regular lon lat)
> >
> > It's more evident looking at the image produced with
plot_data_plane :
> > gray/blue colors should be the Alps in Italy Whereas political
> > boundaries show part of Denmark , Sweden ...
> > Any advice?
> >
> > Thanks
> > Irene
> >
> >
> > +------------------------------------+
> > WE ALL are environmentally responsible. By not printing this
email
> > you've helped save paper, ink, and millions of trees.
> > +------------------------------------+
> >
> > Irene Gallai
> >
> > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > SOC Direzione Tecnica Scientifica
> > CRMA - Centro Regionale di Modellistica Ambientale
> >
> > Via Cairoli 14
> > I-33057 Palmanova (UD) ITALY
> > tel pers.: +39 0432 1918044
> > tel centr.: +39 0432 1918120
> > e-mail :
> > irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> > +------------------------------------+
> >
> >
> > AVVISO DI RISERVATEZZA
> > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > suoi allegati. Se non siete i destinatari indicati nel messaggio,
o
> > responsabili per la sua consegna alla persona, o se avete ricevuto
> > il messaggio per errore, non dovete trascriverlo, copiarlo o
> > inviarlo a nessuno. In tal caso, dovete cancellare/distruggere il
> > messaggio ed i
> suoi allegati.
> > Grazie."
> >
> > CONFIDENTIALITY NOTICE
> > "Confidential information may be contained in this message or in
its
> > attachments. If you are not the addressee indicated in this
message,
> > or responsible for message delivering to that person, or if you
have
> > received this message in error, you may not transcribe, copy or
> > deliver this message to anyone. In that case, you should
> > delete/destroy this message and its attachments. Thank you .
> >
> >
> > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > contenute nel messaggio o nei suoi allegati. Se non siete i
> > destinatari indicati nel messaggio, o responsabili per la sua
> > consegna alla persona, o se avete ricevuto il messaggio per
errore,
> > siete pregati di non trascriverlo, copiarlo o inviarlo ad alcuno.
In
> > tal caso vi invitiamo a cancellare il messaggio ed i suoi
allegati. Grazie.
> > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > this message or in its attachments. If you are not the addressee
> > indicated in this message, or responsible for message delivering
to
> > that person, or if you have received this message in error, you
may
> > not transcribe, copy or deliver this message to anyone. In that
> > case, you should delete this message and its attachments. Thank
you.
> >
> >
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute
> nel messaggio o nei suoi allegati. Se non siete i destinatari
indicati
> nel messaggio, o responsabili per la sua consegna alla persona, o se
> avete ricevuto il messaggio per errore, siete pregati di non
> trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi
invitiamo
> a cancellare il messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
> this message or in its attachments. If you are not the addressee
> indicated in this message, or responsible for message delivering to
> that person, or if you have received this message in error, you may
> not transcribe, copy or deliver this message to anyone. In that
case,
> you should delete this message and its attachments. Thank you.
>
>
>
AVVISO DI RISERVATEZZA Informazioni riservate possono essere contenute
nel messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel messaggio, o responsabili per la sua consegna alla persona, o se
avete ricevuto il messaggio per errore, siete pregati di non
trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo
a cancellare il messaggio ed i suoi allegati. Grazie.
CONFIDENTIALITY NOTICE Confidential information may be contained in
this message or in its attachments. If you are not the addressee
indicated in this message, or responsible for message delivering to
that person, or if you have received this message in error, you may
not transcribe, copy or deliver this message to anyone. In that case,
you should delete this message and its attachments. Thank you.
------------------------------------------------
Subject: pointstat problem
From: John Halley Gotway
Time: Mon May 16 13:22:55 2016
Irene,
After looking very closely, I suspect that the problem lies in how
your
GRIB data is encoded. Specifically, I believe the scanning mode flag
in
octet 28 is set wrong. It is set to 0 but should be set to 64.
Here's some output from wgrib for your sample data:
wgrib -GDS10 CRMA_WRF-4-CRMA_d02_2016050500+000000h_TMP_SFC.GRIB
1:0:d=16050500:TMP:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0:GDS10=
0 0 32 0 255 0 0 135 0 155 0 137 118 0 17 211 136 0 191 127 0 84 81 0
127 0
89 0 0 0 0 0
The 28th entry in the GDS list is set to 0. We use the following
table to
interpret that value of 0:
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table8.html
Since all of the bits are set to 0, that means the data is listed in
the +x
and -y directions. However, I believe your data is actually listed in
the
+x and +y directions.
By way of comparison, I ran both your lat/lon data and some GFS
lat/lon
data through an NCL utility named "ncl_convert2nc". Then I used
ncview to
see the result. Note that both your data and GFS have the scan flag
set to
0. However, the attached ncview images (gfs_globe.png and
d02_italy.png)
show that your data is "upside" down.
When I manually edit the MET code to reset the scanning mode flag to
64,
the resulting images from MET look great.
So at this point, I suspect there isn't a bug in MET, but a problem in
how
your GRIB record in encoded.
Can you please tell me how these GRIB files were created?
Thanks,
John
On Mon, May 16, 2016 at 1:40 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
>
> John,
> many thanks to you,
> I'm happy to contribute to MET development, even if only discovering
> bugs...
> I will be waiting for the fix.
>
> Thank you again
> Irene
>
>
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail : irene.gallai at arpa.fvg.it
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
suoi
> allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili
> per la sua consegna alla persona, o se avete ricevuto il messaggio
per
> errore, non dovete trascriverlo, copiarlo o inviarlo a nessuno. In
tal
> caso, dovete cancellare/distruggere il messaggio ed i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
or
> responsible for message delivering to that person, or if you have
received
> this message in error, you may not transcribe, copy or deliver this
message
> to anyone. In that case, you should delete/destroy this message and
its
> attachments. Thank you .
>
> -----Messaggio originale-----
> Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Inviato: sabato 14 maggio 2016 01:18
> A: Gallai Irene
> Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
>
> Irene,
>
> Thank you for sending me those sample data files. I took a look and
am
> able to reproduce the problem you're seeing. I ran your data file
both
> through MET and IDV and see that IDV plots it correctly while MET
does not.
>
> I strongly suspect that MET has a bug in the processing of the GRIB
grid
> description section. I'll work on narrowing down the problem,
fixing it,
> and posting a bugfix next week.
>
> I will email when a bugfix for this issue has been posted to the MET
> website.
>
> Thanks for letting us know about this problem!
>
> John
>
> On Fri, May 13, 2016 at 3:52 AM, Gallai Irene via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
> >
> > Hello,
> > As requested I posted 3 GRIB files (one for each of our domains)
and 2
> > png images.
> > The one named config_domain show the 3 domains.
> > d02.png instead is been obtained plotting with plot_data_plane
the
> > GRIB file for the second domain obtained with UPP starting from
the
> > same wrf output files but leaving the coordinate in lambert
conformal
> > and not transforming it in regular lat lon: in this case domain
> > localization is correct
> >
> > Concerning the version of MET I'm using is the 5.1 with bug-fix
> > available at the date of 7 january 2016.
> >
> > Thanks
> > Irene
> > +------------------------------------+
> > WE ALL are environmentally responsible. By not printing this
email
> > you've helped save paper, ink, and millions of trees.
> > +------------------------------------+
> >
> > Irene Gallai
> >
> > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > SOC Direzione Tecnica Scientifica
> > CRMA - Centro Regionale di Modellistica Ambientale
> >
> > Via Cairoli 14
> > I-33057 Palmanova (UD) ITALY
> > tel pers.: +39 0432 1918044
> > tel centr.: +39 0432 1918120
> > e-mail : irene.gallai at arpa.fvg.it
> > +------------------------------------+
> >
> >
> > AVVISO DI RISERVATEZZA
> > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > suoi allegati. Se non siete i destinatari indicati nel messaggio,
o
> > responsabili per la sua consegna alla persona, o se avete ricevuto
il
> > messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo
a
> > nessuno. In tal caso, dovete cancellare/distruggere il messaggio
ed i
> suoi allegati.
> > Grazie."
> >
> > CONFIDENTIALITY NOTICE
> > "Confidential information may be contained in this message or in
its
> > attachments. If you are not the addressee indicated in this
message,
> > or responsible for message delivering to that person, or if you
have
> > received this message in error, you may not transcribe, copy or
> > deliver this message to anyone. In that case, you should
> > delete/destroy this message and its attachments. Thank you .
> >
> >
> > -----Messaggio originale-----
> > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Inviato: giovedì 12 maggio 2016 22:40
> > A: Gallai Irene
> > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> >
> > Hello Irene,
> >
> > I see you're having trouble with the grid definition in a GRIB
file.
> > Could you please post a sample GRIB file to our anonymous ftp site
> > following these instructions?
> >
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > I'd like to look at the grid definition in the file, plot is using
> > plot_data_plane (as you've done), and also plot it using Unidata's
IDV
> > tool. Also, do you have a plot showing the "intended" domain?
> >
> > Hopefully that info will help us track down the discrepancy you're
> seeing.
> >
> > Lastly, please tell me the version of MET you're using.
> >
> > Thanks,
> > John
> >
> >
> > On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
> > <met_help at ucar.edu>
> > wrote:
> >
> > >
> > > Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> > > Transaction: Ticket created by irene.gallai at arpa.fvg.it
> > > Queue: met_help
> > > Subject: pointstat problem
> > > Owner: Nobody
> > > Requestors: irene.gallai at arpa.fvg.it
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
> > > >
> > >
> > >
> > >
> > > Dear All,
> > >
> > > I'm facing some troubles with pointstat routine, in detail I
> > > couldn't obtain matching pairs using my data.
> > > Using verbose debugging I see that Grid definition is not
correct
> > > because for example as lower left (ll) latitude is indicated a
> > > value that in effects is my upper level value , as you can see
in
> > > the following screenshots that are taken from pointstat log the
> > > first and a call to cdo libraries the latter.
> > >
> > > DEBUG 2:
> > > DEBUG 2:
> > >
--------------------------------------------------------------------
> > > --
> > > ----------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for TMP/Z2.
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range
match
> > > for VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> > > "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> > > f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
> > > records matching VarInfo "TMP/Z2" in GRIB file
> > > "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> > > 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > DEBUG 3: Use the matching forecast and observation grids.
> > > DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90
lat_ll:
> > > 46.920
> > > lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025 DEBUG 2: For
> > > TMP/Z2 found 1 forecast levels and 0 climatology levels.
> > > DEBUG 2:
> > > DEBUG 2:
> > >
--------------------------------------------------------------------
> > > --
> > > ----------
> > > .....
> > >
> > > Whereas interrogating the same GRIB file with cdo I obtain
> > > Grid coordinates :
> > > 1 : lonlat > size : dim = 8100 nx = 90 ny = 90
> > > lon : first = 11.888 last =
14.149 inc
> =
> > > 0.0254044944 degrees_east
> > > lat : first = 45.322 last = 46.92
inc =
> > > 0.0179550562 degrees_north
> > >
> > >
> > > I'm missing something obvious? GRIB files come from UPP
(remapping
> > > the grid from lambert conformal to regular lon lat)
> > >
> > > It's more evident looking at the image produced with
plot_data_plane :
> > > gray/blue colors should be the Alps in Italy Whereas political
> > > boundaries show part of Denmark , Sweden ...
> > > Any advice?
> > >
> > > Thanks
> > > Irene
> > >
> > >
> > > +------------------------------------+
> > > WE ALL are environmentally responsible. By not printing this
email
> > > you've helped save paper, ink, and millions of trees.
> > > +------------------------------------+
> > >
> > > Irene Gallai
> > >
> > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > > SOC Direzione Tecnica Scientifica
> > > CRMA - Centro Regionale di Modellistica Ambientale
> > >
> > > Via Cairoli 14
> > > I-33057 Palmanova (UD) ITALY
> > > tel pers.: +39 0432 1918044
> > > tel centr.: +39 0432 1918120
> > > e-mail :
> > > irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> > > +------------------------------------+
> > >
> > >
> > > AVVISO DI RISERVATEZZA
> > > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > > suoi allegati. Se non siete i destinatari indicati nel
messaggio, o
> > > responsabili per la sua consegna alla persona, o se avete
ricevuto
> > > il messaggio per errore, non dovete trascriverlo, copiarlo o
> > > inviarlo a nessuno. In tal caso, dovete cancellare/distruggere
il
> > > messaggio ed i
> > suoi allegati.
> > > Grazie."
> > >
> > > CONFIDENTIALITY NOTICE
> > > "Confidential information may be contained in this message or in
its
> > > attachments. If you are not the addressee indicated in this
message,
> > > or responsible for message delivering to that person, or if you
have
> > > received this message in error, you may not transcribe, copy or
> > > deliver this message to anyone. In that case, you should
> > > delete/destroy this message and its attachments. Thank you .
> > >
> > >
> > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > destinatari indicati nel messaggio, o responsabili per la sua
> > > consegna alla persona, o se avete ricevuto il messaggio per
errore,
> > > siete pregati di non trascriverlo, copiarlo o inviarlo ad
alcuno. In
> > > tal caso vi invitiamo a cancellare il messaggio ed i suoi
allegati.
> Grazie.
> > > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > > this message or in its attachments. If you are not the addressee
> > > indicated in this message, or responsible for message delivering
to
> > > that person, or if you have received this message in error, you
may
> > > not transcribe, copy or deliver this message to anyone. In that
> > > case, you should delete this message and its attachments. Thank
you.
> > >
> > >
> >
> >
> > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute
> > nel messaggio o nei suoi allegati. Se non siete i destinatari
indicati
> > nel messaggio, o responsabili per la sua consegna alla persona, o
se
> > avete ricevuto il messaggio per errore, siete pregati di non
> > trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi
invitiamo
> > a cancellare il messaggio ed i suoi allegati. Grazie.
> > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > this message or in its attachments. If you are not the addressee
> > indicated in this message, or responsible for message delivering
to
> > that person, or if you have received this message in error, you
may
> > not transcribe, copy or deliver this message to anyone. In that
case,
> > you should delete this message and its attachments. Thank you.
> >
> >
> >
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute nel
> messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel
> messaggio, o responsabili per la sua consegna alla persona, o se
avete
> ricevuto il messaggio per errore, siete pregati di non trascriverlo,
> copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo a cancellare
il
> messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
this
> message or in its attachments. If you are not the addressee
indicated in
> this message, or responsible for message delivering to that person,
or if
> you have received this message in error, you may not transcribe,
copy or
> deliver this message to anyone. In that case, you should delete this
> message and its attachments. Thank you.
>
>
>
------------------------------------------------
Subject: R: [rt.rap.ucar.edu #76349] pointstat problem
From: Gallai Irene
Time: Wed May 18 05:19:36 2016
Hi Jhon,
We create GRIB files with UPP and then we change map projection using
the copygb utility. Maybe the "copygb step" is the problem.
Unfortunately I'm not the one that make those passages so I don't know
If and which parameters and arguments have been used .
I asked a colleague to check the procedure to see where and how the
scanning mode flag in octet 28 is set.
I'll let you know when I have an answer
For the moment thanks again
Irene
+------------------------------------+
WE ALL are environmentally responsible. By not printing this email
you've helped save paper, ink, and millions of trees.
+------------------------------------+
Irene Gallai
Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
SOC Direzione Tecnica Scientifica
CRMA - Centro Regionale di Modellistica Ambientale
Via Cairoli 14
I-33057 Palmanova (UD) ITALY
tel pers.: +39 0432 1918044
tel centr.: +39 0432 1918120
e-mail : irene.gallai at arpa.fvg.it
+------------------------------------+
AVVISO DI RISERVATEZZA
"Informazioni riservate possono essere contenute nel messaggio o nei
suoi allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili per la sua consegna alla persona, o se avete ricevuto il
messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a
nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed i
suoi allegati. Grazie."
CONFIDENTIALITY NOTICE
"Confidential information may be contained in this message or in its
attachments. If you are not the addressee indicated in this message,
or responsible for message delivering to that person, or if you have
received this message in error, you may not transcribe, copy or
deliver this message to anyone. In that case, you should
delete/destroy this message and its attachments. Thank you .
-----Messaggio originale-----
Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Inviato: lunedì 16 maggio 2016 21:23
A: Gallai Irene
Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
Irene,
After looking very closely, I suspect that the problem lies in how
your GRIB data is encoded. Specifically, I believe the scanning mode
flag in octet 28 is set wrong. It is set to 0 but should be set to
64.
Here's some output from wgrib for your sample data:
wgrib -GDS10 CRMA_WRF-4-CRMA_d02_2016050500+000000h_TMP_SFC.GRIB
1:0:d=16050500:TMP:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0:GDS10=
0 0 32 0 255 0 0 135 0 155 0 137 118 0 17 211 136 0 191 127 0 84 81 0
127 0
89 0 0 0 0 0
The 28th entry in the GDS list is set to 0. We use the following
table to interpret that value of 0:
http://www.nco.ncep.noaa.gov/pmb/docs/on388/table8.html
Since all of the bits are set to 0, that means the data is listed in
the +x and -y directions. However, I believe your data is actually
listed in the
+x and +y directions.
By way of comparison, I ran both your lat/lon data and some GFS
lat/lon data through an NCL utility named "ncl_convert2nc". Then I
used ncview to see the result. Note that both your data and GFS have
the scan flag set to 0. However, the attached ncview images
(gfs_globe.png and d02_italy.png) show that your data is "upside"
down.
When I manually edit the MET code to reset the scanning mode flag to
64, the resulting images from MET look great.
So at this point, I suspect there isn't a bug in MET, but a problem in
how your GRIB record in encoded.
Can you please tell me how these GRIB files were created?
Thanks,
John
On Mon, May 16, 2016 at 1:40 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
>
> John,
> many thanks to you,
> I'm happy to contribute to MET development, even if only discovering
> bugs...
> I will be waiting for the fix.
>
> Thank you again
> Irene
>
>
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail : irene.gallai at arpa.fvg.it
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
> suoi allegati. Se non siete i destinatari indicati nel messaggio, o
> responsabili per la sua consegna alla persona, o se avete ricevuto
il
> messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a
> nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed
i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
> or responsible for message delivering to that person, or if you have
> received this message in error, you may not transcribe, copy or
> deliver this message to anyone. In that case, you should
> delete/destroy this message and its attachments. Thank you .
>
> -----Messaggio originale-----
> Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Inviato: sabato 14 maggio 2016 01:18
> A: Gallai Irene
> Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
>
> Irene,
>
> Thank you for sending me those sample data files. I took a look and
> am able to reproduce the problem you're seeing. I ran your data
file
> both through MET and IDV and see that IDV plots it correctly while
MET does not.
>
> I strongly suspect that MET has a bug in the processing of the GRIB
> grid description section. I'll work on narrowing down the problem,
> fixing it, and posting a bugfix next week.
>
> I will email when a bugfix for this issue has been posted to the MET
> website.
>
> Thanks for letting us know about this problem!
>
> John
>
> On Fri, May 13, 2016 at 3:52 AM, Gallai Irene via RT
> <met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
> >
> > Hello,
> > As requested I posted 3 GRIB files (one for each of our domains)
and
> > 2 png images.
> > The one named config_domain show the 3 domains.
> > d02.png instead is been obtained plotting with plot_data_plane
the
> > GRIB file for the second domain obtained with UPP starting from
the
> > same wrf output files but leaving the coordinate in lambert
> > conformal and not transforming it in regular lat lon: in this case
> > domain localization is correct
> >
> > Concerning the version of MET I'm using is the 5.1 with bug-fix
> > available at the date of 7 january 2016.
> >
> > Thanks
> > Irene
> > +------------------------------------+
> > WE ALL are environmentally responsible. By not printing this
email
> > you've helped save paper, ink, and millions of trees.
> > +------------------------------------+
> >
> > Irene Gallai
> >
> > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > SOC Direzione Tecnica Scientifica
> > CRMA - Centro Regionale di Modellistica Ambientale
> >
> > Via Cairoli 14
> > I-33057 Palmanova (UD) ITALY
> > tel pers.: +39 0432 1918044
> > tel centr.: +39 0432 1918120
> > e-mail : irene.gallai at arpa.fvg.it
> > +------------------------------------+
> >
> >
> > AVVISO DI RISERVATEZZA
> > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > suoi allegati. Se non siete i destinatari indicati nel messaggio,
o
> > responsabili per la sua consegna alla persona, o se avete ricevuto
> > il messaggio per errore, non dovete trascriverlo, copiarlo o
> > inviarlo a nessuno. In tal caso, dovete cancellare/distruggere il
> > messaggio ed i
> suoi allegati.
> > Grazie."
> >
> > CONFIDENTIALITY NOTICE
> > "Confidential information may be contained in this message or in
its
> > attachments. If you are not the addressee indicated in this
message,
> > or responsible for message delivering to that person, or if you
have
> > received this message in error, you may not transcribe, copy or
> > deliver this message to anyone. In that case, you should
> > delete/destroy this message and its attachments. Thank you .
> >
> >
> > -----Messaggio originale-----
> > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Inviato: giovedì 12 maggio 2016 22:40
> > A: Gallai Irene
> > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> >
> > Hello Irene,
> >
> > I see you're having trouble with the grid definition in a GRIB
file.
> > Could you please post a sample GRIB file to our anonymous ftp site
> > following these instructions?
> >
> > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> >
> > I'd like to look at the grid definition in the file, plot is using
> > plot_data_plane (as you've done), and also plot it using Unidata's
> > IDV tool. Also, do you have a plot showing the "intended" domain?
> >
> > Hopefully that info will help us track down the discrepancy you're
> seeing.
> >
> > Lastly, please tell me the version of MET you're using.
> >
> > Thanks,
> > John
> >
> >
> > On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
> > <met_help at ucar.edu>
> > wrote:
> >
> > >
> > > Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> > > Transaction: Ticket created by irene.gallai at arpa.fvg.it
> > > Queue: met_help
> > > Subject: pointstat problem
> > > Owner: Nobody
> > > Requestors: irene.gallai at arpa.fvg.it
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
> > > >
> > >
> > >
> > >
> > > Dear All,
> > >
> > > I'm facing some troubles with pointstat routine, in detail I
> > > couldn't obtain matching pairs using my data.
> > > Using verbose debugging I see that Grid definition is not
correct
> > > because for example as lower left (ll) latitude is indicated a
> > > value that in effects is my upper level value , as you can see
in
> > > the following screenshots that are taken from pointstat log the
> > > first and a call to cdo libraries the latter.
> > >
> > > DEBUG 2:
> > > DEBUG 2:
> > >
------------------------------------------------------------------
> > > --
> > > --
> > > ----------
> > > DEBUG 2:
> > > DEBUG 2: Reading data for TMP/Z2.
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range
match
> > > for VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> > > "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> > > f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
> > > records matching VarInfo "TMP/Z2" in GRIB file
> > > "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> > > 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > DEBUG 3: Use the matching forecast and observation grids.
> > > DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90
lat_ll:
> > > 46.920
> > > lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025 DEBUG 2: For
> > > TMP/Z2 found 1 forecast levels and 0 climatology levels.
> > > DEBUG 2:
> > > DEBUG 2:
> > >
------------------------------------------------------------------
> > > --
> > > --
> > > ----------
> > > .....
> > >
> > > Whereas interrogating the same GRIB file with cdo I obtain
> > > Grid coordinates :
> > > 1 : lonlat > size : dim = 8100 nx = 90 ny = 90
> > > lon : first = 11.888 last =
14.149 inc
> =
> > > 0.0254044944 degrees_east
> > > lat : first = 45.322 last = 46.92
inc =
> > > 0.0179550562 degrees_north
> > >
> > >
> > > I'm missing something obvious? GRIB files come from UPP
(remapping
> > > the grid from lambert conformal to regular lon lat)
> > >
> > > It's more evident looking at the image produced with
plot_data_plane :
> > > gray/blue colors should be the Alps in Italy Whereas political
> > > boundaries show part of Denmark , Sweden ...
> > > Any advice?
> > >
> > > Thanks
> > > Irene
> > >
> > >
> > > +------------------------------------+
> > > WE ALL are environmentally responsible. By not printing this
> > > email you've helped save paper, ink, and millions of trees.
> > > +------------------------------------+
> > >
> > > Irene Gallai
> > >
> > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > > SOC Direzione Tecnica Scientifica
> > > CRMA - Centro Regionale di Modellistica Ambientale
> > >
> > > Via Cairoli 14
> > > I-33057 Palmanova (UD) ITALY
> > > tel pers.: +39 0432 1918044
> > > tel centr.: +39 0432 1918120
> > > e-mail :
> > > irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> > > +------------------------------------+
> > >
> > >
> > > AVVISO DI RISERVATEZZA
> > > "Informazioni riservate possono essere contenute nel messaggio o
> > > nei suoi allegati. Se non siete i destinatari indicati nel
> > > messaggio, o responsabili per la sua consegna alla persona, o se
> > > avete ricevuto il messaggio per errore, non dovete trascriverlo,
> > > copiarlo o inviarlo a nessuno. In tal caso, dovete
> > > cancellare/distruggere il messaggio ed i
> > suoi allegati.
> > > Grazie."
> > >
> > > CONFIDENTIALITY NOTICE
> > > "Confidential information may be contained in this message or in
> > > its attachments. If you are not the addressee indicated in this
> > > message, or responsible for message delivering to that person,
or
> > > if you have received this message in error, you may not
> > > transcribe, copy or deliver this message to anyone. In that
case,
> > > you should delete/destroy this message and its attachments.
Thank you .
> > >
> > >
> > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > destinatari indicati nel messaggio, o responsabili per la sua
> > > consegna alla persona, o se avete ricevuto il messaggio per
> > > errore, siete pregati di non trascriverlo, copiarlo o inviarlo
ad
> > > alcuno. In tal caso vi invitiamo a cancellare il messaggio ed i
suoi allegati.
> Grazie.
> > > CONFIDENTIALITY NOTICE Confidential information may be contained
> > > in this message or in its attachments. If you are not the
> > > addressee indicated in this message, or responsible for message
> > > delivering to that person, or if you have received this message
in
> > > error, you may not transcribe, copy or deliver this message to
> > > anyone. In that case, you should delete this message and its
attachments. Thank you.
> > >
> > >
> >
> >
> > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > contenute nel messaggio o nei suoi allegati. Se non siete i
> > destinatari indicati nel messaggio, o responsabili per la sua
> > consegna alla persona, o se avete ricevuto il messaggio per
errore,
> > siete pregati di non trascriverlo, copiarlo o inviarlo ad alcuno.
In
> > tal caso vi invitiamo a cancellare il messaggio ed i suoi
allegati. Grazie.
> > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > this message or in its attachments. If you are not the addressee
> > indicated in this message, or responsible for message delivering
to
> > that person, or if you have received this message in error, you
may
> > not transcribe, copy or deliver this message to anyone. In that
> > case, you should delete this message and its attachments. Thank
you.
> >
> >
> >
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute
> nel messaggio o nei suoi allegati. Se non siete i destinatari
indicati
> nel messaggio, o responsabili per la sua consegna alla persona, o se
> avete ricevuto il messaggio per errore, siete pregati di non
> trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi
invitiamo
> a cancellare il messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
> this message or in its attachments. If you are not the addressee
> indicated in this message, or responsible for message delivering to
> that person, or if you have received this message in error, you may
> not transcribe, copy or deliver this message to anyone. In that
case,
> you should delete this message and its attachments. Thank you.
>
>
>
AVVISO DI RISERVATEZZA Informazioni riservate possono essere contenute
nel messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel messaggio, o responsabili per la sua consegna alla persona, o se
avete ricevuto il messaggio per errore, siete pregati di non
trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo
a cancellare il messaggio ed i suoi allegati. Grazie.
CONFIDENTIALITY NOTICE Confidential information may be contained in
this message or in its attachments. If you are not the addressee
indicated in this message, or responsible for message delivering to
that person, or if you have received this message in error, you may
not transcribe, copy or deliver this message to anyone. In that case,
you should delete this message and its attachments. Thank you.
------------------------------------------------
Subject: pointstat problem
From: John Halley Gotway
Time: Wed May 18 08:30:58 2016
Irene,
Hmm, that's interesting. Do you be any chance have a sample GRIB
output
file from UPP before you run it through coypgb? Or even a sample WRF
NetCDF output file prior to running through UPP? If you post one of
those
along with the UPP namelist or copygb command you use, I could take a
look
at the process.
The DTC is also responsible for supporting UPP, which includes copygb.
So
I'd like to figure out where the problem is coming in.
And it is still possible that MET isn't processing the grid
information
correctly. It is suspicious that changing the scanning mode flag in
octet
28 to 64 fixes the problem. But it's also suspicious that IDV is able
to
display the data correctly
John
On Wed, May 18, 2016 at 5:19 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
>
> Hi Jhon,
>
> We create GRIB files with UPP and then we change map projection
using the
> copygb utility. Maybe the "copygb step" is the problem.
> Unfortunately I'm not the one that make those passages so I don't
know If
> and which parameters and arguments have been used .
> I asked a colleague to check the procedure to see where and how the
> scanning mode flag in octet 28 is set.
> I'll let you know when I have an answer
> For the moment thanks again
>
> Irene
>
>
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail : irene.gallai at arpa.fvg.it
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
suoi
> allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili
> per la sua consegna alla persona, o se avete ricevuto il messaggio
per
> errore, non dovete trascriverlo, copiarlo o inviarlo a nessuno. In
tal
> caso, dovete cancellare/distruggere il messaggio ed i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
or
> responsible for message delivering to that person, or if you have
received
> this message in error, you may not transcribe, copy or deliver this
message
> to anyone. In that case, you should delete/destroy this message and
its
> attachments. Thank you .
>
>
> -----Messaggio originale-----
> Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Inviato: lunedì 16 maggio 2016 21:23
> A: Gallai Irene
> Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
>
> Irene,
>
> After looking very closely, I suspect that the problem lies in how
your
> GRIB data is encoded. Specifically, I believe the scanning mode
flag in
> octet 28 is set wrong. It is set to 0 but should be set to 64.
>
> Here's some output from wgrib for your sample data:
> wgrib -GDS10 CRMA_WRF-4-CRMA_d02_2016050500+000000h_TMP_SFC.GRIB
>
>
1:0:d=16050500:TMP:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc:anl:NAve=0:GDS10=
> 0 0 32 0 255 0 0 135 0 155 0 137 118 0 17 211 136 0 191 127 0 84 81
0 127 0
> 89 0 0 0 0 0
>
> The 28th entry in the GDS list is set to 0. We use the following
table to
> interpret that value of 0:
> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table8.html
>
> Since all of the bits are set to 0, that means the data is listed in
the
> +x and -y directions. However, I believe your data is actually
listed in
> the
> +x and +y directions.
>
> By way of comparison, I ran both your lat/lon data and some GFS
lat/lon
> data through an NCL utility named "ncl_convert2nc". Then I used
ncview to
> see the result. Note that both your data and GFS have the scan flag
set to
> 0. However, the attached ncview images (gfs_globe.png and
d02_italy.png)
> show that your data is "upside" down.
>
> When I manually edit the MET code to reset the scanning mode flag to
64,
> the resulting images from MET look great.
>
> So at this point, I suspect there isn't a bug in MET, but a problem
in how
> your GRIB record in encoded.
>
> Can you please tell me how these GRIB files were created?
>
> Thanks,
> John
>
>
> On Mon, May 16, 2016 at 1:40 AM, Gallai Irene via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
> >
> > John,
> > many thanks to you,
> > I'm happy to contribute to MET development, even if only
discovering
> > bugs...
> > I will be waiting for the fix.
> >
> > Thank you again
> > Irene
> >
> >
> > +------------------------------------+
> > WE ALL are environmentally responsible. By not printing this
email
> > you've helped save paper, ink, and millions of trees.
> > +------------------------------------+
> >
> > Irene Gallai
> >
> > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > SOC Direzione Tecnica Scientifica
> > CRMA - Centro Regionale di Modellistica Ambientale
> >
> > Via Cairoli 14
> > I-33057 Palmanova (UD) ITALY
> > tel pers.: +39 0432 1918044
> > tel centr.: +39 0432 1918120
> > e-mail : irene.gallai at arpa.fvg.it
> > +------------------------------------+
> >
> >
> > AVVISO DI RISERVATEZZA
> > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > suoi allegati. Se non siete i destinatari indicati nel messaggio,
o
> > responsabili per la sua consegna alla persona, o se avete ricevuto
il
> > messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo
a
> > nessuno. In tal caso, dovete cancellare/distruggere il messaggio
ed i
> suoi allegati.
> > Grazie."
> >
> > CONFIDENTIALITY NOTICE
> > "Confidential information may be contained in this message or in
its
> > attachments. If you are not the addressee indicated in this
message,
> > or responsible for message delivering to that person, or if you
have
> > received this message in error, you may not transcribe, copy or
> > deliver this message to anyone. In that case, you should
> > delete/destroy this message and its attachments. Thank you .
> >
> > -----Messaggio originale-----
> > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Inviato: sabato 14 maggio 2016 01:18
> > A: Gallai Irene
> > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> >
> > Irene,
> >
> > Thank you for sending me those sample data files. I took a look
and
> > am able to reproduce the problem you're seeing. I ran your data
file
> > both through MET and IDV and see that IDV plots it correctly while
MET
> does not.
> >
> > I strongly suspect that MET has a bug in the processing of the
GRIB
> > grid description section. I'll work on narrowing down the
problem,
> > fixing it, and posting a bugfix next week.
> >
> > I will email when a bugfix for this issue has been posted to the
MET
> > website.
> >
> > Thanks for letting us know about this problem!
> >
> > John
> >
> > On Fri, May 13, 2016 at 3:52 AM, Gallai Irene via RT
> > <met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
> > >
> > > Hello,
> > > As requested I posted 3 GRIB files (one for each of our domains)
and
> > > 2 png images.
> > > The one named config_domain show the 3 domains.
> > > d02.png instead is been obtained plotting with plot_data_plane
the
> > > GRIB file for the second domain obtained with UPP starting from
the
> > > same wrf output files but leaving the coordinate in lambert
> > > conformal and not transforming it in regular lat lon: in this
case
> > > domain localization is correct
> > >
> > > Concerning the version of MET I'm using is the 5.1 with bug-fix
> > > available at the date of 7 january 2016.
> > >
> > > Thanks
> > > Irene
> > > +------------------------------------+
> > > WE ALL are environmentally responsible. By not printing this
email
> > > you've helped save paper, ink, and millions of trees.
> > > +------------------------------------+
> > >
> > > Irene Gallai
> > >
> > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > > SOC Direzione Tecnica Scientifica
> > > CRMA - Centro Regionale di Modellistica Ambientale
> > >
> > > Via Cairoli 14
> > > I-33057 Palmanova (UD) ITALY
> > > tel pers.: +39 0432 1918044
> > > tel centr.: +39 0432 1918120
> > > e-mail : irene.gallai at arpa.fvg.it
> > > +------------------------------------+
> > >
> > >
> > > AVVISO DI RISERVATEZZA
> > > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > > suoi allegati. Se non siete i destinatari indicati nel
messaggio, o
> > > responsabili per la sua consegna alla persona, o se avete
ricevuto
> > > il messaggio per errore, non dovete trascriverlo, copiarlo o
> > > inviarlo a nessuno. In tal caso, dovete cancellare/distruggere
il
> > > messaggio ed i
> > suoi allegati.
> > > Grazie."
> > >
> > > CONFIDENTIALITY NOTICE
> > > "Confidential information may be contained in this message or in
its
> > > attachments. If you are not the addressee indicated in this
message,
> > > or responsible for message delivering to that person, or if you
have
> > > received this message in error, you may not transcribe, copy or
> > > deliver this message to anyone. In that case, you should
> > > delete/destroy this message and its attachments. Thank you .
> > >
> > >
> > > -----Messaggio originale-----
> > > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Inviato: giovedì 12 maggio 2016 22:40
> > > A: Gallai Irene
> > > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> > >
> > > Hello Irene,
> > >
> > > I see you're having trouble with the grid definition in a GRIB
file.
> > > Could you please post a sample GRIB file to our anonymous ftp
site
> > > following these instructions?
> > >
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > I'd like to look at the grid definition in the file, plot is
using
> > > plot_data_plane (as you've done), and also plot it using
Unidata's
> > > IDV tool. Also, do you have a plot showing the "intended"
domain?
> > >
> > > Hopefully that info will help us track down the discrepancy
you're
> > seeing.
> > >
> > > Lastly, please tell me the version of MET you're using.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
> > > <met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> > > > Transaction: Ticket created by irene.gallai at arpa.fvg.it
> > > > Queue: met_help
> > > > Subject: pointstat problem
> > > > Owner: Nobody
> > > > Requestors: irene.gallai at arpa.fvg.it
> > > > Status: new
> > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
> > > > >
> > > >
> > > >
> > > >
> > > > Dear All,
> > > >
> > > > I'm facing some troubles with pointstat routine, in detail I
> > > > couldn't obtain matching pairs using my data.
> > > > Using verbose debugging I see that Grid definition is not
correct
> > > > because for example as lower left (ll) latitude is indicated
a
> > > > value that in effects is my upper level value , as you can see
in
> > > > the following screenshots that are taken from pointstat log
the
> > > > first and a call to cdo libraries the latter.
> > > >
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
------------------------------------------------------------------
> > > > --
> > > > --
> > > > ----------
> > > > DEBUG 2:
> > > > DEBUG 2: Reading data for TMP/Z2.
> > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range
match
> > > > for VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> > > > "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> > > > f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
> > > > records matching VarInfo "TMP/Z2" in GRIB file
> > > > "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> > > > 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > > DEBUG 3: Use the matching forecast and observation grids.
> > > > DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90
lat_ll:
> > > > 46.920
> > > > lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025 DEBUG 2: For
> > > > TMP/Z2 found 1 forecast levels and 0 climatology levels.
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
------------------------------------------------------------------
> > > > --
> > > > --
> > > > ----------
> > > > .....
> > > >
> > > > Whereas interrogating the same GRIB file with cdo I obtain
> > > > Grid coordinates :
> > > > 1 : lonlat > size : dim = 8100 nx = 90 ny =
90
> > > > lon : first = 11.888 last =
14.149
> inc
> > =
> > > > 0.0254044944 degrees_east
> > > > lat : first = 45.322 last =
46.92
> inc =
> > > > 0.0179550562 degrees_north
> > > >
> > > >
> > > > I'm missing something obvious? GRIB files come from UPP
(remapping
> > > > the grid from lambert conformal to regular lon lat)
> > > >
> > > > It's more evident looking at the image produced with
plot_data_plane
> :
> > > > gray/blue colors should be the Alps in Italy Whereas political
> > > > boundaries show part of Denmark , Sweden ...
> > > > Any advice?
> > > >
> > > > Thanks
> > > > Irene
> > > >
> > > >
> > > > +------------------------------------+
> > > > WE ALL are environmentally responsible. By not printing this
> > > > email you've helped save paper, ink, and millions of trees.
> > > > +------------------------------------+
> > > >
> > > > Irene Gallai
> > > >
> > > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG
)
> > > > SOC Direzione Tecnica Scientifica
> > > > CRMA - Centro Regionale di Modellistica Ambientale
> > > >
> > > > Via Cairoli 14
> > > > I-33057 Palmanova (UD) ITALY
> > > > tel pers.: +39 0432 1918044
> > > > tel centr.: +39 0432 1918120
> > > > e-mail :
> > > > irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> > > > +------------------------------------+
> > > >
> > > >
> > > > AVVISO DI RISERVATEZZA
> > > > "Informazioni riservate possono essere contenute nel messaggio
o
> > > > nei suoi allegati. Se non siete i destinatari indicati nel
> > > > messaggio, o responsabili per la sua consegna alla persona, o
se
> > > > avete ricevuto il messaggio per errore, non dovete
trascriverlo,
> > > > copiarlo o inviarlo a nessuno. In tal caso, dovete
> > > > cancellare/distruggere il messaggio ed i
> > > suoi allegati.
> > > > Grazie."
> > > >
> > > > CONFIDENTIALITY NOTICE
> > > > "Confidential information may be contained in this message or
in
> > > > its attachments. If you are not the addressee indicated in
this
> > > > message, or responsible for message delivering to that person,
or
> > > > if you have received this message in error, you may not
> > > > transcribe, copy or deliver this message to anyone. In that
case,
> > > > you should delete/destroy this message and its attachments.
Thank
> you .
> > > >
> > > >
> > > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > > destinatari indicati nel messaggio, o responsabili per la sua
> > > > consegna alla persona, o se avete ricevuto il messaggio per
> > > > errore, siete pregati di non trascriverlo, copiarlo o inviarlo
ad
> > > > alcuno. In tal caso vi invitiamo a cancellare il messaggio ed
i suoi
> allegati.
> > Grazie.
> > > > CONFIDENTIALITY NOTICE Confidential information may be
contained
> > > > in this message or in its attachments. If you are not the
> > > > addressee indicated in this message, or responsible for
message
> > > > delivering to that person, or if you have received this
message in
> > > > error, you may not transcribe, copy or deliver this message to
> > > > anyone. In that case, you should delete this message and its
> attachments. Thank you.
> > > >
> > > >
> > >
> > >
> > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > destinatari indicati nel messaggio, o responsabili per la sua
> > > consegna alla persona, o se avete ricevuto il messaggio per
errore,
> > > siete pregati di non trascriverlo, copiarlo o inviarlo ad
alcuno. In
> > > tal caso vi invitiamo a cancellare il messaggio ed i suoi
allegati.
> Grazie.
> > > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > > this message or in its attachments. If you are not the addressee
> > > indicated in this message, or responsible for message delivering
to
> > > that person, or if you have received this message in error, you
may
> > > not transcribe, copy or deliver this message to anyone. In that
> > > case, you should delete this message and its attachments. Thank
you.
> > >
> > >
> > >
> >
> >
> > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute
> > nel messaggio o nei suoi allegati. Se non siete i destinatari
indicati
> > nel messaggio, o responsabili per la sua consegna alla persona, o
se
> > avete ricevuto il messaggio per errore, siete pregati di non
> > trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi
invitiamo
> > a cancellare il messaggio ed i suoi allegati. Grazie.
> > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > this message or in its attachments. If you are not the addressee
> > indicated in this message, or responsible for message delivering
to
> > that person, or if you have received this message in error, you
may
> > not transcribe, copy or deliver this message to anyone. In that
case,
> > you should delete this message and its attachments. Thank you.
> >
> >
> >
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute nel
> messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel
> messaggio, o responsabili per la sua consegna alla persona, o se
avete
> ricevuto il messaggio per errore, siete pregati di non trascriverlo,
> copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo a cancellare
il
> messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
this
> message or in its attachments. If you are not the addressee
indicated in
> this message, or responsible for message delivering to that person,
or if
> you have received this message in error, you may not transcribe,
copy or
> deliver this message to anyone. In that case, you should delete this
> message and its attachments. Thank you.
>
>
>
------------------------------------------------
Subject: R: [rt.rap.ucar.edu #76349] pointstat problem
From: Gallai Irene
Time: Thu May 19 05:13:40 2016
Hello Jhon,
I posted a sample GRIB output file from UPP before we run it through
coypgb, it's named CRMA_WRF-4-CRMA_d02_2016051900+000000h.GRIB.tmp.
My colleague also handed to me the copygb command used and I think
that the problem is in there. If I'm not mistaken the last entry of
the kgds string should be the one that will set the scanning mode
and as you can see there is a 0 and not a 64 as you pointed out.
Here is the command line
copygb.exe -xg"255 0 135 155 35190 4563 136 49023 21585 127 89 0"
CRMA_WRF-4-CRMA_d02_2016051900+000000h.GRIB.tmp CRMA_WRF-4-
CRMA_d02_2016051900+000000h.GRIB
Probably it remained some typo in the script from tests but we didn't
detected it because generally we look at outputs with GrADS after post
processing them with grib2ctl and gribmap and the plots never showed
reversed domains. It's possible that some "degribbing utilities"
don't consider the scanning mode flag? As you say it's suspicious
that IDV or GrADS are able to display the data correctly.
In any case we would correct copygb command
Many thanks again
Irene
+------------------------------------+
WE ALL are environmentally responsible. By not printing this email
you've helped save paper, ink, and millions of trees.
+------------------------------------+
Irene Gallai
Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
SOC Direzione Tecnica Scientifica
CRMA - Centro Regionale di Modellistica Ambientale
Via Cairoli 14
I-33057 Palmanova (UD) ITALY
tel pers.: +39 0432 1918044
tel centr.: +39 0432 1918120
e-mail : irene.gallai at arpa.fvg.it
+------------------------------------+
AVVISO DI RISERVATEZZA
"Informazioni riservate possono essere contenute nel messaggio o nei
suoi allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili per la sua consegna alla persona, o se avete ricevuto il
messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a
nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed i
suoi allegati. Grazie."
CONFIDENTIALITY NOTICE
"Confidential information may be contained in this message or in its
attachments. If you are not the addressee indicated in this message,
or responsible for message delivering to that person, or if you have
received this message in error, you may not transcribe, copy or
deliver this message to anyone. In that case, you should
delete/destroy this message and its attachments. Thank you .
-----Messaggio originale-----
Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
Inviato: mercoledì 18 maggio 2016 16:31
A: Gallai Irene
Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
Irene,
Hmm, that's interesting. Do you be any chance have a sample GRIB
output file from UPP before you run it through coypgb? Or even a
sample WRF NetCDF output file prior to running through UPP? If you
post one of those along with the UPP namelist or copygb command you
use, I could take a look at the process.
The DTC is also responsible for supporting UPP, which includes copygb.
So I'd like to figure out where the problem is coming in.
And it is still possible that MET isn't processing the grid
information correctly. It is suspicious that changing the scanning
mode flag in octet
28 to 64 fixes the problem. But it's also suspicious that IDV is able
to display the data correctly
John
On Wed, May 18, 2016 at 5:19 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
>
> Hi Jhon,
>
> We create GRIB files with UPP and then we change map projection
using
> the copygb utility. Maybe the "copygb step" is the problem.
> Unfortunately I'm not the one that make those passages so I don't
know
> If and which parameters and arguments have been used .
> I asked a colleague to check the procedure to see where and how the
> scanning mode flag in octet 28 is set.
> I'll let you know when I have an answer For the moment thanks again
>
> Irene
>
>
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail : irene.gallai at arpa.fvg.it
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
> suoi allegati. Se non siete i destinatari indicati nel messaggio, o
> responsabili per la sua consegna alla persona, o se avete ricevuto
il
> messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo a
> nessuno. In tal caso, dovete cancellare/distruggere il messaggio ed
i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
> or responsible for message delivering to that person, or if you have
> received this message in error, you may not transcribe, copy or
> deliver this message to anyone. In that case, you should
> delete/destroy this message and its attachments. Thank you .
>
>
> -----Messaggio originale-----
> Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Inviato: lunedì 16 maggio 2016 21:23
> A: Gallai Irene
> Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
>
> Irene,
>
> After looking very closely, I suspect that the problem lies in how
> your GRIB data is encoded. Specifically, I believe the scanning
mode
> flag in octet 28 is set wrong. It is set to 0 but should be set to
64.
>
> Here's some output from wgrib for your sample data:
> wgrib -GDS10 CRMA_WRF-4-CRMA_d02_2016050500+000000h_TMP_SFC.GRIB
>
>
1:0:d=16050500:TMP:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc
> :anl:NAve=0:GDS10=
> 0 0 32 0 255 0 0 135 0 155 0 137 118 0 17 211 136 0 191 127 0 84 81
0
> 127 0
> 89 0 0 0 0 0
>
> The 28th entry in the GDS list is set to 0. We use the following
> table to interpret that value of 0:
> http://www.nco.ncep.noaa.gov/pmb/docs/on388/table8.html
>
> Since all of the bits are set to 0, that means the data is listed in
> the
> +x and -y directions. However, I believe your data is actually
listed
> +in
> the
> +x and +y directions.
>
> By way of comparison, I ran both your lat/lon data and some GFS
> lat/lon data through an NCL utility named "ncl_convert2nc". Then I
> used ncview to see the result. Note that both your data and GFS
have
> the scan flag set to 0. However, the attached ncview images
> (gfs_globe.png and d02_italy.png) show that your data is "upside"
down.
>
> When I manually edit the MET code to reset the scanning mode flag to
> 64, the resulting images from MET look great.
>
> So at this point, I suspect there isn't a bug in MET, but a problem
in
> how your GRIB record in encoded.
>
> Can you please tell me how these GRIB files were created?
>
> Thanks,
> John
>
>
> On Mon, May 16, 2016 at 1:40 AM, Gallai Irene via RT
> <met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
> >
> > John,
> > many thanks to you,
> > I'm happy to contribute to MET development, even if only
discovering
> > bugs...
> > I will be waiting for the fix.
> >
> > Thank you again
> > Irene
> >
> >
> > +------------------------------------+
> > WE ALL are environmentally responsible. By not printing this
email
> > you've helped save paper, ink, and millions of trees.
> > +------------------------------------+
> >
> > Irene Gallai
> >
> > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > SOC Direzione Tecnica Scientifica
> > CRMA - Centro Regionale di Modellistica Ambientale
> >
> > Via Cairoli 14
> > I-33057 Palmanova (UD) ITALY
> > tel pers.: +39 0432 1918044
> > tel centr.: +39 0432 1918120
> > e-mail : irene.gallai at arpa.fvg.it
> > +------------------------------------+
> >
> >
> > AVVISO DI RISERVATEZZA
> > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > suoi allegati. Se non siete i destinatari indicati nel messaggio,
o
> > responsabili per la sua consegna alla persona, o se avete ricevuto
> > il messaggio per errore, non dovete trascriverlo, copiarlo o
> > inviarlo a nessuno. In tal caso, dovete cancellare/distruggere il
> > messaggio ed i
> suoi allegati.
> > Grazie."
> >
> > CONFIDENTIALITY NOTICE
> > "Confidential information may be contained in this message or in
its
> > attachments. If you are not the addressee indicated in this
message,
> > or responsible for message delivering to that person, or if you
have
> > received this message in error, you may not transcribe, copy or
> > deliver this message to anyone. In that case, you should
> > delete/destroy this message and its attachments. Thank you .
> >
> > -----Messaggio originale-----
> > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Inviato: sabato 14 maggio 2016 01:18
> > A: Gallai Irene
> > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> >
> > Irene,
> >
> > Thank you for sending me those sample data files. I took a look
and
> > am able to reproduce the problem you're seeing. I ran your data
> > file both through MET and IDV and see that IDV plots it correctly
> > while MET
> does not.
> >
> > I strongly suspect that MET has a bug in the processing of the
GRIB
> > grid description section. I'll work on narrowing down the
problem,
> > fixing it, and posting a bugfix next week.
> >
> > I will email when a bugfix for this issue has been posted to the
MET
> > website.
> >
> > Thanks for letting us know about this problem!
> >
> > John
> >
> > On Fri, May 13, 2016 at 3:52 AM, Gallai Irene via RT
> > <met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
> > >
> > > Hello,
> > > As requested I posted 3 GRIB files (one for each of our domains)
> > > and
> > > 2 png images.
> > > The one named config_domain show the 3 domains.
> > > d02.png instead is been obtained plotting with plot_data_plane
> > > the GRIB file for the second domain obtained with UPP starting
> > > from the same wrf output files but leaving the coordinate in
> > > lambert conformal and not transforming it in regular lat lon: in
> > > this case domain localization is correct
> > >
> > > Concerning the version of MET I'm using is the 5.1 with bug-fix
> > > available at the date of 7 january 2016.
> > >
> > > Thanks
> > > Irene
> > > +------------------------------------+
> > > WE ALL are environmentally responsible. By not printing this
> > > email you've helped save paper, ink, and millions of trees.
> > > +------------------------------------+
> > >
> > > Irene Gallai
> > >
> > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > > SOC Direzione Tecnica Scientifica
> > > CRMA - Centro Regionale di Modellistica Ambientale
> > >
> > > Via Cairoli 14
> > > I-33057 Palmanova (UD) ITALY
> > > tel pers.: +39 0432 1918044
> > > tel centr.: +39 0432 1918120
> > > e-mail : irene.gallai at arpa.fvg.it
> > > +------------------------------------+
> > >
> > >
> > > AVVISO DI RISERVATEZZA
> > > "Informazioni riservate possono essere contenute nel messaggio o
> > > nei suoi allegati. Se non siete i destinatari indicati nel
> > > messaggio, o responsabili per la sua consegna alla persona, o se
> > > avete ricevuto il messaggio per errore, non dovete trascriverlo,
> > > copiarlo o inviarlo a nessuno. In tal caso, dovete
> > > cancellare/distruggere il messaggio ed i
> > suoi allegati.
> > > Grazie."
> > >
> > > CONFIDENTIALITY NOTICE
> > > "Confidential information may be contained in this message or in
> > > its attachments. If you are not the addressee indicated in this
> > > message, or responsible for message delivering to that person,
or
> > > if you have received this message in error, you may not
> > > transcribe, copy or deliver this message to anyone. In that
case,
> > > you should delete/destroy this message and its attachments.
Thank you .
> > >
> > >
> > > -----Messaggio originale-----
> > > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Inviato: giovedì 12 maggio 2016 22:40
> > > A: Gallai Irene
> > > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> > >
> > > Hello Irene,
> > >
> > > I see you're having trouble with the grid definition in a GRIB
file.
> > > Could you please post a sample GRIB file to our anonymous ftp
site
> > > following these instructions?
> > >
> > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > >
> > > I'd like to look at the grid definition in the file, plot is
using
> > > plot_data_plane (as you've done), and also plot it using
Unidata's
> > > IDV tool. Also, do you have a plot showing the "intended"
domain?
> > >
> > > Hopefully that info will help us track down the discrepancy
you're
> > seeing.
> > >
> > > Lastly, please tell me the version of MET you're using.
> > >
> > > Thanks,
> > > John
> > >
> > >
> > > On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
> > > <met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> > > > Transaction: Ticket created by irene.gallai at arpa.fvg.it
> > > > Queue: met_help
> > > > Subject: pointstat problem
> > > > Owner: Nobody
> > > > Requestors: irene.gallai at arpa.fvg.it
> > > > Status: new
> > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
> > > > >
> > > >
> > > >
> > > >
> > > > Dear All,
> > > >
> > > > I'm facing some troubles with pointstat routine, in detail I
> > > > couldn't obtain matching pairs using my data.
> > > > Using verbose debugging I see that Grid definition is not
> > > > correct because for example as lower left (ll) latitude is
> > > > indicated a value that in effects is my upper level value ,
as
> > > > you can see in the following screenshots that are taken from
> > > > pointstat log the first and a call to cdo libraries the
latter.
> > > >
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
----------------------------------------------------------------
> > > > --
> > > > --
> > > > --
> > > > ----------
> > > > DEBUG 2:
> > > > DEBUG 2: Reading data for TMP/Z2.
> > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range
> > > > match for VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> > > > "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> > > > f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1 GRIB
> > > > records matching VarInfo "TMP/Z2" in GRIB file
> > > > "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> > > > 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > > DEBUG 3: Use the matching forecast and observation grids.
> > > > DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90
lat_ll:
> > > > 46.920
> > > > lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025 DEBUG 2: For
> > > > TMP/Z2 found 1 forecast levels and 0 climatology levels.
> > > > DEBUG 2:
> > > > DEBUG 2:
> > > >
----------------------------------------------------------------
> > > > --
> > > > --
> > > > --
> > > > ----------
> > > > .....
> > > >
> > > > Whereas interrogating the same GRIB file with cdo I obtain
> > > > Grid coordinates :
> > > > 1 : lonlat > size : dim = 8100 nx = 90 ny =
90
> > > > lon : first = 11.888 last =
14.149
> inc
> > =
> > > > 0.0254044944 degrees_east
> > > > lat : first = 45.322 last =
46.92
> inc =
> > > > 0.0179550562 degrees_north
> > > >
> > > >
> > > > I'm missing something obvious? GRIB files come from UPP
> > > > (remapping the grid from lambert conformal to regular lon lat)
> > > >
> > > > It's more evident looking at the image produced with
> > > > plot_data_plane
> :
> > > > gray/blue colors should be the Alps in Italy Whereas political
> > > > boundaries show part of Denmark , Sweden ...
> > > > Any advice?
> > > >
> > > > Thanks
> > > > Irene
> > > >
> > > >
> > > > +------------------------------------+
> > > > WE ALL are environmentally responsible. By not printing this
> > > > email you've helped save paper, ink, and millions of trees.
> > > > +------------------------------------+
> > > >
> > > > Irene Gallai
> > > >
> > > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG
)
> > > > SOC Direzione Tecnica Scientifica
> > > > CRMA - Centro Regionale di Modellistica Ambientale
> > > >
> > > > Via Cairoli 14
> > > > I-33057 Palmanova (UD) ITALY
> > > > tel pers.: +39 0432 1918044
> > > > tel centr.: +39 0432 1918120
> > > > e-mail :
> > > > irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> > > > +------------------------------------+
> > > >
> > > >
> > > > AVVISO DI RISERVATEZZA
> > > > "Informazioni riservate possono essere contenute nel messaggio
o
> > > > nei suoi allegati. Se non siete i destinatari indicati nel
> > > > messaggio, o responsabili per la sua consegna alla persona, o
se
> > > > avete ricevuto il messaggio per errore, non dovete
trascriverlo,
> > > > copiarlo o inviarlo a nessuno. In tal caso, dovete
> > > > cancellare/distruggere il messaggio ed i
> > > suoi allegati.
> > > > Grazie."
> > > >
> > > > CONFIDENTIALITY NOTICE
> > > > "Confidential information may be contained in this message or
in
> > > > its attachments. If you are not the addressee indicated in
this
> > > > message, or responsible for message delivering to that person,
> > > > or if you have received this message in error, you may not
> > > > transcribe, copy or deliver this message to anyone. In that
> > > > case, you should delete/destroy this message and its
> > > > attachments. Thank
> you .
> > > >
> > > >
> > > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > > destinatari indicati nel messaggio, o responsabili per la sua
> > > > consegna alla persona, o se avete ricevuto il messaggio per
> > > > errore, siete pregati di non trascriverlo, copiarlo o inviarlo
> > > > ad alcuno. In tal caso vi invitiamo a cancellare il messaggio
ed
> > > > i suoi
> allegati.
> > Grazie.
> > > > CONFIDENTIALITY NOTICE Confidential information may be
contained
> > > > in this message or in its attachments. If you are not the
> > > > addressee indicated in this message, or responsible for
message
> > > > delivering to that person, or if you have received this
message
> > > > in error, you may not transcribe, copy or deliver this message
> > > > to anyone. In that case, you should delete this message and
its
> attachments. Thank you.
> > > >
> > > >
> > >
> > >
> > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > destinatari indicati nel messaggio, o responsabili per la sua
> > > consegna alla persona, o se avete ricevuto il messaggio per
> > > errore, siete pregati di non trascriverlo, copiarlo o inviarlo
ad
> > > alcuno. In tal caso vi invitiamo a cancellare il messaggio ed i
suoi allegati.
> Grazie.
> > > CONFIDENTIALITY NOTICE Confidential information may be contained
> > > in this message or in its attachments. If you are not the
> > > addressee indicated in this message, or responsible for message
> > > delivering to that person, or if you have received this message
in
> > > error, you may not transcribe, copy or deliver this message to
> > > anyone. In that case, you should delete this message and its
attachments. Thank you.
> > >
> > >
> > >
> >
> >
> > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > contenute nel messaggio o nei suoi allegati. Se non siete i
> > destinatari indicati nel messaggio, o responsabili per la sua
> > consegna alla persona, o se avete ricevuto il messaggio per
errore,
> > siete pregati di non trascriverlo, copiarlo o inviarlo ad alcuno.
In
> > tal caso vi invitiamo a cancellare il messaggio ed i suoi
allegati. Grazie.
> > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > this message or in its attachments. If you are not the addressee
> > indicated in this message, or responsible for message delivering
to
> > that person, or if you have received this message in error, you
may
> > not transcribe, copy or deliver this message to anyone. In that
> > case, you should delete this message and its attachments. Thank
you.
> >
> >
> >
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute
> nel messaggio o nei suoi allegati. Se non siete i destinatari
indicati
> nel messaggio, o responsabili per la sua consegna alla persona, o se
> avete ricevuto il messaggio per errore, siete pregati di non
> trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi
invitiamo
> a cancellare il messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
> this message or in its attachments. If you are not the addressee
> indicated in this message, or responsible for message delivering to
> that person, or if you have received this message in error, you may
> not transcribe, copy or deliver this message to anyone. In that
case,
> you should delete this message and its attachments. Thank you.
>
>
>
AVVISO DI RISERVATEZZA Informazioni riservate possono essere contenute
nel messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel messaggio, o responsabili per la sua consegna alla persona, o se
avete ricevuto il messaggio per errore, siete pregati di non
trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo
a cancellare il messaggio ed i suoi allegati. Grazie.
CONFIDENTIALITY NOTICE Confidential information may be contained in
this message or in its attachments. If you are not the addressee
indicated in this message, or responsible for message delivering to
that person, or if you have received this message in error, you may
not transcribe, copy or deliver this message to anyone. In that case,
you should delete this message and its attachments. Thank you.
------------------------------------------------
Subject: pointstat problem
From: John Halley Gotway
Time: Thu May 19 09:45:43 2016
Irene,
That's great news! I'm glad you were able to track down the problem.
Just
to confirm your findings, I grabbed your original GRIB file, ran it
through
coypgb using both scanning mode 0 and 64, and plotted the 2-m
temperature
output. The resulting images are attached.
As you can see, the original one is oriented correctly and so it the
scanning mode 64 copygb output. But scanning mode 0 is upside down,
as you
originally found.
Hopefully you'll now be able to get matched pairs from your runs of
Point-Stat.
I'll go ahead and resolve this helpdesk ticket.
Thanks,
John
On Thu, May 19, 2016 at 5:13 AM, Gallai Irene via RT
<met_help at ucar.edu>
wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
>
>
> Hello Jhon,
> I posted a sample GRIB output file from UPP before we run it
through
> coypgb, it's named CRMA_WRF-4-CRMA_d02_2016051900+000000h.GRIB.tmp.
> My colleague also handed to me the copygb command used and I think
that
> the problem is in there. If I'm not mistaken the last entry of the
kgds
> string should be the one that will set the scanning mode and as you
can
> see there is a 0 and not a 64 as you pointed out.
>
> Here is the command line
>
> copygb.exe -xg"255 0 135 155 35190 4563 136 49023 21585 127 89 0"
> CRMA_WRF-4-CRMA_d02_2016051900+000000h.GRIB.tmp
> CRMA_WRF-4-CRMA_d02_2016051900+000000h.GRIB
>
> Probably it remained some typo in the script from tests but we
didn't
> detected it because generally we look at outputs with GrADS after
post
> processing them with grib2ctl and gribmap and the plots never
showed
> reversed domains. It's possible that some "degribbing utilities"
don't
> consider the scanning mode flag? As you say it's suspicious that
IDV or
> GrADS are able to display the data correctly.
>
> In any case we would correct copygb command
>
> Many thanks again
> Irene
>
> +------------------------------------+
> WE ALL are environmentally responsible. By not printing this email
> you've helped save paper, ink, and millions of trees.
> +------------------------------------+
>
> Irene Gallai
>
> Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> SOC Direzione Tecnica Scientifica
> CRMA - Centro Regionale di Modellistica Ambientale
>
> Via Cairoli 14
> I-33057 Palmanova (UD) ITALY
> tel pers.: +39 0432 1918044
> tel centr.: +39 0432 1918120
> e-mail : irene.gallai at arpa.fvg.it
> +------------------------------------+
>
>
> AVVISO DI RISERVATEZZA
> "Informazioni riservate possono essere contenute nel messaggio o nei
suoi
> allegati. Se non siete i destinatari indicati nel messaggio, o
responsabili
> per la sua consegna alla persona, o se avete ricevuto il messaggio
per
> errore, non dovete trascriverlo, copiarlo o inviarlo a nessuno. In
tal
> caso, dovete cancellare/distruggere il messaggio ed i suoi allegati.
> Grazie."
>
> CONFIDENTIALITY NOTICE
> "Confidential information may be contained in this message or in its
> attachments. If you are not the addressee indicated in this message,
or
> responsible for message delivering to that person, or if you have
received
> this message in error, you may not transcribe, copy or deliver this
message
> to anyone. In that case, you should delete/destroy this message and
its
> attachments. Thank you .
>
> -----Messaggio originale-----
> Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> Inviato: mercoledì 18 maggio 2016 16:31
> A: Gallai Irene
> Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
>
> Irene,
>
> Hmm, that's interesting. Do you be any chance have a sample GRIB
output
> file from UPP before you run it through coypgb? Or even a sample
WRF
> NetCDF output file prior to running through UPP? If you post one of
those
> along with the UPP namelist or copygb command you use, I could take
a look
> at the process.
>
> The DTC is also responsible for supporting UPP, which includes
copygb. So
> I'd like to figure out where the problem is coming in.
>
> And it is still possible that MET isn't processing the grid
information
> correctly. It is suspicious that changing the scanning mode flag in
octet
> 28 to 64 fixes the problem. But it's also suspicious that IDV is
able to
> display the data correctly
>
> John
>
>
> On Wed, May 18, 2016 at 5:19 AM, Gallai Irene via RT
<met_help at ucar.edu>
> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
> >
> > Hi Jhon,
> >
> > We create GRIB files with UPP and then we change map projection
using
> > the copygb utility. Maybe the "copygb step" is the problem.
> > Unfortunately I'm not the one that make those passages so I don't
know
> > If and which parameters and arguments have been used .
> > I asked a colleague to check the procedure to see where and how
the
> > scanning mode flag in octet 28 is set.
> > I'll let you know when I have an answer For the moment thanks
again
> >
> > Irene
> >
> >
> > +------------------------------------+
> > WE ALL are environmentally responsible. By not printing this
email
> > you've helped save paper, ink, and millions of trees.
> > +------------------------------------+
> >
> > Irene Gallai
> >
> > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > SOC Direzione Tecnica Scientifica
> > CRMA - Centro Regionale di Modellistica Ambientale
> >
> > Via Cairoli 14
> > I-33057 Palmanova (UD) ITALY
> > tel pers.: +39 0432 1918044
> > tel centr.: +39 0432 1918120
> > e-mail : irene.gallai at arpa.fvg.it
> > +------------------------------------+
> >
> >
> > AVVISO DI RISERVATEZZA
> > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > suoi allegati. Se non siete i destinatari indicati nel messaggio,
o
> > responsabili per la sua consegna alla persona, o se avete ricevuto
il
> > messaggio per errore, non dovete trascriverlo, copiarlo o inviarlo
a
> > nessuno. In tal caso, dovete cancellare/distruggere il messaggio
ed i
> suoi allegati.
> > Grazie."
> >
> > CONFIDENTIALITY NOTICE
> > "Confidential information may be contained in this message or in
its
> > attachments. If you are not the addressee indicated in this
message,
> > or responsible for message delivering to that person, or if you
have
> > received this message in error, you may not transcribe, copy or
> > deliver this message to anyone. In that case, you should
> > delete/destroy this message and its attachments. Thank you .
> >
> >
> > -----Messaggio originale-----
> > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > Inviato: lunedì 16 maggio 2016 21:23
> > A: Gallai Irene
> > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> >
> > Irene,
> >
> > After looking very closely, I suspect that the problem lies in how
> > your GRIB data is encoded. Specifically, I believe the scanning
mode
> > flag in octet 28 is set wrong. It is set to 0 but should be set
to 64.
> >
> > Here's some output from wgrib for your sample data:
> > wgrib -GDS10 CRMA_WRF-4-CRMA_d02_2016050500+000000h_TMP_SFC.GRIB
> >
> >
1:0:d=16050500:TMP:kpds5=11:kpds6=1:kpds7=0:TR=0:P1=0:P2=0:TimeU=1:sfc
> > :anl:NAve=0:GDS10=
> > 0 0 32 0 255 0 0 135 0 155 0 137 118 0 17 211 136 0 191 127 0 84
81 0
> > 127 0
> > 89 0 0 0 0 0
> >
> > The 28th entry in the GDS list is set to 0. We use the following
> > table to interpret that value of 0:
> > http://www.nco.ncep.noaa.gov/pmb/docs/on388/table8.html
> >
> > Since all of the bits are set to 0, that means the data is listed
in
> > the
> > +x and -y directions. However, I believe your data is actually
listed
> > +in
> > the
> > +x and +y directions.
> >
> > By way of comparison, I ran both your lat/lon data and some GFS
> > lat/lon data through an NCL utility named "ncl_convert2nc". Then
I
> > used ncview to see the result. Note that both your data and GFS
have
> > the scan flag set to 0. However, the attached ncview images
> > (gfs_globe.png and d02_italy.png) show that your data is "upside"
down.
> >
> > When I manually edit the MET code to reset the scanning mode flag
to
> > 64, the resulting images from MET look great.
> >
> > So at this point, I suspect there isn't a bug in MET, but a
problem in
> > how your GRIB record in encoded.
> >
> > Can you please tell me how these GRIB files were created?
> >
> > Thanks,
> > John
> >
> >
> > On Mon, May 16, 2016 at 1:40 AM, Gallai Irene via RT
> > <met_help at ucar.edu>
> > wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349 >
> > >
> > > John,
> > > many thanks to you,
> > > I'm happy to contribute to MET development, even if only
discovering
> > > bugs...
> > > I will be waiting for the fix.
> > >
> > > Thank you again
> > > Irene
> > >
> > >
> > > +------------------------------------+
> > > WE ALL are environmentally responsible. By not printing this
email
> > > you've helped save paper, ink, and millions of trees.
> > > +------------------------------------+
> > >
> > > Irene Gallai
> > >
> > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG )
> > > SOC Direzione Tecnica Scientifica
> > > CRMA - Centro Regionale di Modellistica Ambientale
> > >
> > > Via Cairoli 14
> > > I-33057 Palmanova (UD) ITALY
> > > tel pers.: +39 0432 1918044
> > > tel centr.: +39 0432 1918120
> > > e-mail : irene.gallai at arpa.fvg.it
> > > +------------------------------------+
> > >
> > >
> > > AVVISO DI RISERVATEZZA
> > > "Informazioni riservate possono essere contenute nel messaggio o
nei
> > > suoi allegati. Se non siete i destinatari indicati nel
messaggio, o
> > > responsabili per la sua consegna alla persona, o se avete
ricevuto
> > > il messaggio per errore, non dovete trascriverlo, copiarlo o
> > > inviarlo a nessuno. In tal caso, dovete cancellare/distruggere
il
> > > messaggio ed i
> > suoi allegati.
> > > Grazie."
> > >
> > > CONFIDENTIALITY NOTICE
> > > "Confidential information may be contained in this message or in
its
> > > attachments. If you are not the addressee indicated in this
message,
> > > or responsible for message delivering to that person, or if you
have
> > > received this message in error, you may not transcribe, copy or
> > > deliver this message to anyone. In that case, you should
> > > delete/destroy this message and its attachments. Thank you .
> > >
> > > -----Messaggio originale-----
> > > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > Inviato: sabato 14 maggio 2016 01:18
> > > A: Gallai Irene
> > > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> > >
> > > Irene,
> > >
> > > Thank you for sending me those sample data files. I took a look
and
> > > am able to reproduce the problem you're seeing. I ran your data
> > > file both through MET and IDV and see that IDV plots it
correctly
> > > while MET
> > does not.
> > >
> > > I strongly suspect that MET has a bug in the processing of the
GRIB
> > > grid description section. I'll work on narrowing down the
problem,
> > > fixing it, and posting a bugfix next week.
> > >
> > > I will email when a bugfix for this issue has been posted to the
MET
> > > website.
> > >
> > > Thanks for letting us know about this problem!
> > >
> > > John
> > >
> > > On Fri, May 13, 2016 at 3:52 AM, Gallai Irene via RT
> > > <met_help at ucar.edu>
> > > wrote:
> > >
> > > >
> > > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
>
> > > >
> > > > Hello,
> > > > As requested I posted 3 GRIB files (one for each of our
domains)
> > > > and
> > > > 2 png images.
> > > > The one named config_domain show the 3 domains.
> > > > d02.png instead is been obtained plotting with
plot_data_plane
> > > > the GRIB file for the second domain obtained with UPP starting
> > > > from the same wrf output files but leaving the coordinate in
> > > > lambert conformal and not transforming it in regular lat lon:
in
> > > > this case domain localization is correct
> > > >
> > > > Concerning the version of MET I'm using is the 5.1 with bug-
fix
> > > > available at the date of 7 january 2016.
> > > >
> > > > Thanks
> > > > Irene
> > > > +------------------------------------+
> > > > WE ALL are environmentally responsible. By not printing this
> > > > email you've helped save paper, ink, and millions of trees.
> > > > +------------------------------------+
> > > >
> > > > Irene Gallai
> > > >
> > > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA FVG
)
> > > > SOC Direzione Tecnica Scientifica
> > > > CRMA - Centro Regionale di Modellistica Ambientale
> > > >
> > > > Via Cairoli 14
> > > > I-33057 Palmanova (UD) ITALY
> > > > tel pers.: +39 0432 1918044
> > > > tel centr.: +39 0432 1918120
> > > > e-mail : irene.gallai at arpa.fvg.it
> > > > +------------------------------------+
> > > >
> > > >
> > > > AVVISO DI RISERVATEZZA
> > > > "Informazioni riservate possono essere contenute nel messaggio
o
> > > > nei suoi allegati. Se non siete i destinatari indicati nel
> > > > messaggio, o responsabili per la sua consegna alla persona, o
se
> > > > avete ricevuto il messaggio per errore, non dovete
trascriverlo,
> > > > copiarlo o inviarlo a nessuno. In tal caso, dovete
> > > > cancellare/distruggere il messaggio ed i
> > > suoi allegati.
> > > > Grazie."
> > > >
> > > > CONFIDENTIALITY NOTICE
> > > > "Confidential information may be contained in this message or
in
> > > > its attachments. If you are not the addressee indicated in
this
> > > > message, or responsible for message delivering to that person,
or
> > > > if you have received this message in error, you may not
> > > > transcribe, copy or deliver this message to anyone. In that
case,
> > > > you should delete/destroy this message and its attachments.
Thank
> you .
> > > >
> > > >
> > > > -----Messaggio originale-----
> > > > Da: John Halley Gotway via RT [mailto:met_help at ucar.edu]
> > > > Inviato: giovedì 12 maggio 2016 22:40
> > > > A: Gallai Irene
> > > > Oggetto: Re: [rt.rap.ucar.edu #76349] pointstat problem
> > > >
> > > > Hello Irene,
> > > >
> > > > I see you're having trouble with the grid definition in a GRIB
file.
> > > > Could you please post a sample GRIB file to our anonymous ftp
site
> > > > following these instructions?
> > > >
> > > > http://www.dtcenter.org/met/users/support/met_help.php#ftp
> > > >
> > > > I'd like to look at the grid definition in the file, plot is
using
> > > > plot_data_plane (as you've done), and also plot it using
Unidata's
> > > > IDV tool. Also, do you have a plot showing the "intended"
domain?
> > > >
> > > > Hopefully that info will help us track down the discrepancy
you're
> > > seeing.
> > > >
> > > > Lastly, please tell me the version of MET you're using.
> > > >
> > > > Thanks,
> > > > John
> > > >
> > > >
> > > > On Thu, May 12, 2016 at 6:23 AM, Gallai Irene via RT
> > > > <met_help at ucar.edu>
> > > > wrote:
> > > >
> > > > >
> > > > > Thu May 12 06:23:59 2016: Request 76349 was acted upon.
> > > > > Transaction: Ticket created by irene.gallai at arpa.fvg.it
> > > > > Queue: met_help
> > > > > Subject: pointstat problem
> > > > > Owner: Nobody
> > > > > Requestors: irene.gallai at arpa.fvg.it
> > > > > Status: new
> > > > > Ticket <URL:
> > > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=76349
> > > > > >
> > > > >
> > > > >
> > > > >
> > > > > Dear All,
> > > > >
> > > > > I'm facing some troubles with pointstat routine, in detail I
> > > > > couldn't obtain matching pairs using my data.
> > > > > Using verbose debugging I see that Grid definition is not
> > > > > correct because for example as lower left (ll) latitude is
> > > > > indicated a value that in effects is my upper level value ,
as
> > > > > you can see in the following screenshots that are taken from
> > > > > pointstat log the first and a call to cdo libraries the
latter.
> > > > >
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
----------------------------------------------------------------
> > > > > --
> > > > > --
> > > > > --
> > > > > ----------
> > > > > DEBUG 2:
> > > > > DEBUG 2: Reading data for TMP/Z2.
> > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found range
> > > > > match for VarInfo "TMP/Z2" in GRIB record 181 of GRIB file
> > > > > "/lustre/arpa/gallaii/scratch/wrf4oper//wr
> > > > > f_verif//d03//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > > > DEBUG 3: MetGrib1DataFile::data_plane_array() -> Found 1
GRIB
> > > > > records matching VarInfo "TMP/Z2" in GRIB file
> > > > > "/lustre/arpa/gallaii/scratch/wrf4oper//wrf_verif//d0
> > > > > 3//CRMA_WRF-4-CRMA_d03_2016050300+000000h.GRIB".
> > > > > DEBUG 3: Use the matching forecast and observation grids.
> > > > > DEBUG 3: Grid Definition: Projection: Lat/Lon Nx: 90 Ny: 90
lat_ll:
> > > > > 46.920
> > > > > lon_ll: -11.888 delta_lat: 0.017 delta_lon: 0.025 DEBUG 2:
For
> > > > > TMP/Z2 found 1 forecast levels and 0 climatology levels.
> > > > > DEBUG 2:
> > > > > DEBUG 2:
> > > > >
----------------------------------------------------------------
> > > > > --
> > > > > --
> > > > > --
> > > > > ----------
> > > > > .....
> > > > >
> > > > > Whereas interrogating the same GRIB file with cdo I obtain
> > > > > Grid coordinates :
> > > > > 1 : lonlat > size : dim = 8100 nx = 90 ny
= 90
> > > > > lon : first = 11.888 last =
14.149
> > inc
> > > =
> > > > > 0.0254044944 degrees_east
> > > > > lat : first = 45.322 last =
46.92
> > inc =
> > > > > 0.0179550562 degrees_north
> > > > >
> > > > >
> > > > > I'm missing something obvious? GRIB files come from UPP
> > > > > (remapping the grid from lambert conformal to regular lon
lat)
> > > > >
> > > > > It's more evident looking at the image produced with
> > > > > plot_data_plane
> > :
> > > > > gray/blue colors should be the Alps in Italy Whereas
political
> > > > > boundaries show part of Denmark , Sweden ...
> > > > > Any advice?
> > > > >
> > > > > Thanks
> > > > > Irene
> > > > >
> > > > >
> > > > > +------------------------------------+
> > > > > WE ALL are environmentally responsible. By not printing
this
> > > > > email you've helped save paper, ink, and millions of trees.
> > > > > +------------------------------------+
> > > > >
> > > > > Irene Gallai
> > > > >
> > > > > Agenzia Regionale per la Protezione dell'Ambiente (ARPA
FVG )
> > > > > SOC Direzione Tecnica Scientifica
> > > > > CRMA - Centro Regionale di Modellistica Ambientale
> > > > >
> > > > > Via Cairoli 14
> > > > > I-33057 Palmanova (UD) ITALY
> > > > > tel pers.: +39 0432 1918044
> > > > > tel centr.: +39 0432 1918120
> > > > > e-mail :
> > > > > irene.gallai at arpa.fvg.it<mailto:irene.gallai at osmer.fvg.it>
> > > > > +------------------------------------+
> > > > >
> > > > >
> > > > > AVVISO DI RISERVATEZZA
> > > > > "Informazioni riservate possono essere contenute nel
messaggio o
> > > > > nei suoi allegati. Se non siete i destinatari indicati nel
> > > > > messaggio, o responsabili per la sua consegna alla persona,
o se
> > > > > avete ricevuto il messaggio per errore, non dovete
trascriverlo,
> > > > > copiarlo o inviarlo a nessuno. In tal caso, dovete
> > > > > cancellare/distruggere il messaggio ed i
> > > > suoi allegati.
> > > > > Grazie."
> > > > >
> > > > > CONFIDENTIALITY NOTICE
> > > > > "Confidential information may be contained in this message
or in
> > > > > its attachments. If you are not the addressee indicated in
this
> > > > > message, or responsible for message delivering to that
person,
> > > > > or if you have received this message in error, you may not
> > > > > transcribe, copy or deliver this message to anyone. In that
> > > > > case, you should delete/destroy this message and its
> > > > > attachments. Thank
> > you .
> > > > >
> > > > >
> > > > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > > > destinatari indicati nel messaggio, o responsabili per la
sua
> > > > > consegna alla persona, o se avete ricevuto il messaggio per
> > > > > errore, siete pregati di non trascriverlo, copiarlo o
inviarlo
> > > > > ad alcuno. In tal caso vi invitiamo a cancellare il
messaggio ed
> > > > > i suoi
> > allegati.
> > > Grazie.
> > > > > CONFIDENTIALITY NOTICE Confidential information may be
contained
> > > > > in this message or in its attachments. If you are not the
> > > > > addressee indicated in this message, or responsible for
message
> > > > > delivering to that person, or if you have received this
message
> > > > > in error, you may not transcribe, copy or deliver this
message
> > > > > to anyone. In that case, you should delete this message and
its
> > attachments. Thank you.
> > > > >
> > > > >
> > > >
> > > >
> > > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > > destinatari indicati nel messaggio, o responsabili per la sua
> > > > consegna alla persona, o se avete ricevuto il messaggio per
> > > > errore, siete pregati di non trascriverlo, copiarlo o inviarlo
ad
> > > > alcuno. In tal caso vi invitiamo a cancellare il messaggio ed
i suoi
> allegati.
> > Grazie.
> > > > CONFIDENTIALITY NOTICE Confidential information may be
contained
> > > > in this message or in its attachments. If you are not the
> > > > addressee indicated in this message, or responsible for
message
> > > > delivering to that person, or if you have received this
message in
> > > > error, you may not transcribe, copy or deliver this message to
> > > > anyone. In that case, you should delete this message and its
> attachments. Thank you.
> > > >
> > > >
> > > >
> > >
> > >
> > > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
> > > contenute nel messaggio o nei suoi allegati. Se non siete i
> > > destinatari indicati nel messaggio, o responsabili per la sua
> > > consegna alla persona, o se avete ricevuto il messaggio per
errore,
> > > siete pregati di non trascriverlo, copiarlo o inviarlo ad
alcuno. In
> > > tal caso vi invitiamo a cancellare il messaggio ed i suoi
allegati.
> Grazie.
> > > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > > this message or in its attachments. If you are not the addressee
> > > indicated in this message, or responsible for message delivering
to
> > > that person, or if you have received this message in error, you
may
> > > not transcribe, copy or deliver this message to anyone. In that
> > > case, you should delete this message and its attachments. Thank
you.
> > >
> > >
> > >
> >
> >
> > AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute
> > nel messaggio o nei suoi allegati. Se non siete i destinatari
indicati
> > nel messaggio, o responsabili per la sua consegna alla persona, o
se
> > avete ricevuto il messaggio per errore, siete pregati di non
> > trascriverlo, copiarlo o inviarlo ad alcuno. In tal caso vi
invitiamo
> > a cancellare il messaggio ed i suoi allegati. Grazie.
> > CONFIDENTIALITY NOTICE Confidential information may be contained
in
> > this message or in its attachments. If you are not the addressee
> > indicated in this message, or responsible for message delivering
to
> > that person, or if you have received this message in error, you
may
> > not transcribe, copy or deliver this message to anyone. In that
case,
> > you should delete this message and its attachments. Thank you.
> >
> >
> >
>
>
> AVVISO DI RISERVATEZZA Informazioni riservate possono essere
contenute nel
> messaggio o nei suoi allegati. Se non siete i destinatari indicati
nel
> messaggio, o responsabili per la sua consegna alla persona, o se
avete
> ricevuto il messaggio per errore, siete pregati di non trascriverlo,
> copiarlo o inviarlo ad alcuno. In tal caso vi invitiamo a cancellare
il
> messaggio ed i suoi allegati. Grazie.
> CONFIDENTIALITY NOTICE Confidential information may be contained in
this
> message or in its attachments. If you are not the addressee
indicated in
> this message, or responsible for message delivering to that person,
or if
> you have received this message in error, you may not transcribe,
copy or
> deliver this message to anyone. In that case, you should delete this
> message and its attachments. Thank you.
>
>
>
------------------------------------------------
More information about the Met_help
mailing list