/[MITgcm]/MITgcm/pkg/ggl90/ggl90_calc.F
ViewVC logotype

Log of /MITgcm/pkg/ggl90/ggl90_calc.F

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

Links to HEAD: (view) (download) (annotate)
Links to MAIN: (view) (download) (annotate)
Sticky Tag:
Revision 1.35 - (view) (download) (annotate) - [select for diffs]
Wed Oct 26 00:49:04 2016 UTC (7 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: HEAD, checkpoint66b, checkpoint66c, checkpoint66d, checkpoint66e, checkpoint66f, checkpoint66g, checkpoint66h, checkpoint66i, checkpoint66j, checkpoint66k, checkpoint66l, checkpoint66m, checkpoint66n, checkpoint66o
Changes since 1.34: +2 -1 lines
Diff to previous 1.34 
- initialise to -1 errCode argument before calling S/R SOLVE_TRIDIAGONAL
  (only useful for SOLVE_DIAGONAL_LOWMEMORY version)
Revision 1.34 - (view) (download) (annotate) - [select for diffs]
Thu Jan 14 17:50:25 2016 UTC (8 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65t, checkpoint65u, checkpoint65v, checkpoint65w, checkpoint65x, checkpoint65y, checkpoint65z, checkpoint66a
Changes since 1.33: +125 -114 lines
Diff to previous 1.33 
- add alternative discretisation of vertical shear, to compute the mean (@ grid
  cell center) of vertical shear compon instead of vert. shear of mean flow.
- add correct parenthesis in ALLOW_GGL90_SMOOTH code to get the same truncated
  results on adjacent faces of cubed-sphere grids.
Revision 1.33 - (view) (download) (annotate) - [select for diffs]
Thu Jul 2 04:42:43 2015 UTC (8 years, 10 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint65n, checkpoint65o, checkpoint65p, checkpoint65q, checkpoint65r, checkpoint65s
Changes since 1.32: +2 -2 lines
Diff to previous 1.32 
Some minor rearrangement to support AD
Revision 1.32 - (view) (download) (annotate) - [select for diffs]
Thu Feb 26 16:45:24 2015 UTC (9 years, 2 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint65k, checkpoint65l, checkpoint65m
Changes since 1.31: +121 -85 lines
Diff to previous 1.31 
  - re-organize code in s/r ggl90_calc a little to move computations that
    are different for idemix into separate loops (changes verification
    experiements at truncation level)
  - requires 2 new 2D fields KappaM, verticalShear
  - make hFacI and recip_hFacI local 3D variables that are passed to
    s/r ggl90_idemix
  - fix and improve some comments
Revision 1.31 - (view) (download) (annotate) - [select for diffs]
Mon Feb 23 21:20:15 2015 UTC (9 years, 2 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65j
Changes since 1.30: +3 -3 lines
Diff to previous 1.30 
- change background vertical diffusivity in vertical mixing pkgs ggl90,
  kl10, my82 and pp81 from temperature diffusivity to salinity diffusivity.
  This makes ptracers default diffusivity (that uses salt diffKr) more
  consistent with vertical mixing schemes.
Revision 1.30 - (view) (download) (annotate) - [select for diffs]
Sat Feb 21 20:11:20 2015 UTC (9 years, 2 months ago) by jmc
Branch: MAIN
Changes since 1.29: +2 -2 lines
Diff to previous 1.29 
add one missing "_d 0" in IDEMIX part
Revision 1.29 - (view) (download) (annotate) - [select for diffs]
Sat Feb 21 17:13:20 2015 UTC (9 years, 2 months ago) by jmc
Branch: MAIN
Changes since 1.28: +42 -44 lines
Diff to previous 1.28 
fix newly added DIAGNOSTICS_FILL calls
Revision 1.28 - (view) (download) (annotate) - [select for diffs]
Sat Feb 21 01:46:39 2015 UTC (9 years, 2 months ago) by mlosch
Branch: MAIN
Changes since 1.27: +17 -2 lines
Diff to previous 1.27 
initialise prTemp = 1. that caused a vermix.ggl90 to fail with
-finit-real=inf (-devel option)

initialise all arrays for safety
ifdefs around variable IDEMIX_RiNumber
move a commented line into the appropriate ifdef-block
Revision 1.27 - (view) (download) (annotate) - [select for diffs]
Thu Feb 19 15:44:12 2015 UTC (9 years, 2 months ago) by mlosch
Branch: MAIN
Changes since 1.26: +103 -17 lines
Diff to previous 1.26 
add IDEMIX (Olbers and Eden, 2013, Eden and Olbers, 2014):
  - code provided by Carsten Eden as an extension of ggl90
  - so far the code is turned on within ggl90 by setting a CPP-flag at
    compile time; a runtime flag implementation is not yet complete
Revision 1.26 - (view) (download) (annotate) - [select for diffs]
Thu Aug 14 16:42:35 2014 UTC (9 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65c, checkpoint65d, checkpoint65e, checkpoint65f, checkpoint65g, checkpoint65h, checkpoint65i
Changes since 1.25: +2 -3 lines
Diff to previous 1.25 
change gTracer argument (drop bi,bj indices) in S/R SOLVE_TRIDIAGONAL
Revision 1.25 - (view) (download) (annotate) - [select for diffs]
Mon Nov 5 13:11:11 2012 UTC (11 years, 6 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint64a, checkpoint64b, checkpoint64c, checkpoint64d, checkpoint64e, checkpoint64f, checkpoint64g, checkpoint64h, checkpoint64i, checkpoint64j, checkpoint64k, checkpoint64l, checkpoint64m, checkpoint64n, checkpoint64o, checkpoint64p, checkpoint64q, checkpoint64r, checkpoint64s, checkpoint64t, checkpoint64u, checkpoint64v, checkpoint64w, checkpoint64x, checkpoint64y, checkpoint64z, checkpoint65, checkpoint65a, checkpoint65b
Changes since 1.24: +1 -4 lines
Diff to previous 1.24 
fix previous modification: remove spurious (and incorrect) masking
Revision 1.24 - (view) (download) (annotate) - [select for diffs]
Mon Nov 5 08:11:45 2012 UTC (11 years, 6 months ago) by mlosch
Branch: MAIN
Changes since 1.23: +22 -11 lines
Diff to previous 1.23 
fix computation of control volumes for the case of
ALLOW_GGL90_HORIZDIFF (probably not used very often): move control
volumes up half a grid cell from C-centered to W-centeredupdating
(pointed out by C. Eden)
Revision 1.23 - (view) (download) (annotate) - [select for diffs]
Wed Aug 8 22:22:42 2012 UTC (11 years, 9 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint63r, checkpoint63s, checkpoint64
Changes since 1.22: +10 -8 lines
Diff to previous 1.22 
- added run time flag mxlSurfFlag to include the code that ensure
  mixing between first and second level (previously included as a comment)
Revision 1.22 - (view) (download) (annotate) - [select for diffs]
Thu Jun 28 15:35:52 2012 UTC (11 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint63o, checkpoint63p, checkpoint63q
Changes since 1.21: +2 -20 lines
Diff to previous 1.21 
use sigmaR for N^2 calculation (this save density computations)
Revision 1.21 - (view) (download) (annotate) - [select for diffs]
Wed Jun 27 22:39:09 2012 UTC (11 years, 10 months ago) by jmc
Branch: MAIN
Changes since 1.20: +27 -22 lines
Diff to previous 1.20 
pass sigmaR as argument (but not yet used); will save several RHO computation
Revision 1.20 - (view) (download) (annotate) - [select for diffs]
Thu Mar 15 15:23:22 2012 UTC (12 years, 1 month ago) by jmc
Branch: MAIN
CVS Tags: checkpoint63k, checkpoint63l, checkpoint63m, checkpoint63n
Changes since 1.19: +19 -14 lines
Diff to previous 1.19 
-rename arrays a,b,c to a3d,b3d,c3d (easier to grep for)
-dirty fix to use vectorized & differentiable solve_tridiagonal.F
Revision 1.19 - (view) (download) (annotate) - [select for diffs]
Thu Aug 19 23:52:37 2010 UTC (13 years, 8 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint62j, checkpoint62k, checkpoint62l, checkpoint62m, checkpoint62n, checkpoint62o, checkpoint62p, checkpoint62q, checkpoint62r, checkpoint62s, checkpoint62t, checkpoint62u, checkpoint62v, checkpoint62w, checkpoint62x, checkpoint62y, checkpoint62z, checkpoint63, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63h, checkpoint63i, checkpoint63j
Changes since 1.18: +224 -216 lines
Diff to previous 1.18 
- add argument "myIter" to S/R GGL90_CALC
- remove 3-D temp array "gTKE" for future TKE (replaced by 2-D tendency
  array if using Horiz. Diffusion)
- rename "ab15,ab05" (not related to AB) to "implDissFac,explDissFac"
Revision 1.18 - (view) (download) (annotate) - [select for diffs]
Wed Aug 11 03:32:29 2010 UTC (13 years, 9 months ago) by gforget
Branch: MAIN
Changes since 1.17: +2 -2 lines
Diff to previous 1.17 
bug fix. gTKE is a local array, whereas
SOLVE_TRIDIAGONAL expects a global array along with bi,bj.
So 1,1 (not bi,bj) needs to be passed as argument to
SOLVE_TRIDIAGONAL in this case.
Revision 1.17 - (view) (download) (annotate) - [select for diffs]
Mon Aug 9 20:34:02 2010 UTC (13 years, 9 months ago) by gforget
Branch: MAIN
Changes since 1.16: +11 -10 lines
Diff to previous 1.16 
Variable initializations and loop indices
Revision 1.16 - (view) (download) (annotate) - [select for diffs]
Fri Aug 6 18:37:05 2010 UTC (13 years, 9 months ago) by gforget
Branch: MAIN
Changes since 1.15: +106 -32 lines
Diff to previous 1.15 
Minor changes in pkg/ggl90:
 o GGL90diffKrS was removed --> always use GGL90diffKr
 o GGL90viscAr was removed --> replaced with GGL90viscArU, GGL90viscArV
 o hack of mxlMaxFlag=2 --> ensure mixing between first and second level (commented out for now)
 o change in max/min operations to ensure that smoothing is ok
 o smoothing of GGL90viscAr was moved to ggl90_calc.F (as done for GGL90diffKr)
 o always use diffKrNrT as background profile (i.e. never use diffKr field)
Revision 1.15 - (view) (download) (annotate) - [select for diffs]
Sat Nov 14 14:27:56 2009 UTC (14 years, 5 months ago) by dfer
Branch: MAIN
CVS Tags: checkpoint61y, checkpoint61z, checkpoint62, checkpoint62a, checkpoint62b, checkpoint62c, checkpoint62d, checkpoint62e, checkpoint62f, checkpoint62g, checkpoint62h, checkpoint62i
Changes since 1.14: +29 -27 lines
Diff to previous 1.14 
Some adjustments (changed the solver --> changes the numerics)
Revision 1.14 - (view) (download) (annotate) - [select for diffs]
Sat Oct 31 03:18:50 2009 UTC (14 years, 6 months ago) by jmc
Branch: MAIN
Changes since 1.13: +21 -12 lines
Diff to previous 1.13 
fix initialisation of SQRTTKE ;
add few comments (marked "C- jmc:) on questionable pieces of code.
Revision 1.13 - (view) (download) (annotate) - [select for diffs]
Tue Oct 27 15:23:22 2009 UTC (14 years, 6 months ago) by dfer
Branch: MAIN
Changes since 1.12: +99 -23 lines
Diff to previous 1.12 
Add option for mixing length min/max
Revision 1.12 - (view) (download) (annotate) - [select for diffs]
Thu Oct 8 20:07:18 2009 UTC (14 years, 7 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint61w, checkpoint61x
Changes since 1.11: +20 -23 lines
Diff to previous 1.11 
modif for vertical profile of background viscosity
Revision 1.11 - (view) (download) (annotate) - [select for diffs]
Fri Jan 30 02:23:56 2009 UTC (15 years, 3 months ago) by dfer
Branch: MAIN
CVS Tags: checkpoint61i, checkpoint61j, checkpoint61k, checkpoint61l, checkpoint61m, checkpoint61n, checkpoint61o, checkpoint61p, checkpoint61q, checkpoint61r, checkpoint61s, checkpoint61t, checkpoint61u, checkpoint61v
Changes since 1.10: +159 -66 lines
Diff to previous 1.10 
A few adjustments
Revision 1.10 - (view) (download) (annotate) - [select for diffs]
Tue Oct 7 19:35:10 2008 UTC (15 years, 7 months ago) by dfer
Branch: MAIN
CVS Tags: checkpoint61e, checkpoint61f, checkpoint61g, checkpoint61h
Changes since 1.9: +60 -33 lines
Diff to previous 1.9 
A bit of rewriting for ggl90 (see tag-index)
Revision 1.9 - (view) (download) (annotate) - [select for diffs]
Mon Sep 29 13:47:23 2008 UTC (15 years, 7 months ago) by dfer
Branch: MAIN
Changes since 1.8: +2 -2 lines
Diff to previous 1.8 
Add missing sqrt(2) factor in computation of mixing length.
Revision 1.8 - (view) (download) (annotate) - [select for diffs]
Mon Aug 11 22:28:06 2008 UTC (15 years, 9 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint61c, checkpoint61d
Changes since 1.7: +43 -42 lines
Diff to previous 1.7 
replace calls to "FIND_RHO" by new version "FIND_RHO_2D".
Revision 1.7 - (view) (download) (annotate) - [select for diffs]
Tue Jun 6 22:15:19 2006 UTC (17 years, 11 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint58g_post, checkpoint58h_post, checkpoint58i_post, checkpoint58j_post, checkpoint58k_post, checkpoint58l_post, checkpoint58m_post, checkpoint58n_post, checkpoint58o_post, checkpoint58p_post, checkpoint58q_post, checkpoint58r_post, checkpoint58s_post, checkpoint58t_post, checkpoint58u_post, checkpoint58v_post, checkpoint58w_post, checkpoint58x_post, checkpoint58y_post, checkpoint59, checkpoint59a, checkpoint59b, checkpoint59c, checkpoint59d, checkpoint59e, checkpoint59f, checkpoint59g, checkpoint59h, checkpoint59i, checkpoint59j, checkpoint59k, checkpoint59l, checkpoint59m, checkpoint59n, checkpoint59o, checkpoint59p, checkpoint59q, checkpoint59r, checkpoint60, checkpoint61, checkpoint61a, checkpoint61b
Changes since 1.6: +8 -4 lines
Diff to previous 1.6 
 fix bug: avoid division by zero
Revision 1.6 - (view) (download) (annotate) - [select for diffs]
Tue Jun 6 16:18:18 2006 UTC (17 years, 11 months ago) by mlosch
Branch: MAIN
Changes since 1.5: +2 -2 lines
Diff to previous 1.5 
 add additional parameter following David's suggestion, interestingly
 ggl90 seems to be quite insensitive to any change in parameters ...
Revision 1.5 - (view) (download) (annotate) - [select for diffs]
Fri Aug 5 23:46:09 2005 UTC (18 years, 9 months ago) by ce107
Branch: MAIN
CVS Tags: checkpint57u_post, checkpoint57q_post, checkpoint57r_post, checkpoint57s_post, checkpoint57t_post, checkpoint57v_post, checkpoint57w_post, checkpoint57x_post, checkpoint57y_post, checkpoint57y_pre, checkpoint57z_post, checkpoint58, checkpoint58a_post, checkpoint58b_post, checkpoint58c_post, checkpoint58d_post, checkpoint58e_post, checkpoint58f_post
Changes since 1.4: +2 -2 lines
Diff to previous 1.4 
Added _d 0 to constants in calls to MIN() to avoid the errors the IBM XL
compilers produce otherwise. This is a recurring error... :-(
Revision 1.4 - (view) (download) (annotate) - [select for diffs]
Sat Dec 4 00:16:49 2004 UTC (19 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint57, checkpoint57a_post, checkpoint57a_pre, checkpoint57b_post, checkpoint57c_post, checkpoint57c_pre, checkpoint57d_post, checkpoint57e_post, checkpoint57f_post, checkpoint57f_pre, checkpoint57g_post, checkpoint57g_pre, checkpoint57h_done, checkpoint57h_post, checkpoint57h_pre, checkpoint57i_post, checkpoint57j_post, checkpoint57k_post, checkpoint57l_post, checkpoint57m_post, checkpoint57n_post, checkpoint57o_post, checkpoint57p_post, eckpoint57e_pre
Changes since 1.3: +2 -2 lines
Diff to previous 1.3 
depth convergence accelerator: replace deltaTtracer by dTtracerLev(k)
Revision 1.3 - (view) (download) (annotate) - [select for diffs]
Sun Oct 17 23:11:59 2004 UTC (19 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint55h_post, checkpoint55i_post, checkpoint55j_post, checkpoint56, checkpoint56a_post, checkpoint56b_post, checkpoint56c_post
Changes since 1.2: +2 -2 lines
Diff to previous 1.2 
allow to set a vertical profile of vertical diffusivity for T & S
Revision 1.2 - (view) (download) (annotate) - [select for diffs]
Mon Sep 27 08:02:04 2004 UTC (19 years, 7 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint55d_post, checkpoint55d_pre, checkpoint55e_post, checkpoint55f_post, checkpoint55g_post
Changes since 1.1: +76 -15 lines
Diff to previous 1.1 
  - add horizontal diffusion of TKE, requires exchanges in
    do_fields_blocking_exchanges, horizontal diffusivity is zero by
    default. In OPA there is no horizontal diffusion of TKE but the
    mixing coefficients are computed from a horizontal average of TKE of
    6 points or so. I think that diffusion has a little more physical
    justification.
  - clean up ggl90_calc in the hope of reducing memory usage (this hope
    was in vain)
  - mask tke-variable in ggl90_init
Revision 1.1 - (view) (download) (annotate) - [select for diffs]
Thu Sep 16 11:27:18 2004 UTC (19 years, 7 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpoint54f_post, checkpoint55, checkpoint55a_post, checkpoint55b_post, checkpoint55c_post 
o initial check-in of TKE-model of Gaspar et al. 1990

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.

  Diffs between and
  Type of Diff should be a
Sort log by:
  ViewVC Help
Powered by ViewVC 1.1.22