two-way nesting with wet-dry treatment

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
klettmann
Posts: 6
Joined: Mon Mar 31, 2008 8:28 pm
Location: ICBM, University of Oldenburg

two-way nesting with wet-dry treatment

#1 Unread post by 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

Post Reply