Stretched Grid Runs Failing with "Error calling DO_WETDEP" about gchp HOT 7 CLOSED

Thanks for writing @anasali1999. We were recently informed about a bug in GCPy cubed-sphere regridding. We will bring this into the development branch shortly but in the meantime you may want to apply the fix in geoschem/gcpy#311 and then create a new stretched-grid restart file.

from gchp.

HI @anasali1999, thanks for reaching out. Have you validated that your run using a stretch factor of 2.0 compares well with a uniform global run? If yes, I think we can rule out problems with generation of the restart file. I wonder if there is an instability running at higher grid resolution at the particular region of interest. Here are a few things to try:

1. Visually inspect your restart file to make sure nothing looks wrong
2. Turn off wet deposition in setCommonRunSettings.sh and rerun to see where it crashes
3. Reduce the timesteps in setCommonRunSettings.sh. These are automatically changed if the stretched grid high resolution area is above c180. You are approaching that with c172. You can change the threshold in the config to halve both chemistry and dynamic timestep.

from gchp.

Hello everyone, thank you for the suggestions. Sorry I seem to have closed this issue earlier by accident.

@lizziel I reduced all the time steps in setCommonRunSettings.sh by half, this lead to the same DO_WETDEP error.

I then returned the time steps to normal and turned off wet deposition; this time the run failed with "Error calling DO_CONVECTION". Here are all the files for the turned-off wet dep run: wet dep off files.zip

I have yet to inspect the restart files, but I will report back once I compare the working unstretched, working 2.0 factor stretched, and the failing stretched files.

@yantosca This fix looks easy enough for me to try, I will report back once I try it.

from gchp.

Hello

@yantosca I applied the changes to GCPy and generated a new restart file, but the run failed again with the same DO_WETDEP error.

from gchp.

Hi @anasali1999, have you looked through the GEOS-Chem debugging guide? It has a specific section on the issue you are encountering in DO_WETDEP. See if there are any suggestions on that page (found here) that you have not yet tried.

from gchp.

Hello, I have seen that section of the guide but I've been following some other leads in the meantime. I'll try the debugging strategy shown in the guide

from gchp.

Hello, I am closing this issue, my GCHP setup overall was incorrect.

from gchp.