555 |
|
|
556 |
C- Compute effective gW_[n+1/2] terms (including Adams-Bashforth weights) |
C- Compute effective gW_[n+1/2] terms (including Adams-Bashforth weights) |
557 |
C and save gW_[n] into gwNm1 for the next time step. |
C and save gW_[n] into gwNm1 for the next time step. |
558 |
c#ifdef ALLOW_ADAMSBASHFORTH_3 |
#ifdef ALLOW_ADAMSBASHFORTH_3 |
559 |
c CALL ADAMS_BASHFORTH3( |
CALL ADAMS_BASHFORTH3( |
560 |
c I bi, bj, k, |
I bi, bj, k, |
561 |
c U gW, gwNm, |
U gW, gwNm, |
562 |
c I nHydStartAB, myIter, myThid ) |
I nHydStartAB, myIter, myThid ) |
563 |
c#else /* ALLOW_ADAMSBASHFORTH_3 */ |
#else /* ALLOW_ADAMSBASHFORTH_3 */ |
564 |
CALL ADAMS_BASHFORTH2( |
CALL ADAMS_BASHFORTH2( |
565 |
I bi, bj, k, |
I bi, bj, k, |
566 |
U gW, gwNm1, |
U gW, gwNm1, |
567 |
I nHydStartAB, myIter, myThid ) |
I nHydStartAB, myIter, myThid ) |
568 |
c#endif /* ALLOW_ADAMSBASHFORTH_3 */ |
#endif /* ALLOW_ADAMSBASHFORTH_3 */ |
569 |
|
|
570 |
C-- Dissipation term outside the Adams-Bashforth: |
C-- Dissipation term outside the Adams-Bashforth: |
571 |
IF ( momViscosity .AND. .NOT.momDissip_In_AB ) THEN |
IF ( momViscosity .AND. .NOT.momDissip_In_AB ) THEN |