--- ray/src/cv/bsdf2rad.cal 2017/04/09 21:39:26 2.1 +++ ray/src/cv/bsdf2rad.cal 2017/04/10 06:09:14 2.2 @@ -1,4 +1,4 @@ -{ RCSid $Id: bsdf2rad.cal,v 2.1 2017/04/09 21:39:26 greg Exp $ } +{ RCSid $Id: bsdf2rad.cal,v 2.2 2017/04/10 06:09:14 greg Exp $ } { Calculation of false color scale for bsdf2rad. @@ -6,11 +6,15 @@ } { Compute lines of latitude and longitude } LATSTEP : 15*DEGREE; -LONGSTEP : 30*DEGREE; -LWIDTH : 1*DEGREE; +LONGSTEP : 30*DEGREE; { must divide evenly into 90 degrees } +LATWIDTH : 0.7*DEGREE; +LONGWIDTH : 1.*DEGREE; -latlong = if(LWIDTH-mod(Acos(Nz)+LWIDTH/2,LATSTEP), .01, - if(LWIDTH-mod(atan2(Ny,Nx)+LWIDTH/2,LONGSTEP), .01, 1)); +mylong = atan2(Ny,Nx); + +latlong = if(LATWIDTH-mod(Acos(Nz)+ LATWIDTH/2,LATSTEP), .05, + if(LONGWIDTH-mod(mylong+LONGWIDTH/2,90*DEGREE), .01, + if(LONGWIDTH-mod(mylong+LONGWIDTH/2,LONGSTEP), .2, 1))); 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));