From d5d5c7c1b1c6fd88ac3f09525b40ca9d0f2e6a35 Mon Sep 17 00:00:00 2001
From: Howard Soh <hsoh@seneca.rap.ucar.edu>
Date: Thu, 2 Feb 2023 12:17:54 -0700
Subject: [PATCH] #2218 Initialize polar stereographics parametres for
 ellipsoial earth

---
 src/libcode/vx_data2d_grib/grib_utils.cc             |  6 ++++++
 src/libcode/vx_data2d_grib2/data2d_grib2.cc          |  5 +++++
 src/libcode/vx_data2d_nc_met/get_met_grid.cc         | 12 ++++++++++++
 src/libcode/vx_data2d_nc_pinterp/get_pinterp_grid.cc |  6 ++++++
 4 files changed, 29 insertions(+)

diff --git a/src/libcode/vx_data2d_grib/grib_utils.cc b/src/libcode/vx_data2d_grib/grib_utils.cc
index c1ef70f775..71b3ee8dc2 100644
--- a/src/libcode/vx_data2d_grib/grib_utils.cc
+++ b/src/libcode/vx_data2d_grib/grib_utils.cc
@@ -572,6 +572,12 @@ data.nx = char2_to_int(gds.nx);
    // Number of points in the y-direction
 data.ny = char2_to_int(gds.ny);
 
+data.eccentricity = 0.;
+data.false_east = 0.;
+data.false_north = 0.;
+data.scale_factor = 1.0;
+data.dy_km = data.d_km;
+
 data.dump();
 
    //
diff --git a/src/libcode/vx_data2d_grib2/data2d_grib2.cc b/src/libcode/vx_data2d_grib2/data2d_grib2.cc
index 42582b781e..45d931b9c2 100644
--- a/src/libcode/vx_data2d_grib2/data2d_grib2.cc
+++ b/src/libcode/vx_data2d_grib2/data2d_grib2.cc
@@ -1172,6 +1172,11 @@ void MetGrib2DataFile::read_grib2_grid( gribfield *gfld) {
       data.r_km         = r_km;
       data.nx           = gfld->igdtmpl[7];
       data.ny           = gfld->igdtmpl[8];
+      data.eccentricity = 0.;
+      data.false_east = 0.;
+      data.false_north = 0.;
+      data.scale_factor = 1.0;
+      data.dy_km = data.d_km;
 
       //  check for dx != dy
       if( !is_eq((double)gfld->igdtmpl[14] / 1000000.0,
diff --git a/src/libcode/vx_data2d_nc_met/get_met_grid.cc b/src/libcode/vx_data2d_nc_met/get_met_grid.cc
index 6016ab4dd2..67ddad0fa8 100644
--- a/src/libcode/vx_data2d_nc_met/get_met_grid.cc
+++ b/src/libcode/vx_data2d_nc_met/get_met_grid.cc
@@ -447,6 +447,12 @@ get_global_att(ncfile, string("nx"), data.nx);
    // Number of points in the y-direction
 get_global_att(ncfile, string("ny"), data.ny);
 
+data.eccentricity = 0.;
+data.false_east = 0.;
+data.false_north = 0.;
+data.scale_factor = 1.0;
+data.dy_km = data.d_km;
+
 data.dump();
 
    //
@@ -643,6 +649,12 @@ get_global_att(ncfile, string("nx"), data.nx);
    // Number of points in the y-direction
 get_global_att(ncfile, string("ny"), data.ny);
 
+data.eccentricity = 0.;
+data.false_east = 0.;
+data.false_north = 0.;
+data.scale_factor = 1.0;
+data.dy_km = data.d_km;
+
 data.dump();
 
    //
diff --git a/src/libcode/vx_data2d_nc_pinterp/get_pinterp_grid.cc b/src/libcode/vx_data2d_nc_pinterp/get_pinterp_grid.cc
index b5a78b109a..d93f0d3701 100644
--- a/src/libcode/vx_data2d_nc_pinterp/get_pinterp_grid.cc
+++ b/src/libcode/vx_data2d_nc_pinterp/get_pinterp_grid.cc
@@ -185,6 +185,12 @@ data.d_km *= 0.001;
 
 data.r_km = default_grib_radius_km;
 
+data.eccentricity = 0.;
+data.false_east = 0.;
+data.false_north = 0.;
+data.scale_factor = 1.0;
+data.dy_km = data.d_km;
+
    //
    //  done
    //