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	c $Header: /u/gcmpack/MITgcm_contrib/mlosch/optim_m1qn3/optim_main.F,v 1.2 2012/04/27 09:48:53 mlosch Exp $
2	c $Name: $
3	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	nn=1
69	ff=-9999.
70	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