MITgcm/lsopt/instore.F

4cee17c1be Patr*0001       subroutine instore( n, fc, gnorm0, is, m, jmin, jmax,
                0002      &                    cold, ifail )
                0003 c
                0004 c         arguments
                0005 c
f7e8e52ac8 Patr*0006       integer n, is, m, jmin, jmax, ifail
4cee17c1be Patr*0007       double precision    fc, gnorm0
                0008       logical cold
f7e8e52ac8 Patr*0009       integer*8 rectmp, tmp1, tmp2
4cee17c1be Patr*0010 c----
                0011 c
                0012       common /xxstore/ itape, ntape, isize
                0013       integer itape, ntape, isize
                0014
                0015       ifail = 0
                0016       itape = 91
                0017       ntape = 92
                0018       isize = is
                0019       cold  = .true.
                0020
60c4531875 Patr*0021       print *, 'pathei: in instore '
                0022
4cee17c1be Patr*0023 cc one record = two arrays of real*(isize)
                0024
                0025       open(  itape
                0026      $    , file   = 'OPWARMI'
                0027      $    , status = 'unknown'
                0028      $    , form   = 'formatted'
                0029      $    , access = 'sequential'
                0030      $    )
                0031
                0032       read( itape, *, end=800, err=900 ) n, fc, gnorm0,
                0033      $                                   isize, m, jmin, jmax
                0034       cold  = .false.
                0035
                0036  800  continue
                0037       close(itape)
                0038
                0039       print*
                0040       write(*,'(a,i10,a,i2)')
                0041      $        '  opening direct access with recl = ',n,' * ',isize
                0042       print*
                0043
f7e8e52ac8 Patr*0044       tmp1 = n
                0045       tmp2 = isize
                0046       rectmp = tmp1*tmp2
60c4531875 Patr*0047 cph(
                0048       print *, 'pathei rectmp ', rectmp, n, isize
                0049 cph)
                0050
4cee17c1be Patr*0051       open( ntape
                0052      $    , file   = 'OPWARMD'
                0053      $    , status = 'unknown'
                0054      $    , form   = 'unformatted'
                0055      $    , access = 'direct'
60c4531875 Patr*0056      $    , recl   = rectmp
4cee17c1be Patr*0057      $    )
                0058
                0059       return
                0060
                0061  900  continue
                0062       close(itape)
                0063       ifail = 1
                0064
                0065       return
                0066       end