optim_m1qn3/testbed/model.F

C $Header: /u/gcmpack/MITgcm_contrib/mlosch/optim_m1qn3/testbed/model.F,v 1.1 2012/05/02 12:27:42 mlosch Exp $
C $Name:  $
      program model
      implicit none
      integer indic, n, nn
      parameter (nn=20)
      double precision f,x(nn),g(nn)
      integer izs(1)
      real    rzs(1)
      double precision dzs(1)

      integer io, i
      character*(80) fmt

      io = 99
      open(io,file='model_in',status='old',form='unformatted')
      read(io) indic,n
      if ( n.gt.nn ) stop 'n > nn'
      read(io) (x(i),i=1,n)
      read(io) izs,rzs,dzs
      close(io)
      do i=1,n
       g(i)=0.
      enddo
      call mysimul(indic,n,f,x,g,izs,rzs,dzs)
      write(fmt,'(A,I2,A)') '(A,',n+1,'E12.4)'
      write(*,fmt) 'model finished f    = ', f, (g(i),i=1,n)
      io = 99
      open(io,file='model_out',status='unknown',form='unformatted')
      write(io) indic,n
      write(io) f
      write(io) (g(i),i=1,n)
      write(io) izs,rzs,dzs
      close(io)
      end
1	mlosch	1.2	C $Header: /u/gcmpack/MITgcm_contrib/mlosch/optim_m1qn3/testbed/model.F,v 1.1 2012/05/02 12:27:42 mlosch Exp $
2			C $Name: $
3	mlosch	1.1	program model
4			implicit none
5			integer indic, n, nn
6			parameter (nn=20)
7			double precision f,x(nn),g(nn)
8			integer izs(1)
9			real rzs(1)
10			double precision dzs(1)
11
12			integer io, i
13			character*(80) fmt
14
15			io = 99
16			open(io,file='model_in',status='old',form='unformatted')
17			read(io) indic,n
18			if ( n.gt.nn ) stop 'n > nn'
19			read(io) (x(i),i=1,n)
20			read(io) izs,rzs,dzs
21			close(io)
22			do i=1,n
23			g(i)=0.
24			enddo
25			call mysimul(indic,n,f,x,g,izs,rzs,dzs)
26			write(fmt,'(A,I2,A)') '(A,',n+1,'E12.4)'
27			write(*,fmt) 'model finished f = ', f, (g(i),i=1,n)
28			io = 99
29			open(io,file='model_out',status='unknown',form='unformatted')
30			write(io) indic,n
31			write(io) f
32			write(io) (g(i),i=1,n)
33			write(io) izs,rzs,dzs
34			close(io)
35			end