diff --git a/astro.F b/astro.F index a0338dc00..b97c9cc0b 100644 --- a/astro.F +++ b/astro.F @@ -37,8 +37,6 @@ C NB: may want to smooth the Tsky map to 10 degrees or so. open(13,file=fname,status='old',access='stream',err=10) read(13) nsky close(13) -! call rfile2(fname,nsky,129600,nr) -! if(nr.ne.129600) go to 10 #endif ltsky=.true. first=.false. diff --git a/getfile2.F90 b/getfile2.F90 index 914d6f013..ad2ede349 100644 --- a/getfile2.F90 +++ b/getfile2.F90 @@ -33,22 +33,12 @@ subroutine getfile2(fname,len) ndecoding=4 monitoring=0 kbuf=1 -!#ifdef CVF -! open(10,file=fname,form='binary',status='old',err=998) + call rfile3a(fname,id,n,ierr) if(ierr.ne.0) then print*,'Error opening or reading file: ',fname,ierr go to 999 endif -!#else -! call rfile2(fname,id,n,nr) -! if(nr.ne.n) then -! print*,'Error opening or reading file: ',fname,n,nr -! ierr=1002 -! go to 999 -! endif -! -!#endif sq=0. ka=0.1*NSMAX diff --git a/map65.py b/map65.py index 631af4923..ff2489229 100644 --- a/map65.py +++ b/map65.py @@ -1,4 +1,4 @@ -#--------------------------------------------------------------------- MAP65 +#-------------------------------------------------------------------- MAP65 # $Date$ $Revision$ # from Tkinter import * diff --git a/sec_midn.F90 b/sec_midn.F90 index 2a919d000..3af19984c 100644 --- a/sec_midn.F90 +++ b/sec_midn.F90 @@ -1,14 +1,5 @@ real function sec_midn() - -real*8 sec8,hrtime - -#ifdef CVF sec_midn=secnds(0.0) -#else - sec8=hrtime() - sec_midn=mod(sec8,86400.d0) -#endif - return end function sec_midn diff --git a/spec.f90 b/spec.f90 index 40f399453..db5f34e9c 100644 --- a/spec.f90 +++ b/spec.f90 @@ -89,9 +89,13 @@ subroutine spec(brightness,contrast,ngain,nspeed,a,a2) a2(i)=a2(i-newpts) enddo - gain=40*sqrt(nstep(nspeed)/5.0) * 5.0**(0.01*contrast) + logmap=1 gamma=1.3 + 0.01*contrast offset=(brightness+64.0)/2 + if(logmap.eq.1) then + gain=40*sqrt(nstep(nspeed)/5.0) * 5.0**(0.01*contrast) + offset=brightness/2 + 10 + endif fac=20.0/nadd fac=fac*0.065/base ! fac=fac*(0.1537/base) @@ -120,14 +124,16 @@ subroutine spec(brightness,contrast,ngain,nspeed,a,a2) x=max(x,s(ii,j)) enddo x=fac*x - if(x.gt.0.0) n=(2.0*x)**gamma + offset + if(x.gt.0.0 .and. logmap.eq.0) n=(2.0*x)**gamma + offset + if(x.gt.0.0 .and. logmap.eq.1) n=gain*log10(1.0*x) + offset n=min(252,max(0,n)) a(k)=n ! Now do the lower (zoomed) waterfall with one FFT bin per pixel. n=0 x=fac*s(i0+i-1,j) - if(x.gt.0.0) n=(3.0*x)**gamma + offset + if(x.gt.0.0 .and. logmap.eq.0) n=(3.0*x)**gamma + offset + if(x.gt.0.0 .and. logmap.eq.1) n=1.2*gain*log10(1.0*x) + offset n=min(252,max(0,n)) a2(k)=n diff --git a/wsjtgen.F90 b/wsjtgen.F90 index 225d5eddd..8e8b6fcbb 100644 --- a/wsjtgen.F90 +++ b/wsjtgen.F90 @@ -64,6 +64,9 @@ subroutine wsjtgen testfile=msg(2:) #ifdef CVF open(18,file=testfile,form='binary',status='old',err=12) +#else + open(18,file=testfile,access='stream',status='old',err=12) +#endif go to 14 12 print*,'Cannot open test file ',msg(2:) go to 999 @@ -73,17 +76,17 @@ subroutine wsjtgen close(18) if(ierr.ne.0) print*,'Error reading test file ',msg(2:) -#else - tfile2=testfile - call rfile2(tfile2,hdr,44+2*661500,nr) - if(nr.le.0) then - print*,'Error reading ',testfile - stop - endif - do i=1,ndata/2 - iwave(i)=jwave(i) - enddo -#endif +!#else +! tfile2=testfile +! call rfile2(tfile2,hdr,44+2*661500,nr) +! if(nr.le.0) then +! print*,'Error reading ',testfile +! stop +! endif +! do i=1,ndata/2 +! iwave(i)=jwave(i) +! enddo +!#endif nwave=ndata/2 do i=nwave,NTXMAX iwave(i)=0