blowing up on nested grid

[jitraporn13]

Dear All,

I’m trying to run nested grid, It blows out since the beginning
I noticed on DZmax and courant number on Grid02:
Minimum X-grid spacing, DXmin = 1.00000000E-40 km
Maximum X-grid spacing, DXmax = 5.12482819E+00 km
Minimum Y-grid spacing, DYmin = 1.00000000E-40 km
Maximum Y-grid spacing, DYmax = 5.01111003E+00 km
Minimum Z-grid spacing, DZmin = 5.60947682E-02 m
Maximum Z-grid spacing, DZmax = 4.64003949E+35 m

Minimum barotropic Courant Number = 1.23681706E-03
Maximum barotropic Courant Number = 1.40071410E+56
Maximum Coriolis Courant Number = 1.00000000E+37

When I checked Grid02 by itself. It works well. Please give me some suggestion.

Bests,
Jitraporn

[arango]

The information is in the front of you. If you don't understand that information, you have homework to do before running ROMS. Look the large and a small numbers. Check your grid files. (DXmin, DYmin) cannot be zero!!! It is a trillion times smaller than the scales of electrons... I am sure that it is not what you are modeling.

[wilkin]

When I checked Grid02 by itself. It works well. Please give me some suggestion.

By this do you mean you get the same extreme DXmin, DYmin and DZmax yet it runs?

If so, then the issues you are having (this is a guess) may be associated with a poor match of land/sea mask in the contact points region. If you have "wet" points in grid 2 that map to "dry" points in grid 1, the initial value of zeta in the contact points region may be a special value (1E+35) from the initial conditions file which could explain your extreme DZmax of 4E+35 meters.

My advice is to have the land/sea mask of parent and child grid be exactly the same in the few grid points adjacent to the child perimeter.

[arango]

That can be also a problem. However, the main issue here is that there are one or more grid cells with a value of 1E+40 1/km for both inverse metric factors pm and pn. We cannot multiply those variables with the mask because of division by zero in ROMS. The grid spacing (DX=1/pm and DY=1/pn) is a fundamental physical quantity of the grid. The value 1E+40 and its inverse is quite peculiar and we don't use that special value in ROMS or contact point scripts. Something went wrong in the grid generation. We don't have any information if it is a refined grid or not.

[jitraporn13]

Dear All

Thank you so much for your replying. After your suggestion, I rechecked all grid files and contact file.
The problem was from my contact file.
So, I fixed it already.

Thank you so much,
Jitraporn

[Guotong_Deng]

Hi jitraporn13,

Could you please tell me the problem in your contact file? It seems we have the same problem.

Best

Guotong