I have a grid where most of the land is unmasked. After switching to COAWST version 3.1 I started getting spurious wet cells right after a few steps in to the simulation. They are spurious because the cells that get wet have elevations much higher than the tidal amplitude (e.g. -3 m of elevation getting wet with 0.5 m water surface level). Digging into the code I found that a piece of code that applied wet-dry mask to the right hand side momentum equations was missing in the newer version. After I changed it things looked much better. I had this issue with COAWST, however, it looks like Rutgers trunk is missing the same block of code in rev741, too. I have just looked at the cod on Trac, I did not test it. I wonder if anybody else has had a similar issue (either with ROMS or COAWST). I'd happy to listen and learn from their experiences...
Zafer
Code: Select all
# ifdef WET_DRY
DO j=Jstr,Jend
DO i=IstrU,Iend
cff5=ABS(ABS(umask_wet(i,j))-1.0_r8)
cff6=0.5_r8+DSIGN(0.5_r8,rhs_ubar(i,j))*umask_wet(i,j)
cff7=0.5_r8*umask_wet(i,j)*cff5+cff6*(1.0_r8-cff5)
rhs_ubar(i,j)=rhs_ubar(i,j)*cff7
END DO
END DO
DO j=JstrV,Jend
DO i=Istr,Iend
cff5=ABS(ABS(vmask_wet(i,j))-1.0_r8)
cff6=0.5_r8+DSIGN(0.5_r8,rhs_vbar(i,j))*vmask_wet(i,j)
cff7=0.5_r8*vmask_wet(i,j)*cff5+cff6*(1.0_r8-cff5)
rhs_vbar(i,j)=rhs_vbar(i,j)*cff7
END DO
END DO
# endif