[MITgcm-support] tutorial_global_oce_optim optimisation failed
Andrew McRae
andrew.mcrae at physics.ox.ac.uk
Fri Apr 27 11:25:28 EDT 2018
Just separating this from the other thread
<http://mailman.mitgcm.org/pipermail/mitgcm-support/2018-April/011521.html>,
I got the bundled MITgcm optim routine built (having made these changes
<https://github.com/MITgcm/MITgcm/compare/master...dorugeber:optim_fix>,
based on this
<http://mailman.mitgcm.org/pipermail/mitgcm-support/2010-September/006825.html>
thread from 2010 and this
<http://mailman.mitgcm.org/pipermail/mitgcm-support/2016-July/010527.html>
one from 2016).
I use OpenAD to create the adjoint.
My steps are:
1) in the build directory, run ../../../tools/genmake2 -oad
-mods=../code_oad
2) run make depend and make adAll
3) copy input_oad/ into a new folder scratch/
4) within scratch/, run ./prepare_run
5) copy mitgcmuv_ad from build/ into scratch/, copy optim.x into
scratch/OPTIM/
6) run ./mitgcmuv_ad
7) in scratch/OPTIM, create symlinks to ../data.optim and ../data.ctrl
8) copy the files ecco_cost_MIT_CE_000.opt0000 and
ecco_ctrl_MIT_CE_000.opt0000 into the OPTIM subdirectory
9) run ./optim.x within the subdirectory
The full output is attached, but I assume the optimisation failed since the
last lines are
optimization stopped because :
ifail = 4 the search direction is not a descent one
Any ideas? (I guess this isn't something that is tested in the daily
builds?)
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.
Andrew
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20180427/699c3f5a/attachment-0001.html>
-------------- next part --------------
==================================================
Large Scale Optimization with off-line capability.
==================================================
Version 2.1.0
OPTIM_NUMBMOD: Control options have been read.
OPTIM_NUMBMOD: Minimization options have been read.
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading control vector
for optimization cycle: 0
opened file ecco_ctrl_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
OPTIM_NUMBMOD: Iteration number = 0
number of control variables = 2315
Data will be read from the following files:
pathei: vor optim_sub
OPTIMUM: Calling lsopt for iteration: 0
OPTIMUM: with nn, REAL_BYTE = 2315 4
===============================================
=== O P T I M I Z A T I O N ===
===============================================
number of control variables....... 2315
precision on x.................... 0.10E-05
precision on g.................... 0.10E-05
expected optimal function value... 0.57E+01
maximal number of iterations...... 1
maximal number of simulations..... 3
information level................. 10
number of updates................. 8
size of used memory............... 6945
pathei: vor instore
pathei: in instore
opening direct access with recl = 2315 * 4
pathei rectmp 9260 2315 4
pathei-lsopt: cold start
pathei-lsopt vor simul 2315
pathei-lsopt xx(1), gg(1) 0.0000000000000000 0.0000000000000000
pathei-lsopt in simul
pathei-lsopt vor optim_readdata
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading control vector
for optimization cycle: 0
opened file ecco_ctrl_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading cost function and
gradient of cost function
for optimization cycle: 0
opened file ecco_cost_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
leaving simul with nn, objf = 2315 6.2002322818232924
leaving simul with xx, adxx = 0.0000000000000000 -6.7879882408306003E-005
pathei: nach simul: nn, ff = 2315 6.2002322818232924
pathei: nach simul: xx(1), gg(1) = 0.0000000000000000 -6.7879882408306003E-005
pathei: vor dostore
pathei: vor lswri
pathei: vor gnorm0
pathei: gnorm0 = 1.2730926712322933E-002
pathei-lsopt: cold; first call simul: ff = 6.2002322818232924
Itn Step Nfun Objective Norm G Norm X Norm (X(k-1)-X(k))
0 1.0E+00 0 6.2002E+00 1.3E-02 0.0E+00 0.0E+00
pathei-lsopt: ++++++++++++++++++++++++
pathei-lsopt: entering iiter = 1
pathei-lsopt: using fold = fmin = 5.7400000000000002
pathei-lsopt: first estimate of dd via fold - ff
pathei-lsopt: 1 simul.
pathei-lsopt in simul
pathei-lsopt vor optim_readdata
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading control vector
for optimization cycle: 0
opened file ecco_ctrl_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading cost function and
gradient of cost function
for optimization cycle: 0
opened file ecco_cost_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
leaving simul with nn, objf = 2315 6.2002322818232924
leaving simul with xx, adxx = 0.0000000000000000 -6.7879882408306003E-005
pathei-lsopt: interpolation: td, tg, tact, ifail = 1.0000000000000000 0.0000000000000000 0.21132486482402846 0
pathei-lsopt: 2 simul.
pathei-lsopt in simul
pathei-lsopt vor optim_readdata
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading control vector
for optimization cycle: 0
opened file ecco_ctrl_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading cost function and
gradient of cost function
for optimization cycle: 0
opened file ecco_cost_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
leaving simul with nn, objf = 2315 6.2002322818232924
leaving simul with xx, adxx = 0.0000000000000000 -6.7879882408306003E-005
pathei-lsopt: interpolation: td, tg, tact, ifail = 0.21132486482402846 0.0000000000000000 0.20921161622302295 0
pathei-lsopt: 3 simul.
pathei-lsopt in simul
pathei-lsopt vor optim_readdata
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading control vector
for optimization cycle: 0
opened file ecco_ctrl_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
pathei-lsopt in optim_readdata
OPTIM_READDATA: Reading cost function and
gradient of cost function
for optimization cycle: 0
opened file ecco_cost_MIT_CE_000.opt0000
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: filenopt 0
pathei: fileff 6.2002322818232924
pathei: fileiG 1
pathei: filejG 1
pathei: filensx 2
pathei: filensy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
leaving simul with nn, objf = 2315 6.2002322818232924
leaving simul with xx, adxx = 0.0000000000000000 -6.7879882408306003E-005
pathei-lsopt: interpolation: td, tg, tact, ifail = 0.20921161622302295 0.0000000000000000 0.19875103564568400 0
pathei-lsopt: back from optline; ifail = 6
pathei-lsopt: dostore 1,2 at iiter 1
1 2.0E-01 4 6.2002E+00 1.3E-02 0.0E+00 1.4E+01
pathei-lsopt: ys < 0; ifail = 4
pathei-lsopt: using fold = fmin = 5.7400000000000002
pathei-lsopt: first estimate of dd via fold - ff
pathei: writing 2315 sized control to fileecco_ctrl
OPTIM_WRITEDATA: Writing new control vector to file(s)
for optimization cycle: 1
pathei: nvartype 1
pathei: nvarlength 2315
pathei: yctrlid MIT_CE_000
pathei: nopt 1
pathei: ff -9999.0000000000000
pathei: iG 1
pathei: jG 1
pathei: nsx 2
pathei: nsy 2
pathei: nWetcGlobal 2315 2315 2254 2215 2178 2142 2114 2076 2048 1999 1948 1850 1655 1372 828
pathei: nWetsGlobal 2149 2149 2102 2058 2027 1987 1959 1918 1887 1831 1771 1653 1404 1118 623
pathei: nWetwGlobal 2206 2206 2146 2104 2070 2029 2004 1959 1925 1869 1808 1705 1458 1164 671
pathei: nWetvGlobal 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarindex -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 124 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1 -1
pathei: ncvarrecs 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarxmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 45 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarymax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 20 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvarnrmax 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
pathei: ncvargrd ???????????????????????c????????????????????????????????????
in owd: icvoffset 2315
number of iterations.............. 1
number of simultations............ 4
relative precision on g........... 0.10E+01
cost function............... 0.62002323E+01
norm of x................... 0.00000000E+00
norm of g................... 0.12730927E-01
optimization stopped because :
ifail = 4 the search direction is not a descent one
======================================
Large Scale Optimization run finished.
======================================
More information about the MITgcm-support
mailing list