sponge layer

[kobl1201]

hi i have some questions

i saw ROMS code about sponge layer

#if defined ADRIA02
!
! Adriatic Sea southern sponge areas.
!
fac=4.0_r8

# if defined UV_VIS2
DO i=IstrT,IendT
DO j=JstrT,MIN(6,JendT)
cff=visc2(ng)+REAL(6-j,r8)*(fac*visc2(ng)-visc2(ng))/6.0_r8
MIXING(ng) % visc2_r(i,j)=cff
MIXING(ng) % visc2_p(i,j)=cff
END DO
DO j=MAX(JstrT,7),JendT
MIXING(ng) % visc2_r(i,j)=0.0_r8
MIXING(ng) % visc2_p(i,j)=0.0_r8
END DO
END DO
# endif

Here is code, i don't know what role is fac variable in the roms model.

I mean if i change fac = 10_r8m what difference between fac = 4_r8 and fac = r_8 in roms model

Thank you~

[susonic]

Hi,
The name in the code 'fac' means factor.

Therefore fac determines that the largest value of the viscosity(or diffusivity) in your model boundary.

For example, if you set 4.0_r8 and your viscosity is 10 m2/s(in ocean.in), you will have maximum 40 m2/s along your model boundary.

-JH