[MITgcm-devel] offline_exf_seaice

[Martin Losch]

Hi Jean-Michel,

I have a specific testreport question. I want to create a test for the ITD code with advection etc. Ideally something like the offline_exf_seaice, but rather than coming up with a new directory, my plan is to modify the code and offline_exf_seaice, that we can use it for both testing the code with SEAICE_multDim = 1 and larger, should be possible.

The first hurdle (after removing code that resets useThSice in seaice_readparms, defining SEAICE_ITD and setting nITD=1 in SEAICE_SIZE.h) is that testreport for thermo gives only 5 and 4 digits of agreement for mean and s.d. of U, but really it happens in seaice_area_mean, seaice_area_sd, seaice_area_del2 only once in timestep 30. Before and after that everything is perfect except for exf_hflux_del2 and exf_sflux_del2 in timestep 0. All other experiments look perfect.

Two questions: 
1. Am I right that the table header that testreport generates cannot be trusted entirely and that I should really have a look at tr_checklist instead? 
2. Are these differences described above typical? I mean I find it hard to believe that that results are precise to 16 digits except for one variable at one timestep in the middle of the test (not the beginning or the end). Should I just ignore that and update the output file?

Martin