/[MITgcm]/MITgcm/verification/lab_sea/code_ad/GMREDI_OPTIONS.h
ViewVC logotype

Contents of /MITgcm/verification/lab_sea/code_ad/GMREDI_OPTIONS.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Fri Dec 15 18:11:58 2006 UTC (17 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint58u_post, checkpoint58w_post, checkpoint58x_post, checkpoint58t_post, checkpoint59d, checkpoint59a, checkpoint59c, checkpoint59b, checkpoint59, checkpoint58y_post, checkpoint58v_post
File MIME type: text/plain
Initial check-in of lab_sea adjoint
(very similar to natl_box_adjoint/code_seaice)
* used C-grid and EVP solver
* used KPP, GM/Redi
* uses inAdExact=.TRUE.
* gradient checks fail for both ADM and TLM
  due to sensitive snow melt code in seaice_grwoth (look for "JZ")
  f.d. gradients actually seem more problematic than ADM, TLM gradients;
  issue is related to forward mismatch between g77 an ifort

1 C $Header: /u/gcmpack/MITgcm/verification/lab_sea/code_ad_seaice_ocean/GMREDI_OPTIONS.h,v 1.1 2005/09/10 14:17:07 heimbach Exp $
2 C $Name: $
3
4 C CPP options file for GM/Redi package
5 C
6 C Use this file for selecting options within the GM/Redi package
7 C
8 #ifndef GMREDI_OPTIONS_H
9 #define GMREDI_OPTIONS_H
10 #include "PACKAGES_CONFIG.h"
11
12 #include "CPP_OPTIONS.h"
13
14 #ifdef ALLOW_GMREDI
15
16 C Designed to simplify the Ajoint code:
17 C exclude the clipping/tapering part of the code that is not used
18 #define GM_EXCLUDE_CLIPPING
19 #define GM_EXCLUDE_AC02_TAP
20 #undef GM_EXCLUDE_TAPERING
21
22 C This allows to use Visbeck et al formulation to compute K_GM+Redi
23 #undef GM_VISBECK_VARIABLE_K
24
25 C This allows the leading diagonal (top two rows) to be non-unity
26 C (a feature required when tapering adiabatically).
27 #define GM_NON_UNITY_DIAGONAL
28
29 C Allows to use different values of K_GM and K_Redi ; also to
30 C be used with the advective form (Bolus velocity) of GM
31 #undef GM_EXTRA_DIAGONAL
32
33 C Allows to use the advective form (Bolus velocity) of GM
34 C instead of the Skew-Flux form (=default)
35 #undef GM_BOLUS_ADVEC
36
37 C Following option avoids specific recomputation in adjoint
38 C routines of gmredi_x/y/rtransport
39 C It's not needed, only for tests, and very memory-consuming
40 #undef GM_AUTODIFF_EXCESSIVE_STORE
41
42 #endif /* ALLOW_GMREDI */
43 #endif /* GMREDI_OPTIONS_H */
44

  ViewVC Help
Powered by ViewVC 1.1.22