[MITgcm-support] Tides in OBCS
Bahman Alavi
bahman.al at gmail.com
Tue Apr 18 16:50:18 EDT 2017
Dear all
I am trying to use tidal force on open boundaries in my model. The file
"data.obcs" is like:
# Open-boundaries
&OBCS_PARM01
OB_Ieast=160*-1,
OB_Iwest=288*1,
OB_Jsouth=114*1
&
&OBCS_PARM02
OBCSfixTopo=.TRUE.,
useOBCSprescribe=.TRUE.,
OBWtFile='sstW.bin',
OBWs=40.8,
OBEtFile='sstE.bin',
OBEsFile='sssE.bin',
OBStFile='sstS.bin',
OBSsFile='sssS.bin',
useOBCStides = .TRUE.,
tidalPeriod =
44714.16,43200.,45569.88,43081.92,86164.2,92949.48,86637.24,96726.24,1180295.64,2380716.,
OBSamFile='OBSam.bin',
OBEamFile='OBEam.bin',
OBWamFile='OBWam.bin',
OBSphFile='OBSph.bin',
OBEphFile='OBEph.bin',
OBWphFile='OBWph.bin',
&
and change the "OBCS_OPTIONS.h" as below:
#define ALLOW_OBCS_PRESCRIBE
C Enable OB conditions following Stevens (1990)
#undef ALLOW_OBCS_STEVENS
C This includes hooks to sponge layer treatment of uvel, vvel
#undef ALLOW_OBCS_SPONGE
C Include hooks to sponge layer treatment of pkg/seaice variables
#define ALLOW_OBCS_SEAICE_SPONGE
C balance barotropic velocity
#undef ALLOW_OBCS_BALANCE
C add tidal contributions to normal OB flow
#define ALLOW_OBCS_TIDES
but in after run, the output of OBCS_CHECK is
(PID.TID 0000.0001) == Packages configuration : Check & print summary ==
(PID.TID 0000.0001)
(PID.TID 0000.0001) OBCS_CHECK: #define ALLOW_OBCS
(PID.TID 0000.0001) OBCS_CHECK: start summary:
(PID.TID 0000.0001) useOBCSprescribe = /* prescribe OB values */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) useOBCSbalance = /* balance the flow through OB */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) OBCS_uvApplyFac = /* Factor to apply to U,V 2nd
column/row */
(PID.TID 0000.0001) 1.000000000000000E+00
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) OBCS_u1_adv_T = /* Temp uses upwind adv-scheme @ OB */
(PID.TID 0000.0001) 0
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) OBCS_u1_adv_S = /* Salt uses upwind adv-scheme @ OB */
(PID.TID 0000.0001) 0
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) OBCS_monitorFreq = /* monitor output frequency [s] */
(PID.TID 0000.0001) 1.200000000000000E+03
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) OBCS_monSelect = /* select group of variables to
monitor */
(PID.TID 0000.0001) 0
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) useOBCStides = /* apply tidal forcing through OB */
(PID.TID 0000.0001) F
(PID.TID 0000.0001) ;
(PID.TID 0000.0001) tidalPeriod = /* (s) */
(PID.TID 0000.0001) 10 @ 0.000000000000000E+00 /* I = 1:
10 */
that means the "useOBCStide" and "usePRESCRIBE" do not changed to true. I
can't understand what or where went wrong?
Thanks for your help
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20170419/e2d4e2f9/attachment.htm>
More information about the MITgcm-support
mailing list