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. |