1 |
C $Header: /u/gcmpack/MITgcm/pkg/shap_filt/SHAP_FILT_OPTIONS.h,v 1.5 2003/10/09 04:19:20 edhill Exp $ |
2 |
C $Name: $ |
3 |
|
4 |
C Header file for package pkg/shap_filt |
5 |
|
6 |
#ifndef SHAP_FILT_OPTIONS_H |
7 |
#define SHAP_FILT_OPTIONS_H |
8 |
#include "PACKAGES_CONFIG.h" |
9 |
#ifdef ALLOW_SHAP_FILT |
10 |
|
11 |
#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 |
|
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 |
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 |
#endif /* ALLOW_SHAP_FILT */ |
41 |
#endif /* SHAP_FILT_OPTIONS_H */ |