445 |
ENDIF |
ENDIF |
446 |
#endif /* ALLOW_AIM */ |
#endif /* ALLOW_AIM */ |
447 |
|
|
448 |
|
#ifndef DISABLE_MULTIDIM_ADVECTION |
449 |
|
C The CPP flag DISABLE_MULTIDIM_ADVECTION is currently set in GAD_OPTIONS.h |
450 |
|
C and is currently set only for auto-differentiating. Ideally this would |
451 |
|
C not be necessary. If you need mutli-dimensional advection in a |
452 |
|
C differentiated code, it does work but costs more), you can edit |
453 |
|
C GAD_OPTIONS.h to re-enable it. |
454 |
|
|
455 |
C-- Some advection schemes are better calculated using a multi-dimensional |
C-- Some advection schemes are better calculated using a multi-dimensional |
456 |
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. |
|
|
|
|
#ifdef ALLOW_MULTIDIM_ADVECTION |
|
457 |
IF (multiDimAdvection) THEN |
IF (multiDimAdvection) THEN |
458 |
IF (tempStepping .AND. |
IF (tempStepping .AND. |
459 |
& tempAdvScheme.NE.ENUM_CENTERED_2ND .AND. |
& tempAdvScheme.NE.ENUM_CENTERED_2ND .AND. |
472 |
I myTime,myIter,myThid) |
I myTime,myIter,myThid) |
473 |
ENDIF |
ENDIF |
474 |
ENDIF |
ENDIF |
475 |
#endif /* ALLOW_MULTIDIM_ADVECTION */ |
#endif /* DISABLE_MULTIDIM_ADVECTION */ |
476 |
|
|
477 |
C-- Start of thermodynamics loop |
C-- Start of thermodynamics loop |
478 |
DO k=Nr,1,-1 |
DO k=Nr,1,-1 |