MITgcm/lsopt/instore.F

      subroutine instore( n, fc, gnorm0, is, m, jmin, jmax,
     &                    cold, ifail )
c
c         arguments
c
      integer n, is, m, jmin, jmax, ifail
      double precision    fc, gnorm0
      logical cold
c----
c
      common /xxstore/ itape, ntape, isize
      integer itape, ntape, isize

      ifail = 0
      itape = 91
      ntape = 92
      isize = is
      cold  = .true.

cc one record = two arrays of real*(isize)

      open(  itape
     $    , file   = 'OPWARMI'
     $    , status = 'unknown'
     $    , form   = 'formatted'
     $    , access = 'sequential'
     $    )

      read( itape, *, end=800, err=900 ) n, fc, gnorm0,
     $                                   isize, m, jmin, jmax
      cold  = .false.

 800  continue
      close(itape)

      print*
      write(*,'(a,i10,a,i2)')
     $        '  opening direct access with recl = ',n,' * ',isize
      print*

      open( ntape
     $    , file   = 'OPWARMD'
     $    , status = 'unknown'
     $    , form   = 'unformatted'
     $    , access = 'direct'
     $    , recl   = n * isize
     $    )

      return

 900  continue
      close(itape)
      ifail = 1

      return
      end
1	heimbach	1.1.2.1	subroutine instore( n, fc, gnorm0, is, m, jmin, jmax,
2			& cold, ifail )
3			c
4			c arguments
5			c
6			integer n, is, m, jmin, jmax, ifail
7			double precision fc, gnorm0
8			logical cold
9			c----
10			c
11			common /xxstore/ itape, ntape, isize
12			integer itape, ntape, isize
13
14			ifail = 0
15			itape = 91
16			ntape = 92
17			isize = is
18			cold = .true.
19
20			cc one record = two arrays of real*(isize)
21
22			open( itape
23			$ , file = 'OPWARMI'
24			$ , status = 'unknown'
25			$ , form = 'formatted'
26			$ , access = 'sequential'
27			$ )
28
29			read( itape, *, end=800, err=900 ) n, fc, gnorm0,
30			$ isize, m, jmin, jmax
31			cold = .false.
32
33			800 continue
34			close(itape)
35
36			print*
37			write(*,'(a,i10,a,i2)')
38			$ ' opening direct access with recl = ',n,' * ',isize
39			print*
40
41			open( ntape
42			$ , file = 'OPWARMD'
43			$ , status = 'unknown'
44			$ , form = 'unformatted'
45			$ , access = 'direct'
46			$ , recl = n * isize
47			$ )
48
49			return
50
51			900 continue
52			close(itape)
53			ifail = 1
54
55			return
56			end