Add boundary call to scale_omega

[jcwarner]

Can we add a boundary call in scale_omega so that the output of omega is consistent during periodic simulations? When used now, the variable omega (W) is never filled on the outer RHO-points (since it is not used during computations). But for output purposes, especially periodic simulations, it would be nice to see omega be uniform in the along-periodic direction, instead of having zeros along that boundary. This can be accomplished by adding to scale_omega:

1) 
      USE bc_3d_mod, ONLY : bc_w3d_tile

2)

!
!  Set lateral boundary conditions.
!
      CALL bc_w3d_tile (ng, tile,                       &
     &                  LBi, UBi, LBj, UBj, 0, N(ng),   &
     &                  Wscl)

thanks, -j

[arango]

Yes, good catch. This is only needed in periodic boundary condition since the IstrR, IendR, JstrR, JendR ranges are defined differently. So instead, we need:

# if defined EW_PERIODIC || defined NS_PERIODIC
!
      USE exchange_3d_mod, ONLY : exchange_w3d_tile
# endif

      ...

# if defined EW_PERIODIC || defined NS_PERIODIC
!
!  Exchange boundary data.
!
      CALL exchange_w3d_tile (ng, tile,                                 &
     &                        LBi, UBi, LBj, UBj, LBk, UBk,             &
     &                        Wscl)
# endif

In non-periodic applications, the values in the outer RHO-points were assigned in routine omega.F. They are rescaled for output purposes in the history NetCDF file.