| 1 |
jmc |
1.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. |
Parent Directory
| 
Revision Log
| 
Revision Graph