--- ray/src/px/falsecolor.csh 2005/11/12 06:14:51 2.13 +++ ray/src/px/falsecolor.csh 2005/11/18 22:17:16 2.14 @@ -1,5 +1,5 @@ #!/bin/csh -fe -# RCSid: $Id: falsecolor.csh,v 2.13 2005/11/12 06:14:51 greg Exp $ +# RCSid: $Id: falsecolor.csh,v 2.14 2005/11/18 22:17:16 greg Exp $ # # Create false color image with legend # @@ -115,22 +115,22 @@ boundary(a,b) : neq(floor(ndivs*a+.5),floor(ndivs*b+.5 interp_arr2(i,x,f):(i+1-x)*f(i)+(x-i)*f(i+1); interp_arr(x,f):if(x-1,if(f(0)-x,interp_arr2(floor(x),x,f),f(f(0))),f(1)); def_redp(i):select(i,0.18848,0.05468174, -0.00103547,0,7.449763e-06,0.0004390987,0.001367254,0.003076, -0.01376382,0.06085327,0.1712073,0.2889212,0.3404111,0.3803934, -0.4135148,0.4461151,0.4966744,0.6343959,0.7709373,0.9048474,1, -1,0.9863); +0.00103547,8.311144e-08,7.449763e-06,0.0004390987,0.001389302, +0.003306035,0.01381413,0.0563861,0.1534137,0.2581243,0.3059902, +0.3449827,0.3791056,0.4138921,0.4678283,0.624231,0.7732902, +0.9067587,1,1,0.9863); def_red(x):interp_arr(x/0.0454545+1,def_redp); -def_grnp(i):select(i,0.0009766,0,0.0008966244, -0.0264977,0.1256843,0.2865799,0.4251992,0.4775091,0.4406337, -0.3546771,0.2374256,0.1424751,0.09136063,0.0474073,0.0181899, -0.004181549,9.587855e-05,0,0.008718861,0.1025931,0.3106831, -0.6447838,0.9707); +def_grnp(i):select(i,0.0009766,2.35501e-05, +0.0008966244,0.0264977,0.1256843,0.2865799,0.424867,0.4750973, +0.4403898,0.3568169,0.244569,0.1643434,0.1297663,0.1001566, +0.07522873,0.05465921,0.03394894,0.001412969,0.00681687, +0.0974651,0.3106831,0.6447838,0.9707); def_grn(x):interp_arr(x/0.0454545+1,def_grnp); def_blup(i):select(i,0.2666,0.3638662,0.4770437, 0.5131397,0.5363797,0.5193677,0.4091852,0.1732726,0.05331749, -0.0457587,0.06656252,0.06819438,0.04973373,0.02796644, -0.01131437,0.00280772,0.0002246617,0.0003480934,0.001044561, -0.001006845,0.002297005,0.01128842,0.02539); +0.04184916,0.05768565,0.070331,0.0738812,0.07805368,0.08024318, +0.07775604,0.06449649,0.01028498,0.0003090336,0.0004275006, +0.002297005,0.01128842,0.02539); def_blu(x):interp_arr(x/0.0454545+1,def_blup); isconta = if(btwn(0,v,1),or(boundary(vleft,vright),boundary(vabove,vbelow)),-1);