[MITgcm-support] bug in gad_init.F ?
samar khatiwala
spk at ldeo.columbia.edu
Fri Sep 17 11:48:13 EDT 2004
Hello
I may be missing something here, but it seems to me, that if
tempAdvScheme & saltAdvScheme = linear (2,3,4) useMultiDimAdvec will
always be false regardless of what MultiDimAdvection is set to in
'data'.
gad_init.F has the following lines:
useMultiDimAdvec = .FALSE. <-- ?????
useMultiDimAdvec = useMultiDimAdvec.OR.tempMultiDimAdvec
useMultiDimAdvec = useMultiDimAdvec.OR.saltMultiDimAdvec
...
#ifdef ALLOW_PASSIVE_TRACER
tr1_MultiDimAdvec = multiDimAdvection
tr1_AdamsBashforth = .FALSE.
IF ( tracerAdvScheme.EQ.ENUM_CENTERED_2ND
& .OR.tracerAdvScheme.EQ.ENUM_UPWIND_3RD
& .OR.tracerAdvScheme.EQ.ENUM_CENTERED_4TH ) THEN
tr1_MultiDimAdvec = .FALSE.
tr1_AdamsBashforth = .TRUE.
ENDIF
useMultiDimAdvec = useMultiDimAdvec.OR.tr1_MultiDimAdvec
#endif
Without "useMultiDimAdvec = .FALSE.", the code makes sense. Where have I
screwed up?
Thanks, Samar
More information about the MITgcm-support
mailing list