89 |
# undef SEAICE_ZETA_SMOOTHREG |
# undef SEAICE_ZETA_SMOOTHREG |
90 |
C allow the truncated ellipse rheology (runtime flag SEAICEuseTEM) |
C allow the truncated ellipse rheology (runtime flag SEAICEuseTEM) |
91 |
# undef SEAICE_ALLOW_TEM |
# undef SEAICE_ALLOW_TEM |
92 |
|
C Use LSR vector code; not useful on non-vector machines, because it |
93 |
|
C slows down convergence considerably, but the extra iterations are |
94 |
|
C more than made up by the much faster code on vector machines. For |
95 |
|
C the only regularly test vector machine these flags a specified |
96 |
|
C in the build options file SUPER-UX_SX-8_sxf90_awi, so that we comment |
97 |
|
C them out here. |
98 |
|
C# define SEAICE_VECTORIZE_LSR |
99 |
|
C# ifdef SEAICE_VECTORIZE_LSR |
100 |
|
C Use modified LSR vector code that splits vector loop into two with |
101 |
|
C step size 2. This modification improves the convergence of the vector |
102 |
|
C code dramatically, so that is may actually be useful in general, but |
103 |
|
C that needs to be tested. |
104 |
|
C# define SEAICE_VECTORIZE_LSR_ZEBRA |
105 |
|
C# endif |
106 |
#else /* not SEAICE_CGRID, but old B-grid */ |
#else /* not SEAICE_CGRID, but old B-grid */ |
107 |
C-- By default for B-grid dynamics solver wind stress under sea-ice is |
C-- By default for B-grid dynamics solver wind stress under sea-ice is |
108 |
C set to the same value as it would be if there was no sea-ice. |
C set to the same value as it would be if there was no sea-ice. |