Hi,
Can any one check this get_bry.F issues:
zetabry_west(j) = zetabry_west(j) - (patm2d(Istr,j)-paref)/(g*rho0)
possible change:
zetabry_west(j) = zetabry_west(j) - (patm2d(IstrR,j)-paref)/(g*rho0)
east, north, and south is correct but west incorrect.
Thank you.
best regards