25 implicit double precision(a-h,o-z)
27 integer hrec1(9999),hrec2(9999)
28 real*4 zrec1(9999),zrec2(9999)
29 equivalence(hrec1(1),zrec1(1))
30 equivalence(hrec2(1),zrec2(1))
36 write(6,*)
'program subtrc' 39 3
format(
'Enter "plus" input file: ',$)
42 open(lin1,file=fname,status=
'old',form=
'unformatted')
45 4
format(
'Enter "minus" input file: ',$)
47 open(lin2,file=fname,status=
'old',form=
'unformatted')
50 2
format(
'Enter "a-b" output file: ',$)
52 open(lout,file=fname,status=
'new',form=
'unformatted')
54 read(lin1) glamn1,glomn1,dgla1,dglo1,nla1,nlo1,ikind1
55 read(lin2) glamn2,glomn2,dgla2,dglo2,nla2,nlo2,ikind2
59 if(dabs(glomn1-glomn2).gt.1.d-7) stop 1
60 if(dabs(glamn1-glamn2).gt.1.d-7) stop 6
61 if(dabs(dgla1 -dgla2 ).gt.1.d-7) stop 2
62 if(dabs(dglo1 -dglo2 ).gt.1.d-7) stop 3
63 if(nlo1 .ne.nlo2 ) stop 4
64 if(nla1 .ne.nla2 ) stop 7
65 if(ikind1.eq.0.and.ikind2.ne.0.or.
66 * ikind2.eq.0.and.ikind1.ne.0) stop 5
68 write(lout) glamn1,glomn1,dgla1,dglo1,nla1,nlo1,ikind1
71 read (lin1) (hrec1(i),i=1,nlo1)
72 read (lin2) (hrec2(i),i=1,nlo1)
74 11 hrec1(i)=hrec1(i)-hrec2(i)
75 10
write(lout) (hrec1(i),i=1,nlo1)
78 read (lin1) (zrec1(i),i=1,nlo1)
79 read (lin2) (zrec2(i),i=1,nlo1)
81 21 zrec1(i)=zrec1(i)-zrec2(i)
82 20
write(lout) (zrec1(i),i=1,nlo1)
program subtrc
Part of the NADCON5 NADCON5 Core Library , Subtract one grid from another.