--- ray/src/px/falsecolor.csh 2005/11/18 22:17:16 2.14 +++ ray/src/px/falsecolor.csh 2005/11/18 22:50:21 2.15 @@ -1,5 +1,5 @@ #!/bin/csh -fe -# RCSid: $Id: falsecolor.csh,v 2.14 2005/11/18 22:17:16 greg Exp $ +# RCSid: $Id: falsecolor.csh,v 2.15 2005/11/18 22:50:21 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,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); +0.00103547,8.311144e-08,7.449763e-06,0.0004390987,0.001367254, +0.003076,0.01376382,0.06170773,0.1739422,0.2881156,0.3299725, +0.3552663,0.372552,0.3921184,0.4363976,0.6102754,0.7757267, +0.9087369,1,1,0.9863); def_red(x):interp_arr(x/0.0454545+1,def_redp); 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); +0.0008966244,0.0264977,0.1256843,0.2865799,0.4247083,0.4739468, +0.4402732,0.3671876,0.2629843,0.1725325,0.1206819,0.07316644, +0.03761026,0.01612362,0.004773749,6.830967e-06,0.00803605, +0.1008085,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.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); +0.5131397,0.5363797,0.5193677,0.4085123,0.1702815,0.05314236, +0.05194055,0.08564082,0.09881395,0.08324373,0.06072902, +0.0391076,0.02315354,0.01284458,0.005184709,0.001691774, +2.432735e-05,1.212949e-05,0.006659406,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);