linse/sics
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

new versions

Browse Source
This commit is contained in:
cvs
2002-01-24 13:24:43 +00:00
parent c7347e4129
commit e3c9ba56e6
4 changed files with 303 additions and 154 deletions
Download Patch File Download Diff File Expand all files Collapse all files

35
tecs/tecs_plot.f90
View File

@@ -171,7 +171,7 @@ subroutine tecs_plot(auxpar)
sel2=sel
endif
do rl=1,2
if (mode==zoom) then
if (mode==zoom .and. rl==1) then
ymin(1)=y1
ymax(1)=y2
else
@@ -361,7 +361,8 @@ subroutine tecs_plot(auxpar)
numl=0
numb=' '
7 if (mode==live) then
7 ex=undef
if (mode==live) then
! if (device(1:1)=='X') then
! call pgmtxt('T', 0.5, 0.0, 0.0, 'LIVE MODE (click on text window before pressing any further key)')
! endif
@@ -369,7 +370,7 @@ subroutine tecs_plot(auxpar)
i=chartperiod-mod(myc_now(), chartperiod)
call get_cursor(ex, ey, key, -i)
do while (key == char(0)) ! no key pressed
do while (key == char(0) .or. key == ' ') ! no key or space key pressed
iret=tecs_get_mult(pars, t, nset, yy1)
if (iret<0) goto 99
tim1=t-tbase
@@ -428,7 +429,7 @@ subroutine tecs_plot(auxpar)
if (mode==zoom) then
if (ex==undef) then
ex=(x1+x2)/2
ey=(min(y1,ymax(1))+max(y1,ymin(1)))/2
ey=(min(y2,ymax(1))+max(y1,ymin(1)))/2
end if
fx=max(winmin,x2-x1)
fy=max(y2-y1,1e-3,y2*1e-5)
@@ -532,6 +533,7 @@ subroutine tecs_plot(auxpar)
elseif (key == 'S') then
sel=sel+1
if (sel > tmax) sel=0
if (mode==zoom) mode=0
elseif (key == 'C') then
auxsel=auxsel+1
if (auxsel > amax) auxsel=1
@@ -553,18 +555,16 @@ subroutine tecs_plot(auxpar)
read(*,'(a)') filnam
open(lund, file=filnam, status='unknown', carriagecontrol='list')
gap=.false.
do j=1,nset
if (j == 1) then
write(lund, *) ' time [h]',char(9),' Tm [K]'
elseif (j == 2) then
write(lund,*)
write(lund, *) ' time [h]',char(9),' Ts [K]'
l=0
i2=0
do i1=1,nset
if (i1 > nset-naux) then
j=i1
else
write(lund,*)
write(lund, *) ' time [h]',char(9), ' ',trim(parnam(j))
j=nset-i1+(1-naux)
endif
gap=.false.
do i=1,retLen(j)
if (yd(i,j)==undef) then
if (gap) then
@@ -572,13 +572,20 @@ subroutine tecs_plot(auxpar)
gap=.false.
endif
else
if (i2/=i1) then
if (l > 0) write(lund, *)
write(lund, *) ' time [h]',char(9), ' ',trim(parnam(j))
l=l+1
i2=i1
endif
write(lund, '(f9.4,a,f9.4)') xd(i,j)/3600., char(9), max(-999.,min(9999.,yd(i,j)))
l=l+1
gap=.true.
endif
enddo
enddo
close(lund)
print *, retLen(1)+1, ' lines written to ',filnam(1:48)
print *, l, ' lines written to ',filnam(1:48)
endif
end subroutine