MITgcm/optim/simul.F


      subroutine simul(
     I                  indic,
     I                  nn,
     I                  xx,
     O                  objf,
     O                  adxx
     &                )

c     ==================================================================
c     SUBROUTINE simul
c     ==================================================================
c
c     o This routine is called by the large-scale optimization lsopt.
c
c       Input  : indic - Parameter (not used here).
c                nn    - Number of control variables.
c                xx    - Array  of control variables.
c
c       Output : objf  - Value of objective function.
c                adxx  - Gradients of objective function with respect
c                        to the control variables.
c
c
c     started: Christian Eckert eckert@mit.edu 15-Feb-2000
c
c     changed: Christian Eckert eckert@mit.edu 10-Mar-2000
c
c              - Added ECCO layout.
c
c     changed:  Patrick Heimbach heimbach@mit.edu 19-Jun-2000
c               - finished, revised and debugged
c
c     ==================================================================
c     SUBROUTINE simul
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "ctrl.h"

c     == routine arguments ==

      integer indic
      integer nn
      _RL   xx(nn)
      _RL   objf
      _RL   adxx(nn)

c     == local variables ==

      integer i
      _RL   adobjf

      logical lheaderonly

c     == end of interface ==

      print *, 'pathei-lsopt in simul'

c--   Call the combined modified forward model and the adjoint model.
      do i = 1,nn
        adxx(i) = 0.
      enddo

      adobjf = 1.
c
      lheaderonly = .false.

      print *, 'pathei-lsopt vor optim_readdata'

      call optim_readdata( nn, ctrlname, lheaderonly, objf,   xx )
      call optim_readdata( nn, costname, lheaderonly, objf, adxx )

      indic = indic

cph(
      print *, ' leaving simul with nn, objf = ', nn, objf
      print *, ' leaving simul with xx, adxx = ', xx(1), adxx(1)
      do i=1,nn
c         if (xx(i).EQ.'NaN') then
c            print *, 'pathei - out: i = ', i
c         end if
      end do
cph)

      return
      end
1
2	subroutine simul(
3	I indic,
4	I nn,
5	I xx,
6	O objf,
7	O adxx
8	& )
9
10	c ==================================================================
11	c SUBROUTINE simul
12	c ==================================================================
13	c
14	c o This routine is called by the large-scale optimization lsopt.
15	c
16	c Input : indic - Parameter (not used here).
17	c nn - Number of control variables.
18	c xx - Array of control variables.
19	c
20	c Output : objf - Value of objective function.
21	c adxx - Gradients of objective function with respect
22	c to the control variables.
23	c
24	c
25	c started: Christian Eckert eckert@mit.edu 15-Feb-2000
26	c
27	c changed: Christian Eckert eckert@mit.edu 10-Mar-2000
28	c
29	c - Added ECCO layout.
30	c
31	c changed: Patrick Heimbach heimbach@mit.edu 19-Jun-2000
32	c - finished, revised and debugged
33	c
34	c ==================================================================
35	c SUBROUTINE simul
36	c ==================================================================
37
38	implicit none
39
40	c == global variables ==
41
42	#include "EEPARAMS.h"
43	#include "SIZE.h"
44	#include "ctrl.h"
45
46	c == routine arguments ==
47
48	integer indic
49	integer nn
50	_RL xx(nn)
51	_RL objf
52	_RL adxx(nn)
53
54	c == local variables ==
55
56	integer i
57	_RL adobjf
58
59	logical lheaderonly
60
61	c == end of interface ==
62
63	print *, 'pathei-lsopt in simul'
64
65	c-- Call the combined modified forward model and the adjoint model.
66	do i = 1,nn
67	adxx(i) = 0.
68	enddo
69
70	adobjf = 1.
71	c
72	lheaderonly = .false.
73
74	print *, 'pathei-lsopt vor optim_readdata'
75
76	call optim_readdata( nn, ctrlname, lheaderonly, objf, xx )
77	call optim_readdata( nn, costname, lheaderonly, objf, adxx )
78
79	indic = indic
80
81	cph(
82	print *, ' leaving simul with nn, objf = ', nn, objf
83	print *, ' leaving simul with xx, adxx = ', xx(1), adxx(1)
84	do i=1,nn
85	c if (xx(i).EQ.'NaN') then
86	c print *, 'pathei - out: i = ', i
87	c end if
88	end do
89	cph)
90
91	return
92	end