[Met_help] [rt.rap.ucar.edu #91162] History for Subtracting Rotated Grids in PCP Combine
John Halley Gotway via RT
met_help at ucar.edu
Mon Aug 5 14:12:53 MDT 2019
----------------------------------------------------------------
Initial Request
----------------------------------------------------------------
John, the problem I sent you was from a program I used to subtract two grid parameters in the same file. Because it was a rotated grid, I was placing the wrong grid header data in the NetCDF file and then running point_stat. I was doing this because of a problem in PCP Combine. When I do a subtraction of two grid parameters on PCP combine for a rotated grid I get the following error:
/h/MET/bin/pcp_combine -subtract /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2 'name="CEIL"; level="L0"; ' /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2 'name="MTERH"; level="L0"; ' /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-a8ab4b677452 -v 3
DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file: /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file: /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
DEBUG 2: Since the "-field" command line option was not used, parsing the command line arguments a list of files, each followed by a configuration string.
DEBUG 1: Reading input file: /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input file: /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
DEBUG 3:
DEBUG 3: Multiple GRIB2 table entries match lookup criteria (parm_name = CEIL):
DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6, index_c = 13
DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14, grib2_cntr = 57, grib2_ltab = 1, index_b = 6, index_c = 11
DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6, index_c = 13
DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 6, index_c = 13
DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14, grib2_cntr = 57, grib2_ltab = 0, index_b = 6, index_c = 13
DEBUG 3: Using the first match found: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6, index_c = 13
DEBUG 3:
DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for 'CEIL/L0' in GRIB2 record 1 field 1 of GRIB2 file '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
DEBUG 1: Reading input file: /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input file: /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
DEBUG 3:
DEBUG 3: Multiple GRIB2 table entries match lookup criteria (parm_name = MTERH):
DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0, index_c = 7
DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14, grib2_cntr = 57, grib2_ltab = 0, index_b = 0, index_c = 7
DEBUG 3: Using the first match found: parm_name: MTERH, index_a = 2, grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0, index_c = 7
DEBUG 3:
DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
ERROR :
ERROR : do_sub_command() -> the input fields must be on the same grid.
ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll: -7.542 rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036 true_lat_south_pole: -34.515 true_lon_south_pole: -22.179 aux_rotation: 0.000
ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll: -7.542 rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036 true_lat_south_pole: -34.515 true_lon_south_pole: -22.179 aux_rotation: 0.000
ERROR :
This was the real problem I need to get past.
Thanks
Bob
----------------------------------------------------------------
Complete Ticket History
----------------------------------------------------------------
Subject: Subtracting Rotated Grids in PCP Combine
From: John Halley Gotway
Time: Fri Jul 19 14:18:05 2019
Bob,
Are you able to send those two sample rotated lat/lon files? If I can
replicate the problem here, I should be able to debug it.
Thanks,
John
On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
> Queue: met_help
> Subject: Subtracting Rotated Grids in PCP Combine
> Owner: Nobody
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
>
>
> John, the problem I sent you was from a program I used to subtract
two
> grid parameters in the same file. Because it was a rotated grid, I
was
> placing the wrong grid header data in the NetCDF file and then
running
> point_stat. I was doing this because of a problem in PCP Combine.
When I
> do a subtraction of two grid parameters on PCP combine for a rotated
grid I
> get the following error:
>
> /h/MET/bin/pcp_combine -subtract
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> 'name="CEIL"; level="L0"; '
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> 'name="MTERH"; level="L0"; '
> /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b677452 -v 3
> DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> DEBUG 2: Since the "-field" command line option was not used,
parsing the
> command line arguments a list of files, each followed by a
configuration
> string.
> DEBUG 1: Reading input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 3:
> DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name =
> CEIL):
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 11
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 6, index_c = 13
> DEBUG 3: Using the first match found: parm_name: CEIL, index_a =
0,
> grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
index_c = 13
> DEBUG 3:
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
'CEIL/L0'
> in GRIB2 record 1 field 1 of GRIB2 file
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> DEBUG 1: Reading input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 3:
> DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name =
> MTERH):
> DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 0, index_c = 7
> DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
grib2_cntr =
> 57, grib2_ltab = 0, index_b = 0, index_c = 7
> DEBUG 3: Using the first match found: parm_name: MTERH, index_a =
2,
> grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0, index_c
= 7
> DEBUG 3:
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> ERROR :
> ERROR : do_sub_command() -> the input fields must be on the same
grid.
> ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
-7.542
> rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> 0.000
> ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
-7.542
> rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> 0.000
> ERROR :
>
> This was the real problem I need to get past.
>
> Thanks
> Bob
>
>
------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting Rotated Grids in PCP Combine
From: robert.craig.2 at us.af.mil
Time: Fri Jul 19 14:23:00 2019
The grid was sent in AMRDEC with the last problem report. Let me know
if I need to resend.
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Friday, July 19, 2019 3:18 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
Rotated Grids in PCP Combine
Bob,
Are you able to send those two sample rotated lat/lon files? If I can
replicate the problem here, I should be able to debug it.
Thanks,
John
On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
> Queue: met_help
> Subject: Subtracting Rotated Grids in PCP Combine
> Owner: Nobody
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162
> >
>
>
> John, the problem I sent you was from a program I used to subtract
two
> grid parameters in the same file. Because it was a rotated grid, I
> was placing the wrong grid header data in the NetCDF file and then
> running point_stat. I was doing this because of a problem in PCP
> Combine. When I do a subtraction of two grid parameters on PCP
> combine for a rotated grid I get the following error:
>
> /h/MET/bin/pcp_combine -subtract
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> 'name="CEIL"; level="L0"; '
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> 'name="MTERH"; level="L0"; '
> /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b677452
> -v 3 DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> DEBUG 2: Since the "-field" command line option was not used,
parsing
> the command line arguments a list of files, each followed by a
> configuration string.
> DEBUG 1: Reading input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 3:
> DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name
> =
> CEIL):
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 11
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 6, index_c = 13
> DEBUG 3: Using the first match found: parm_name: CEIL, index_a =
0,
> grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
index_c
> = 13 DEBUG 3:
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
'CEIL/L0'
> in GRIB2 record 1 field 1 of GRIB2 file
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> DEBUG 1: Reading input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 3:
> DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name
> =
> MTERH):
> DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 0, index_c = 7
> DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
grib2_cntr =
> 57, grib2_ltab = 0, index_b = 0, index_c = 7
> DEBUG 3: Using the first match found: parm_name: MTERH, index_a =
2,
> grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0, index_c
=
> 7 DEBUG 3:
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> ERROR :
> ERROR : do_sub_command() -> the input fields must be on the same
grid.
> ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> -7.542
> rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> 0.000
> ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> -7.542
> rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> 0.000
> ERROR :
>
> This was the real problem I need to get past.
>
> Thanks
> Bob
>
>
------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting Rotated Grids in PCP Combine
From: robert.craig.2 at us.af.mil
Time: Fri Jul 19 14:24:38 2019
There is only one rotated lat lon file, I am subtracting to fields in
the file.
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Friday, July 19, 2019 3:18 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
Rotated Grids in PCP Combine
Bob,
Are you able to send those two sample rotated lat/lon files? If I can
replicate the problem here, I should be able to debug it.
Thanks,
John
On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> Transaction: Ticket created by robert.craig.2 at us.af.mil
> Queue: met_help
> Subject: Subtracting Rotated Grids in PCP Combine
> Owner: Nobody
> Requestors: robert.craig.2 at us.af.mil
> Status: new
> Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162
> >
>
>
> John, the problem I sent you was from a program I used to subtract
two
> grid parameters in the same file. Because it was a rotated grid, I
> was placing the wrong grid header data in the NetCDF file and then
> running point_stat. I was doing this because of a problem in PCP
> Combine. When I do a subtraction of two grid parameters on PCP
> combine for a rotated grid I get the following error:
>
> /h/MET/bin/pcp_combine -subtract
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> 'name="CEIL"; level="L0"; '
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> 'name="MTERH"; level="L0"; '
> /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b677452
> -v 3 DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> DEBUG 2: Since the "-field" command line option was not used,
parsing
> the command line arguments a list of files, each followed by a
> configuration string.
> DEBUG 1: Reading input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 3:
> DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name
> =
> CEIL):
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 11
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 57,
> grib2_ltab = 1, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 6, index_c = 13
> DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14, grib2_cntr
= 57,
> grib2_ltab = 0, index_b = 6, index_c = 13
> DEBUG 3: Using the first match found: parm_name: CEIL, index_a =
0,
> grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
index_c
> = 13 DEBUG 3:
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
'CEIL/L0'
> in GRIB2 record 1 field 1 of GRIB2 file
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> DEBUG 1: Reading input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input file:
> /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> DEBUG 3:
> DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name
> =
> MTERH):
> DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1, grib2_cntr
= 0,
> grib2_ltab = 0, index_b = 0, index_c = 7
> DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
grib2_cntr =
> 57, grib2_ltab = 0, index_b = 0, index_c = 7
> DEBUG 3: Using the first match found: parm_name: MTERH, index_a =
2,
> grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0, index_c
=
> 7 DEBUG 3:
> DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> ERROR :
> ERROR : do_sub_command() -> the input fields must be on the same
grid.
> ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> -7.542
> rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> 0.000
> ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> -7.542
> rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> 0.000
> ERROR :
>
> This was the real problem I need to get past.
>
> Thanks
> Bob
>
>
------------------------------------------------
Subject: Subtracting Rotated Grids in PCP Combine
From: John Halley Gotway
Time: Fri Jul 19 14:56:33 2019
Bob,
Got it thanks. Unfortunately, this is a bug in met-8.1. I documented
it
in this GitHub issue:
https://github.com/NCAR/MET/issues/1170
I was able to fix this by updating these two files:
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.h
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.cc
These updates will be included in the met-8.1.2 bugfix release.
However,
you're welcome to make whatever changes you'd like. If you drop in
the new
versions of these files and recompile, it should fix the problem.
Thanks,
John
On Fri, Jul 19, 2019 at 2:24 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
>
> There is only one rotated lat lon file, I am subtracting to fields
in the
> file.
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Friday, July 19, 2019 3:18 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
> Rotated Grids in PCP Combine
>
> Bob,
>
> Are you able to send those two sample rotated lat/lon files? If I
can
> replicate the problem here, I should be able to debug it.
>
> Thanks,
> John
>
> On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> > Transaction: Ticket created by robert.craig.2 at us.af.mil
> > Queue: met_help
> > Subject: Subtracting Rotated Grids in PCP Combine
> > Owner: Nobody
> > Requestors: robert.craig.2 at us.af.mil
> > Status: new
> > Ticket <URL:
https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162
> > >
> >
> >
> > John, the problem I sent you was from a program I used to subtract
two
> > grid parameters in the same file. Because it was a rotated grid,
I
> > was placing the wrong grid header data in the NetCDF file and then
> > running point_stat. I was doing this because of a problem in PCP
> > Combine. When I do a subtraction of two grid parameters on PCP
> > combine for a rotated grid I get the following error:
> >
> > /h/MET/bin/pcp_combine -subtract
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > 'name="CEIL"; level="L0"; '
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > 'name="MTERH"; level="L0"; '
> > /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b677452
> > -v 3 DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> > /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> > DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> > /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> > DEBUG 2: Since the "-field" command line option was not used,
parsing
> > the command line arguments a list of files, each followed by a
> > configuration string.
> > DEBUG 1: Reading input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 3:
> > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name
> > =
> > CEIL):
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 57,
> > grib2_ltab = 1, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> 57,
> > grib2_ltab = 1, index_b = 6, index_c = 11
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 57,
> > grib2_ltab = 1, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 0,
> > grib2_ltab = 0, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> 57,
> > grib2_ltab = 0, index_b = 6, index_c = 13
> > DEBUG 3: Using the first match found: parm_name: CEIL, index_a =
0,
> > grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
index_c
> > = 13 DEBUG 3:
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'CEIL/L0'
> > in GRIB2 record 1 field 1 of GRIB2 file
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > DEBUG 1: Reading input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input
file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 3:
> > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
(parm_name
> > =
> > MTERH):
> > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1,
grib2_cntr = 0,
> > grib2_ltab = 0, index_b = 0, index_c = 7
> > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
grib2_cntr =
> > 57, grib2_ltab = 0, index_b = 0, index_c = 7
> > DEBUG 3: Using the first match found: parm_name: MTERH, index_a
= 2,
> > grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0,
index_c =
> > 7 DEBUG 3:
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > ERROR :
> > ERROR : do_sub_command() -> the input fields must be on the same
grid.
> > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > -7.542
> > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > 0.000
> > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > -7.542
> > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > 0.000
> > ERROR :
> >
> > This was the real problem I need to get past.
> >
> > Thanks
> > Bob
> >
> >
>
>
>
>
------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting Rotated Grids in PCP Combine
From: robert.craig.2 at us.af.mil
Time: Fri Jul 19 15:05:47 2019
Okay, I will try that.
Thanks
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Friday, July 19, 2019 3:57 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
Rotated Grids in PCP Combine
Bob,
Got it thanks. Unfortunately, this is a bug in met-8.1. I documented
it in this GitHub issue:
https://github.com/NCAR/MET/issues/1170
I was able to fix this by updating these two files:
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.h
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.cc
These updates will be included in the met-8.1.2 bugfix release.
However, you're welcome to make whatever changes you'd like. If you
drop in the new versions of these files and recompile, it should fix
the problem.
Thanks,
John
On Fri, Jul 19, 2019 at 2:24 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
>
> There is only one rotated lat lon file, I am subtracting to fields
in
> the file.
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Friday, July 19, 2019 3:18 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> <robert.craig.2 at us.af.mil>
> Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
> Rotated Grids in PCP Combine
>
> Bob,
>
> Are you able to send those two sample rotated lat/lon files? If I
can
> replicate the problem here, I should be able to debug it.
>
> Thanks,
> John
>
> On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> > Transaction: Ticket created by robert.craig.2 at us.af.mil
> > Queue: met_help
> > Subject: Subtracting Rotated Grids in PCP Combine
> > Owner: Nobody
> > Requestors: robert.craig.2 at us.af.mil
> > Status: new
> > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162
> > >
> >
> >
> > John, the problem I sent you was from a program I used to subtract
> > two grid parameters in the same file. Because it was a rotated
> > grid, I was placing the wrong grid header data in the NetCDF file
> > and then running point_stat. I was doing this because of a
problem
> > in PCP Combine. When I do a subtraction of two grid parameters on
> > PCP combine for a rotated grid I get the following error:
> >
> > /h/MET/bin/pcp_combine -subtract
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > 'name="CEIL"; level="L0"; '
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > 'name="MTERH"; level="L0"; '
> > /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b67745
> > 2 -v 3 DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> > /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> > DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> > /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> > DEBUG 2: Since the "-field" command line option was not used,
> > parsing the command line arguments a list of files, each followed
by
> > a configuration string.
> > DEBUG 1: Reading input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 3:
> > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > (parm_name =
> > CEIL):
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 57,
> > grib2_ltab = 1, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> 57,
> > grib2_ltab = 1, index_b = 6, index_c = 11
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 57,
> > grib2_ltab = 1, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 0,
> > grib2_ltab = 0, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> 57,
> > grib2_ltab = 0, index_b = 6, index_c = 13
> > DEBUG 3: Using the first match found: parm_name: CEIL, index_a =
0,
> > grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
> > index_c = 13 DEBUG 3:
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'CEIL/L0'
> > in GRIB2 record 1 field 1 of GRIB2 file
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > DEBUG 1: Reading input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input
file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 3:
> > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > (parm_name =
> > MTERH):
> > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1,
grib2_cntr = 0,
> > grib2_ltab = 0, index_b = 0, index_c = 7
> > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
grib2_cntr =
> > 57, grib2_ltab = 0, index_b = 0, index_c = 7
> > DEBUG 3: Using the first match found: parm_name: MTERH, index_a
= 2,
> > grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0,
index_c
> > =
> > 7 DEBUG 3:
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > ERROR :
> > ERROR : do_sub_command() -> the input fields must be on the same
grid.
> > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > -7.542
> > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > 0.000
> > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > -7.542
> > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > 0.000
> > ERROR :
> >
> > This was the real problem I need to get past.
> >
> > Thanks
> > Bob
> >
> >
>
>
>
>
------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting Rotated Grids in PCP Combine
From: robert.craig.2 at us.af.mil
Time: Mon Jul 22 11:43:39 2019
John, I tried to recompile with your changes but get the following
error:
/vx_gnomon -I../../../src/libcode/vx_grid
-I../../../src/libcode/vx_gsl_prob -I../../../src/libcode/vx_nav
-I../../../src/libcode/vx_nc_obs -I../../../src/libcode/vx_nc_util
-I../../../src/libcode/vx_pb_util -I../../../src/libcode/vx_plot_util
-I../../../src/libcode/vx_ps -I../../../src/libcode/vx_pxm
-I../../../src/libcode/vx_render -I../../../src/libcode/vx_regrid
-I../../../src/libcode/vx_shapedata -I../../../src/libcode/vx_solar
-I../../../src/libcode/vx_statistics
-I../../../src/libcode/vx_stat_out -I../../../src/libcode/vx_summary
-I../../../src/libcode/vx_time_series
-I../../../src/libcode/vx_tc_util -I/h/WXQC/met-
8.1_fixes/external_libs/netcdf-c-4.4.1.1/include -I/h/WXQC/met-
8.1_fixes/external_libs/hdf5-1.8.18/include -I/h/WXQC/met-
8.1_fixes/external_libs/g2clib-1.4.0 -I/h/WXQC/met-
8.1_fixes/external_libs/gsl-1.15/include -I/h/WXQC/met-
8.1_fixes/external_libs/HDF4.2.12/include -I/h/WXQC/met-
8.1_fixes/external_libs/hdfeos/include -DBLOCK4 -DWITH_GRIB2
-DMET_BASE="\"/h/WXQC/met-8.1_fixes/share/met\"" -g -MT libvx_grid_a-
grid_base.o -MD -MP -MF .deps/libvx_grid_a-grid_base.Tpo -c -o
libvx_grid_a-grid_base.o `test -f 'grid_base.cc' || echo
'./'`grid_base.cc
grid_base.cc: In function ‘bool is_eq(const RotatedLatLonData*, const
RotatedLatLonData*)’:
grid_base.cc:1290: error: ‘aux_rotation’ was not declared in this
scope
make[3]: *** [libvx_grid_a-grid_base.o] Error 1
Was there any other files that needed to be changed?
Thanks
Bob
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Friday, July 19, 2019 3:57 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
Rotated Grids in PCP Combine
Bob,
Got it thanks. Unfortunately, this is a bug in met-8.1. I documented
it in this GitHub issue:
https://github.com/NCAR/MET/issues/1170
I was able to fix this by updating these two files:
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.h
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.cc
These updates will be included in the met-8.1.2 bugfix release.
However, you're welcome to make whatever changes you'd like. If you
drop in the new versions of these files and recompile, it should fix
the problem.
Thanks,
John
On Fri, Jul 19, 2019 at 2:24 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
>
> There is only one rotated lat lon file, I am subtracting to fields
in
> the file.
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Friday, July 19, 2019 3:18 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> <robert.craig.2 at us.af.mil>
> Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
> Rotated Grids in PCP Combine
>
> Bob,
>
> Are you able to send those two sample rotated lat/lon files? If I
can
> replicate the problem here, I should be able to debug it.
>
> Thanks,
> John
>
> On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> > Transaction: Ticket created by robert.craig.2 at us.af.mil
> > Queue: met_help
> > Subject: Subtracting Rotated Grids in PCP Combine
> > Owner: Nobody
> > Requestors: robert.craig.2 at us.af.mil
> > Status: new
> > Ticket <URL:
> > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162
> > >
> >
> >
> > John, the problem I sent you was from a program I used to subtract
> > two grid parameters in the same file. Because it was a rotated
> > grid, I was placing the wrong grid header data in the NetCDF file
> > and then running point_stat. I was doing this because of a
problem
> > in PCP Combine. When I do a subtraction of two grid parameters on
> > PCP combine for a rotated grid I get the following error:
> >
> > /h/MET/bin/pcp_combine -subtract
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > 'name="CEIL"; level="L0"; '
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > 'name="MTERH"; level="L0"; '
> > /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b67745
> > 2 -v 3 DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> > /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> > DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> > /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> > DEBUG 2: Since the "-field" command line option was not used,
> > parsing the command line arguments a list of files, each followed
by
> > a configuration string.
> > DEBUG 1: Reading input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 3:
> > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > (parm_name =
> > CEIL):
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 57,
> > grib2_ltab = 1, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> 57,
> > grib2_ltab = 1, index_b = 6, index_c = 11
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 57,
> > grib2_ltab = 1, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr = 0,
> > grib2_ltab = 0, index_b = 6, index_c = 13
> > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> 57,
> > grib2_ltab = 0, index_b = 6, index_c = 13
> > DEBUG 3: Using the first match found: parm_name: CEIL, index_a =
0,
> > grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
> > index_c = 13 DEBUG 3:
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> 'CEIL/L0'
> > in GRIB2 record 1 field 1 of GRIB2 file
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > DEBUG 1: Reading input file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input
file:
> > /home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_
> > GR
> > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > DEBUG 3:
> > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > (parm_name =
> > MTERH):
> > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1,
grib2_cntr = 0,
> > grib2_ltab = 0, index_b = 0, index_c = 7
> > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
grib2_cntr =
> > 57, grib2_ltab = 0, index_b = 0, index_c = 7
> > DEBUG 3: Using the first match found: parm_name: MTERH, index_a
= 2,
> > grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0,
index_c
> > =
> > 7 DEBUG 3:
> > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > ERROR :
> > ERROR : do_sub_command() -> the input fields must be on the same
grid.
> > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > -7.542
> > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > 0.000
> > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > -7.542
> > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > 0.000
> > ERROR :
> >
> > This was the real problem I need to get past.
> >
> > Thanks
> > Bob
> >
> >
>
>
>
>
------------------------------------------------
Subject: Subtracting Rotated Grids in PCP Combine
From: John Halley Gotway
Time: Mon Jul 22 12:06:29 2019
Bob,
You should update grid_base.h and grid_base.cc. I just tested this
myself
to make sure.
(1) Untarred the met-8.1 release tarball.
(2) Copied the updated version of grid_base.h:
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.h
(3) Copied the updated version of grid_base.cc:
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.cc
(4) Ran configure.
(5) Ran make install.
And it all compiled well.
Thanks,
John
On Mon, Jul 22, 2019 at 11:43 AM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
>
> John, I tried to recompile with your changes but get the following
error:
>
> /vx_gnomon -I../../../src/libcode/vx_grid
> -I../../../src/libcode/vx_gsl_prob -I../../../src/libcode/vx_nav
> -I../../../src/libcode/vx_nc_obs -I../../../src/libcode/vx_nc_util
> -I../../../src/libcode/vx_pb_util
-I../../../src/libcode/vx_plot_util
> -I../../../src/libcode/vx_ps -I../../../src/libcode/vx_pxm
> -I../../../src/libcode/vx_render -I../../../src/libcode/vx_regrid
> -I../../../src/libcode/vx_shapedata -I../../../src/libcode/vx_solar
> -I../../../src/libcode/vx_statistics
-I../../../src/libcode/vx_stat_out
> -I../../../src/libcode/vx_summary
-I../../../src/libcode/vx_time_series
> -I../../../src/libcode/vx_tc_util
> -I/h/WXQC/met-8.1_fixes/external_libs/netcdf-c-4.4.1.1/include
> -I/h/WXQC/met-8.1_fixes/external_libs/hdf5-1.8.18/include
> -I/h/WXQC/met-8.1_fixes/external_libs/g2clib-1.4.0
> -I/h/WXQC/met-8.1_fixes/external_libs/gsl-1.15/include
> -I/h/WXQC/met-8.1_fixes/external_libs/HDF4.2.12/include
> -I/h/WXQC/met-8.1_fixes/external_libs/hdfeos/include -DBLOCK4
-DWITH_GRIB2
> -DMET_BASE="\"/h!
> /WXQC/met-8.1_fixes/share/met\"" -g -MT libvx_grid_a-grid_base.o
-MD -MP
> -MF .deps/libvx_grid_a-grid_base.Tpo -c -o libvx_grid_a-grid_base.o
`test
> -f 'grid_base.cc' || echo './'`grid_base.cc
> grid_base.cc: In function ‘bool is_eq(const RotatedLatLonData*,
const
> RotatedLatLonData*)’:
> grid_base.cc:1290: error: ‘aux_rotation’ was not declared in this
scope
> make[3]: *** [libvx_grid_a-grid_base.o] Error 1
>
> Was there any other files that needed to be changed?
>
> Thanks
> Bob
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Friday, July 19, 2019 3:57 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162]
Subtracting
> Rotated Grids in PCP Combine
>
> Bob,
>
> Got it thanks. Unfortunately, this is a bug in met-8.1. I
documented it
> in this GitHub issue:
> https://github.com/NCAR/MET/issues/1170
>
> I was able to fix this by updating these two files:
>
>
>
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.h
>
>
>
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.cc
>
> These updates will be included in the met-8.1.2 bugfix release.
However,
> you're welcome to make whatever changes you'd like. If you drop in
the new
> versions of these files and recompile, it should fix the problem.
>
> Thanks,
> John
>
>
> On Fri, Jul 19, 2019 at 2:24 PM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
> >
> > There is only one rotated lat lon file, I am subtracting to fields
in
> > the file.
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT <met_help at ucar.edu>
> > Sent: Friday, July 19, 2019 3:18 PM
> > To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> > <robert.craig.2 at us.af.mil>
> > Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
> > Rotated Grids in PCP Combine
> >
> > Bob,
> >
> > Are you able to send those two sample rotated lat/lon files? If I
can
> > replicate the problem here, I should be able to debug it.
> >
> > Thanks,
> > John
> >
> > On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> > > Transaction: Ticket created by robert.craig.2 at us.af.mil
> > > Queue: met_help
> > > Subject: Subtracting Rotated Grids in PCP Combine
> > > Owner: Nobody
> > > Requestors: robert.craig.2 at us.af.mil
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162
> > > >
> > >
> > >
> > > John, the problem I sent you was from a program I used to
subtract
> > > two grid parameters in the same file. Because it was a rotated
> > > grid, I was placing the wrong grid header data in the NetCDF
file
> > > and then running point_stat. I was doing this because of a
problem
> > > in PCP Combine. When I do a subtraction of two grid parameters
on
> > > PCP combine for a rotated grid I get the following error:
> > >
> > > /h/MET/bin/pcp_combine -subtract
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > 'name="CEIL"; level="L0"; '
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > 'name="MTERH"; level="L0"; '
> > > /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b67745
> > > 2 -v 3 DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> > > /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> > > DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> > > /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> > > DEBUG 2: Since the "-field" command line option was not used,
> > > parsing the command line arguments a list of files, each
followed by
> > > a configuration string.
> > > DEBUG 1: Reading input file:
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input
file:
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > DEBUG 3:
> > > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > > (parm_name =
> > > CEIL):
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> 57,
> > > grib2_ltab = 1, index_b = 6, index_c = 13
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> > 57,
> > > grib2_ltab = 1, index_b = 6, index_c = 11
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> 57,
> > > grib2_ltab = 1, index_b = 6, index_c = 13
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> 0,
> > > grib2_ltab = 0, index_b = 6, index_c = 13
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> > 57,
> > > grib2_ltab = 0, index_b = 6, index_c = 13
> > > DEBUG 3: Using the first match found: parm_name: CEIL, index_a
= 0,
> > > grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
> > > index_c = 13 DEBUG 3:
> > > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'CEIL/L0'
> > > in GRIB2 record 1 field 1 of GRIB2 file
> > >
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > > DEBUG 1: Reading input file:
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input
file:
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-RD_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > DEBUG 3:
> > > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > > (parm_name =
> > > MTERH):
> > > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1,
grib2_cntr =
> 0,
> > > grib2_ltab = 0, index_b = 0, index_c = 7
> > > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
grib2_cntr =
> > > 57, grib2_ltab = 0, index_b = 0, index_c = 7
> > > DEBUG 3: Using the first match found: parm_name: MTERH,
index_a = 2,
> > > grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0,
index_c
> > > =
> > > 7 DEBUG 3:
> > > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > > 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> > >
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > > ERROR :
> > > ERROR : do_sub_command() -> the input fields must be on the
same grid.
> > > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > > -7.542
> > > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > > 0.000
> > > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > > -7.542
> > > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > > 0.000
> > > ERROR :
> > >
> > > This was the real problem I need to get past.
> > >
> > > Thanks
> > > Bob
> > >
> > >
> >
> >
> >
> >
>
>
>
>
------------------------------------------------
Subject: RE: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting Rotated Grids in PCP Combine
From: robert.craig.2 at us.af.mil
Time: Mon Jul 22 12:24:11 2019
John, could you email me those two files. GITHUB is not letting me
get them.
Thanks
Bob
-----Original Message-----
From: John Halley Gotway via RT <met_help at ucar.edu>
Sent: Monday, July 22, 2019 1:06 PM
To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
Rotated Grids in PCP Combine
Bob,
You should update grid_base.h and grid_base.cc. I just tested this
myself to make sure.
(1) Untarred the met-8.1 release tarball.
(2) Copied the updated version of grid_base.h:
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.h
(3) Copied the updated version of grid_base.cc:
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.cc
(4) Ran configure.
(5) Ran make install.
And it all compiled well.
Thanks,
John
On Mon, Jul 22, 2019 at 11:43 AM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
>
> John, I tried to recompile with your changes but get the following
error:
>
> /vx_gnomon -I../../../src/libcode/vx_grid
> -I../../../src/libcode/vx_gsl_prob -I../../../src/libcode/vx_nav
> -I../../../src/libcode/vx_nc_obs -I../../../src/libcode/vx_nc_util
> -I../../../src/libcode/vx_pb_util
-I../../../src/libcode/vx_plot_util
> -I../../../src/libcode/vx_ps -I../../../src/libcode/vx_pxm
> -I../../../src/libcode/vx_render -I../../../src/libcode/vx_regrid
> -I../../../src/libcode/vx_shapedata -I../../../src/libcode/vx_solar
> -I../../../src/libcode/vx_statistics
> -I../../../src/libcode/vx_stat_out
> -I../../../src/libcode/vx_summary
> -I../../../src/libcode/vx_time_series
> -I../../../src/libcode/vx_tc_util
> -I/h/WXQC/met-8.1_fixes/external_libs/netcdf-c-4.4.1.1/include
> -I/h/WXQC/met-8.1_fixes/external_libs/hdf5-1.8.18/include
> -I/h/WXQC/met-8.1_fixes/external_libs/g2clib-1.4.0
> -I/h/WXQC/met-8.1_fixes/external_libs/gsl-1.15/include
> -I/h/WXQC/met-8.1_fixes/external_libs/HDF4.2.12/include
> -I/h/WXQC/met-8.1_fixes/external_libs/hdfeos/include -DBLOCK4
> -DWITH_GRIB2 -DMET_BASE="\"/h!
> /WXQC/met-8.1_fixes/share/met\"" -g -MT libvx_grid_a-grid_base.o
-MD
> -MP -MF .deps/libvx_grid_a-grid_base.Tpo -c -o
> libvx_grid_a-grid_base.o `test -f 'grid_base.cc' || echo
> './'`grid_base.cc
> grid_base.cc: In function ‘bool is_eq(const RotatedLatLonData*,
const
> RotatedLatLonData*)’:
> grid_base.cc:1290: error: ‘aux_rotation’ was not declared in this
> scope
> make[3]: *** [libvx_grid_a-grid_base.o] Error 1
>
> Was there any other files that needed to be changed?
>
> Thanks
> Bob
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Friday, July 19, 2019 3:57 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> <robert.craig.2 at us.af.mil>
> Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162]
Subtracting
> Rotated Grids in PCP Combine
>
> Bob,
>
> Got it thanks. Unfortunately, this is a bug in met-8.1. I
documented
> it in this GitHub issue:
> https://github.com/NCAR/MET/issues/1170
>
> I was able to fix this by updating these two files:
>
>
>
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/g
> rid_base.h
>
>
>
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/g
> rid_base.cc
>
> These updates will be included in the met-8.1.2 bugfix release.
> However, you're welcome to make whatever changes you'd like. If you
> drop in the new versions of these files and recompile, it should fix
the problem.
>
> Thanks,
> John
>
>
> On Fri, Jul 19, 2019 at 2:24 PM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
> >
> > There is only one rotated lat lon file, I am subtracting to fields
> > in the file.
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT <met_help at ucar.edu>
> > Sent: Friday, July 19, 2019 3:18 PM
> > To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> > <robert.craig.2 at us.af.mil>
> > Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162] Subtracting
> > Rotated Grids in PCP Combine
> >
> > Bob,
> >
> > Are you able to send those two sample rotated lat/lon files? If I
> > can replicate the problem here, I should be able to debug it.
> >
> > Thanks,
> > John
> >
> > On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> > > Transaction: Ticket created by robert.craig.2 at us.af.mil
> > > Queue: met_help
> > > Subject: Subtracting Rotated Grids in PCP Combine
> > > Owner: Nobody
> > > Requestors: robert.craig.2 at us.af.mil
> > > Status: new
> > > Ticket <URL:
> > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162
> > > >
> > >
> > >
> > > John, the problem I sent you was from a program I used to
subtract
> > > two grid parameters in the same file. Because it was a rotated
> > > grid, I was placing the wrong grid header data in the NetCDF
file
> > > and then running point_stat. I was doing this because of a
> > > problem in PCP Combine. When I do a subtraction of two grid
> > > parameters on PCP combine for a rotated grid I get the following
error:
> > >
> > > /h/MET/bin/pcp_combine -subtract
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > D_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > 'name="CEIL"; level="L0"; '
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > D_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > 'name="MTERH"; level="L0"; '
> > > /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b677
> > > 45
> > > 2 -v 3 DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES file:
> > > /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> > > DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> > > /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> > > DEBUG 2: Since the "-field" command line option was not used,
> > > parsing the command line arguments a list of files, each
followed
> > > by a configuration string.
> > > DEBUG 1: Reading input file:
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > D_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input
file:
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > D_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > DEBUG 3:
> > > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > > (parm_name =
> > > CEIL):
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> 57,
> > > grib2_ltab = 1, index_b = 6, index_c = 13
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> > 57,
> > > grib2_ltab = 1, index_b = 6, index_c = 11
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> 57,
> > > grib2_ltab = 1, index_b = 6, index_c = 13
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> 0,
> > > grib2_ltab = 0, index_b = 6, index_c = 13
> > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr =
> > 57,
> > > grib2_ltab = 0, index_b = 6, index_c = 13
> > > DEBUG 3: Using the first match found: parm_name: CEIL, index_a
= 0,
> > > grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
> > > index_c = 13 DEBUG 3:
> > > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > 'CEIL/L0'
> > > in GRIB2 record 1 field 1 of GRIB2 file
> > >
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > > DEBUG 1: Reading input file:
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > D_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input
file:
> > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > D_
> > > GR
> > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > DEBUG 3:
> > > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > > (parm_name =
> > > MTERH):
> > > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1,
grib2_cntr =
> 0,
> > > grib2_ltab = 0, index_b = 0, index_c = 7
> > > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
grib2_cntr =
> > > 57, grib2_ltab = 0, index_b = 0, index_c = 7
> > > DEBUG 3: Using the first match found: parm_name: MTERH,
index_a = 2,
> > > grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0,
> > > index_c =
> > > 7 DEBUG 3:
> > > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match for
> > > 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> > >
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > > ERROR :
> > > ERROR : do_sub_command() -> the input fields must be on the
same grid.
> > > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > > -7.542
> > > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > > 0.000
> > > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420 rot_lat_ll:
> > > -7.542
> > > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
aux_rotation:
> > > 0.000
> > > ERROR :
> > >
> > > This was the real problem I need to get past.
> > >
> > > Thanks
> > > Bob
> > >
> > >
> >
> >
> >
> >
>
>
>
>
------------------------------------------------
Subject: Subtracting Rotated Grids in PCP Combine
From: John Halley Gotway
Time: Mon Jul 22 13:03:33 2019
They are attached.
Thanks,
John
On Mon, Jul 22, 2019 at 12:24 PM robert.craig.2 at us.af.mil via RT <
met_help at ucar.edu> wrote:
>
> <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
>
> John, could you email me those two files. GITHUB is not letting me
get
> them.
>
> Thanks
> Bob
>
> -----Original Message-----
> From: John Halley Gotway via RT <met_help at ucar.edu>
> Sent: Monday, July 22, 2019 1:06 PM
> To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
<robert.craig.2 at us.af.mil>
> Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162]
Subtracting
> Rotated Grids in PCP Combine
>
> Bob,
>
> You should update grid_base.h and grid_base.cc. I just tested this
myself
> to make sure.
>
> (1) Untarred the met-8.1 release tarball.
> (2) Copied the updated version of grid_base.h:
>
>
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.h
> (3) Copied the updated version of grid_base.cc:
>
>
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/grid_base.cc
> (4) Ran configure.
> (5) Ran make install.
>
> And it all compiled well.
>
> Thanks,
> John
>
> On Mon, Jul 22, 2019 at 11:43 AM robert.craig.2 at us.af.mil via RT <
> met_help at ucar.edu> wrote:
>
> >
> > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
> >
> > John, I tried to recompile with your changes but get the following
error:
> >
> > /vx_gnomon -I../../../src/libcode/vx_grid
> > -I../../../src/libcode/vx_gsl_prob -I../../../src/libcode/vx_nav
> > -I../../../src/libcode/vx_nc_obs -I../../../src/libcode/vx_nc_util
> > -I../../../src/libcode/vx_pb_util
-I../../../src/libcode/vx_plot_util
> > -I../../../src/libcode/vx_ps -I../../../src/libcode/vx_pxm
> > -I../../../src/libcode/vx_render -I../../../src/libcode/vx_regrid
> > -I../../../src/libcode/vx_shapedata
-I../../../src/libcode/vx_solar
> > -I../../../src/libcode/vx_statistics
> > -I../../../src/libcode/vx_stat_out
> > -I../../../src/libcode/vx_summary
> > -I../../../src/libcode/vx_time_series
> > -I../../../src/libcode/vx_tc_util
> > -I/h/WXQC/met-8.1_fixes/external_libs/netcdf-c-4.4.1.1/include
> > -I/h/WXQC/met-8.1_fixes/external_libs/hdf5-1.8.18/include
> > -I/h/WXQC/met-8.1_fixes/external_libs/g2clib-1.4.0
> > -I/h/WXQC/met-8.1_fixes/external_libs/gsl-1.15/include
> > -I/h/WXQC/met-8.1_fixes/external_libs/HDF4.2.12/include
> > -I/h/WXQC/met-8.1_fixes/external_libs/hdfeos/include -DBLOCK4
> > -DWITH_GRIB2 -DMET_BASE="\"/h!
> > /WXQC/met-8.1_fixes/share/met\"" -g -MT libvx_grid_a-grid_base.o
-MD
> > -MP -MF .deps/libvx_grid_a-grid_base.Tpo -c -o
> > libvx_grid_a-grid_base.o `test -f 'grid_base.cc' || echo
> > './'`grid_base.cc
> > grid_base.cc: In function ‘bool is_eq(const RotatedLatLonData*,
const
> > RotatedLatLonData*)’:
> > grid_base.cc:1290: error: ‘aux_rotation’ was not declared in this
> > scope
> > make[3]: *** [libvx_grid_a-grid_base.o] Error 1
> >
> > Was there any other files that needed to be changed?
> >
> > Thanks
> > Bob
> >
> > -----Original Message-----
> > From: John Halley Gotway via RT <met_help at ucar.edu>
> > Sent: Friday, July 19, 2019 3:57 PM
> > To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> > <robert.craig.2 at us.af.mil>
> > Subject: Re: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162]
Subtracting
> > Rotated Grids in PCP Combine
> >
> > Bob,
> >
> > Got it thanks. Unfortunately, this is a bug in met-8.1. I
documented
> > it in this GitHub issue:
> > https://github.com/NCAR/MET/issues/1170
> >
> > I was able to fix this by updating these two files:
> >
> >
> >
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/g
> > rid_base.h
> >
> >
> >
https://github.com/NCAR/MET/blob/master_v8.1/met/src/libcode/vx_grid/g
> > rid_base.cc
> >
> > These updates will be included in the met-8.1.2 bugfix release.
> > However, you're welcome to make whatever changes you'd like. If
you
> > drop in the new versions of these files and recompile, it should
fix the
> problem.
> >
> > Thanks,
> > John
> >
> >
> > On Fri, Jul 19, 2019 at 2:24 PM robert.craig.2 at us.af.mil via RT <
> > met_help at ucar.edu> wrote:
> >
> > >
> > > <URL: https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162 >
> > >
> > > There is only one rotated lat lon file, I am subtracting to
fields
> > > in the file.
> > >
> > > -----Original Message-----
> > > From: John Halley Gotway via RT <met_help at ucar.edu>
> > > Sent: Friday, July 19, 2019 3:18 PM
> > > To: CRAIG, ROBERT J GS-12 USAF ACC 16 WS/WXN
> > > <robert.craig.2 at us.af.mil>
> > > Subject: [Non-DoD Source] Re: [rt.rap.ucar.edu #91162]
Subtracting
> > > Rotated Grids in PCP Combine
> > >
> > > Bob,
> > >
> > > Are you able to send those two sample rotated lat/lon files? If
I
> > > can replicate the problem here, I should be able to debug it.
> > >
> > > Thanks,
> > > John
> > >
> > > On Fri, Jul 19, 2019 at 2:05 PM robert.craig.2 at us.af.mil via RT
<
> > > met_help at ucar.edu> wrote:
> > >
> > > >
> > > > Fri Jul 19 14:05:15 2019: Request 91162 was acted upon.
> > > > Transaction: Ticket created by robert.craig.2 at us.af.mil
> > > > Queue: met_help
> > > > Subject: Subtracting Rotated Grids in PCP Combine
> > > > Owner: Nobody
> > > > Requestors: robert.craig.2 at us.af.mil
> > > > Status: new
> > > > Ticket <URL:
> > > > https://rt.rap.ucar.edu/rt/Ticket/Display.html?id=91162
> > > > >
> > > >
> > > >
> > > > John, the problem I sent you was from a program I used to
subtract
> > > > two grid parameters in the same file. Because it was a
rotated
> > > > grid, I was placing the wrong grid header data in the NetCDF
file
> > > > and then running point_stat. I was doing this because of a
> > > > problem in PCP Combine. When I do a subtraction of two grid
> > > > parameters on PCP combine for a rotated grid I get the
following
> error:
> > > >
> > > > /h/MET/bin/pcp_combine -subtract
> > > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > > D_
> > > > GR
> > > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > > 'name="CEIL"; level="L0"; '
> > > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > > D_
> > > > GR
> > > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > > 'name="MTERH"; level="L0"; '
> > > > /h/data/global/WXQC/data/met/tmp/2dc3122b-4f20-42c1-9d3c-
a8ab4b677
> > > > 45
> > > > 2 -v 3 DEBUG 1: Reading user-defined grib1 MET_GRIB_TABLES
file:
> > > > /h/data/global/WXQC/data/met/config/grib2/grib1_af_131_57.txt
> > > > DEBUG 1: Reading user-defined grib2 MET_GRIB_TABLES file:
> > > > /h/data/global/WXQC/data/met/config/grib2/grib2_geps.txt
> > > > DEBUG 2: Since the "-field" command line option was not used,
> > > > parsing the command line arguments a list of files, each
followed
> > > > by a configuration string.
> > > > DEBUG 1: Reading input file:
> > > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > > D_
> > > > GR
> > > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > > DEBUG 1: Reading data (name="CEIL"; level="L0"; ) from input
file:
> > > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > > D_
> > > > GR
> > > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > > DEBUG 3:
> > > > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > > > (parm_name =
> > > > CEIL):
> > > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> > 57,
> > > > grib2_ltab = 1, index_b = 6, index_c = 13
> > > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr
> =
> > > 57,
> > > > grib2_ltab = 1, index_b = 6, index_c = 11
> > > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> > 57,
> > > > grib2_ltab = 1, index_b = 6, index_c = 13
> > > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 1,
grib2_cntr =
> > 0,
> > > > grib2_ltab = 0, index_b = 6, index_c = 13
> > > > DEBUG 3: parm_name: CEIL, index_a = 0, grib2_mtab = 14,
grib2_cntr
> =
> > > 57,
> > > > grib2_ltab = 0, index_b = 6, index_c = 13
> > > > DEBUG 3: Using the first match found: parm_name: CEIL,
index_a = 0,
> > > > grib2_mtab = 1, grib2_cntr = 57, grib2_ltab = 1, index_b = 6,
> > > > index_c = 13 DEBUG 3:
> > > > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match
for
> > > 'CEIL/L0'
> > > > in GRIB2 record 1 field 1 of GRIB2 file
> > > >
> > >
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > > > DEBUG 1: Reading input file:
> > > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > > D_
> > > > GR
> > > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > > DEBUG 1: Reading data (name="MTERH"; level="L0"; ) from input
file:
> > > >
/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-R
> > > > D_
> > > > GR
> > > > .C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2
> > > > DEBUG 3:
> > > > DEBUG 3: Multiple GRIB2 table entries match lookup criteria
> > > > (parm_name =
> > > > MTERH):
> > > > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 1,
grib2_cntr
> =
> > 0,
> > > > grib2_ltab = 0, index_b = 0, index_c = 7
> > > > DEBUG 3: parm_name: MTERH, index_a = 2, grib2_mtab = 14,
> grib2_cntr =
> > > > 57, grib2_ltab = 0, index_b = 0, index_c = 7
> > > > DEBUG 3: Using the first match found: parm_name: MTERH,
index_a =
> 2,
> > > > grib2_mtab = 1, grib2_cntr = 0, grib2_ltab = 0, index_b = 0,
> > > > index_c =
> > > > 7 DEBUG 3:
> > > > DEBUG 3: MetGrib2DataFile::data_plane() -> Found exact match
for
> > > > 'MTERH/L0' in GRIB2 record 157 field 1 of GRIB2 file
> > > >
> > >
> >
> '/home/specmdev/scratch/test/baltic4/PS.557WW_SC.U_DI.C_GP.GALWEM-
RD_GR.C4KM_AR.BALTIC_DD.20180102_CY.00_FH.000_DF.GR2'
> > > > ERROR :
> > > > ERROR : do_sub_command() -> the input fields must be on the
same
> grid.
> > > > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420
rot_lat_ll:
> > > > -7.542
> > > > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > > > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
> aux_rotation:
> > > > 0.000
> > > > ERROR : Projection: Rotated Lat/Lon Nx: 408 Ny: 420
rot_lat_ll:
> > > > -7.542
> > > > rot_lon_ll: 8.326 delta_rot_lat: 0.036 delta_rot_lon: 0.036
> > > > true_lat_south_pole: -34.515 true_lon_south_pole: -22.179
> aux_rotation:
> > > > 0.000
> > > > ERROR :
> > > >
> > > > This was the real problem I need to get past.
> > > >
> > > > Thanks
> > > > Bob
> > > >
> > > >
> > >
> > >
> > >
> > >
> >
> >
> >
> >
>
>
>
>
------------------------------------------------
More information about the Met_help
mailing list