pkg/ctrl/optim_readparms.F

C $Header:  $
C $Name:  $

#include "CTRL_CPPOPTIONS.h"


      subroutine optim_readparms(
     I                            mythid
     &                          )

c     ==================================================================
c     SUBROUTINE optim_readparms
c     ==================================================================
c
c     o Initialise the optimization part of the ECCO release.
c
c     started: Christian Eckert eckert@mit.edu 03-Mar-2000
c
c     changed:
c
c     ==================================================================
c     SUBROUTINE optim_readparms
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"

#include "optim.h"

c     == routine arguments ==

      integer mythid

c     == local variables ==

      integer errio
      integer il
      integer iUnit

      character*(max_len_prec) record
      character*(max_len_mbuf) msgbuf

c     == external ==

      integer  ilnblnk
      external ilnblnk

c     == end of interface ==

c--   Optimization parameters.
      namelist /optim/
     &                 optimcycle, nvars,
     &                 nondimcontrol,
     &                 numiter, nfunc, fmin, iprint,
     &                 epsf, epsx, epsg,
     &                 nupdate, eps

      _BEGIN_MASTER( mythid )

c--     Set default values.
        optimcycle    =   0
        nvars         =  10
        nondimcontrol =  .false.
        numiter       =   0
        nfunc         =   1
        fmin          =   0.0 _d 0
        iprint        =  10
        epsx          =   1.e-6
        epsg          =   1.e-6
        eps           =  -1.e-6
        nupdate       =   0

        WRITE(msgBuf,'(A)') 'OPTIM_READPARMS: opening data.optim'
        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                    SQUEEZE_RIGHT , 1)

        CALL OPEN_COPY_DATA_FILE(
     I                          'data.optim', 'OPTIM_READPARMS',
     O                          iUnit,
     I                          myThid )

        READ(unit = iUnit, nml = optim)

        WRITE(msgBuf,'(A)')
     &       'OPTIM_READPARMS: finished reading data.optim'
        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
     &                  SQUEEZE_RIGHT , 1)

        CLOSE( iUnit )

      _END_MASTER( mythid )

      _BARRIER

      return
      end

1	jmc	1.3	C $Header: $
2			C $Name: $
3	heimbach	1.1
4			#include "CTRL_CPPOPTIONS.h"
5
6
7			subroutine optim_readparms(
8			I mythid
9			& )
10
11			c ==================================================================
12			c SUBROUTINE optim_readparms
13			c ==================================================================
14			c
15			c o Initialise the optimization part of the ECCO release.
16			c
17			c started: Christian Eckert eckert@mit.edu 03-Mar-2000
18			c
19	jmc	1.3	c changed:
20	heimbach	1.1	c
21			c ==================================================================
22			c SUBROUTINE optim_readparms
23			c ==================================================================
24
25			implicit none
26
27			c == global variables ==
28
29			#include "EEPARAMS.h"
30
31			#include "optim.h"
32
33			c == routine arguments ==
34
35	jmc	1.3	integer mythid
36	heimbach	1.1
37			c == local variables ==
38
39			integer errio
40			integer il
41	heimbach	1.2	integer iUnit
42	heimbach	1.1
43			character*(max_len_prec) record
44	heimbach	1.2	character*(max_len_mbuf) msgbuf
45	heimbach	1.1
46			c == external ==
47
48			integer ilnblnk
49			external ilnblnk
50
51			c == end of interface ==
52
53			c-- Optimization parameters.
54			namelist /optim/
55			& optimcycle, nvars,
56			& nondimcontrol,
57			& numiter, nfunc, fmin, iprint,
58			& epsf, epsx, epsg,
59			& nupdate, eps
60
61			_BEGIN_MASTER( mythid )
62
63			c-- Set default values.
64			optimcycle = 0
65			nvars = 10
66			nondimcontrol = .false.
67			numiter = 0
68			nfunc = 1
69			fmin = 0.0 _d 0
70			iprint = 10
71			epsx = 1.e-6
72			epsg = 1.e-6
73			eps = -1.e-6
74			nupdate = 0
75
76	heimbach	1.2	WRITE(msgBuf,'(A)') 'OPTIM_READPARMS: opening data.optim'
77			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
78			& SQUEEZE_RIGHT , 1)
79
80			CALL OPEN_COPY_DATA_FILE(
81			I 'data.optim', 'OPTIM_READPARMS',
82			O iUnit,
83			I myThid )
84
85			READ(unit = iUnit, nml = optim)
86
87	jmc	1.3	WRITE(msgBuf,'(A)')
88	heimbach	1.2	& 'OPTIM_READPARMS: finished reading data.optim'
89			CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
90			& SQUEEZE_RIGHT , 1)
91	heimbach	1.1
92	heimbach	1.2	CLOSE( iUnit )
93	heimbach	1.1
94			_END_MASTER( mythid )
95
96			_BARRIER
97
98			return
99			end
100