2 |
C $Name$ |
C $Name$ |
3 |
|
|
4 |
#include "THSICE_OPTIONS.h" |
#include "THSICE_OPTIONS.h" |
5 |
|
#ifdef ALLOW_GENERIC_ADVDIFF |
6 |
|
# include "GAD_OPTIONS.h" |
7 |
|
#endif |
8 |
|
|
9 |
CBOP |
CBOP |
10 |
C !ROUTINE: THSICE_ADVDIFF |
C !ROUTINE: THSICE_ADVDIFF |
35 |
#include "EEPARAMS.h" |
#include "EEPARAMS.h" |
36 |
#include "PARAMS.h" |
#include "PARAMS.h" |
37 |
#include "GRID.h" |
#include "GRID.h" |
|
#include "GAD.h" |
|
38 |
#include "THSICE_SIZE.h" |
#include "THSICE_SIZE.h" |
39 |
#include "THSICE_PARAMS.h" |
#include "THSICE_PARAMS.h" |
40 |
#include "THSICE_VARS.h" |
#include "THSICE_VARS.h" |
41 |
#include "THSICE_2DYN.h" |
#include "THSICE_2DYN.h" |
42 |
|
#ifdef ALLOW_GENERIC_ADVDIFF |
43 |
|
# include "GAD.h" |
44 |
|
#endif |
45 |
|
|
46 |
C !INPUT PARAMETERS: =================================================== |
C !INPUT PARAMETERS: =================================================== |
47 |
C === Routine arguments === |
C === Routine arguments === |
126 |
IF ( minIcArea.GT.0. _d 0 ) r_minArea = 1. _d 0 / minIcArea |
IF ( minIcArea.GT.0. _d 0 ) r_minArea = 1. _d 0 / minIcArea |
127 |
|
|
128 |
thSIce_multiDimAdv = .TRUE. |
thSIce_multiDimAdv = .TRUE. |
129 |
|
#ifdef ALLOW_GENERIC_ADVDIFF |
130 |
IF ( thSIceAdvScheme.EQ.ENUM_CENTERED_2ND |
IF ( thSIceAdvScheme.EQ.ENUM_CENTERED_2ND |
131 |
& .OR.thSIceAdvScheme.EQ.ENUM_UPWIND_3RD |
& .OR.thSIceAdvScheme.EQ.ENUM_UPWIND_3RD |
132 |
& .OR.thSIceAdvScheme.EQ.ENUM_CENTERED_4TH ) THEN |
& .OR.thSIceAdvScheme.EQ.ENUM_CENTERED_4TH ) THEN |
133 |
thSIce_multiDimAdv = .FALSE. |
thSIce_multiDimAdv = .FALSE. |
134 |
ENDIF |
ENDIF |
135 |
|
#endif /* ALLOW_GENERIC_ADVDIFF */ |
136 |
|
|
137 |
C-- Initialisation (+ build oceanic mask) |
C-- Initialisation (+ build oceanic mask) |
138 |
DO j=1-OLy,sNy+OLy |
DO j=1-OLy,sNy+OLy |