lower temp, salt value at open east bounday move in domain

General scientific issues regarding ROMS

Moderators: arango, robertson

Post Reply
Message
Author
poppai
Posts: 9
Joined: Wed Sep 30, 2009 6:05 pm
Location: IMHEN
Contact:

lower temp, salt value at open east bounday move in domain

#1 Unread post by poppai »

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.
Attachments
snapshot2.png
snapshot2.png (70.51 KiB) Viewed 2356 times

User avatar
susonic
Posts: 167
Joined: Tue Aug 21, 2007 5:44 pm
Location: UST21 / Korea
Contact:

Re: lower temp, salt value at open east bounday move in doma

#2 Unread post by susonic »

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
Joonho Lee

Post Reply