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

Annotation 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 - (hide annotations) (download)
Fri Dec 15 18:11:58 2006 UTC (17 years, 5 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 heimbach 1.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