at this step the intermediate 'topog_raw.nc' file already has vertica…

…l coordinate named elevation. Threw an error when vcoord had different name

regional_mom6/regional_mom6.py

@@ -1022,7 +1022,7 @@ def bathymetry(
 
         ## Ensure correct encoding
         topog = xr.Dataset(
-            {"depth": (["ny", "nx"], topog[varnames["elevation"]].values)}
+            {"depth": (["ny", "nx"], topog["elevation"].values)}
         )
         topog.attrs["depth"] = "meters"
         topog.attrs["standard_name"] = "topographic depth at T-cell centers"