[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