1 |
jmc |
1.4 |
C $Header: /u/gcmpack/MITgcm/pkg/shap_filt/SHAP_FILT_OPTIONS.h,v 1.3 2001/02/04 14:38:50 cnh Exp $ |
2 |
|
|
C $Name: $ |
3 |
adcroft |
1.2 |
|
4 |
|
|
C Header file for package pkg/shap_filt |
5 |
|
|
|
6 |
|
|
#include "CPP_OPTIONS.h" |
7 |
|
|
|
8 |
|
|
C Use [1-d_yy^n)(1-d_xx^n] instead of [1-d_xx^n-d_yy^n] |
9 |
|
|
C This changes the spectral response function dramatically. |
10 |
|
|
C You need to do some analysis before changing this option. ;^) |
11 |
|
|
#define SEQUENTIAL_2D_SHAP |
12 |
|
|
|
13 |
|
|
C This invokes the older code which produces the same numbers |
14 |
|
|
C and might even be faster because it does not do multiple |
15 |
|
|
C exchanges within the filter. However, this requires the |
16 |
|
|
C overlap to be sufficiently wide and also does not work |
17 |
|
|
C for arbitrarily arranged tiles (ie. as in cubed-sphere). |
18 |
|
|
C *DO NOT USE THIS OPTION UNLESS YOU REALLY WANT TO* :-( |
19 |
|
|
#undef USE_OLD_SHAPIRO_FILTERS |
20 |
jmc |
1.4 |
|
21 |
|
|
#ifdef USE_OLD_SHAPIRO_FILTERS |
22 |
|
|
C Boundary conditions for U,V Shapiro code |
23 |
|
|
C Horizontal shear is calculated as if the boundaries are no-slip |
24 |
|
|
C Note: option NO_SLIP_SHAP only used in OLD_SHAPIRO_FILTERS ; |
25 |
|
|
C it is replaced by parameter "Shap_noSlip=1." in new S/R. |
26 |
|
|
#undef NO_SLIP_SHAP |
27 |
|
|
#endif |
28 |
|
|
|