about m3obc.h Question

[pip39]

l modify 'm2obc.h' variable
'analytical.f90' successfully produce the expected variables

but
l modify 'm3obc.h' variable
'analytical.f90' no produce the expected variables


i already #define ANA_M3OBC in 'estuary_test.h'


Help me to explain




---'m2obc.h'------
!-----------------------------------------------------------------------
! 2D momentum open boundary conditions.
!-----------------------------------------------------------------------!

#if defined ESTUARY_TEST
# ifdef WEST_M2OBC
IF (WESTERN_EDGE) THEN
DO j=JstrR,JendR
BOUNDARY(ng)%ubar_west(j)=3.0_r8
END DO
DO j=Jstr,JendR
BOUNDARY(ng)%vbar_west(j)=0.0_r8
END DO



---analytical.f90-----
!-----------------------------------------------------------------------
! 2D momentum open boundary conditions.
!-----------------------------------------------------------------------!

IF (Istr.eq.1) THEN
DO j=JstrR,JendR
BOUNDARY(ng)%ubar_west(j)=3.0_r8
END DO
DO j=Jstr,JendR
BOUNDARY(ng)%vbar_west(j)=0.0_r8
END DO
END IF
RETURN
END SUBROUTINE ana_m2obc_tile
SUBROUTINE ana_nudgcoef (ng, tile, model)





---'m3obc.h'----------
!-----------------------------------------------------------------------
! 3D momentum open boundary conditions.
!-----------------------------------------------------------------------#

if defined ESTUARY_TEST
# ifdef WEST_M3OBC

IF (WESTERN_EDGE) THEN
DO k=1,N(ng)
DO j=JstrR,JendR
BOUNDARY(ng)%u_west(j,k)=2.0_r8
END DO
DO j=Jstr,JendR
BOUNDARY(ng)%v_west(j,k)=0.0_r8
END DO
END DO



---analytical.f90-----
!-----------------------------------------------------------------------
! 3D momentum open boundary conditions.
!-----------------------------------------------------------------------

RETURN
END SUBROUTINE ana_m3obc_tile
SUBROUTINE ana_nudgcoef (ng, tile, model)

[Entrapmen]

What are the boundary conditions for your domain?

[pip39]

What are the boundary conditions for your domain?

#define UV_ADV
#define UV_LOGDRAG
#define TS_U3HADVECTION
#define SALINITY
#define SOLVE3D
#define SPLINES
#define SEDIMENT
#ifdef SEDIMENT
# define SUSPLOAD
#endif
#define AVERAGES
#define AVERAGES_AKV
#define AVERAGES_AKS
#define NORTHERN_WALL
#define SOUTHERN_WALL
#define EAST_FSGRADIENT
#define EAST_M2CLAMPED
#define EAST_M3GRADIENT
#define EAST_TCLAMPED
#define WEST_FSCHAPMAN
#define WEST_M2CLAMPED
#undef WEST_M2REDUCED
#undef FSOBC_REDUCED
#define WEST_M3GRADIENT
#define WEST_TRADIATION
#define WEST_TNUDGING
#define GLS_MIXING
#undef MY25_MIXING
#if defined GLS_MIXING || defined MY25_MIXING
# define KANTHA_CLAYSON
# undef CANUTO_A
# define N2S2_HORAVG
#endif
#define ANA_GRID
#define ANA_INITIAL
#define ANA_SEDIMENT
#define ANA_SMFLUX
#define ANA_STFLUX
#define ANA_BTFLUX
#define ANA_SSFLUX
#define ANA_BSFLUX
#define ANA_SPFLUX
#define ANA_BPFLUX
#define ANA_FSOBC
#define ANA_M2OBC
#define ANA_M3OBC
#define ANA_TOBC

[kate]

Perhaps it's because the gradient condition doesn't need a boundary value. The boundary chunk is surrounded by '# ifdef WEST_M3OBS' which is defined in globaldefs.h:

Code: Select all

#if (defined WEST_M3RADIATION  && defined WEST_M3NUDGING)  || \
     defined WEST_M3CLAMPED
# define WEST_M3OBC
#endif