[MITgcm-devel] problem with argument list of adthe_main_loop

[Martin Losch]

Hi Patrick,


sorry for the inaccurate description of the problem. I was just looking for a general explanation/direction. Both you have provided, thanks. 


I found that with ecco, you need at least one "time-dependent" costfunction contribution (e.g. ALLOW_THETA_COST_CONTRIBUTION) for the "mytime" to appear in the argument list of adthe_main_loop. This corresponds to your description.


In the end I/we want to use the ecco package in a "standard" form so that this problem will not affect use, but right now I want to compute the sensitivity to meridional overturning in a cs32 configuration with seaice. We (Takasumi and I) tried with and without the ecco package.
1. without ecco: I put my cost function computations into cost_test (just as in bottom_ctrl_5x5) and get many recomputations, most likely due to the seaice (at least they all go away when I turn off seaice). Since I did not have time then (and do not have it now either, unfortunately) to dig into that, I recommended to Takasumi to try this:
2. with the ecco package, I put my cost function into cost_usercost_all.F and I do not get recomputations (however so far I also did not get a cost contributions from cost_usercost, I have just started with that, so give me some more time), but we have this "mytime, mythid" problem. I can now circumvent that by having fake theta-cost functions (with mult_theta=0.)


So my question is, according to your experience, which route will be faster to the goal (sensitivity of the overturning with respect to surface forcing), the ecco or the non-ecco route. I have the feeling the non-ecco route, but only if you point me to where I have insert store directives to fix the recompuations. If you find the time (today), I'd love to show you my configuration.


Martin


----- Original Message -----
From: Patrick Heimbach <heimbach at MIT.EDU>
Date: Tuesday, April 27, 2010 5:07
Subject: Re: [MITgcm-devel] problem with argument list of adthe_main_loop
To: MITgcm-devel at mitgcm.org

> 
> Hi Martin,
> 
> I haven't looked into this in a long time.
> As you see, the "broad" behavior is that
> mycurrenttime is kept for ECCO, but not otherwise.
> 
> If I remember correctly, I tried a number of ways to "recompute" 
> mycurrenttimewithin the_main_loop and make it "self-sufficient" 
> (e.g. use starttime), but didn't help much.
> I think I also tried the other way round, i.e. "force" TAF to 
> keep mycurrenttime always,
> but didn't help much either.
> My hunch is it has to do with the use (or non-use) of the 
> cost_averages... routines.
> 
> But would have to see your setup in detail ;)
> in particular, what's the difference between the two cases.
> 
> I had not pursued this further, since it seems stable for most 
> applications,and can be easily changed in case of unexpected behavior.
> 
> Cheers
> -p.
> 
> On Apr 26, 2010, at 9:08 AM, Martin Losch wrote:
> 
> >Hi there (Patrick in particular),
> >
> >I found this in the archives (the search functions seems to 
> work better):
> >http://forge.csail.mit.edu/pipermail/mitgcm-devel/2008-
> June/003425.html>and I seem to have a similar problem in a 
> particular configuration:
> >in ad_taf_output.f I find
> >subroutine adthe_main_loop(mythid)
> >but it's called as adthe_main_loop(mycurrenttime, mythid) from 
> the_model_main.F>
> >I have other configurations where taf properly produces 
> "subroutine adthe_main_loop(mytime,mythid)", but I cannot see 
> why. Both configurations use the ecco package (so ALLOW_ECCO is 
> defined, that's why the_model_main requires mytime and mythid) 
> and pretty similar ECCO_CPPOPTIONS.h.
> >
> >Can you point me to the critical parameters that let taf 
> produce the argument list with (mytime, mythid)? Is the type of 
> cost function?
> >
> >Martin
> >
> >
> >
> >_______________________________________________
> >MITgcm-devel mailing list
> >MITgcm-devel at mitgcm.org
> >http://mitgcm.org/mailman/listinfo/mitgcm-devel
> 
> ---
> Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
> MIT | EAPS 54-1518 | 77 Massachusetts Ave | Cambridge MA 02139 USA
> FON +1-617-253-5259 | FAX +1-617-253-4464 | SKYPE patrick.heimbach
> 
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel

Martin Losch
Alfred Wegener Institute 
Postfach 120161, 27515 Bremerhaven, Germany; 
Tel./Fax: ++49(0471)4831-1872/1797




-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20100427/87559003/attachment-0001.htm>