Setting of Sponge Layers in Ocean Modeling

Report or discuss software problems and other woes

Moderators: arango, robertson

Post Reply
Message
Author
hyc006
Posts: 24
Joined: Thu Nov 30, 2023 2:12 am
Location: UCSD

Setting of Sponge Layers in Ocean Modeling

#1 Unread post by hyc006 »

I'm currently attempting to implement a model with a sponge layer, but I'm encountering confusion regarding the specification of the sponge layer's width. Below is the relevant excerpt from my code in the ana_sponge.h file:

Code: Select all

#elif defined SHALLOW
!
!  Shallow Ocean With An Island, Island Wake be reduce by sponge Layer at Boundary.
!
      width=10.0_r8
!
!  Set horizontal mixing factor over 10 grid points in the eastern boundary.
!
        MIXING(ng) % visc2_r(IstrT:IendT,JstrT:JendT)=0.0_r8
        MIXING(ng) % visc2_p(IstrT:IendT,JstrT:JendT)=0.0_r8
!
!  Eastern edge.
!
      DO i=MAX(IstrT,IendT-INT(width)),IendT
        DO j=JstrT,JendT
          MIXING(ng) % visc2_r(i,j)=400.0_r8
          MIXING(ng) % visc2_p(i,j)=400.0_r8
        END DO
      END DO

#endif 
Initially, I interpreted the width parameter as indicating that the sponge layer would be applied to the easternmost 10 grid points. However, the resulting simulation suggests that the sponge layer extends further than intended.
width10.jpg
Furthermore, I attempted to adjust the width by setting it to width = 3, but this adjustment did not seem to significantly alter the behavior of the sponge layer.
width3.jpg
Could someone offer insight into what might be amiss in my implementation? Any assistance would be greatly appreciated.

Post Reply