mlosch/optim_m1qn3/optim_main.F

c $Header: /u/gcmpack/MITgcm_contrib/mlosch/optim_m1qn3/optim_main.F,v 1.2 2012/04/27 09:48:53 mlosch Exp $
c $Name:  $
c     ==================================================================
c
c     optim2: Routines for doing an off-line optimization with m1qn3
c             after the ECCO forward and adjoint model have been run.
c
c     optim_main - driver routine
c     - optim_readparms   :: read namelist files and number of control 
c                            parameters
c     - optim_sub         :: driver for m1qn3
c       + optim_readdata  :: read control vector or gradient and cost 
c                           function
c       + optim_writedata :: write control vector
c       + optim_store_m1qn3 :: store/restore the state of the optimizer
c       + m1qn3_offline     :: optimization routine
c         > ddot            :: LAPACK routine, required by m1qn3_offline
c
c     Documentation:
c
c     m1qn3_offline is based on m1qn3 v3.3, downloaded here:
c     https://who.rocq.inria.fr/Jean-Charles.Gilbert/modulopt/
c             optimization-routines/m1qn3/m1qn3.html
c     To make work offline, the "save" statement have been replace
c     by common blocks, by which variables are passed to
c     optim_store_m1qn3 where they are stored on disk. In order to
c     reduce the number of variables to be saved, some have been move
c     down to where they are first used, e.g. barmul, tesf, etc.
c     
c     The driver routines optime_*.F are adapted from the original
c     optim-pkg of the MITgcm.
c
c     Martin.Losch@awi.de 24-Apr-2012
c
c     ==================================================================

      program optim_main

c     ==================================================================
c     program optim_main
c     ==================================================================
c
c     o Driver routine for the ECCO optimization package.
c     ==================================================================

      implicit none

c     == global variables ==

c     == routine arguments ==

c     == local variables ==

      integer nn
      _RL     ff

c     == end of interface ==

c--   Headline.
      print*
      print*,' =================================================='
      print*,' Large Scale Optimization with off-line capability.'
      print*,' =================================================='
      print*

c--   Get the number of control variables and other namelist 
c     parameters 
      nn=1
      ff=-9999.
      call optim_readparms( nn, ff )

c--   Call the subroutine that will set up the optimization routine
      call optim_sub( nn, ff )

c--   Succesful termination.
      print*
      print*,' ======================================'
      print*,' Large Scale Optimization run finished.'
      print*,' ======================================'
      print*
      end
1	mlosch	1.3	c $Header: /u/gcmpack/MITgcm_contrib/mlosch/optim_m1qn3/optim_main.F,v 1.2 2012/04/27 09:48:53 mlosch Exp $
2	mlosch	1.2	c $Name: $
3	mlosch	1.1	c ==================================================================
4			c
5			c optim2: Routines for doing an off-line optimization with m1qn3
6			c after the ECCO forward and adjoint model have been run.
7			c
8			c optim_main - driver routine
9			c - optim_readparms :: read namelist files and number of control
10			c parameters
11			c - optim_sub :: driver for m1qn3
12			c + optim_readdata :: read control vector or gradient and cost
13			c function
14			c + optim_writedata :: write control vector
15			c + optim_store_m1qn3 :: store/restore the state of the optimizer
16			c + m1qn3_offline :: optimization routine
17			c > ddot :: LAPACK routine, required by m1qn3_offline
18			c
19			c Documentation:
20			c
21			c m1qn3_offline is based on m1qn3 v3.3, downloaded here:
22			c https://who.rocq.inria.fr/Jean-Charles.Gilbert/modulopt/
23			c optimization-routines/m1qn3/m1qn3.html
24			c To make work offline, the "save" statement have been replace
25			c by common blocks, by which variables are passed to
26			c optim_store_m1qn3 where they are stored on disk. In order to
27			c reduce the number of variables to be saved, some have been move
28			c down to where they are first used, e.g. barmul, tesf, etc.
29			c
30			c The driver routines optime_*.F are adapted from the original
31			c optim-pkg of the MITgcm.
32			c
33			c Martin.Losch@awi.de 24-Apr-2012
34			c
35			c ==================================================================
36
37			program optim_main
38
39			c ==================================================================
40			c program optim_main
41			c ==================================================================
42			c
43			c o Driver routine for the ECCO optimization package.
44			c ==================================================================
45
46			implicit none
47
48			c == global variables ==
49
50			c == routine arguments ==
51
52			c == local variables ==
53
54			integer nn
55			_RL ff
56
57			c == end of interface ==
58
59			c-- Headline.
60			print*
61			print*,' =================================================='
62			print*,' Large Scale Optimization with off-line capability.'
63			print*,' =================================================='
64			print*
65
66			c-- Get the number of control variables and other namelist
67			c parameters
68	mlosch	1.2	nn=1
69			ff=-9999.
70	mlosch	1.1	call optim_readparms( nn, ff )
71
72			c-- Call the subroutine that will set up the optimization routine
73			call optim_sub( nn, ff )
74
75			c-- Succesful termination.
76			print*
77			print*,' ======================================'
78			print*,' Large Scale Optimization run finished.'
79			print*,' ======================================'
80			print*
81			end