This may or may not be similar to a previous posting (Nov 8 : scale_region strange behavior), but I, too, am having difficulty with scale_region.py. To summarize, I used grid-rotate and MPAS-Limited-Area-master on the 46km-12km variable mesh to center and regionalize my domain over the Gulf (Helene.grid.nc). Then I ran two simulations, one using Helene.grid.nc, and another after scaling from 46km-12km to 15km-4km using scale_region.py (./scale_region.py Helene.grid.nc Helene.rgrid.nc 3 22 -88). Aside from changing config_dt based on resolution, the namelist and streams files are exactly the same for each simulation (graph.info as well).
The Helele.grid and Helene.rgrid mesh is attached (orig_mesh and scaled_mesh, respectively). I don't believe there are issues here, and the grid, init, and static files look ok. Problems arise when simulations are performed. Figures Helene_unscaled_ts1, ts5, and ts10 show MSLP and 850 hPa winds at the initial timestep, hour 5, and hour 10, respectively. It looks fine. However, when using Helene.rgrid.nc, the hurricane quickly dissipates (Helene_scaled_ts1, ts5, ts10). The log file shows no errors, and the simulation does not terminate prematurely.
Am I missing an important step or setting? Any suggestions are greatly appreciated.
I am using MPAS version 8.3.1
The Helele.grid and Helene.rgrid mesh is attached (orig_mesh and scaled_mesh, respectively). I don't believe there are issues here, and the grid, init, and static files look ok. Problems arise when simulations are performed. Figures Helene_unscaled_ts1, ts5, and ts10 show MSLP and 850 hPa winds at the initial timestep, hour 5, and hour 10, respectively. It looks fine. However, when using Helene.rgrid.nc, the hurricane quickly dissipates (Helene_scaled_ts1, ts5, ts10). The log file shows no errors, and the simulation does not terminate prematurely.
Am I missing an important step or setting? Any suggestions are greatly appreciated.
I am using MPAS version 8.3.1
