444 |
|
|
445 |
C-- Some advection schemes are better calculated using a multi-dimensional |
C-- Some advection schemes are better calculated using a multi-dimensional |
446 |
C method in the absence of any other terms and, if used, is done here. |
C method in the absence of any other terms and, if used, is done here. |
447 |
cph IF (multiDimAdvection) THEN |
IF (multiDimAdvection) THEN |
448 |
cph IF (tempStepping .AND. |
IF (tempStepping .AND. |
449 |
cph & tempAdvScheme.NE.ENUM_CENTERED_2ND .AND. |
& tempAdvScheme.NE.ENUM_CENTERED_2ND .AND. |
450 |
cph & tempAdvScheme.NE.ENUM_UPWIND_3RD .AND. |
& tempAdvScheme.NE.ENUM_UPWIND_3RD .AND. |
451 |
cph & tempAdvScheme.NE.ENUM_CENTERED_4TH ) |
& tempAdvScheme.NE.ENUM_CENTERED_4TH ) |
452 |
cph & CALL GAD_ADVECTION(bi,bj,tempAdvScheme,GAD_TEMPERATURE,theta, |
& CALL GAD_ADVECTION(bi,bj,tempAdvScheme,GAD_TEMPERATURE,theta, |
453 |
cph U gT, |
U gT, |
454 |
cph I myTime,myIter,myThid) |
I myTime,myIter,myThid) |
455 |
cph IF (saltStepping .AND. |
IF (saltStepping .AND. |
456 |
cph & saltAdvScheme.NE.ENUM_CENTERED_2ND .AND. |
& saltAdvScheme.NE.ENUM_CENTERED_2ND .AND. |
457 |
cph & saltAdvScheme.NE.ENUM_UPWIND_3RD .AND. |
& saltAdvScheme.NE.ENUM_UPWIND_3RD .AND. |
458 |
cph & saltAdvScheme.NE.ENUM_CENTERED_4TH ) |
& saltAdvScheme.NE.ENUM_CENTERED_4TH ) |
459 |
cph & CALL GAD_ADVECTION(bi,bj,saltAdvScheme,GAD_SALINITY,salt, |
& CALL GAD_ADVECTION(bi,bj,saltAdvScheme,GAD_SALINITY,salt, |
460 |
cph U gS, |
U gS, |
461 |
cph I myTime,myIter,myThid) |
I myTime,myIter,myThid) |
462 |
cph ENDIF |
ENDIF |
463 |
|
|
464 |
|
|
465 |
C-- Start of thermodynamics loop |
C-- Start of thermodynamics loop |