Save and restore the "fac0" scaling value.

This commit is contained in:
Joe Taylor 2024-01-05 07:01:55 -05:00
parent c733c0c123
commit 9e0425ef40
2 changed files with 10 additions and 8 deletions

View File

@ -5,17 +5,19 @@ subroutine read_qm(fname)
character*(*) fname
character prog_id*24,mycall*12,mygrid*6
real*8 fcenter
integer nxtra(16) !For possible future additions
integer nxtra(15) !For possible future additions
integer*1 id1(2,NMAX)
common/datcom/dd(2,5760000),ss(400,NFFT),savg(NFFT),fcenter,nutc,junk(NJUNK)
open(28,file=trim(fname),status='old',access='stream',err=900)
read(28) prog_id,mycall,mygrid,fcenter,nutc,ntx30a,ntx30b,ndop00,ndop58, &
ia,ib,nxtra
ia,ib,fac0,nxtra
fac=1.0
if(fac0.gt.0.0) fac=1.0/fac0
id1=0
read(28) id1(1:2,ia:ib)
dd=0.
dd(1:2,ia:ib)=2.4*id1(1:2,ia:ib) !### Why the boost by 2.4 ???
dd(1:2,ia:ib)=fac*id1(1:2,ia:ib) !Boost back to previous level
! write(*,3001) prog_id,mycall(1:6),mygrid,fcenter,nutc,ntx30a,ntx30b, &
! ndop00,ndop58,ia,ib

View File

@ -6,7 +6,7 @@ subroutine save_qm(fname,prog_id,mycall,mygrid,dd,ntx30a,ntx30b,fcenter, &
character prog_id_24*24,mycall_12*12,mygrid_6*6
real*4 dd(2,NMAX)
real*8 fcenter
integer nxtra(16) !For possible future additions
integer nxtra(15) !For possible future additions
integer*1 id1(2,NMAX)
ia=1
@ -25,10 +25,10 @@ subroutine save_qm(fname,prog_id,mycall,mygrid,dd,ntx30a,ntx30b,fcenter, &
nbad=0
dmax=0.
fac=10.0/rms
fac0=10.0/rms
do i=ia,ib
x=fac*dd(1,i)
y=fac*dd(2,i)
x=fac0*dd(1,i)
y=fac0*dd(2,i)
ax=abs(x)
ay=abs(y)
dmax=max(dmax,ax,ay)
@ -57,7 +57,7 @@ subroutine save_qm(fname,prog_id,mycall,mygrid,dd,ntx30a,ntx30b,fcenter, &
mycall_12=mycall
mygrid_6=mygrid
write(29) prog_id_24,mycall_12,mygrid_6,fcenter,nutc,ntx30a,ntx30b, &
ndop00,ndop58,ia,ib,nxtra !Write header to disk
ndop00,ndop58,ia,ib,fac0,nxtra !Write header to disk
write(29) id1(1:2,ia:ib) !Write 8-bit data to disk
close(29)