1 |
jmc |
1.6 |
C $Header: /u/gcmpack/MITgcm/pkg/shap_filt/SHAP_FILT_OPTIONS.h,v 1.5 2003/10/09 04:19:20 edhill Exp $ |
2 |
jmc |
1.4 |
C $Name: $ |
3 |
adcroft |
1.2 |
|
4 |
|
|
C Header file for package pkg/shap_filt |
5 |
|
|
|
6 |
edhill |
1.5 |
#ifndef SHAP_FILT_OPTIONS_H |
7 |
|
|
#define SHAP_FILT_OPTIONS_H |
8 |
|
|
#include "PACKAGES_CONFIG.h" |
9 |
|
|
#ifdef ALLOW_SHAP_FILT |
10 |
|
|
|
11 |
adcroft |
1.2 |
#include "CPP_OPTIONS.h" |
12 |
|
|
|
13 |
|
|
C Use [1-d_yy^n)(1-d_xx^n] instead of [1-d_xx^n-d_yy^n] |
14 |
|
|
C This changes the spectral response function dramatically. |
15 |
|
|
C You need to do some analysis before changing this option. ;^) |
16 |
|
|
#define SEQUENTIAL_2D_SHAP |
17 |
|
|
|
18 |
|
|
C This invokes the older code which produces the same numbers |
19 |
|
|
C and might even be faster because it does not do multiple |
20 |
|
|
C exchanges within the filter. However, this requires the |
21 |
|
|
C overlap to be sufficiently wide and also does not work |
22 |
|
|
C for arbitrarily arranged tiles (ie. as in cubed-sphere). |
23 |
|
|
C *DO NOT USE THIS OPTION UNLESS YOU REALLY WANT TO* :-( |
24 |
|
|
#undef USE_OLD_SHAPIRO_FILTERS |
25 |
jmc |
1.4 |
|
26 |
|
|
#ifdef USE_OLD_SHAPIRO_FILTERS |
27 |
|
|
C Boundary conditions for U,V Shapiro code |
28 |
|
|
C Horizontal shear is calculated as if the boundaries are no-slip |
29 |
|
|
C Note: option NO_SLIP_SHAP only used in OLD_SHAPIRO_FILTERS ; |
30 |
|
|
C it is replaced by parameter "Shap_noSlip=1." in new S/R. |
31 |
|
|
#undef NO_SLIP_SHAP |
32 |
|
|
#endif |
33 |
|
|
|
34 |
jmc |
1.6 |
C Masking of vorticity in dissipation terms (and Shap-Filt) |
35 |
|
|
C can be different from the masking applied for momentum advection; |
36 |
|
|
C This option allows to use the local S/R: SHAP_FILT_RELVORT3 to |
37 |
|
|
C compute vorticity, instead of pkg/mom_common S/R: MOM_CALC_RELVORT3 |
38 |
|
|
#undef USE_SHAP_CALC_VORTICITY |
39 |
|
|
|
40 |
edhill |
1.5 |
#endif /* ALLOW_SHAP_FILT */ |
41 |
|
|
#endif /* SHAP_FILT_OPTIONS_H */ |