macro give_specs * * vec/crea hmm(2) r 0.001 1.2 ! can change it if you like - it is max and min of spectra vec/crea ncalc(1) r 4 ncalcs = $SIGMA(ncalc(1)) vec/crea calcs([ncalcs]) r 1 7 8 9 * * following are the x- and y-coordinates and angles of the labels "pion" (index=1) and "proton" (index=2) vec/crea xlab(2) r 1 1.4 vec/crea ylab(2) r 0.03 0.2 vec/crea alab(2) r -40 -30 * vec/crea vykey(1) r 0.0045 return * *=================================================== * macro make_key icalc lt ykey = $SIGMA(vykey(1)) * vec/rea params '../../../DataFiles/Calc'//[icalc]//'/input.inp' *** index below should be 1=Ry 2=Rx 3=T 4=rho0 5=rho2 6=tau0 7=a_s 8=DeltaTau pval1 = $SIGMA(params(7)) pval2 = $SIGMA(params(4)) if ([pval1].lt.0.05) then key 0.23 [ykey] [lt] 'a?S!=0.0 [r]?0! = '//[pval2] 1.2 L else key 0.23 [ykey] [lt] 'a?S!='//[pval1]//' [r]?0! = '//[pval2] 1.2 L endif * vec/in vykey(1) $SIGMA(vykey(1)/1.45) return *===================================================