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
      integer*8 rectmp, tmp1, tmp2
c----
c
      common /xxstore/ itape, ntape, isize
      integer itape, ntape, isize

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

      print *, 'pathei: in instore '

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*

      tmp1 = n
      tmp2 = isize
      rectmp = tmp1*tmp2
cph(
      print *, 'pathei rectmp ', rectmp, n, isize
cph)

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

      return

 900  continue
      close(itape)
      ifail = 1

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