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, rectmp
      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.

      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*

      rectmp = n * isize
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	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, rectmp
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	print *, 'pathei: in instore '
21
22	cc one record = two arrays of real*(isize)
23
24	open( itape
25	$ , file = 'OPWARMI'
26	$ , status = 'unknown'
27	$ , form = 'formatted'
28	$ , access = 'sequential'
29	$ )
30
31	read( itape, *, end=800, err=900 ) n, fc, gnorm0,
32	$ isize, m, jmin, jmax
33	cold = .false.
34
35	800 continue
36	close(itape)
37
38	print*
39	write(*,'(a,i10,a,i2)')
40	$ ' opening direct access with recl = ',n,' * ',isize
41	print*
42
43	rectmp = n * isize
44	cph(
45	print *, 'pathei rectmp ', rectmp, n, isize
46	cph)
47
48	open( ntape
49	$ , file = 'OPWARMD'
50	$ , status = 'unknown'
51	$ , form = 'unformatted'
52	$ , access = 'direct'
53	$ , recl = rectmp
54	$ )
55
56	return
57
58	900 continue
59	close(itape)
60	ifail = 1
61
62	return
63	end