[MITgcm-devel] changes in thermodynamics.F (tomorrow)

Jean-Michel Campin jmc at ocean.mit.edu
Mon Nov 18 13:28:12 EST 2013


Hi,

Tomorrow, I am going to check-in significant changes in thermodynamics.F
(+ few other routines, including new ones and removing few others) to 
move the k-loop inside the tracer loop (doing T, S, and all ptracers)
similar to what is currently done for the multidimension advection part.

I am not sure if this will speed up the code, but to be fair, there
is a second set of changes which involves calling some small subroutines
(ADAMS_BASHFORTH, FREESURF_RESCALE_G, TIMESTEP_TRACER) outside the k-loop
that could make it faster; will need to check later on if this the case.

With the little changes in gmredi_calc_uvflow.F & gmredi_calc_wflow.F
(define GM_CALC_TRANSP_FROM_VEL) that I made this morning, the forward
experiment results will not change (with -ieee). Some machine truncation
changes are expected for adjoint experiments; also I am not sure that all
the new TAF store directives are strictly necessary (might be worth to check),
or some could be missing when using some untested features.

Cheers,
Jean-Michel



More information about the MITgcm-devel mailing list