1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
|
subroutine pltcog
c******************************************************************************
c This subroutine controls the decisions that are made around the
c plots of curves-of-growth
c******************************************************************************
implicit real*8 (a-h,o-z)
include 'Atmos.com'
include 'Pstuff.com'
real*8 xnum
real*4 shortnum
c call up the curve-of-growth plot
if (plotopt == 0) return
10 choice = 'y'
plotroutine = 'term_land_cog '
lscreen = 12
whichwin = '1of1'
call makeplot (lscreen)
c make a hardcopy, write to a postscript file, point at a place on the
c screen, try a new model atmosphere, or replot?
1 array = 'WHAT TO DO NEXT ([n]/h/f/p/v/m/r)? '
lscreen = 12
nchars = 37
call getasci (nchars,lscreen)
choice = chinfo(1:1)
if (choice=='n' .or. nchars<=0) then
return
elseif (choice == 'h') then
plotroutine = 'hard_land_cog '
call makeplot (lscreen)
go to 10
elseif (choice == 'v') then
write (array,*) 'What is the new microturbulence (km/s)? '
nchars = 41
lscreen = lscreen + 2
call getnum (nchars,lscreen,xnum,shortnum)
do i=1,ntau
vturb(i) = xnum*1.0e5
enddo
write (moditle(57:64),1010) xnum
rewind nf1out
rewind nf2out
rewind nfmodel
rewind nflines
return
elseif (choice == 'm') then
return
elseif (choice == 'r') then
go to 10
elseif (choice == 'p') then
array = 'MARK THE POSITION WITH THE CURSOR'
istat=ivcleof(21,1)
istat=ivwrite(13,3,array,34)
call drawcurs
go to 1
elseif (choice == 'f') then
plotroutine = 'file_land_cog '
call makeplot (lscreen)
go to 10
endif
return
c*****format statements
1010 format ('vt=',f5.2)
end
|
