MACRO viewmaps weighted=0 scale=1 * displays the reconstruction maps GRAPHICS/MISC/CLR SET '*' OPTION '*' ZONE 1 1 SET YGTI 1.22 EXEC bluered 32 OPTION LINY OPTION LINX OPTION DATE OPTION ZFL1 dir1 = '//LUN1' READ dir1 'Enter first hbook-dir' IF [ncols] = 0 THEN plottitle = 'G?en!-RC Coeff."J# ' ELSE plottitle = 'G?en!-RC Coeff.(weighted)"J# ' ENDIF READ num1 'Enter theta [deg]' IF [num1] <> '' THEN plottitle = [plottitle]//' [q]='//[num1]//'^o!]' ENDIF READ num2 'Enter x [cm]' IF [num2] <> '' THEN plottitle = [plottitle]//' x='//[num2]//'cm' ENDIF READ num3 'Enter y [cm]' IF [num3] <> '' THEN plottitle = [plottitle]//' y='//[num3]//'cm' ENDIF TITLE [plottitle] * plot the matrix coeffs EXEC viewmap [dir1]//'/10' [weighted] [scale] PICTURE/PRINT paw1.ps READ conf 'Press CR (1= Ok to all)' EXEC viewmap [dir1]//'/20' [weighted] [scale] PICTURE/PRINT paw2.ps IF [conf]<>1 THEN READ conf 'Press CR' ENDIF EXEC viewmap [dir1]//'/30' [weighted] [scale] PICTURE/PRINT paw3.ps IF [conf]<>1 THEN READ conf 'Press CR' ENDIF EXEC viewmap [dir1]//'/40' [weighted] [scale] PICTURE/PRINT paw4.ps RETURN MACRO viewmap id weighted=0 scale=1 * display a single map IF $HEXIST(100) = 1 THEN H/DEL 100 ENDIF IF [weighted] = 0 THEN H/COP [id] 100 ELSE H/OP/MUL [id] 1 100 1. 1. ENDIF a = $HINFO(100,'max')*[scale] b = $HINFO(100,'min')*[scale] IF ([a] < -[b]) THEN a = -[b] ENDIF MIN 100 -[a] MAX 100 [a] HISTO/PLOT 100 COLZ RETURN MACRO bluered ncols=16 * create a palette from blue->white->red k = [ncols]*2+8 SET NCOL [k] PALETTE 1 k = 8+[ncols] COLOR [k] 1 1 1 DO i=1,[ncols] i1 = [k]-[i] i2 = [k]+[i] rb = (1-[i]/[ncols])*0.9 COLOR [i1] [rb] [rb] 1 COLOR [i2] 1 [rb] [rb] ENDDO RETURN