Arbitrary wetting over the land cells

Bug reports, work arounds and fixes

Moderators: arango, robertson

Post Reply
Message
Author
nacholibre
Posts: 81
Joined: Thu Dec 07, 2006 3:14 pm
Location: USGS
Contact:

Arbitrary wetting over the land cells

#1 Unread post by nacholibre »

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
Attachments
with_rhs_wdmask.png
with_rhs_wdmask.png (9.82 KiB) Viewed 3776 times
spurious_wet_cells.png
spurious_wet_cells.png (11.15 KiB) Viewed 3776 times

lanerolle
Posts: 157
Joined: Mon Apr 28, 2003 5:12 pm
Location: NOAA

Re: Arbitrary wetting over the land cells

#2 Unread post by lanerolle »

Hi Zafer,

Which routine(s) in ROMS did you modify to multiply by the wet/dry mask?

Lyon.

Post Reply