Issues with scale_region.py simulations

[skIANJ86]

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

 

[Ming Chen]

Hi,

Thank you for the detailed description of your issue.

In your namelist.atmophere, I found that you turned off convection, i.e.,

config_physics_suite = 'convection_permitting'
config_convection_scheme = 'off'

Is there any special reason that you don't want to activate convection? Did you also turn off convection in your 46-12km run?

 

[skIANJ86]

Hi Ming,

Thank you for looking into this. Both the 46-12km and the 15-4km were run with the default "convection_permitting" suite. They were also run using the "mesoscale_reference" suite, but the results were qualitatively similar.

The convection off setting was from the follow-up calculations. Sorry for the confusion. We implemented this setting to determine if the scale-aware aspect of the Grell-Frietas was the source of the issue, as areas at/near the tropical cyclone clearly exhibited something unusual. This, however, was not the issue, as convection off runs were also dissipating the TC rapidly.

We note that similar issues arise from another TC case in the Western Pacific.

 

[Ming Chen]

Hi,

Thank you for the kind information. Please upload the following files for me to take a look:

(1) namelist.input for grid_rotate
(2) points.txt for regronal grid rotate
(3) namelist.init_atmopshere and streams.init_atmopshere for static, initial and lbc conditions
(4) namelist.atmosphere and streams.atmosphere for mpas run
(5) Helene.grid.nc and Helene.rgrid.nc
(6) Helene.static.nc and Helene.rstatic.nc (static data corresponding to 46km-12km and 15-4km meshes, respectively)

Please upload these files to upload to our Box storage and let me know after you are done. Thank you.

I will look into this case and get back to you once I figure out what is wrong.

 

[skIANJ86]

Hi,

I have uploaded the requested files. Thank you for looking into this so close to Thanksgiving.

 

[skIANJ86]

Hello,

I was wondering if there were any further updates for this issue.
Thank you in advance.

 

[skIANJ86]

So sorry to bother, but are there any updates on this issue?

 

[Ming Chen]

Hi,
Thank you for checking. I apologize for the delayed answer due to the holiday season and the tutorial we just finished. Sorry that I haven't found a solution to your issue. I will get back to you once I know for sure what is wrong.