[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