Bugfix #2115 develop Rotated LatLon (#2117)

dtcenter · Apr 1, 2022 · f75e27c · f75e27c
1 parent f1e7e9f
commit f75e27c
Show file tree

Hide file tree

Showing 2 changed files with 18 additions and 2 deletions.
diff --git a/met/src/libcode/vx_data2d_nccf/nccf_file.cc b/met/src/libcode/vx_data2d_nccf/nccf_file.cc
@@ -2430,10 +2430,12 @@ void NcCfFile::get_grid_mapping_rotated_latitude_longitude(const NcVar *grid_map
 
   data.name = rotated_latlon_proj_type;
 
-  // Derive south pole location from the north pole
+  // Derive south pole location from the north pole:
+  // - Reverse the sign of the latitude
+  // - Add 180 to the longitude and switch from degrees east to west
   data.true_lat_south_pole = -1.0 * get_att_value_double(grid_np_lat_att);
   double np_lon = rescale_lon(get_att_value_double(grid_np_lon_att));
-  data.true_lon_south_pole = rescale_lon(-1.0 * (180.0 - fabs(np_lon)));
+  data.true_lon_south_pole = rescale_lon(-1.0 * (np_lon + 180.0));
 
   // Copied from the LatLon data structure
   data.rot_lat_ll = ll_data.lat_ll;

diff --git a/test/xml/unit_plot_data_plane.xml b/test/xml/unit_plot_data_plane.xml
@@ -271,6 +271,20 @@
     </output>
   </test>
 
+  <test name="plot_data_plane_NCCF_rotlatlon">
+    <exec>&MET_BIN;/plot_data_plane</exec>
+    <param> \
+      &DATA_DIR_MODEL;/nccf/prods_op_ukv_rotated_latlon_20220301_00_000.nc \
+      &OUTPUT_DIR;/plot_data_plane/plot_data_plane_NCCF_rotlatlon.ps \
+      'name="surface_altitude"; level="(*,*)";' \
+      -title "Rotated Lat/Lon Surface Altitude" \
+      -v 1
+    </param>
+    <output>
+      <ps>&OUTPUT_DIR;/plot_data_plane/plot_data_plane_NCCF_rotlatlon.ps</ps>
+    </output>
+  </test>
+
   <test name="plot_data_plane_GRIB2_PROB">
     <exec>&MET_BIN;/plot_data_plane</exec>
     <param> \