Welcome to our forums
In my fortran program its showing error as ( INCREMENTRUNNUM - in file onechr.f at line 858 [+02f5] MAIN# -)? What is that..????
Show us line 858 of your code together with a few lines before it. Also show us the full error report.
c Main for New Chromosome Simulation Program c-----------------------------------------------------------------------
include 'mc.inc'
include 'mcds.inc'
integer i,jobno,ii,j
logical firstloop
character*60 fncomeq
equivalence(fncomeq,fncom(1))
character*12 pofn
character*1 epofn(12)
equivalence(epofn(1),pofn)
bowie=.false.
deskside=.true.
strucread=0
thischr=1
useonewidth=0
c-----------------------------------------------------------------------
if(bowie) then
poffset=7
bos=22
else
poffset=0
c For ANA1 and ANA2
if(deskside) then
bos=15
else
c123456 cctron:cseq.mmm:mmm.rules.raw
bos=6
endif
endif
c-----------------------------------------------------------------------
c For Bowie work
firstloop=.true.
call incrementrunnum
inullchar=0
c write(wunit,1101) c1101 format('Job number =',$) c read(runit,1102) jobno c1102 format(i5) c if(jobno .eq. 0) jobno=1
c 0 = functional names not generated c 1 = functional names will be generated
funnames=1
rikenflg=0
permflg=1
timeseed=59823
group5=.true.
block5=.true.
group6=.false.
group7=.true.
group8=.false.
stopprog=.false.
separatepo2=.false.
ltotalseq=0
selfflg=0
useexon=1
chrnum=21
rw2flg=0
ws4override=0
usemfa=.false.
rawnflg=0
zaplf=.false.
fmissing=.true.
ucset=0
highereukaryote=0
rmaxmachinenum=22
maxmachinenum=rmaxmachinenum
maxsspt=maxss+10
numzf=7
minint=21
forcepercent=0
if(smallgenome .eq. 1) then !a
numzf=5
minint=15
endif !a
showfilenames=.true.
ihour=0
iminute=0
isecond=0
oldseconds=0
numchr=22
c numchr=1
write(wunit,102)
102 format('Find the clusters of all chromosomes')
c-----------------------------------------------------------------------
c For the smaller eukaryotes, use 5-finger rules
if( ((thisfn(1) .eq. 'a') .and.
1 (thisfn(2) .eq. 't') .and. (thisfn(3) .eq. 'h')) .or.
1 ((thisfn(1) .eq. 'w') .and.
1 (thisfn(2) .eq. 'r') .and. (thisfn(3) .eq. 'm')) .or.
1 ((thisfn(1) .eq. 'y') .and.
1 (thisfn(2) .eq. 's') .and. (thisfn(3) .eq. 't')) .or.
1 ((thisfn(1) .eq. 'd') .and.
1 (thisfn(2) .eq. 'm') .and. (thisfn(3) .eq. 'e')) ) then !c
c numzf=5 c minint=15
c 11/26/04 all the smaller eukaryotes are - for the moment - c considered to be 8-finger genomes
numzf=7
minint=21
zaplf=.true.
endif !c
numzbases=numzf*3
c-----------------------------------------------------------------------
if ( ((thisfn(1) .eq. 'h') .and.
1 (thisfn(2) .eq. 's')) .or.
1 ((thisfn(1) .eq. 'a') .and.
1 (thisfn(2) .eq. 't') .and. (thisfn(3) .eq. 'h')) .or.
1 ((thisfn(1) .eq. 'w') .and.
1 (thisfn(2) .eq. 'r') .and. (thisfn(3) .eq. 'm')) .or.
1 ((thisfn(1) .eq. 'd') .and.
1 (thisfn(2) .eq. 'm') .and. (thisfn(3) .eq. 'e')) .or.
1 ((thisfn(1) .eq. 'm') .and.
1 (thisfn(2) .eq. 'm') .and.
1 ((thisfn(3) .eq. 'm') .or. (thisfn(3) .eq. 'u'))) .or.
1 ((thisfn(1) .eq. 'r') .and. (thisfn(2) .eq. 'n') .and.
1 (thisfn(3) .eq. 'o')) ) then !x
if(smallgenome .eq. 1) then !b
write(wunit,7265)
7265 format('This is a large eukaryotic genome'/ 1 'The variable smallgenome must = 0 - Exiting')
goto 11
endif !b
endif !x
c-----------------------------------------------------------------------
if ( ((thisfn(1) .eq. 'h') .and.
1 (thisfn(2) .eq. 's')) .or.
1 ((thisfn(1) .eq. 'm') .and.
1 (thisfn(2) .eq. 'm') .and.
1 ((thisfn(3) .eq. 'm') .or. (thisfn(3) .eq. 'u'))) .or.
1 ((thisfn(1) .eq. 'r') .and. (thisfn(2) .eq. 'n') .and.
1 (thisfn(3) .eq. 'o')) ) then !x
highereukaryote=1
if((thisfn(1) .eq. 'm') .and.
1 (thisfn(2) .eq. 'm') .and.
1 ((thisfn(3) .eq. 'm') .or. (thisfn(3) .eq. 'u'))) then !s
rikenflg=1
write(wunit,7356)
7356 format('RIKEN flag set')
endif !s
endif !x
c-----------------------------------------------------------------------
c 1 2 3 4
c12345678901234567890123456789012345678901234567 - 47
c/mnt/disk3/allgseq/mmm/003/allgroups.000.pv.000
fncomeq=
1 '/mnt/disk3/allgseq/mmm/003/allgroups.000.pv.000'
fncom(20)=thisfn(1)
fncom(21)=thisfn(2)
fncom(22)=thisfn(3)
fncom(1)=fchar
fncom(5)=barrier
fncom(11)=barrier
fncom(19)=barrier
fncom(23)=barrier
fncom(27)=barrier
efnpos=45
c-----------------------------------------------------------------------
c 1 2 3 4 5
c12345678901234567890123456789012345678901234567890123 - 53
c/home/rjfeldma/allgseq/mmm/003/group1/local.group1.01
cctron:rjfeldma:allgseq:mmm:000:group1:local.group1.01
rfile=
1 'ctron:rjfeldma:allgseq:mmm:000:group1:local.group1.01'
erfile(24)=thisfn(1)
erfile(25)=thisfn(2)
erfile(26)=thisfn(3)
erfile(1)=fchar
erfile(6)=barrier
erfile(15)=barrier
erfile(23)=barrier
erfile(27)=barrier
erfile(31)=barrier
erfile(38)=barrier
erfile(54)=nullchar
erfile(55)=nullchar
if(wunit .ne. 9) then !s
c 1 2 3 4 5
c123456789012345678901234567890123456789012345678901 - 51
c/mnt/disk2/allgseq/mmm/01.41/group1/local.group1.01
rfile=
1 '/mnt/disk2/allgseq/mmm/01.41/group1/local.group1.01'
erfile(20)=thisfn(1)
erfile(21)=thisfn(2)
erfile(22)=thisfn(3)
erfile(50)=nullchar
erfile(51)=nullchar
erfile(52)=nullchar
erfile(53)=nullchar
endif !s
c-----------------------------------------------------------------------
c 1 2 3 4
c12345678901234567890123456789012345678901234567 - 44
c/home/rjfeldma/allgseq/mmm/003/allgroups.000.p0
cctron:rjfeldma:allgseq:mmm:000:allgroups.000.p0
pfile=
1 'ctron:rjfeldma:allgseq:mmm:000:allgroups.000'
c For Bowie work
c 1 2 3 4
c123456789012345678901234567890123456789012345678901 - 51
c/share/homes/feldmann:allgseq:mmm:000:allgroups.000
if(poffset .ne. 0) pfile=
1 '/share/homes/feldmann:allgseq:mmm:000:allgroups.000'
epfile(24+poffset)=thisfn(1)
epfile(25+poffset)=thisfn(2)
epfile(26+poffset)=thisfn(3)
epfile(1)=fchar
epfile(6+poffset)=barrier
epfile(15+poffset)=barrier
epfile(23+poffset)=barrier
epfile(27+poffset)=barrier
epfile(31+poffset)=barrier
epfile(47+poffset)=nullchar
epfile(48+poffset)=nullchar
c 1 2 3 4 5 c12345678901234567890123456789012345678901234567890 - 47 c/home/rjfeldma/allgseq/mmm/003/allsequences.000.p0 cctron:rjfeldma:allgseq:mmm:000:allsequences.000.p0
pqfile=
1 'ctron:rjfeldma:allgseq:mmm:000:allsequences.000.p0'
epqfile(24)=thisfn(1)
epqfile(25)=thisfn(2)
epqfile(26)=thisfn(3)
epqfile(1)=fchar
epqfile(6+poffset)=barrier
epqfile(15+poffset)=barrier
epqfile(23+poffset)=barrier
epqfile(27+poffset)=barrier
epqfile(31+poffset)=barrier
c-----------------------------------------------------------------------
c 1 2 3 4 5 c12345678901234567890123456789012345678901234567890 - 43 c/home/rjfeldma/sseq/eco/eco01.fna cctron:rjfeldma:sseq:eco:eco01.fna
sfile=
1 'ctron:rjfeldma:sseq:eco:eco01.fna'
c For Bowie work
c 1 2 3 4
c1234567890123456789012345678901234567890 - 40
c/share/homes/feldmann:sseq:eco:eco01.fna
if(poffset .ne. 0) then !a
sfile=
1 '/share/homes/feldmann/sseq/eco/eco01.fna'
else !a
esfile(1)=fchar
esfile(6)=barrier
esfile(15)=barrier
esfile(20)=barrier
esfile(24)=barrier
esfile(34)=nullchar
esfile(35)=nullchar
endif !a
esfile(21+poffset)=thisfn(1)
esfile(22+poffset)=thisfn(2)
esfile(23+poffset)=thisfn(3)
esfile(25+poffset)=thisfn(1)
esfile(26+poffset)=thisfn(2)
esfile(27+poffset)=thisfn(3)
c-----------------------------------------------------------------------
c 1 2 3 4 5 6 c12345678901234567890123456789012345678901234567890123456789012 - 49 c/home/rjfeldma/gseq/mmu/mmun.map.group1.XXXX.YYYY cctron:rjfeldma:gseq:mmu:mmun.map.group1.XXXX.YYYY
sgroup=
1 'ctron:rjfeldma:gseq:mmu:mmun.map.group1.XXXX.YYYY'
esgroup(1)=fchar
esgroup(6)=barrier
esgroup(15)=barrier
esgroup(20)=barrier
esgroup(24)=barrier
esgroup(21)=thisfn(1)
esgroup(22)=thisfn(2)
esgroup(23)=thisfn(3)
esgroup(25)=thisfn(1)
esgroup(26)=thisfn(2)
esgroup(27)=thisfn(3)
esgroup(50)=nullchar
esgroup(51)=nullchar
egnfptr=41
c-----------------------------------------------------------------------
c 1 2 3 4 5 6 c12345678901234567890123456789012345678901234567890123456789012345678 - 52 c/home/rjfeldma/gseq/mmu/mmun.map.sequences.XXXX.YYYY cctron:rjfeldma:gseq:mmu:mmun.map.sequences.XXXX.YYYY
ssequence=
1 'ctron:rjfeldma:gseq:mmu:mmun.map.sequences.XXXX.YYYY'
essequence(1)=fchar
essequence(6)=barrier
essequence(15)=barrier
essequence(20)=barrier
essequence(24)=barrier
essequence(21)=thisfn(1)
essequence(22)=thisfn(2)
essequence(23)=thisfn(3)
essequence(25)=thisfn(1)
essequence(26)=thisfn(2)
essequence(27)=thisfn(3)
essequence(53)=nullchar
essequence(54)=nullchar
esnfptr=44
c-----------------------------------------------------------------------
if((poffset .eq. 0) .and. (wunit .ne. 9)) then !a
esgroup(2)='h'
esgroup(3)='o'
esgroup(4)='m'
esgroup(5)='e'
essequence(2)='h'
essequence(3)='o'
essequence(4)='m'
essequence(5)='e'
esfile(2)='h'
esfile(3)='o'
esfile(4)='m'
esfile(5)='e'
c if(smallgenome .eq. 1) then !s
c erfile(2)='h' c erfile(3)='o' c erfile(4)='m' c erfile(5)='e'
c endif !s
epfile(2)='h'
epfile(3)='o'
epfile(4)='m'
epfile(5)='e'
epqfile(2)='h'
epqfile(3)='o'
epqfile(4)='m'
epqfile(5)='e'
endif !a
c-----------------------------------------------------------------------
do 10 i=1,100
c write(wunit,103) 103 format(/'1 = Simplify and prepare files'/ 1 '2 = Mage format'/ 1 '3 = Run genome simulation'/ 1 '4 = Onesort for genome simulation'/ 1 '5 = Run chromosome simulation'/ 1 '6 = Onesort for chromosome simulation'/ 1 '7 = Run unique endpoint chromosome simulation'/ 1 '8 = Determine topological number of the qj files'/ 1 '9 = Simplify the chromosome files'/ 1 '10 = Determine topological number of the qm files'/ 1 '11 = Reduce the local features in the ct files'/ 1 '12 = All Connectrons in each chromosome'/ 1 '13 = Generate Initial States'/ 1 '14 = Calculate State Matrix'/ 1 '15 = Run Experiments'/ 1 '16 = Make C and T sequence files for Kishore'/ 1 '17 = Determined C1/C2 lifetimes'/ 1 '18 = Merge sorted c.sequence') write(wunit,303) 303 format(/ 1 '1 = Specify chromosome'/ 1 '2 = Chromosomal plot (R0)'/ 1 '3 = Chromosomal plot (R1)'/ 1 '4 = Chromosomal plot (R2)'/ 1 '5 = Chromosomal plot (R3)'/ 1 '6 = Chromosomal plot (R4)'/ 1 '-----------------------------------------------------'/ 1 '22 = Make state difference A - B → C'/ 1 '23 = Wide lines'/ 1 '24 = Thin lines'/ 1 '-----------------------------------------------------') write(wunit,403) 403 format( 1 '25 = Cluster relationships'/ 1 '26 = Plot Simulation (R9)'/ 1 '27 = Plot Cross Chromosome - OneShot Oriented (R10)'/ 1 '28 = Find OS & Driven by OS Connectrons'/ 1 '29 = Plot Cross Chromosome - New/Dead Oriented (R11)'/ 1 '30 = Do 26, 27, 29 and 31'/ 1 '31 = Plot GC/NC Tedrad (R12)'/ 1 '-----------------------------------------------------'/ 1 '40 = Scan mammalian genome'/ 1 '41 = New Split'/ 1 '42 = Make Split file'// 1 'cmd = ',$)
read(runit,104,err=10) cc
104 format(i3)
write(wunit,104) cc
c-----------------------------------------------------------------------
if(cc .eq. 0) goto 11
c-----------------------------------------------------------------------
if(firstloop) then !f
c 1234567 pofn='ms.poXX'
epofn(6)=chonear(cc/10)
epofn(7)=chonear(mod(cc,10))
write(wunit,203) pofn
203 format('W-',a7)
open(unit=2,file=pofn,status='unknown',err=205)
goto 204
205 continue
close(unit=2,status='delete')
j=2
open(unit=j,file=pofn,status='new')
204 continue
write(2,102)
write(2,202) runnum
202 format('Run=',i9,' ')
c write(2,103) write(2,303) write(2,104) cc
endif !f
call rjob
firstloop=.false.
c-----------------------------------------------------------------------
call rjfsince
if(smallgenome .eq. 0) then !d
write(wunit,3373)
write(2,3373)
3373 format('This is a big genome')
else !d
write(wunit,3374)
write(2,3374)
3374 format('This is a small genome')
endif !d
write(2,110) '*',sfile
c write(wunit,110) '',sgroup write(2,110) '',sgroup 110 format(a1,'-',a62)
c write(wunit,111) '',ssequence write(2,111) '',ssequence 111 format(a1,'-',a68)
write(2,109)
109 format('-------------------------------------------------')
c-----------------------------------------------------------------------
if(cc .eq. 1) then !1
write(wunit,115)
115 format('Chromosome =',$)
read(runit,104) thischr
if(thischr .le. 0) thischr=1
write(wunit,115)
write(wunit,104) thischr
write(2,115)
write(2,104) thischr
endif !1
c-----------------------------------------------------------------------
if(cc .eq. 2) then !2
ioord='i'
call mageformat
endif !2
c-----------------------------------------------------------------------
if(cc .eq. 3) then !3
ioord='o'
call mageformat
endif !3
c-----------------------------------------------------------------------
if(cc .eq. 4) then !4
ioord='l'
call mageformat
endif !4
c-----------------------------------------------------------------------
if(cc .eq. 5) then !5
ioord='t'
call mageformat
endif !5
c-----------------------------------------------------------------------
if(cc .eq. 6) then !6
ioord='6'
call mageformat
endif !6
c-----------------------------------------------------------------------
if(cc .eq. 22) then !6
call msdifference
if(stopprog) goto 11
endif !6
c-----------------------------------------------------------------------
if(cc .eq. 23) useonewidth=0
if(cc .eq. 24) useonewidth=1
c-----------------------------------------------------------------------
if(cc .eq. 25) then !25
call aftgraphstructure
endif !25
c-----------------------------------------------------------------------
if(cc .eq. 26) then !26
call plotsimulation
if(stopprog) goto 11
endif !26
c-----------------------------------------------------------------------
if(cc .eq. 27) then !27
style=1
call plotcrosschromosome
if(stopprog) goto 11
endif !27
c-----------------------------------------------------------------------
if(cc .eq. 28) then !28
call finddrivenbyos
if(stopprog) goto 11
endif !28
c-----------------------------------------------------------------------
if(cc .eq. 29) then !29
style=2
call plotcrosschromosome
if(stopprog) goto 11
endif !29
c-----------------------------------------------------------------------
if(cc .eq. 30) then !30
call plotsimulation
if(stopprog) goto 11
style=1
call plotcrosschromosome
if(stopprog) goto 11
style=2
call plotcrosschromosome
if(stopprog) goto 11
call plotgcnc
if(stopprog) goto 11
endif !30
c-----------------------------------------------------------------------
if(cc .eq. 31) then !31
call plotgcnc
if(stopprog) goto 11
endif !31
c-----------------------------------------------------------------------
if(cc .eq. 40) then !40
call scanmammaliangenome
if(stopprog) goto 11
endif !40
c-----------------------------------------------------------------------
if(cc .eq. 41) then !41
call nsplit
if(stopprog) goto 11
endif !41
c-----------------------------------------------------------------------
if(cc .eq. 42) then !42
call msplitfile
if(stopprog) goto 11
endif !42
c-----------------------------------------------------------------------
10 continue
11 continue
c-----------------------------------------------------------------------
if(wunit .eq. 9) then !b
write(wunit,101)
101 format('at pause')
pause
endif !b
end
subroutine incrementrunnum
c----------------------------------------------------------------------- c c Increment the run counter c
include 'mc.inc'
integer jj,startgene,tnumzf
startgene=0
open(unit=1,file='run.number',status='old',err=121)
read(1,103) thisfn(1),thisfn(2),thisfn(3),thisfn(4),
1 runnum,startgene,tnumzf,jj
close(unit=1)
runnum=runnum+1
if(jj .eq. 9) then
c Macintosh runit=9 wunit=9
barrier=':'
fchar='c'
else
c SGI wunit=6 runit=5
barrier='/'
fchar='/'
endif
write(wunit,104) 'R','run.number'
104 format(a1,'-',a10)
open(unit=1,file='run.number',status='old')
rewind(unit=1)
write(1,103) thisfn(1),thisfn(2),thisfn(3),thisfn(4),
1 runnum,startgene,numzf,jj
103 format(4a1,5i5)
write(wunit,102) runnum
102 format('Run=',i9,' ')
close(unit=1)
goto 122
121 continue write(wunit,123) 123 format('RUN.NUMBER file missing') stopprog=.true.
122 continue
end
block data
c----------------------------------------------------------------------- c c Initialize the data statements c
include 'mc.inc'
data chonear /'0','1','2','3','4','5','6','7','8','9',
1 ' ','-',' ','*','o','|',
1 'f','g','h','i','j',
1 'k','l','m','n','o',
1 'p','q','r','s','t',
1 'u','v','w','x','y',
1 'z','A','B','C','D',
1 'E','F','G','H','I',
1 'J','K','L','M','N',
1 'O','P','Q','R','S'/
data schrlimit /10,10,10,10,10,
1 10,10,10,10,10,
1 10,10,10,10,10,
1 10,10,10,10,10,
1 10,10/
end
Error Line : ** write(wunit,104) 'R','run.number'**
Any One there to help me...???
Quoted from mmaharish Any One there to help me...???
Not possible.
The source code that you provided was truncated. The 'include' files are missing. We do not know how you compiled the program, and what data you ran it with.
The line with the error, 'write(*,104)..' is not in the listing, nor is the Format statement 104.
I'm having INCLUDE FILE and full coding..it's a big project. this project consists lot of FORTRAN files. My boss asked me to search for main class file and run it. I found this(the program i pasted) was the main class file. Plz help me to run the program. i pasted full program but it's not pasting fully. How can i paste it friend?
c Main for New Chromosome Simulation Program c-----------------------------------------------------------------------
include 'mc.inc'
include 'mcds.inc'
integer i,jobno,ii,j
logical firstloop
character*60 fncomeq
equivalence(fncomeq,fncom(1))
character*12 pofn
character*1 epofn(12)
equivalence(epofn(1),pofn)
bowie=.false.
deskside=.true.
strucread=0
thischr=1
useonewidth=0
c-----------------------------------------------------------------------
if(bowie) then
poffset=7
bos=22
else
poffset=0
c For ANA1 and ANA2
if(deskside) then
bos=15
else
c123456 cctron:cseq.mmm:mmm.rules.raw
bos=6
endif
endif
c-----------------------------------------------------------------------
c For Bowie work
firstloop=.true.
call incrementrunnum
inullchar=0
c write(wunit,1101) c1101 format('Job number =',$) c read(runit,1102) jobno c1102 format(i5) c if(jobno .eq. 0) jobno=1
c 0 = functional names not generated c 1 = functional names will be generated
funnames=1
rikenflg=0
permflg=1
timeseed=59823
group5=.true.
block5=.true.
group6=.false.
group7=.true.
group8=.false.
stopprog=.false.
separatepo2=.false.
ltotalseq=0
selfflg=0
useexon=1
chrnum=21
rw2flg=0
ws4override=0
usemfa=.false.
rawnflg=0
zaplf=.false.
fmissing=.true.
ucset=0
highereukaryote=0
rmaxmachinenum=22
maxmachinenum=rmaxmachinenum
maxsspt=maxss+10
numzf=7
minint=21
forcepercent=0
if(smallgenome .eq. 1) then !a
numzf=5
minint=15
endif !a
showfilenames=.true.
ihour=0
iminute=0
isecond=0
oldseconds=0
numchr=22
c numchr=1
write(wunit,102)
102 format('Find the clusters of all chromosomes')
c-----------------------------------------------------------------------
c For the smaller eukaryotes, use 5-finger rules
if( ((thisfn(1) .eq. 'a') .and.
1 (thisfn(2) .eq. 't') .and. (thisfn(3) .eq. 'h')) .or.
1 ((thisfn(1) .eq. 'w') .and.
1 (thisfn(2) .eq. 'r') .and. (thisfn(3) .eq. 'm')) .or.
1 ((thisfn(1) .eq. 'y') .and.
1 (thisfn(2) .eq. 's') .and. (thisfn(3) .eq. 't')) .or.
1 ((thisfn(1) .eq. 'd') .and.
1 (thisfn(2) .eq. 'm') .and. (thisfn(3) .eq. 'e')) ) then !c
c numzf=5 c minint=15
c 11/26/04 all the smaller eukaryotes are - for the moment - c considered to be 8-finger genomes
numzf=7
minint=21
zaplf=.true.
endif !c
numzbases=numzf*3
c-----------------------------------------------------------------------
if ( ((thisfn(1) .eq. 'h') .and.
1 (thisfn(2) .eq. 's')) .or.
1 ((thisfn(1) .eq. 'a') .and.
1 (thisfn(2) .eq. 't') .and. (thisfn(3) .eq. 'h')) .or.
1 ((thisfn(1) .eq. 'w') .and.
1 (thisfn(2) .eq. 'r') .and. (thisfn(3) .eq. 'm')) .or.
1 ((thisfn(1) .eq. 'd') .and.
1 (thisfn(2) .eq. 'm') .and. (thisfn(3) .eq. 'e')) .or.
1 ((thisfn(1) .eq. 'm') .and.
1 (thisfn(2) .eq. 'm') .and.
1 ((thisfn(3) .eq. 'm') .or. (thisfn(3) .eq. 'u'))) .or.
1 ((thisfn(1) .eq. 'r') .and. (thisfn(2) .eq. 'n') .and.
1 (thisfn(3) .eq. 'o')) ) then !x
if(smallgenome .eq. 1) then !b
write(wunit,7265)
7265 format('This is a large eukaryotic genome'/ 1 'The variable smallgenome must = 0 - Exiting')
goto 11
endif !b
endif !x
c-----------------------------------------------------------------------
if ( ((thisfn(1) .eq. 'h') .and.
1 (thisfn(2) .eq. 's')) .or.
1 ((thisfn(1) .eq. 'm') .and.
1 (thisfn(2) .eq. 'm') .and.
1 ((thisfn(3) .eq. 'm') .or. (thisfn(3) .eq. 'u'))) .or.
1 ((thisfn(1) .eq. 'r') .and. (thisfn(2) .eq. 'n') .and.
1 (thisfn(3) .eq. 'o')) ) then !x
highereukaryote=1
if((thisfn(1) .eq. 'm') .and.
1 (thisfn(2) .eq. 'm') .and.
1 ((thisfn(3) .eq. 'm') .or. (thisfn(3) .eq. 'u'))) then !s
rikenflg=1
write(wunit,7356)
7356 format('RIKEN flag set')
endif !s
endif !x
c-----------------------------------------------------------------------
c 1 2 3 4
c12345678901234567890123456789012345678901234567 - 47
c/mnt/disk3/allgseq/mmm/003/allgroups.000.pv.000
fncomeq=
1 '/mnt/disk3/allgseq/mmm/003/allgroups.000.pv.000'
fncom(20)=thisfn(1)
fncom(21)=thisfn(2)
fncom(22)=thisfn(3)
fncom(1)=fchar
fncom(5)=barrier
fncom(11)=barrier
fncom(19)=barrier
fncom(23)=barrier
fncom(27)=barrier
efnpos=45
c-----------------------------------------------------------------------
c 1 2 3 4 5
c12345678901234567890123456789012345678901234567890123 - 53
c/home/rjfeldma/allgseq/mmm/003/group1/local.group1.01
cctron:rjfeldma:allgseq:mmm:000:group1:local.group1.01
rfile=
1 'ctron:rjfeldma:allgseq:mmm:000:group1:local.group1.01'
erfile(24)=thisfn(1)
erfile(25)=thisfn(2)
erfile(26)=thisfn(3)
erfile(1)=fchar
erfile(6)=barrier
erfile(15)=barrier
erfile(23)=barrier
erfile(27)=barrier
erfile(31)=barrier
erfile(38)=barrier
erfile(54)=nullchar
erfile(55)=nullchar
if(wunit .ne. 9) then !s
c 1 2 3 4 5
c123456789012345678901234567890123456789012345678901 - 51
c/mnt/disk2/allgseq/mmm/01.41/group1/local.group1.01
rfile=
1 '/mnt/disk2/allgseq/mmm/01.41/group1/local.group1.01'
erfile(20)=thisfn(1)
erfile(21)=thisfn(2)
erfile(22)=thisfn(3)
erfile(50)=nullchar
erfile(51)=nullchar
erfile(52)=nullchar
erfile(53)=nullchar
endif !s
c-----------------------------------------------------------------------
c 1 2 3 4
c12345678901234567890123456789012345678901234567 - 44
c/home/rjfeldma/allgseq/mmm/003/allgroups.000.p0
cctron:rjfeldma:allgseq:mmm:000:allgroups.000.p0
pfile=
1 'ctron:rjfeldma:allgseq:mmm:000:allgroups.000'
c For Bowie work
c 1 2 3 4
c123456789012345678901234567890123456789012345678901 - 51
c/share/homes/feldmann:allgseq:mmm:000:allgroups.000
if(poffset .ne. 0) pfile=
1 '/share/homes/feldmann:allgseq:mmm:000:allgroups.000'
epfile(24+poffset)=thisfn(1)
epfile(25+poffset)=thisfn(2)
epfile(26+poffset)=thisfn(3)
epfile(1)=fchar
epfile(6+poffset)=barrier
epfile(15+poffset)=barrier
epfile(23+poffset)=barrier
epfile(27+poffset)=barrier
epfile(31+poffset)=barrier
epfile(47+poffset)=nullchar
epfile(48+poffset)=nullchar
c 1 2 3 4 5 c12345678901234567890123456789012345678901234567890 - 47 c/home/rjfeldma/allgseq/mmm/003/allsequences.000.p0 cctron:rjfeldma:allgseq:mmm:000:allsequences.000.p0
pqfile=
1 'ctron:rjfeldma:allgseq:mmm:000:allsequences.000.p0'
epqfile(24)=thisfn(1)
epqfile(25)=thisfn(2)
epqfile(26)=thisfn(3)
epqfile(1)=fchar
epqfile(6+poffset)=barrier
epqfile(15+poffset)=barrier
epqfile(23+poffset)=barrier
epqfile(27+poffset)=barrier
epqfile(31+poffset)=barrier
c-----------------------------------------------------------------------
c 1 2 3 4 5 c12345678901234567890123456789012345678901234567890 - 43 c/home/rjfeldma/sseq/eco/eco01.fna cctron:rjfeldma:sseq:eco:eco01.fna
sfile=
1 'ctron:rjfeldma:sseq:eco:eco01.fna'
c For Bowie work
c 1 2 3 4
c1234567890123456789012345678901234567890 - 40
c/share/homes/feldmann:sseq:eco:eco01.fna
if(poffset .ne. 0) then !a
sfile=
1 '/share/homes/feldmann/sseq/eco/eco01.fna'
else !a
esfile(1)=fchar
esfile(6)=barrier
esfile(15)=barrier
esfile(20)=barrier
esfile(24)=barrier
esfile(34)=nullchar
esfile(35)=nullchar
endif !a
esfile(21+poffset)=thisfn(1)
esfile(22+poffset)=thisfn(2)
esfile(23+poffset)=thisfn(3)
esfile(25+poffset)=thisfn(1)
esfile(26+poffset)=thisfn(2)
esfile(27+poffset)=thisfn(3)
c-----------------------------------------------------------------------
c 1 2 3 4 5 6 c12345678901234567890123456789012345678901234567890123456789012 - 49 c/home/rjfeldma/gseq/mmu/mmun.map.group1.XXXX.YYYY cctron:rjfeldma:gseq:mmu:mmun.map.group1.XXXX.YYYY
sgroup=
1 'ctron:rjfeldma:gseq:mmu:mmun.map.group1.XXXX.YYYY'
esgroup(1)=fchar
esgroup(6)=barrier
esgroup(15)=barrier
esgroup(20)=barrier
esgroup(24)=barrier
esgroup(21)=thisfn(1)
esgroup(22)=thisfn(2)
esgroup(23)=thisfn(3)
esgroup(25)=thisfn(1)
esgroup(26)=thisfn(2)
esgroup(27)=thisfn(3)
esgroup(50)=nullchar
esgroup(51)=nullchar
egnfptr=41
c-----------------------------------------------------------------------
c 1 2 3 4 5 6 c12345678901234567890123456789012345678901234567890123456789012345678 - 52 c/home/rjfeldma/gseq/mmu/mmun.map.sequences.XXXX.YYYY cctron:rjfeldma:gseq:mmu:mmun.map.sequences.XXXX.YYYY
ssequence=
1 'ctron:rjfeldma:gseq:mmu:mmun.map.sequences.XXXX.YYYY'
essequence(1)=fchar
essequence(6)=barrier
essequence(15)=barrier
essequence(20)=barrier
essequence(24)=barrier
essequence(21)=thisfn(1)
essequence(22)=thisfn(2)
essequence(23)=thisfn(3)
essequence(25)=thisfn(1)
essequence(26)=thisfn(2)
essequence(27)=thisfn(3)
essequence(53)=nullchar
essequence(54)=nullchar
esnfptr=44
c-----------------------------------------------------------------------
if((poffset .eq. 0) .and. (wunit .ne. 9)) then !a
esgroup(2)='h'
esgroup(3)='o'
esgroup(4)='m'
esgroup(5)='e'
essequence(2)='h'
essequence(3)='o'
essequence(4)='m'
essequence(5)='e'
esfile(2)='h'
esfile(3)='o'
esfile(4)='m'
esfile(5)='e'
c if(smallgenome .eq. 1) then !s
c erfile(2)='h' c erfile(3)='o' c erfile(4)='m' c erfile(5)='e'
c endif !s
epfile(2)='h'
epfile(3)='o'
epfile(4)='m'
epfile(5)='e'
epqfile(2)='h'
epqfile(3)='o'
epqfile(4)='m'
epqfile(5)='e'
endif !a
c-----------------------------------------------------------------------
do 10 i=1,100
c write(wunit,103) 103 format(/'1 = Simplify and prepare files'/ 1 '2 = Mage format'/ 1 '3 = Run genome simulation'/ 1 '4 = Onesort for genome simulation'/ 1 '5 = Run chromosome simulation'/ 1 '6 = Onesort for chromosome simulation'/ 1 '7 = Run unique endpoint chromosome simulation'/ 1 '8 = Determine topological number of the qj files'/ 1 '9 = Simplify the chromosome files'/ 1 '10 = Determine topological number of the qm files'/ 1 '11 = Reduce the local features in the ct files'/ 1 '12 = All Connectrons in each chromosome'/ 1 '13 = Generate Initial States'/ 1 '14 = Calculate State Matrix'/ 1 '15 = Run Experiments'/ 1 '16 = Make C and T sequence files for Kishore'/ 1 '17 = Determined C1/C2 lifetimes'/ 1 '18 = Merge sorted c.sequence') write(wunit,303) 303 format(/ 1 '1 = Specify chromosome'/ 1 '2 = Chromosomal plot (R0)'/ 1 '3 = Chromosomal plot (R1)'/ 1 '4 = Chromosomal plot (R2)'/ 1 '5 = Chromosomal plot (R3)'/ 1 '6 = Chromosomal plot (R4)'/ 1 '-----------------------------------------------------'/ 1 '22 = Make state difference A - B → C'/ 1 '23 = Wide lines'/ 1 '24 = Thin lines'/ 1 '-----------------------------------------------------') write(wunit,403) 403 format( 1 '25 = Cluster relationships'/ 1 '26 = Plot Simulation (R9)'/ 1 '27 = Plot Cross Chromosome - OneShot Oriented (R10)'/ 1 '28 = Find OS & Driven by OS Connectrons'/ 1 '29 = Plot Cross Chromosome - New/Dead Oriented (R11)'/ 1 '30 = Do 26, 27, 29 and 31'/ 1 '31 = Plot GC/NC Tedrad (R12)'/ 1 '-----------------------------------------------------'/ 1 '40 = Scan mammalian genome'/ 1 '41 = New Split'/ 1 '42 = Make Split file'// 1 'cmd = ',$)
read(runit,104,err=10) cc
104 format(i3)
write(wunit,104) cc
c-----------------------------------------------------------------------
if(cc .eq. 0) goto 11
c-----------------------------------------------------------------------
if(firstloop) then !f
c 1234567 pofn='ms.poXX'
epofn(6)=chonear(cc/10)
epofn(7)=chonear(mod(cc,10))
write(wunit,203) pofn
203 format('W-',a7)
open(unit=2,file=pofn,status='unknown',err=205)
goto 204
205 continue
close(unit=2,status='delete')
j=2
open(unit=j,file=pofn,status='new')
204 continue
write(2,102)
write(2,202) runnum
202 format('Run=',i9,' ')
c write(2,103) write(2,303) write(2,104) cc
endif !f
call rjob
firstloop=.false.
c-----------------------------------------------------------------------
call rjfsince
if(smallgenome .eq. 0) then !d
write(wunit,3373)
write(2,3373)
3373 format('This is a big genome')
else !d
write(wunit,3374)
write(2,3374)
3374 format('This is a small genome')
endif !d
write(2,110) '*',sfile
c write(wunit,110) '',sgroup write(2,110) '',sgroup 110 format(a1,'-',a62)
c write(wunit,111) '',ssequence write(2,111) '',ssequence 111 format(a1,'-',a68)
write(2,109)
109 format('-------------------------------------------------')
c-----------------------------------------------------------------------
if(cc .eq. 1) then !1
write(wunit,115)
115 format('Chromosome =',$)
read(runit,104) thischr
if(thischr .le. 0) thischr=1
write(wunit,115)
write(wunit,104) thischr
write(2,115)
write(2,104) thischr
endif !1
c-----------------------------------------------------------------------
if(cc .eq. 2) then !2
ioord='i'
call mageformat
endif !2
c-----------------------------------------------------------------------
if(cc .eq. 3) then !3
ioord='o'
call mageformat
endif !3
c-----------------------------------------------------------------------
if(cc .eq. 4) then !4
ioord='l'
call mageformat
endif !4
c-----------------------------------------------------------------------
if(cc .eq. 5) then !5
ioord='t'
call mageformat
endif !5
c-----------------------------------------------------------------------
if(cc .eq. 6) then !6
ioord='6'
call mageformat
endif !6
c-----------------------------------------------------------------------
if(cc .eq. 22) then !6
call msdifference
if(stopprog) goto 11
endif !6
c-----------------------------------------------------------------------
if(cc .eq. 23) useonewidth=0
if(cc .eq. 24) useonewidth=1
c-----------------------------------------------------------------------
if(cc .eq. 25) then !25
call aftgraphstructure
endif !25
c-----------------------------------------------------------------------
if(cc .eq. 26) then !26
call plotsimulation
if(stopprog) goto 11
endif !26
c-----------------------------------------------------------------------
if(cc .eq. 27) then !27
style=1
call plotcrosschromosome
if(stopprog) goto 11
endif !27
c-----------------------------------------------------------------------
if(cc .eq. 28) then !28
call finddrivenbyos
if(stopprog) goto 11
endif !28
c-----------------------------------------------------------------------
if(cc .eq. 29) then !29
style=2
call plotcrosschromosome
if(stopprog) goto 11
endif !29
c-----------------------------------------------------------------------
if(cc .eq. 30) then !30
call plotsimulation
if(stopprog) goto 11
style=1
call plotcrosschromosome
if(stopprog) goto 11
style=2
call plotcrosschromosome
if(stopprog) goto 11
call plotgcnc
if(stopprog) goto 11
endif !30
c-----------------------------------------------------------------------
if(cc .eq. 31) then !31
call plotgcnc
if(stopprog) goto 11
endif !31
c-----------------------------------------------------------------------
if(cc .eq. 40) then !40
call scanmammaliangenome
if(stopprog) goto 11
endif !40
c-----------------------------------------------------------------------
if(cc .eq. 41) then !41
call nsplit
if(stopprog) goto 11
endif !41
c-----------------------------------------------------------------------
if(cc .eq. 42) then !42
call msplitfile
if(stopprog) goto 11
endif !42
c-----------------------------------------------------------------------
10 continue
11 continue
c-----------------------------------------------------------------------
if(wunit .eq. 9) then !b
write(wunit,101)
101 format('at pause')
pause
endif !b
end
subroutine incrementrunnum
c----------------------------------------------------------------------- c c Increment the run counter c
include 'mc.inc'
integer jj,startgene,tnumzf
startgene=0
open(unit=1,file='run.number',status='old',err=121)
read(1,103) thisfn(1),thisfn(2),thisfn(3),thisfn(4),
1 runnum,startgene,tnumzf,jj
close(unit=1)
runnum=runnum+1
if(jj .eq. 9) then
c Macintosh runit=9 wunit=9
barrier=':'
fchar='c'
else
c SGI wunit=6 runit=5
barrier='/'
fchar='/'
endif
write(wunit,104) 'R','run.number'
104 format(a1,'-',a10)
open(unit=1,file='run.number',status='old')
rewind(unit=1)
write(1,103) thisfn(1),thisfn(2),thisfn(3),thisfn(4),
1 runnum,startgene,numzf,jj
103 format(4a1,5i5)
write(wunit,102) runnum
102 format('Run=',i9,' ')
close(unit=1)
goto 122
121 continue write(wunit,123) 123 format('RUN.NUMBER file missing') stopprog=.true.
122 continue
end
block data
c----------------------------------------------------------------------- c c Initialize the data statements c
include 'mc.inc'
data chonear /'0','1','2','3','4','5','6','7','8','9',
1 ' ','-',' ','*','o','|',
1 'f','g','h','i','j',
1 'k','l','m','n','o',
1 'p','q','r','s','t',
1 'u','v','w','x','y',
1 'z','A','B','C','D',
1 'E','F','G','H','I',
1 'J','K','L','M','N',
1 'O','P','Q','R','S'/
data schrlimit /10,10,10,10,10,
1 10,10,10,10,10,
1 10,10,10,10,10,
1 10,10,10,10,10,
1 10,10/
end
i pasted full program but only half of the program is displaying. Please say any other way i want a solution to my problem. please help me....im having al files related to this project.
It would help if you can find the line given by line number 858. Show us that line only with a few lines before it. Also give us the full error report.
What is the value of WUNIT? It must be either a standard value (2 or 6 or *) or correspond to an OPEN statement in your code.
this is the some line . . . .
wunit=6
runit=5
barrier='/'
fchar='/'
endif
***ERROR LINE IS: * ** write(wunit,104) 'R','run.number'
ERROR REPORT :
INCREMENTRUNNUM - in file harish.f at line 859 [+02f5] MAIN# - in file harish.f at line 63 [+00fa]
I cannot see anything wrong with this.
Try using CHECKMATE and use the debugger to run the program.
Quoted from mmaharish
INCREMENTRUNNUM - in file harish.f at line 859 [+02f5] MAIN# - in file harish.f at line 63 [+00fa]
That information tells us [u:60257a99e6]where [/u:60257a99e6]the error occurred, but not [u:60257a99e6]what [/u:60257a99e6]the error message was.
What is the debugger..?? I have run many files but every file is telling that, ERROR 94, UNIT NEITHER OPENED NOR PRECONNECTED
Im having many fortran files, i have to run the project. i dont know which is the first file i have to run. please help me to run.
anyone there???????????
Quoted from mmaharish anyone there???????????
Yes, I am here, as are others.
However, we cannot help when you ask questions without providing context. Nor is this the place to come before reaching a minimal level of competence in Fortran program development.
The error message is clear enough. Your program is attempting to read or write from a file which has not been opened. You have not shown any OPEN statements in the code fragments that you showed.
Reading the Salford compiler documentation would have informed you of this fact.
Im new to fortran so i don't know anything, but im having codings. My boss asked me to run the project. it consists of many fortran files, so i need help. i run one program that too is showing error. What i have to do now.??
Quoted from mmaharish What i have to do now.??
Well you could ask your boss how to run the program!
As there don't seem to be any open statements in your program, perhaps it relies on the user re-directing 'standard input' to read from your input file.
Assuming your program is called PROGRAM and your input file is called FILE try typing the following from a command window.
PROGRAM < FILE