pkg/grdchk/grdchk_readparms.F

C $Header: /u/gcmpack/MITgcm/pkg/grdchk/grdchk_readparms.F,v 1.6 2003/06/24 16:08:45 heimbach Exp $

#include "CPP_OPTIONS.h"


      subroutine grdchk_readparms( mythid )

c     ==================================================================
c     SUBROUTINE grdchk_readparms
c     ==================================================================
c
c     o Initialize the ECCO gradient check.
c
c     started: Christian Eckert eckert@mit.edu 03-Mar-2000
c     continued: heimbach@mit.edu: 13-Jun-2001
c
c     ==================================================================
c     SUBROUTINE grdchk_readparms
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"

#include "ctrl.h"
#include "grdchk.h"

c     == routine arguments ==

      integer mythid

#ifdef ALLOW_GRDCHK
c     == local variables ==

      integer i,j,k
      integer bi,bj
      integer itlo,ithi
      integer jtlo,jthi
      integer jmin,jmax
      integer imin,imax

      integer errio
      integer il

      character*(max_len_prec) record

c     == external ==

      integer  ilnblnk
      external ilnblnk

c     == end of interface ==

c--   Optimization parameters.
      namelist /grdchk_nml/
     &                  grdchk_eps,
     &                  nbeg,
     &                  nstep,
     &                  nend,
     &                  grdchkvarindex,
     &                  useCentralDiff,
     &                  grdchkwhichproc

      jtlo = mybylo(mythid)
      jthi = mybyhi(mythid)
      itlo = mybxlo(mythid)
      ithi = mybxhi(mythid)
      jmin = 1-oly
      jmax = sny+oly
      imin = 1-olx
      imax = snx+olx

      _BEGIN_MASTER( mythid )

c--     Set default values.
        grdchk_eps      = 1. _d 0
        nbeg            = 1
        nend            = 1
        nstep           = 1
        useCentralDiff  = .TRUE.
        grdchkwhichproc = 0

c--     Check versions.

        open(unit=scrunit1,status='scratch')

c--     Next, read the forcing data file.
        open(unit = modeldataunit,file = 'data.grdchk',
     &       status = 'old', iostat = errio)
        if ( errio .lt. 0 ) then
          stop ' ... stopped in grdchk_readparms.'
        endif

        do while ( .true. )
          read(modeldataunit, fmt='(a)', end=1001) record
          il = max(ilnblnk(record),1)
          if ( record(1:1) .ne. commentcharacter )
     &        write(unit=scrunit1, fmt='(a)') record(:il)
        enddo
 1001   continue

        close( modeldataunit )

        rewind( scrunit1 )
        read(unit = scrunit1, nml = grdchk_nml)
        close( scrunit1 )

c--     Summarize the gradient check's setup.
        call grdchk_Summary( mythid )
      _END_MASTER( mythid )

      _BARRIER

#endif /* ALLOW_GRDCHK */

      end


1	C $Header: /u/gcmpack/MITgcm/pkg/grdchk/grdchk_readparms.F,v 1.6 2003/06/24 16:08:45 heimbach Exp $
2
3	#include "CPP_OPTIONS.h"
4
5
6	subroutine grdchk_readparms( mythid )
7
8	c ==================================================================
9	c SUBROUTINE grdchk_readparms
10	c ==================================================================
11	c
12	c o Initialize the ECCO gradient check.
13	c
14	c started: Christian Eckert eckert@mit.edu 03-Mar-2000
15	c continued: heimbach@mit.edu: 13-Jun-2001
16	c
17	c ==================================================================
18	c SUBROUTINE grdchk_readparms
19	c ==================================================================
20
21	implicit none
22
23	c == global variables ==
24
25	#include "EEPARAMS.h"
26	#include "SIZE.h"
27
28	#include "ctrl.h"
29	#include "grdchk.h"
30
31	c == routine arguments ==
32
33	integer mythid
34
35	#ifdef ALLOW_GRDCHK
36	c == local variables ==
37
38	integer i,j,k
39	integer bi,bj
40	integer itlo,ithi
41	integer jtlo,jthi
42	integer jmin,jmax
43	integer imin,imax
44
45	integer errio
46	integer il
47
48	character*(max_len_prec) record
49
50	c == external ==
51
52	integer ilnblnk
53	external ilnblnk
54
55	c == end of interface ==
56
57	c-- Optimization parameters.
58	namelist /grdchk_nml/
59	& grdchk_eps,
60	& nbeg,
61	& nstep,
62	& nend,
63	& grdchkvarindex,
64	& useCentralDiff,
65	& grdchkwhichproc
66
67	jtlo = mybylo(mythid)
68	jthi = mybyhi(mythid)
69	itlo = mybxlo(mythid)
70	ithi = mybxhi(mythid)
71	jmin = 1-oly
72	jmax = sny+oly
73	imin = 1-olx
74	imax = snx+olx
75
76	_BEGIN_MASTER( mythid )
77
78	c-- Set default values.
79	grdchk_eps = 1. _d 0
80	nbeg = 1
81	nend = 1
82	nstep = 1
83	useCentralDiff = .TRUE.
84	grdchkwhichproc = 0
85
86	c-- Check versions.
87
88	open(unit=scrunit1,status='scratch')
89
90	c-- Next, read the forcing data file.
91	open(unit = modeldataunit,file = 'data.grdchk',
92	& status = 'old', iostat = errio)
93	if ( errio .lt. 0 ) then
94	stop ' ... stopped in grdchk_readparms.'
95	endif
96
97	do while ( .true. )
98	read(modeldataunit, fmt='(a)', end=1001) record
99	il = max(ilnblnk(record),1)
100	if ( record(1:1) .ne. commentcharacter )
101	& write(unit=scrunit1, fmt='(a)') record(:il)
102	enddo
103	1001 continue
104
105	close( modeldataunit )
106
107	rewind( scrunit1 )
108	read(unit = scrunit1, nml = grdchk_nml)
109	close( scrunit1 )
110
111	c-- Summarize the gradient check's setup.
112	call grdchk_Summary( mythid )
113	_END_MASTER( mythid )
114
115	_BARRIER
116
117	#endif /* ALLOW_GRDCHK */
118
119	end
120
121