two-way nesting with wet-dry treatment

[klettmann]

Dear all,

I hope I did not miss the answer to the following question inside the forum.

We are using two-way nesting for tidal-flat applications.

When comparing the free surface in coarse and fine grid within the same region, one can see a difference in dry areas. Whereas in fine grid, zeta meets the D_crit criterion, zeta in coarse grid is just zero.

This setting is done in fine2coarse, where the interpolation weights from fine grid in dry regions sets zeta-coarse to zero.

In "normal" settings (without nesting), there are code lines that guaranty zeta to meet the D_crit criterion like these in step_2d... :

Code: Select all

# if defined WET_DRY && defined MASKING
          zeta(i,j,knew)=zeta(i,j,knew)+                                &
     &                   (Dcrit(ng)-h(i,j))*(1.0_r8-rmask(i,j))
# endif 

This step does not seem to be applied to zeta values in coarse grid after calling the routine fine2coarse.

My question is now, if this inconsistency might cause any harm?

Best, Karsten