[MITgcm] Log of /MITgcm/pkg/gmredi/GMREDI.h

Revision 1.29 - (view) (download) (annotate) - [select for diffs]
Fri Mar 28 04:22:11 2014 UTC (10 years, 1 month ago) by m_bates
Branch: MAIN
CVS Tags: checkpoint64v, checkpoint64w
Changes since 1.28: +7 -2 lines
Diff to previous 1.28

K3D: Added upper and lower bounds to renormalisation factor as namelist variables. Fixed small bug with interpolation to U & V grid. Also added extra documentation and a diagnostic.

Revision 1.28 - (view) (download) (annotate) - [select for diffs]
Wed Jan 1 23:20:48 2014 UTC (10 years, 4 months ago) by m_bates
Branch: MAIN
CVS Tags: checkpoint64s, checkpoint64t, checkpoint64u
Changes since 1.27: +7 -2 lines
Diff to previous 1.27

K3d: Added option to keep diffusivity constant in the mixed layer. Also added some diagnostics.

Revision 1.27 - (view) (download) (annotate) - [select for diffs]
Mon Oct 21 18:46:05 2013 UTC (10 years, 7 months ago) by m_bates
Branch: MAIN
CVS Tags: checkpoint64q, checkpoint64r
Changes since 1.26: +14 -13 lines
Diff to previous 1.26

1) cleaned up and renamed namelist parameters to better reflect their purpose, 2) added some new diagnostics 3) imposed a maximum length for urms length scale and 4) imposed a minimum length for eddy length scale

Revision 1.26 - (view) (download) (annotate) - [select for diffs]
Sat Sep 28 17:59:30 2013 UTC (10 years, 7 months ago) by m_bates
Branch: MAIN
CVS Tags: checkpoint64p
Changes since 1.25: +4 -2 lines
Diff to previous 1.25

1/Introduced namelist logical to turn the PV sheet on and off. 2/ Fixed bug when GM_K3D_likeGM=.TRUE. which gave the wrong diffusivity to the isoneutral diffusion tensor + added some helpful comments.

Revision 1.25 - (view) (download) (annotate) - [select for diffs]
Fri Sep 27 22:34:35 2013 UTC (10 years, 7 months ago) by m_bates
Branch: MAIN
Changes since 1.24: +4 -2 lines
Diff to previous 1.24

In K3d (for gmredi) 1/facility for a minimum bound on urms length scale introduced, 2/ extra diagnostics, 3/ changes to the Eady growth rate calculation

Revision 1.24 - (view) (download) (annotate) - [select for diffs]
Thu Aug 8 22:39:36 2013 UTC (10 years, 9 months ago) by m_bates
Branch: MAIN
CVS Tags: checkpoint64m, checkpoint64n, checkpoint64o
Changes since 1.23: +5 -2 lines
Diff to previous 1.23

Namelist option to ignore beta in the calculation of grad(q) for GM_K3D

Revision 1.23 - (view) (download) (annotate) - [select for diffs]
Thu Jul 11 14:33:23 2013 UTC (10 years, 10 months ago) by m_bates
Branch: MAIN
CVS Tags: checkpoint64k, checkpoint64l
Changes since 1.22: +14 -5 lines
Diff to previous 1.22

Changes associated with the PV eddy closure (GM_K3D) include:
o To improve efficiency the call to solve for the eigenvectors only happen once every GM_K3D_vecFreq seconds.  This required the following changes:
  - read and write pickup files for the eigenvectors and deformation radius (gmredi_read_pickup.F and gmredi_write_pickup.F)
  - making the number of modes (GM_K3D_NModes) a parameter which must be specified at compile time in GMREDI.h
  - A new namelist variable, GM_K3D_vecFreq
  - Added modesC, modesW, modesS and Rdef to the common block
o If the CPP option use_lapack is undefined, then a WKB approximation to the eigenvectors and deformation radius is now used (although, it seems unstable; so for the moment an error is raised in gmredi_check if GM_K3D is defined but use_lapack is not).
o Changed gmredi_calc_eigs returns the deformation radius rather than the deformation wavenumber
o Fixed bug in calculation of tfluxX and tfluxY for the instance where the surface layer is the depth of the water column.
o Added warning messages if there are problems with calculating eigenmodes and eigenvectors
o Cleaned up code
o Improved documentation
o Rationalised diagnostics
o Added some extra startup checks (gmredi_check)

Revision 1.22 - (view) (download) (annotate) - [select for diffs]
Thu Jun 27 14:51:40 2013 UTC (10 years, 10 months ago) by m_bates
Branch: MAIN
CVS Tags: checkpoint64j
Changes since 1.21: +4 -2 lines
Diff to previous 1.21

1/ use calc_oce_mxlayer to calculate the mixed layer depth for the PV eddy closure in gmredi_k3d.  2/ added smoothing option for constant K when using the PV closure (gmredi_k3d).  3/ added some extra checks for when GM_useK3D 4/ cleaned up gmredi_k3d a bit

Revision 1.21 - (view) (download) (annotate) - [select for diffs]
Fri Jun 21 21:56:18 2013 UTC (10 years, 11 months ago) by m_bates
Branch: MAIN
Changes since 1.20: +7 -6 lines
Diff to previous 1.20

added a debugging option to run the eddy PV closure in a GM limit (GM_K3D_likeGM)

Revision 1.20 - (view) (download) (annotate) - [select for diffs]
Fri Jun 21 17:23:30 2013 UTC (10 years, 11 months ago) by m_bates
Branch: MAIN
Changes since 1.19: +45 -3 lines
Diff to previous 1.19

Added a new eddy diffusivity parameterisation pkg/gmredi.  More detailed description in tag-index.

Revision 1.18 - (view) (download) (annotate) - [select for diffs]
Thu Feb 10 21:24:19 2011 UTC (13 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62s, checkpoint62t, checkpoint62u, checkpoint62v, checkpoint62w, checkpoint62x, checkpoint62y, checkpoint62z, checkpoint63
Changes since 1.17: +18 -4 lines
Diff to previous 1.17

add Boundary-Value-Problem bolus transport (Ferrari etal, OM 2010)
 code from Christopher Wolfe (UCSD).

Revision 1.17 - (view) (download) (annotate) - [select for diffs]
Tue Jan 11 00:54:45 2011 UTC (13 years, 4 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62r
Changes since 1.16: +44 -15 lines
Diff to previous 1.16

allow to specify a 2-D horizontal map and a 1-D vertical profile
 as scaling factor for Isopycnal diffusivity (Redi) and GM diffusivity

Revision 1.15 - (view) (download) (annotate) - [select for diffs]
Tue Oct 7 20:15:27 2008 UTC (15 years, 7 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint61e
Changes since 1.14: +52 -32 lines
Diff to previous 1.14

add params for transition layer (fm07)

Revision 1.14 - (view) (download) (annotate) - [select for diffs]
Fri May 30 02:50:16 2008 UTC (15 years, 11 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint60, checkpoint61, checkpoint61a, checkpoint61b, checkpoint61c, checkpoint61d
Changes since 1.13: +4 -2 lines
Diff to previous 1.13

o bridging the gap between eddy stress and GM.
  -> eddyTau is replaced with eddyPsi (eddyTau = f x rho0 x eddyPsi)
      along with a change in CPP option (now ALLOW_EDDYPSI).
  -> when using GM w/ GM_AdvForm:
      The total eddy streamfunction (Psi = eddyPsi + K x Slope)
      is applied either in the tracer Eq. or in momentum Eq.
      depending on data.gmredi (intro. GM_InMomAsStress).
  -> ALLOW_EDDYPSI_CONTROL for estimation purpose.
  The key modifications are in model/src/taueddy_external_forcing.F
  pkg/gmredi/gmredi_calc_*F pkg/gmredi/gmredi_*transport.F

Revision 1.9 - (view) (download) (annotate) - [select for diffs]
Fri Jan 10 00:48:39 2003 UTC (21 years, 4 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint47g_post, checkpoint47h_post
Changes since 1.8: +8 -6 lines
Diff to previous 1.8

Here they are, before they get changed/lost/stolen.
Mostly modif.'s to fix numerical sensitivities.
Gradient checks OK for
- GM_taper_scheme:
  * clipping
  * ac02
  * linear
  * glw91
  * dm95
  * ldd97
- GMREDI_OPTIONS:
  * GM_VISBECK_VARIABLE_K
  * GM_NON_UNITY_DIAGONAL
  * GM_EXTRA_DIAGONAL
  * GM_BOLUS_ADVEC
in conjunction with data.gmredi parameters to be checked in
in a few minutes under verification/carbon/code/

Revision 1.8 - (view) (download) (annotate) - [select for diffs]
Thu Nov 14 22:43:49 2002 UTC (21 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: branch-exfmods-tag, checkpoint47, checkpoint47a_post, checkpoint47b_post, checkpoint47c_post, checkpoint47d_post, checkpoint47d_pre, checkpoint47e_post, checkpoint47f_post
Branch point for: branch-exfmods-curt
Changes since 1.7: +6 -2 lines
Diff to previous 1.7

o * "clean" adjoint code (in terms of extensive recomputations)
    can now be obtained for all GMREDI options (i.e. for
    - GM_VISBECK_VARIABLE_K
    - GM_NON_UNITY_DIAGONAL
    - GM_EXTRA_DIAGONAL
    - GM_BOLUS_ADVEC )
  * However, wrong gradient check problem remains unsolved.
  * New CPP options have been introduced for different
    tapering schemes

Revision 1.6 - (view) (download) (annotate) - [select for diffs]
Sun Dec 16 18:54:49 2001 UTC (22 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint43a-release1mods, checkpoint44, checkpoint44b_post, checkpoint44b_pre, checkpoint44e_post, checkpoint44e_pre, checkpoint44f_post, checkpoint44f_pre, checkpoint44g_post, chkpt44a_post, chkpt44a_pre, chkpt44c_post, chkpt44c_pre, chkpt44d_post, release1-branch-end, release1-branch_branchpoint, release1-branch_tutorials, release1_final_v1
Branch point for: release1-branch, release1_final
Changes since 1.5: +35 -8 lines
Diff to previous 1.5

Modification to the GMREDI package :
 change units of tensor-K arrays, scale now like diffusivity
 initialise all common block arrays in S/R gmredi_init
 add option to use different isopycnal(Redi) & GM diffusivity
 add option to use the advective GM form or the skew-flux form (=default)
 bug in non_unity_diagonal part fixed.

Revision 1.4 - (view) (download) (annotate) - [select for diffs]
Fri Feb 2 21:36:29 2001 UTC (23 years, 3 months ago) by adcroft
Branch: MAIN
Changes since 1.3: +3 -2 lines
Diff to previous 1.3

Merged changes from branch "branch-atmos-merge" into MAIN (checkpoint34)
 - substantial modifications to algorithm sequence (dynamics.F)
 - packaged OBCS, Shapiro filter, Zonal filter, Atmospheric Physics

Revision 1.3 - (view) (download) (annotate) - [select for diffs]
Mon Jan 29 20:07:39 2001 UTC (23 years, 3 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint34
Changes since 1.2: +7 -7 lines
Diff to previous 1.2

Corrected array dimensions.
Temporary scalars introduced to avoid storing by TAMC.

Revision 1.1 - (view) (download) (annotate) - [select for diffs]
Wed Jun 21 19:42:54 2000 UTC (23 years, 11 months ago) by adcroft
Branch: MAIN
CVS Tags: checkpoint29, checkpoint30, checkpoint31

Packaged GM/Redi routines.

This form allows you to request diffs between any two revisions of this file. For each of the two "sides" of the diff, select a symbolic revision name using the selection box, or choose 'Use Text Field' and enter a numeric revision.

Links to HEAD:	(view) (download) (annotate)
Links to checkpoint65b:	(view) (download) (annotate)
Sticky Tag:

	Diffs between and
	Type of Diff should be a

	ViewVC Help
Powered by ViewVC 1.1.22