Warning, /verification/advect_xz/input.nlfs/check_conserve_TS.txt is written in an unsupported language. File is not indexed.
view on githubraw file Latest commit fc2c9620 on 2011-12-06 18:56:16 UTC
fc2c9620b4 Jean*0001 To check tracer conservation in advect_xz/input.nlfs set-up,
0002 and in general when
0003 - using NonLin FreeSurf
0004 - using or not Adams-Bashforth (used here for Salinity)
0005 - without thermal/salt/fresh-water forcing
0006
0007 From a standard run with input.nlfs parameter files:
0008 > cd advect_xz/tr_run.nlfs
0009
0010 Download (using cvs):
0011 > cvs co -d tmp MITgcm_contrib/jmc_script/extract_StD
0012 > mv tmp/extract_StD .
0013 > rm -rf tmp
0014 > cvs co -d tmp MITgcm_contrib/jmc_script/read_StD.m
0015 > mv tmp/read_StD.m .
0016 > rm -rf tmp
0017
0018 Split multi-variable stats-diags file into variable specific files:
0019 > ./extract_StD dynStDiag.0000000000.txt dynStD std
0020
0021 Plot global stats evolution using matlab:
0022 >> grph_StD_AB
0023
0024 figure-3, top-panel: "Avr T"
0025 evolution of global averaged (potential) Temperature.
0026 (should be flat if conserving)
0027
0028 figure-4, top-panel: "Avr S"
0029 evolution of global averaged Salinity
0030
0031 Since Adams-Bashforth is used (alphAB= 0.6 = 0.5+abEps) for Salt
0032 with NonLin-FreeSurf, we don't expect to exactly conserve Salt,
0033 but should not drift and exactly conserve a slightly modified
0034 quantity: S_ab = S - alphAB*deltaT*gSnm1
0035
0036 After editing grph_StD_AB.m (commenting out the 2nd spec of alphAB):
0037 >> grph_StD_AB
0038
0039 figure-3 is unchanged but top-panel on figure-4, now plot
0040 the evolution of global averaged S_ab.
0041 (should be flat if conserving)
0042 Note a problem of conservation at the very 1rst iteration,
0043 but later on the curve is flat. This problem is not seen
0044 after a restart.