Dear All
Someone help me determine my error.
I run solve3D on south china sea domain with open boundary conditions in header file like following:
#undef BULK_FLUXES
#define LONGWAVE_OUT
#define ANA_BSFLUX
#define ANA_BTFLUX /*Analytical kinematic bottom temperature flux.*/
#define ANA_BHFLUX
#define ASSUMED_SHAPE /*Using assumed-shape arrays.*/
#define AVERAGES /*Writing out time-averaged fields.*/
#define AVERAGES_FLUXES /*Writing out time-averaged surface fluxes.*/
/*CURVGRID Orthogonal curvilinear grid.*/
#define DIAGNOSTICS_TS /*Computing and writing tracer diagnostic terms.*/
#define DIAGNOSTICS_UV /*Computing and writing momentum diagnostic terms.*/
#define DJ_GRADPS /*Parabolic Splines density Jacobian (Shchepetkin, 2002).*/
#define DOUBLE_PRECISION /*Double precision arithmetic.*/
#define KANTHA_CLAYSON /*Kantha and Clayson stability function formulation.*/
#define MASKING /*Land/Sea masking.*/
#define MIX_GEO_TS /*Mixing of tracers along geopotential surfaces.*/
#define MIX_S_UV /*Mixing of momentum along constant S-surfaces.*/
#define MY25_MIXING /*Mellor/Yamada Level-2.5 mixing closure.*/
#define NONLINEAR /*Nonlinear Model.*/
#define NONLIN_EOS /*Nonlinear Equation of State for seawater.*/
#define N2S2_HORAVG /*Horizontal smoothing of buoyancy and shear.*/
#define POWER_LAW /*Power-law shape time-averaging barotropic filter.*/
#define PROFILE /*Time profiling activated .*/
#define QCORRECTION /*Surface net heat flux correction.*/
#define K_GSCHEME /*Third-order upstream advection of TKE fields.*/
#define SALINITY /*Using salinity.*/
#define SCORRECTION /*Surface salinity flux correction.*/
#define SOLVE3D /*Solving 3D Primitive Equations.*/
#define SPLINES /*Conservative parabolic spline reconstruction.*/
#define TS_U3HADVECTION /*Third-order upstream horizontal advection of tracers.*/
#define TS_SVADVECTION /*Parabolic splines vertical advection of tracers.*/
#define TS_DIF2 /*Harmonic mixing of tracers.*/
#undef TS_DIF4
#define UV_ADV /*Advection of momentum.*/
#define UV_COR /*Coriolis term.*/
#define UV_U3HADVECTION /*Third-order*/
#define UV_C4VADVECTION /*Fourth-order centered vertical advection of momentum.*/
#define UV_QDRAG /*Quadratic bottom stress.*/
#define UV_VIS2 /*Harmonic mixing of momentum.*/
#define VAR_RHO_2D /*Variable density barotropic mode.*/
/*
** Options for eastern edge open boundary conditions:
*/
#define EAST_FSCHAPMAN /*– Free Surface Gradient Condition*/
#define EAST_M2FLATHER /*– 2D Momentum Radiation Condition*/
#define EAST_M3RADIATION /*– 3D Momentum Radiation Condition*/
#define EAST_TRADIATION /*– Tracers Radiation Condition*/
#define EAST_M2NUDGING /*Eastern edge, 2D momentum, passive/active outflow/inflow.*/
#define EAST_M3NUDGING /*Eastern edge, 3D momentum, passive/active outflow/inflow.*/
#define EAST_TNUDGING /*Eastern edge, tracers, passive/active outflow/inflow.*/
#define EAST_VOLCONS /*Eastern edge, enforce mass conservation.*/
/*
** Options for western edge open boundary conditions:
*/
#define WESTERN_WALL
/*
** Options for northern edge open boundary conditions:
*/
#define NORTH_FSCHAPMAN
#define NORTH_M2FLATHER
#define NORTH_M3RADIATION
#define NORTH_TRADIATION
#define NORTH_M2NUDGING
#define NORTH_M3NUDGING
#define NORTH_TNUDGING
#define NORTH_VOLCONS
/*
** Options for southern edge open boundary conditions:
*/
#define SOUTH_FSCHAPMAN
#define SOUTH_M2FLATHER
#define SOUTH_M3RADIATION
#define SOUTH_TRADIATION
#define SOUTH_M2NUDGING
#define SOUTH_M3NUDGING
#define SOUTH_TNUDGING
#define SOUTH_VOLCONS
#undef SSH_TIDES /*is used to impose tidal elevation.*/
#undef UV_TIDES /*is used to impose tidal currents.*/
#undef RAMP_TIDES /*is used to ramping (over one day) tidal forcing.*/
#undef FSOBC_REDUCED /*is intended to cope with the case where tidal velocity is not available.*/
#undef ADD_FSOBC /*is used to add tidal elevations to processed OBC data.*/
#undef ADD_M2OBC /*is used to add tidal currents to processed OBC data.*/
My result alway have lower temp, salt value at open east bounday move in domain just after some time step. Althout, I think i get bry and clim file from ECCO data not wrong.
lower temp, salt value at open east bounday move in domain
lower temp, salt value at open east bounday move in domain
- Attachments
-
- snapshot2.png (70.51 KiB) Viewed 2699 times
Re: lower temp, salt value at open east bounday move in doma
Why don't you switch off the M2NUDGING?
How did you set the TNUD, M3NUD and OBCFACTOR in ocean.in file?
This one is little different issue but does it work LONGWAVE_OUT with undefined BULK_FLUX?
It may not harm the lateral boundary value but I'm just curious.
-JH
How did you set the TNUD, M3NUD and OBCFACTOR in ocean.in file?
This one is little different issue but does it work LONGWAVE_OUT with undefined BULK_FLUX?
It may not harm the lateral boundary value but I'm just curious.
-JH
Joonho Lee