[MITgcm-devel] PLOT_FIELD_XY in seaice_do_diags.F

[Jean-Michel Campin]

Hi Martin,

I was going to simplify the MDS part of seaice_do_diags.F
(by the way, can I rename it "seaice_output.F, to clarify that
it has nothing to do with the diagnostics pkg ?) removing
the temp arr copy (since we have pkg/rw routine that
allow to write 1 level among several).

Do you really use the PLOT_FIELD_XY output ? Personally, most of the
time, when I have a problem, it's always to track a relative diff < 10^-3
and the alphabet is not helpful. The only reason why I could look
at the PLOT_FIELD_XY output would be for the overlap, but in
seaice_do_diags.F they are reset by EXCH call just before the plot.
I can either remove all the PLOT_FIELD_XY from seaice_do_diags.F
(so that I don't need the temp array anymore, except for MNC), or
keep them (much more coding), but then I will not check that the 
changes are right for those PLOT_FIELD output.

other unrelated question:
Why all the 2D seaice time-ave fields are stored in 5 indices arrays ?
(FUtave, EmPmRtave, QNETtave ...)

Cheers,
Jean-Michel