Using "scale_region" tool for the 60 to 3 km mesh

[anias]

Hello,
I am using the "scale_region" code to convert a 60–3 km mesh to a 20–1 km mesh, but the result looks strange. Even when I use a reduction factor of 2 instead of 3, the resulting mesh still appears incorrect, as shown in the attached images.
Am I using "scale_region" correctly?
Ania Sindermann

 

[mgduda]

What arguments are you providing to the scale_region tool? Also, how large is your regional cut-out from the 60-3 km mesh?

 

[anias]

The arguments are:
python scale_region.py panama.static.nc panama_reescalat.static.nc 2 8.5 -80.0
The regional mesh was done with the .pts file:

Name: panama
Type: circle
Point: 8.5, -80.0
radius: 4500000.0

I attached an image with its domain

 

[mgduda]

Apologies for the slow reply! Everything in your workflow -- the definition of a limited-area and the arguments to the scale_region.py script -- look correct to me. Just in case this is an issue in plotting the domain resolution, would you be able to attach the script that you've used to plot me contours of mesh resolution?

 

[anias]

Hi, I reviewed my script for plotting the mesh contours, and it turns out the issue was in the script itself.
My code wasn’t accounting for the change from m² to sr in the AreaCell units when the regional mesh is scaled.

 

[mgduda]

Thanks for following up, and glad to hear that you've found the issue.