/[MITgcm]/MITgcm/verification/advect_xz/input.nlfs/check_conserve_TS.txt
ViewVC logotype

Contents of /MITgcm/verification/advect_xz/input.nlfs/check_conserve_TS.txt

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


Revision 1.1 - (show annotations) (download)
Tue Dec 6 18:56:16 2011 UTC (12 years, 3 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint64y, checkpoint64x, checkpoint64z, checkpoint64o, checkpoint64a, checkpoint63r, checkpoint64q, checkpoint64p, checkpoint64s, checkpoint64r, checkpoint64u, checkpoint64t, checkpoint64w, checkpoint64v, checkpoint66g, checkpoint66f, checkpoint66e, checkpoint66d, checkpoint66c, checkpoint66b, checkpoint66a, checkpoint64n, checkpoint66o, checkpoint66n, checkpoint66m, checkpoint66l, checkpoint66k, checkpoint66j, checkpoint66i, checkpoint66h, checkpoint64b, checkpoint63h, checkpoint63m, checkpoint64e, checkpoint63i, checkpoint63q, checkpoint64d, checkpoint64c, checkpoint64g, checkpoint64f, checkpoint63j, checkpoint63l, checkpoint65z, checkpoint65x, checkpoint65y, checkpoint63n, checkpoint65r, checkpoint65s, checkpoint65p, checkpoint65q, checkpoint65v, checkpoint65w, checkpoint65t, checkpoint65u, checkpoint65j, checkpoint65k, checkpoint65h, checkpoint65i, checkpoint65n, checkpoint65o, checkpoint65l, checkpoint65m, checkpoint65b, checkpoint65c, checkpoint63k, checkpoint65a, checkpoint65f, checkpoint65g, checkpoint65d, checkpoint65e, checkpoint64i, checkpoint63o, checkpoint63p, checkpoint64h, checkpoint63s, checkpoint64k, checkpoint64, checkpoint65, checkpoint64j, checkpoint64m, checkpoint63g, checkpoint64l, HEAD
File MIME type: text/plain
add note (check_conserve_TS.txt) and matlab script to check global
conservation of T & S with NonLin-FreeSurf and Adams-Bashforth.

1 To check tracer conservation in advect_xz/input.nlfs set-up,
2 and in general when
3 - using NonLin FreeSurf
4 - using or not Adams-Bashforth (used here for Salinity)
5 - without thermal/salt/fresh-water forcing
6
7 From a standard run with input.nlfs parameter files:
8 > cd advect_xz/tr_run.nlfs
9
10 Download (using cvs):
11 > cvs co -d tmp MITgcm_contrib/jmc_script/extract_StD
12 > mv tmp/extract_StD .
13 > rm -rf tmp
14 > cvs co -d tmp MITgcm_contrib/jmc_script/read_StD.m
15 > mv tmp/read_StD.m .
16 > rm -rf tmp
17
18 Split multi-variable stats-diags file into variable specific files:
19 > ./extract_StD dynStDiag.0000000000.txt dynStD std
20
21 Plot global stats evolution using matlab:
22 >> grph_StD_AB
23
24 figure-3, top-panel: "Avr T"
25 evolution of global averaged (potential) Temperature.
26 (should be flat if conserving)
27
28 figure-4, top-panel: "Avr S"
29 evolution of global averaged Salinity
30
31 Since Adams-Bashforth is used (alphAB= 0.6 = 0.5+abEps) for Salt
32 with NonLin-FreeSurf, we don't expect to exactly conserve Salt,
33 but should not drift and exactly conserve a slightly modified
34 quantity: S_ab = S - alphAB*deltaT*gSnm1
35
36 After editing grph_StD_AB.m (commenting out the 2nd spec of alphAB):
37 >> grph_StD_AB
38
39 figure-3 is unchanged but top-panel on figure-4, now plot
40 the evolution of global averaged S_ab.
41 (should be flat if conserving)
42 Note a problem of conservation at the very 1rst iteration,
43 but later on the curve is flat. This problem is not seen
44 after a restart.

  ViewVC Help
Powered by ViewVC 1.1.22