[Mitgcm-support] Re: sea-ice model adjoint

[mitgcm-support at dev.mitgcm.org]

Hi Dimitris,

it seems to me that it's OK to move EXTERNAL_FORCING_SURF
to the place you mentioned.
However, there are two issues to watch:

1. in conjunction with KPP the surfacetendency? fields
    need to be stored (cf. CADJ STORE block in thermodynamics).
2. When using open boundaries, and SST/SSS relaxation
    it might be that e.g. the THETA-SST terms contain
    wrong values at the boundaries
    (not sure, I'd have to check).

Point 2. may, in fact, give a NOT OK to your question.

Patrick



Dimitris Menemenlis wrote:

Hi Patrick, I am working on implementing sea-ice model on top of
ecco_c44_e22.  Quick question: does EXTERNAL_FORCING_SURF need
to remain in THERMODYNAMICS now that THE_CORRECTION_STEP has been
moved out of there?  Would it be correct to instead call
EXTERNAL_FORCING_SURF after EXF_GETFORCING inside ECCO_THE_MAIN_LOOP
(and changing bi, bj behavior accordingly)?

Cheers, D.




-- 
_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/
Patrick Heimbach     Massachusetts Institute of Technology
FON: +1/617/253-5259                    EAPS, Room 54-1518
FAX: +1/617/253-4464               77 Massachusetts Avenue
mailto:heimbach at mit.edu                 Cambridge MA 02139
http://www.mit.edu/~heimbach/                       U.S.A.


-- 
_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/_/
Patrick Heimbach     Massachusetts Institute of Technology
FON: +1/617/253-5259                    EAPS, Room 54-1518
FAX: +1/617/253-4464               77 Massachusetts Avenue
mailto:heimbach at mit.edu                 Cambridge MA 02139
http://www.mit.edu/~heimbach/                       U.S.A.