Possible Bugs in Radiation OBC in t3dbc_im.F

[Barbara]

The calcuation of tau in t3dbc_im.F of the latest version is writen as:

Code: Select all

                  tau=Tobc_out(itrc,ng,iwest)
                  IF ((dTdt*dTdx).lt.0.0_r8) THEN
                    tau=Tobc_in(itrc,ng,iwest)
                  ELSE
                    tau=tau*dt(ng)
                  END IF

and it should be:

Code: Select all

                  IF ((dTdt*dTdx).lt.0.0_r8) THEN
                    tau=Tobc_in(ng,iwest)
                  ELSE
                    tau=Tobc_out(ng,iwest)
                  END IF
                  tau=tau*dt(ng)

[arango]

Wow, great catch You should get the award for the most significant bug discovered by a regular ROMS user. I updated the code. See the following trac ticket for details.

Thank you for reporting this important bug. It is possible that it has very little impact in some applications. It just implies that boundary tracer value was not nudged as strongly as we specified.

[patrickm]

Hi Hernan, could that explain the "inconsistency" between surface bulk formulations and OBCs in the simulations that you presented in Rio last year?