[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