<div dir="ltr"><div>Dear community,</div><div><br></div><div>I am trying to implement a sponge layer denoted a(x) varying between 0 and 1 in the x-direction so that my velocity forcing relaxes at the targeted value (Uf,Vf)=(0.025,0) m/s at the two edges of my domain after a certain time (defined as some multiples of 3600 s) and (Uf,Vf)=(0,0) i.e. no forcing in the area of interest. My domain has an horizontal extent of Lx = sNx*NPx*dx = 50*40*30 m = 60 km long and I have added the mask function by modifying the right hand side of the Navier Stokes equation gU, gV and gT as followed (see lines implemented in external_forcing.F) to get an exponential decay after a certain time. It seems however that the code is not taking into account the sponge layer at all and keeps resolving Navier Stokes without any mask restrictions (which results in oscillations of U,V and theta in the longer times instead of the expected relaxation). Any idea about this issue ?</div><div><br></div><div>Thank you in advance</div><div><br></div><div>Kind regards</div><div><br></div><div>Kevin.</div><div><br></div><div><br></div></div>