Log of /MITgcm/model/src/temp_integrate.F

---

Revision 1.11 - (view) (download) (annotate) - [select for diffs]
Tue Jul 22 12:04:09 2014 UTC (9 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65a
Changes since 1.10: +6 -18 lines
Diff to previous 1.10

change implemetation of AB on tracer (doAB_onGtGs=F, not AB on tendency):
 store AB(Tr)=Tr^(n+1/2) into gtNm(m2) instead of directly in theta,salt.
This fix some forcing-term calculation (e.g., relaxation) that should not
 be computed from AB(Tr) but instead directly from Tr.

---

Revision 1.10 - (view) (download) (annotate) - [select for diffs]
Sat Jul 19 16:49:59 2014 UTC (9 years, 10 months ago) by jmc
Branch: MAIN
Changes since 1.9: +54 -12 lines
Diff to previous 1.9

- collect forcing tendency into local 2-D array (passed as argument to
  apply_forcing S/R); add 3-D diagnostics for forcing tendency.
  affects results at machine truncation level (e.g., ocean: if using
  SHORTWAVE_HEATING and atmos:, if using FrictionHeating or (atm_Rq<>0)
  moist correction due to vertical advection).

---

Revision 1.9 - (view) (download) (annotate) - [select for diffs]
Mon Jul 14 22:47:39 2014 UTC (9 years, 10 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint65
Changes since 1.8: +3 -3 lines
Diff to previous 1.8

- drop the bi,bj dimension of argument recip_hFac in S/R DWNSLP_APPLY;
  for NLFS case, use new recip_hFac value in order to conserve tracer.

---

Revision 1.8 - (view) (download) (annotate) - [select for diffs]
Fri Jul 11 18:43:41 2014 UTC (9 years, 10 months ago) by jmc
Branch: MAIN
Changes since 1.7: +17 -14 lines
Diff to previous 1.7

- new file "apply_forcing.F" containing all the code previously in
  external_forcing.F, but with new argument list: pass, as new argument,
  the current level tendency array to update (instead of a direct update
  of the common bloc array). Change the corresponding calls.

---

Revision 1.7 - (view) (download) (annotate) - [select for diffs]
Fri Apr 4 20:54:11 2014 UTC (10 years, 1 month ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64v, checkpoint64w, checkpoint64x, checkpoint64y, checkpoint64z
Changes since 1.6: +12 -9 lines
Diff to previous 1.6

- Start to include explicitly AUTODIFF_OPTIONS.h, COST_OPTIONS.h,
  and CTRL_OPTIONS.h in src files (to enable to skip the ECCO_CPPOPTIONS.h)
  For now, only in pkgs used in verification/hs94.1x64x5.
- Replace ALLOW_AUTODIFF_TAMC by ALLOW_AUTODIFF (except for tape/storage
  which are specific to TAF/TAMC).

---

Revision 1.6 - (view) (download) (annotate) - [select for diffs]
Mon Jan 13 20:49:54 2014 UTC (10 years, 4 months ago) by gforget
Branch: MAIN
CVS Tags: checkpoint64t, checkpoint64u
Changes since 1.5: +8 -3 lines
Diff to previous 1.5

- re-include store directives that were removed on Dec 27th.
- note : this only changes the location of the calc_3d_diffusivity
  recomputation call, which does not change results unless packages
  used in forward (gmredi) are omitted in adjoint (common practice).

---

Revision 1.5 - (view) (download) (annotate) - [select for diffs]
Fri Dec 27 15:46:46 2013 UTC (10 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64s
Changes since 1.4: +120 -39 lines
Diff to previous 1.4

- move calls to GAD_IMPLICIT_R & IMPLDIFF + DWNSLP_APPLY from thermodynamics.F
  to inside temp_integrate.F, salt_integrate.F & ptracers_integrate
- add fZon & fMer as output argument of S/R GAD_CALC_RHS

---

Revision 1.4 - (view) (download) (annotate) - [select for diffs]
Fri Dec 6 01:55:42 2013 UTC (10 years, 5 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64r
Changes since 1.3: +99 -16 lines
Diff to previous 1.3

move calls to GAD_ADVECTION (Multi-Dim advection) inside temp_integrate.F
 and salt_integrate.F ; and similar move for PTRACERS_ADVECTION.

---

Revision 1.3 - (view) (download) (annotate) - [select for diffs]
Wed Nov 27 23:58:25 2013 UTC (10 years, 6 months ago) by jmc
Branch: MAIN
Changes since 1.2: +1 -6 lines
Diff to previous 1.2

AUTODIFF only: in calc_adv_flow.F, recompute rTransKp (instead of copying from
rTrans) to simplify dependencies (no need to store rTrans anymore in calling S/R)

---

Revision 1.2 - (view) (download) (annotate) - [select for diffs]
Thu Nov 21 00:16:55 2013 UTC (10 years, 6 months ago) by jmc
Branch: MAIN
Changes since 1.1: +4 -1 lines
Diff to previous 1.1

add initialisation of rTrans (only for AD)

---

Revision 1.1 - (view) (download) (annotate) - [select for diffs]
Tue Nov 19 16:58:38 2013 UTC (10 years, 6 months ago) by jmc
Branch: MAIN

- rename calc_gt.F to temp_integrate.F (includes call to TIMESTEP_TRACER)
 - rename calc_gs.F to salt_integrate.F (includes call to TIMESTEP_TRACER)
 - remove k from thermodynamics.F and move it to temp_integrate.F
   salt_integrate.F and ptracers_integrate.F (now done inside the tracer
   loop).
 - compute locally (in thermodynamics.F) 3-D velocity field that is used to
   advect tracers; pass it as argument to GAD_ADVECTION, GAD_SOM_ADVECT,
   PTRACERS_ADVECTION, TEMP_INTEGRATE, SALT_INTEGRATE, PTRACERS_INTEGRATE,
   GAD_IMPLICIT_R and PTRACERS_IMPLICIT