Tracer Open Boundary Conditions Not Taking Effect

[tasha]

I'm doing something wrong and your help in figuring out what is very much appreciated!

I am trying to set an open boundary condition for temperature along the southern boundary, but the warm waters aren't showing up in the model results. The model is running fine, not blowing up, just not having the expected outcome. My cpp defs are below. And ana_tobc.h contains the following (bogus values for now, just want to see that it's working. It is a high latitude model, so temperature of 20 is very high for the model domain.) I have confirmed that it is executing this code. I have similar code in ana_m2obc.h and ana_m3obc.h:

ana_tobc.h

Code: Select all

# ifdef SOUTH_TOBC
      IF (SOUTHERN_EDGE) THEN
        DO k=1,N(ng)
          DO i=istrR,iendR
            BOUNDARY(ng)%t_south(i,k,itemp)=20.0_r8
            BOUNDARY(ng)%t_south(i,k,isalt)=0.0_r8
          END DO
        END DO
      END IF
# endif

my_app.h:

Code: Select all

#define ANA_SSFLUX 
#define ANA_BSFLUX
#define ANA_BTFLUX
#define ANA_BPFLUX

#define ANA_RAIN   
#define ALBEDO
#define LONGWAVE
#define ANA_CLOUD
#define SPONGE

#define BULK_FLUXES 

#define UV_ADV        
#define UV_COR              
#define UV_C4VADVECTION    

#define TS_U3HADVECTION    
#define TS_C4VADVECTION   
#define SALINITY          
#define NONLIN_EOS       
#define TS_PSOURCE        
#define UV_PSOURCE        
#define ANA_PSOURCE

#define SOLVE3D           
#define CURVGRID           
#define MASKING           
#define DJ_GRADPS

#define ICE_MODEL
#ifdef ICE_MODEL
# define ICE_BULK_FLUXES
# define ICE_EVP
# define ICE_MOMENTUM
# define ICE_MOM_BULK
# define ICE_ADVECT
# define ICE_SMOLAR
# define ICE_UPWIND
# define ICE_THERMO
# define ICE_MK
# define ICE_SMOOTH
#endif

#define UV_QDRAG
#define BVF_MIXING

#define ANA_TOBC
#define ANA_M2OBC
#define ANA_M3OBC
#define ANA_FSOBC

#define RADIATION_2D

#define SOUTH_TOBC
#define SOUTH_FSCHAPMAN
#define SOUTH_M2FLATHER
#define SOUTH_M3RADIATION
#define SOUTH_M3NUDGING
#define SOUTH_TRADIATION

[kate]

Code: Select all

#define SOUTH_M3RADIATION
#define SOUTH_M3NUDGING
#define SOUTH_TRADIATION

You're missing the #define SOUTH_TNUDGING

[tasha]

Yes, I was missing SOUTH_TNUDGING. Oops.

Thanks so much for your reply Kate!