<html><head><meta http-equiv="Content-Type" content="text/html; charset=utf-8"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hello<div class=""><br class=""></div><div class="">Its been awhile, but I am pretty sure that is the normal output. It says “fail", but it did give you a new and ecco_ctrl_MIT_CE_000.opt0001 (correct?) and if you unpack and run likely the cost will descend.</div><div class=""><br class=""></div><div class="">I think it worked correctly. lsopt/optim are just confusing…but I think its working. I think all is good!</div><div class=""><br class=""></div><div class="">Matt</div><div class=""><br class=""></div><div class=""><br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On Apr 27, 2018, at 8:25 AM, Andrew McRae <<a href="mailto:andrew.mcrae@physics.ox.ac.uk" class="">andrew.mcrae@physics.ox.ac.uk</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class=""><div class=""><div class="">Just separating this from the <a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/2018-April/011521.html" class="">other thread</a>, I got the bundled MITgcm optim routine built (having made <a href="https://github.com/MITgcm/MITgcm/compare/master...dorugeber:optim_fix" class="">these changes</a>, based on <a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/2010-September/006825.html" class="">this</a> thread from 2010 and <a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/2016-July/010527.html" class="">this</a> one from 2016).<br class=""><br class=""></div>I use OpenAD to create the adjoint.<br class=""><br class=""></div>My steps are:<br class=""></div><div class="">1) in the build directory, run ../../../tools/genmake2 -oad -mods=../code_oad<br class=""></div><div class="">2) run make depend and make adAll<br class="">3) copy input_oad/ into a new folder scratch/<br class=""></div><div class="">4) within scratch/, run ./prepare_run<br class=""></div><div class="">5) copy mitgcmuv_ad from build/ into scratch/, copy optim.x into scratch/OPTIM/<br class=""></div><div class="">6) run ./mitgcmuv_ad<br class=""></div><div class="">7) in scratch/OPTIM, create symlinks to ../data.optim and ../data.ctrl<br class=""></div><div class="">8) copy the files ecco_cost_MIT_CE_000.opt0000 and ecco_ctrl_MIT_CE_000.opt0000 into the OPTIM subdirectory<br class=""></div><div class="">9) run ./optim.x within the subdirectory<br class=""><br class=""></div><div class="">The full output is attached, but I assume the optimisation failed since the last lines are<br class=""></div><div class=""><div class=""><div class=""><br class=""> optimization stopped because :<br class=""> ifail = 4 the search direction is not a descent one<br class=""><br class=""></div><div class="">Any ideas? (I guess this isn't something that is tested in the daily builds?)<br class=""><br class=""></div><div class="">In the meantime, I'll try the m1qn3 routine as in the other thread, which should help distinguish between a problem with the optimisation routine or the gradient generated by mitgcmuv_ad.<br class=""><br class=""></div><div class="">Andrew<br class=""></div></div></div></div>
<span id="cid:9D33FD64-1548-462F-99BA-72879C96702C@ucsd.edu"><out.txt></span>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></div></blockquote></div><br class=""></div></body></html>