Also, there is a typo in the computation of the radiation stress component derivative for the vbar forcing:
Code:
>> git diff master mybranch ROMS/Nonlinear/nearshore_mellor08.h
...
@@ -1116,9 +1116,9 @@
DO j=JstrV,Jend
DO i=Istr,Iend
cff=(Syyl(i,j )- &
- & Syyl(i,j-1))*on_v(i,j)+ &
- & (Sxy_psi(i+1,j)- &
- & Sxy_psi(i ,j))*on_u(i,j)
+ & Syyl(i,j-1))*om_v(i,j)+ &
+ & (Sxyl_psi(i+1,j)- &
+ & Sxyl_psi(i ,j))*on_v(i,j)
rvstr2d(i,j)=cff
!
! Convert units to m2/s2 for output purposes.
Only affects 2d setups with non-equidistant grids. I noticed it while configuring a 2d version of the 3d "Plane beachâ€”spectral waves oblique incident" test case of Haas & Warner (2009) (also used in Uchiyama et al. (2010)). I'll compare it to the the vortex force method later. Uchiyama et al. (2010) do use a 2d vortex-force version of that experiment in addition to the 3-d cases.