[MITgcm-devel] thsice_get_exf

[Patrick Heimbach]

Hi Martin,

that sounds good in principle, and if done 
properly should be fully OK w.r.t. adjoint. 
Since this routine has become a crucial piece
for the OSE, it would be good to test it carefully
before making it the standard.

Cheers
-p.



On Thu, 2006-06-08 at 12:00, Martin Losch wrote:
> Hi Jean-Michel/Patrick,
> 
> I had another look at thsice_get_exf and compared to  
> bulkf_formula_lay it is a mess, even after Jean-Michel has cleaned up  
> all of the obvious mistakes (such as factors of 1000).
> 
> I would like to have a routine that is called both from  
> thsice_solve4temp and exf_bulkformulae (so a pendant to  
> bulkf_formula_lay). That would make it much easier to switch between  
> different bulk formulae, Large and Pond, vs. Large and Yeager, or the  
> NCEP bulkformulae, whatever they are. I think it is probably a fact  
> overlooked by many "uninitiated oceanographers" (such as myself),  
> that the bulkformulae for computing fluxes have to match the ones  
> used in the atmospheric models that produced the atmospheric fields.
> 
> The problem is that exf_bulkformulae does not need all of the stuff  
> that's needed in thsice_get_exf.F
> Also, for the adjoint it is probably problematic to call subroutines  
> for every gridpoint, so maybe it is possible to have code snippet  
> that is included into these routines from a "header"-type file, such  
> as #include "exf_bulkformulae_core.h"
> or
> #include "exf_bulkformulae_ncep.h"
> The extra stuff that is required for thsice_get_exf.F can be computed  
> outside of this snippet.
> 
> What do you think about this?
> 
> Martin
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel
-- 
--------------------------------------------------------
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/                        USA