subroutine getinp include 'scatter7' include 'partblk7' c namelist/input/ header, n1,iz1,n2,iz2,ndt,nprint,nruns,ref,elgev, c a inttyp,ieos,b,dt,sepfp,sepft,switch,cutfct,iseed, c b idtcl namelist/input/ jeader, n1,iz1,n2,iz2,ndt,nprint,nruns,jef,elgev, a inttyp,ieos,b,dt,sepfp,sepft,jwitch,cutfct,iseed, b idtcl c.fqw character*4 mype,dpet3e external dpet3e mype=dpet3e() c open (unit=5,file='file5',status='old') c open (unit=6,file='file6',status='old',form='formatted') c open (unit=7,form='unformatted',file='file7', c a status='unknown') c open (unit=8,file='file8',status='unknown') c open (unit=9,file='file9',status='old',form='formatted') c open (unit=10,file='file10',status='unknown') c open (unit=11,form='unformatted',file='file11', c a status='unknown') c open (unit=12,form='unformatted',file='file12', c a status='unknown') c open (unit=13,file='file13',status='unknown') c open (unit=14,file='file14',status='old') c open (unit=15,file='file15',status='unknown') open (unit=5,file=mype//'file5',status='old') open (unit=6,file=mype//'file6',status='old',form='formatted') open (unit=7,form='unformatted',file=mype//'file7', a status='unknown') open (unit=8,file=mype//'file8',status='unknown') open (unit=9,file=mype//'file9',status='old',form='formatted') open (unit=10,file=mype//'file10',status='unknown') open (unit=11,form='unformatted',file=mype//'file11', a status='unknown') open (unit=12,form='unformatted',file=mype//'file12', a status='unknown') open (unit=13,file=mype//'file13',status='unknown') open (unit=14,file=mype//'file14',status='old') open (unit=15,file=mype//'file15',status='unknown') c^fqw c.fqw jeader=ichar(header(1:1))+ichar(header(2:2))*256+ a ichar(header(3:3))*256**2+ichar(header(4:4))*256**3 if (ref.eq.'cld') then jef=1 elseif (ref.eq.'lab') then jef=2 else jef=0 endif do 100 iii=1,numswi if (switch(iii)) then jwitch(iii)=1 else jwitch(iii)=0 endif 100 continue c^fqw read(unit=5,nml=input) c.fqw header(1:1)=char(mod(jeader,256)) header(2:2)=char(mod(jeader/256,256)) header(3:3)=char(mod(jeader/256**2,256)) header(4:4)=char(mod(jeader/256**3,256)) if (jef.eq.1) then ref='cld' elseif (jef.eq.2) then ref='lab' else ref='nn ' endif do 101 iii=1,numswi if (jwitch(iii).eq.0) then switch(iii)=.false. else switch(iii)=.true. endif 101 continue c^fqw return end