Fixing definition of surface pressure and upper height variables

radixio/examples/radixsnd2arl.cc

@@ -451,7 +451,6 @@ int main(int argc, char **argv)
                                          "VWND"};
 
     thisIndexHeader.var_names.push_back(surfaceVarNames);
-
     for (size_t level = 1; level < inputPressures.size(); ++level)
     {
       thisIndexHeader.var_names.push_back(varNames);
@@ -480,9 +479,8 @@ int main(int argc, char **argv)
         if (level == 0)
         {
           thisRecordHeader.kvar = "PRSS";
-          // Add ground elevation to height (default is 0)
-          float thisHeight      = inputHeights[level] + groundElevation;
-          thisRecordHeader.var1 = thisHeight;
+          float groundPressure  = inputPressures[level];
+          thisRecordHeader.var1 = groundPressure;
           std::vector<std::vector<float>> thisData(
               numberGridCells, std::vector<float>(numberGridCells, thisHeight));
           std::cout << "      pressure...";
@@ -494,7 +492,7 @@ int main(int argc, char **argv)
         else
         {
           thisRecordHeader.kvar = "HGTS";
-          thisRecordHeader.var1 = inputHeights[level];
+          thisRecordHeader.var1 = inputHeights[level] + groundElevation;
           std::vector<std::vector<float>> thisData(
               numberGridCells,
               std::vector<float>(numberGridCells, inputHeights[level]));