35
36
37
39
40
41
42#include "implicit_f.inc"
43
44
45
46 INTEGER IKAD(0:*),
47 . KFUNC,NFCT
48 CHARACTER KEY0(*)*5
49
50
51
52#include "units_c.inc"
53
54
55
56C
57
58
59
60 INTEGER I, N, NBC, K, IKEY, , NP
62 . t, f
63
64 k=0
65 ikey=kfunc
66 DO n=1,nfct
67 READ(iusc1,rec=ikad(ikey)+k,fmt='(49X,I10)',
68 . err=9990)nbc
69 k=k+1
70 CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
71 k=k+1
72 READ(iusc2,*,err=9990,
END=9990)
if
73 np=nbc-1
74 WRITE(iin,'(2I10)')IF,np
75 DO i=1,np
76 CALL wriusc2(ikad(ikey)+k,1,key0(ikey))
77 k=k+1
78 READ(iusc2,*,err=9990,END=9990)T,f
79 WRITE(iin,'(2E16.9)')t,f
80 ENDDO
81 ENDDO
82
83 RETURN
84
85 9990 CONTINUE
86 CALL ancmsg(msgid=73,anmode=aninfo,
87 . c1=key0(ikey))
if(complex_arithmetic) id
subroutine ancmsg(msgid, msgtype, anmode, i1, i2, i3, i4, i5, i6, i7, i8, i9, i10, i11, i12, i13, i14, i15, i16, i17, i18, i19, i20, r1, r2, r3, r4, r5, r6, r7, r8, r9, c1, c2, c3, c4, c5, c6, c7, c8, c9, prmode)
subroutine wriusc2(irec, nbc, key0)