[MITgcm-devel] vectorization of layers

[Martin Losch]

Hi Ryan,

did you get a chance to look at the changed layers_fluxcalc.F?
I also tried to further the vectorization and I did succeed technically, but at the cost of nearly infinite execution time, so even if my idea gives the correct results, there is no computational gain (although nicely vectorized). I feel that, what we have now (what I sent with my last email), is the best I can do right now. I am happy to check that in, if you approve of it.

I have a general question. Does it make sense to compute (and save) the "mean depth" of the layers? I would need that to compute, e.g. the overturning stream function averaged on isopycnals but plotted as a function of latitude and depth. Reconstructing the depths after the run does work, but it becomes a little noisy (as does the post-mortem computation of the overtuning on density surface, that's why your layers package is so useful).

Martin