optim_m1qn3/testbed/model.F

C $Header: /u/gcmpack/MITgcm_contrib/mlosch/optim_m1qn3/testbed/model.F,v 1.2 2012/05/02 19:41:44 mlosch Exp $
C $Name:  $
      program model
      implicit none
      integer indic, n, nn
      parameter (nn=7)
      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,I4,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	C $Header: /u/gcmpack/MITgcm_contrib/mlosch/optim_m1qn3/testbed/model.F,v 1.2 2012/05/02 19:41:44 mlosch Exp $
2	C $Name: $
3	program model
4	implicit none
5	integer indic, n, nn
6	parameter (nn=7)
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,I4,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