volume
KTlensflareV
(
string     fog="off";
float       distance=1500;
color      background=color(0.94,0.92,0.97);

/**   big flare paramater  **/
/*point       SunP=point(-0.6,0.6,3.464);*/
point       SunP=point(-0.42,0.73,3.464);
float       fov=60;
float       ringS=0.3,ringL1=0.2,ringL2=0.4,ringL3=0.1,ringL4=0.1,ringL5=0.8,ringL6=0.15;
float       ringS_s=0.27,ringL1_s=0.13,ringL2_s=0.2;
float       max_ring1=0.98,max_ring2=0.98,max_ring3=0.98,max_ring4=0.98,max_ring5=0.98,max_ring6=0.98;
float       max_ringS_s=0.6,max_ring1_s=0.98,max_ring2_s=0.98;
float       max_mixS=0.1,max_mix1=0.025,max_mix2=0.025,max_mix3=0.05,max_mix4=0.05;
float       max_mix5=0.015,max_mix6=0.1;
float       max_mixS_s=0.3,max_mix1_s=0.1,max_mix2_s=0.05;
string     typeS="clamp",type1="step",type2="step",type3="clamp",type4="clamp",type5="clamp",type6="clamp";
string     typeS_s="step",type1_s="step",type2_s="step";
color      csS=color(1,0.25,0),cs1=color(0.74,1,0.87),cs2=color(1,1,0),cs3=color(1,1,0.5);
color      cs4=color(1,1,0.5),cs5=color(1,1,1),cs6=color(0.74,1,0.87);
color      csS_s=color(0,0.5,1),cs1_s=color(1,0.9,0.9),cs2_s=color(1,1,0);

/**   small flare paramater  **/
float       length1=0.02,length3=0.2,length4=0.05,length5=0.02,length6=0.04;
float       ring_sL1_1=0.005,ring_sL1_2=0.02,ring_sL2=0.15;
float       ring_sL3_1=0.01,ring_sL3_2=0.01;
float       ring_sL4_1=0.06,ring_sL4_2=0.05,ring_sL4_3=0.06;
float       ring_sL5_1=0.03,ring_sL5_2=0.005;
float       ring_sL6_1=0.02,ring_sL6_2_1=0.03,ring_sL6_2_2=0.02,ring_sL6_3=0.02,ring_sL7=0.08;
float       max_ring_sL1_1=0.98,max_ring_sL1_2=0.98,max_ring_sL2=0.98;
float       max_ring_sL3_1=0.25,max_ring_sL3_2=0.25;
float       max_ring_sL4_1=0.98,max_ring_sL4_2=1,max_ring_sL4_3=0.98;
float       max_ring_sL5_1=0.8,max_ring_sL5_2=0.98;
float       max_ring_sL6_1=0.9,max_ring_sL6_2_1=1,max_ring_sL6_2_2=0.9,max_ring_sL6_3=0.9,max_ring_sL7=0.9;
string     type_s1_1="step",type_s1_2="step",type_s2="clamp";
string     type_s3_1="step",type_s3_2="step";
string     type_s4_1="clamp",type_s4_2="specular",type_s4_3="clamp";
string     type_s5_1="step",type_s5_2="step";
string     type_s6_1="step",type_s6_2_1="specular",type_s6_2_2="step",type_s6_3="step",type_s7="clamp";
float       max_mix_sL1_1=1.5,max_mix_sL1_2=0.69,max_mix_sL2=0.05;
float       max_mix_sL3_1=0.8,max_mix_sL3_2=0.5;
float       max_mix_sL4_1=0.05,max_mix_sL4_2=2,max_mix_sL4_3=0.05;
float       max_mix_sL5_1=0.69,max_mix_sL5_2=1.5;
float       max_mix_sL6_1=0.4,max_mix_sL6_2_1=4,max_mix_sL6_2_2=0.4,max_mix_sL6_3=0.4,max_mix_sL7=0.05;
color      cs_s1_1=color(1,1,1),cs_s1_2=color(0.74,1,0.87),cs_s2=color(0.74,0.9,1);
color      cs_s3_1=color(0.48,0.82,1),cs_s3_2=color(0.82,1,0);
color      cs_s4_1=color(0.7,1,0),cs_s4_2=color(0.85,1,0.8),cs_s4_3=color(0.7,1,0);
color      cs_s5_1=color(1,0.9,0.72),cs_s5_2=color(1,1,1);
color      cs_s6_1=color(0.6,0.82,1),cs_s6_2_1=color(0.8,0.85,1),cs_s6_2_2=color(0.6,0.82,1),cs_s6_3=color(0.6,0.82,1),cs_s7=color(0.92,1,0.92);
float       s_roughness=-5;

/*   highlight paramater  */
/*float       Highlight_roughness=10,Highlight_intensity=1.2;*/
float       Highlight_roughness=3,Highlight_intensity=1.05;
color      Highlightcolor=color(1,1,1);
float       lightscale=10,light_min=0.6,light_max=0.7;
float       light_roughness=12,light_intensity=0.4;
color      lightcolor=color(1,1,1)
)

{
float     fog_d;

point     norP,norSunP,norCenterP;
point     LensP1,LensP2,LensP3,LensP4,LensP5,LensP6;
point     DSunP,DLensP5;
point     SWP,SWSunP,SWLensP1,SWLensP2,SWLensP3,SWLensP4,SWLensP5,SWLensP6;
float     Pz,norSunPz,LensP1z,LensP2z,LensP3z,LensP4z,LensP5z,LensP6z;
float     judgeS,judgeLP1,judgeLP2,judgeLP3,judgeLP4,judgeLP5,judgeLP6;
float     r_unit,SWunit;
float     Dfloat;
color    ciS,ci1,ci2,ci3,ci4,ci5,ci6;
color    ciS_s,ci1_s,ci2_s;
float     r,dis;

point      SWsLensP1,SWsLensP2,SWsLensP3,SWsLensP4,SWsLensP5,SWsLensP6,SWsLensP7;
point      SWsLensP1_1,SWsLensP1_2,SWsLensP3_1,SWsLensP3_2,SWsLensP4_1,SWsLensP4_2,SWsLensP4_3;
point      SWsLensP5_1,SWsLensP5_2,SWsLensP6_1,SWsLensP6_2,SWsLensP6_3;
point      moveP_unit;
float      moveL_unit;
float      judge_sLP1_1,judge_sLP1_2,judge_sLP2,judge_sLP3_1,judge_sLP3_2;
float      judge_sLP4_1,judge_sLP4_2,judge_sLP4_3,judge_sLP5_1,judge_sLP5_2;
float      judge_sLP6_1,judge_sLP6_2,judge_sLP6_3,judge_sLP7;
color     ci_s1_1,ci_s1_2,ci_s2,ci_s3_1,ci_s3_2,ci_s4_1,ci_s4_2,ci_s4_3;
color     ci_s5_1,ci_s5_2,ci_s6_1,ci_s6_2_1,ci_s6_2_2,ci_s6_3,ci_s7;
color     ci_s1,ci_s3,ci_s4,ci_s5,ci_s6;

float     Highlight;
point     LightP,norLightP;
float     PorM;
float     Cosin;
float     lightdis,lighttheta;


if (fog=="on") {
     fog_d=1-exp(-length(I)/distance);
     Ci=Oi*mix(Ci,background,fog_d);
}

norP=normalize(P);
/*norSunP=transform("world",SunP);
*norSunP=normalize(norSunP);*/
norSunP=normalize(SunP);
norCenterP=point(0,0,1);
norCenterP=normalize(norCenterP);


/*  flare  BEGIN  */
LensP1=reflect(-norSunP,norCenterP);
LensP1=normalize(reflect(-norCenterP,LensP1));
/*printf ("LensP1 = %p\n",LensP1);*/
LensP2=normalize(LensP1+norSunP);
LensP3=normalize(LensP1+LensP2);
LensP4=normalize(norSunP+LensP2);
LensP5=normalize(reflect(-norCenterP,norSunP));
DSunP=norSunP;
DLensP5=LensP5;
Dfloat=(1/3)/zcomp(DSunP);
DSunP=Dfloat*DSunP;
Dfloat=(2/3)/zcomp(DLensP5);
DLensP5=Dfloat*DLensP5;
LensP6=normalize(DSunP+DLensP5);

norSunPz=zcomp(norSunP);
SWunit=(1/norSunPz);
SWSunP=SWunit*norSunP;

LensP1z=zcomp(LensP1);
SWunit=(1/LensP1z);
SWLensP1=SWunit*LensP1;

LensP2z=zcomp(LensP2);
SWunit=(1/LensP2z);
SWLensP2=SWunit*LensP2;

LensP3z=zcomp(LensP3);
SWunit=(1/LensP3z);
SWLensP3=SWunit*LensP3;

LensP4z=zcomp(LensP4);
SWunit=(1/LensP4z);
SWLensP4=SWunit*LensP4;

LensP5z=zcomp(LensP5);
SWunit=(1/LensP5z);
SWLensP5=SWunit*LensP5;

LensP6z=zcomp(LensP6);
SWunit=(1/LensP6z);
SWLensP6=SWunit*LensP6;

Pz=zcomp(P);
SWunit=(1/Pz);
SWP=SWunit*P;

r_unit=tan(radians(fov/2));

judgeS=distance(SWP,SWSunP)/r_unit;
judgeLP1=distance(SWP,SWLensP1)/r_unit;
judgeLP2=distance(SWP,SWLensP2)/r_unit;
judgeLP3=distance(SWP,SWLensP3)/r_unit;
judgeLP4=distance(SWP,SWLensP4)/r_unit;
judgeLP5=distance(SWP,SWLensP5)/r_unit;
judgeLP6=distance(SWP,SWLensP6)/r_unit;



if (judgeS<ringS) {
    r=judgeS/ringS;
    dis=r;
    r=smoothstep(ringS_s/ringS,(((ringS-ringS_s)/3)+ringS_s)/ringS,dis)-smoothstep((ringS-((ringS-ringS_s)/3))/ringS,1,dis);
    r*=max_mixS;
    ciS=mix(Ci,csS,r);
    Ci=ciS;
}
     
     if (judgeS<ringS_s) {
          r=judgeS/ringS_s;
          dis=r;
          r=clamp(r,0,max_ringS_s);
          r/=max_ringS_s;
          if (typeS_s=="clamp") {
                  r=r-smoothstep(max_ringS_s,1,dis);
          }
          if (typeS_s=="sine") {
                  r=sin(r*(PI/2));
                  r=r-smoothstep(max_ringS_s,1,dis);
          }
          if (typeS_s=="step") {
                  r=1-smoothstep(max_ringS_s,1,dis);
          }
          r*=max_mixS_s;
          ciS_s=mix(Ci,csS_s,r);
          ciS=(ciS+ciS_s)/2;
          Ci=ciS;
     }     

if (judgeLP1<ringL1) {
    r=judgeLP1/ringL1;
    dis=r;
    r=clamp(r,0,max_ring1);
    r/=max_ring1;
    if (type1=="clamp") {
             r=r-smoothstep(max_ring1,1,dis);
     }
     if (type1=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring1,1,dis);
     }
     if (type1=="step") {
             r=1-smoothstep(max_ring1,1,dis);
     }
     r*=max_mix1;
     ci1=mix(Ci,cs1,r);
     Ci=ci1;
     
     if (judgeLP1<ringL1_s) {
          r=judgeLP1/ringL1_s;
          dis=r;
          r=clamp(r,0,max_ring1_s);
          r/=max_ring1_s;
          if (type1_s=="clamp") {
                  r=r-smoothstep(max_ring1_s,1,dis);
          }
          if (type1_s=="sine") {
                  r=sin(r*(PI/2));
                  r=r-smoothstep(max_ring1_s,1,dis);
          }
          if (type1_s=="step") {
                  r=1-smoothstep(max_ring1_s,1,dis);
          }
          r*=max_mix1_s;
          ci1_s=mix(Ci,cs1_s,r);
/*          ci1=(ci1+ci1_s)/2;*/
         Ci=ci1_s;
     }     
}

if (judgeLP2<ringL2) {
    r=judgeLP2/ringL2;
    dis=r;
    r=clamp(r,0,max_ring2);
    r/=max_ring2;
    if (type2=="clamp") {
             r=r-smoothstep(max_ring2,1,dis);
     }
     if (type2=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring2,1,dis);
     }
     if (type2=="step") {
             r=1-smoothstep(max_ring2,1,dis);
     }
     r*=max_mix2;
     ci2=mix(Ci,cs2,r);
     Ci=ci2;

     if (judgeLP2<ringL2_s) {
          r=judgeLP2/ringL2_s;
          dis=r;
          r=clamp(r,0,max_ring2_s);
          r/=max_ring2_s;
          if (type2_s=="clamp") {
                  r=r-smoothstep(max_ring2_s,1,dis);
          }
          if (type2_s=="sine") {
                  r=sin(r*(PI/2));
                  r=r-smoothstep(max_ring2_s,1,dis);
          }
          if (type2_s=="step") {
                  r=1-smoothstep(max_ring2_s,1,dis);
          }
          r*=max_mix2_s;
          ci2_s=mix(Ci,cs2_s,r);
/*          ci2=(ci2+ci2_s)/2;*/
          Ci=ci2_s;
     }     
}

if (judgeLP3<ringL3) {
    r=judgeLP3/ringL3;
    dis=r;
    r=clamp(r,0,max_ring3);
    r/=max_ring3;
    if (type3=="clamp") {
             r=r-smoothstep(max_ring3,1,dis);
     }
     if (type3=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring3,1,dis);
     }
     if (type3=="step") {
             r=1-smoothstep(max_ring3,1,dis);
     }
     r*=max_mix3;
     ci3=mix(Ci,cs3,r);
     Ci=ci3;
}

if (judgeLP4<ringL4) {
    r=judgeLP4/ringL4;
    dis=r;
    r=clamp(r,0,max_ring4);
    r/=max_ring4;
    if (type4=="clamp") {
             r=r-smoothstep(max_ring4,1,dis);
     }
     if (type4=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring4,1,dis);
     }
     if (type4=="step") {
             r=1-smoothstep(max_ring4,1,dis);
     }
     r*=max_mix4;
     ci4=mix(Ci,cs4,r);
     Ci=ci4;
}

if (judgeLP5<ringL5) {
    r=judgeLP5/ringL5;
    dis=r;
    r=clamp(r,0,max_ring5);
    r/=max_ring5;
    if (type5=="clamp") {
             r=r-smoothstep(max_ring5,1,dis);
     }
     if (type5=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring5,1,dis);
     }
     if (type5=="step") {
             r=1-smoothstep(max_ring5,1,dis);
     }
     r*=max_mix5;
     ci5=mix(Ci,cs5,r);
     Ci=ci5;
}

if (judgeLP6<ringL6) {
    r=judgeLP6/ringL6;
    dis=r;
    r=clamp(r,0,max_ring6);
    r/=max_ring6;
    if (type6=="clamp") {
             r=r-smoothstep(max_ring6,1,dis);
     }
     if (type6=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring6,1,dis);
     }
     if (type6=="step") {
             r=1-smoothstep(max_ring6,1,dis);
     }
     r*=max_mix6;
     ci6=mix(Ci,cs6,r);
     Ci=ci6;
}

/*  flare END  */




/*  small flare  BEGIN  */
SWsLensP1=(SWSunP-(3/5*point(0,0,1)))*(5/2);
SWsLensP2=(1/4)*point(0,0,1)+(3/4)*SWSunP;
SWsLensP3=(3/4)*point(0,0,1)+(1/4)*SWSunP;
SWsLensP5=(point(0,0,1)-(1/2)*SWSunP)*2;
SWsLensP4=(1/2)*point(0,0,1)+(1/2)*SWsLensP5;
SWsLensP7=(point(0,0,1)-(2/3)*SWSunP)*3;
SWsLensP6=(1/4)*SWsLensP7+(3/4)*SWsLensP5;

/*moveP_unit=SWSunP-point(0,0,1);*/
moveP_unit=point(0,0,1)-SWSunP;
moveL_unit=length(moveP_unit);
moveP_unit=normalize(moveP_unit);

SWsLensP1_1=SWsLensP1;
SWsLensP1_2=SWsLensP1+moveL_unit*length1*moveP_unit;
SWsLensP3_1=SWsLensP3;
SWsLensP3_2=SWsLensP3+moveL_unit*length3*moveP_unit;
SWsLensP4_1=SWsLensP4;
SWsLensP4_2=SWsLensP4+moveL_unit*length4*moveP_unit;
SWsLensP4_3=SWsLensP4+2*moveL_unit*length4*moveP_unit;
SWsLensP5_1=SWsLensP5;
SWsLensP5_2=SWsLensP5+moveL_unit*length5*moveP_unit;
SWsLensP6_1=SWsLensP6;
SWsLensP6_2=SWsLensP6+moveL_unit*length6*moveP_unit;
SWsLensP6_3=SWsLensP6+2*moveL_unit*length6*moveP_unit;

judge_sLP1_1=distance(SWP,SWsLensP1_1)/r_unit;
judge_sLP1_2=distance(SWP,SWsLensP1_2)/r_unit;
judge_sLP2=distance(SWP,SWsLensP2)/r_unit;
judge_sLP3_1=distance(SWP,SWsLensP3_1)/r_unit;
judge_sLP3_2=distance(SWP,SWsLensP3_2)/r_unit;
judge_sLP4_1=distance(SWP,SWsLensP4_1)/r_unit;
judge_sLP4_2=distance(SWP,SWsLensP4_2)/r_unit;
judge_sLP4_3=distance(SWP,SWsLensP4_3)/r_unit;
judge_sLP5_1=distance(SWP,SWsLensP5_1)/r_unit;
judge_sLP5_2=distance(SWP,SWsLensP5_2)/r_unit;
judge_sLP6_1=distance(SWP,SWsLensP6_1)/r_unit;
judge_sLP6_2=distance(SWP,SWsLensP6_2)/r_unit;
judge_sLP6_3=distance(SWP,SWsLensP6_3)/r_unit;
judge_sLP7=distance(SWP,SWsLensP7)/r_unit;



if (judge_sLP1_1<ring_sL1_1) {
    r=judge_sLP1_1/ring_sL1_1;
    dis=r;
    r=clamp(r,0,max_ring_sL1_1);
    r/=max_ring_sL1_1;
    if (type_s1_1=="clamp") {
             r=r-smoothstep(max_ring_sL1_1,1,dis);
     }
     if (type_s1_1=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL1_1,1,dis);
     }
     if (type_s1_1=="step") {
             r=1-smoothstep(max_ring_sL1_1,1,dis);
     }
     if (type_s1_1=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL1_1;
     ci_s1_1=mix(Ci,cs_s1_1,r);
     
     Ci=ci_s1_1;
}     

if (judge_sLP1_2<ring_sL1_2) {
    r=judge_sLP1_2/ring_sL1_2;
    dis=r;
    r=clamp(r,0,max_ring_sL1_2);
    r/=max_ring_sL1_2;
    if (type_s1_2=="clamp") {
             r=r-smoothstep(max_ring_sL1_2,1,dis);
     }
     if (type_s1_2=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL1_2,1,dis);
     }
     if (type_s1_2=="step") {
             r=1-smoothstep(max_ring_sL1_2,1,dis);
     }
     if (type_s1_2=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL1_2;
     ci_s1_2=mix(Ci,cs_s1_2,r);
     
     Ci=ci_s1_2;
}


if (judge_sLP2<ring_sL2) {
    r=judge_sLP2/ring_sL2;
    dis=r;
    r=clamp(r,0,max_ring_sL2);
    r/=max_ring_sL2;
    if (type_s2=="clamp") {
             r=r-smoothstep(max_ring_sL2,1,dis);
     }
     if (type_s2=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL2,1,dis);
     }
     if (type_s2=="step") {
             r=1-smoothstep(max_ring_sL2,1,dis);
     }
     if (type_s2=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL2;
     ci_s2=mix(Ci,cs_s2,r);

     Ci=ci_s2;
}


if (judge_sLP3_1<ring_sL3_1) {
    r=judge_sLP3_1/ring_sL3_1;
    dis=r;
    r=clamp(r,0,max_ring_sL3_1);
    r/=max_ring_sL3_1;
    if (type_s3_1=="clamp") {
             r=r-smoothstep(max_ring_sL3_1,1,dis);
     }
     if (type_s3_1=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL3_1,1,dis);
     }
     if (type_s3_1=="step") {
             r=1-smoothstep(max_ring_sL3_1,1,dis);
     }
     if (type_s3_1=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL3_1;
     ci_s3_1=mix(Ci,cs_s3_1,r);
     
     Ci=ci_s3_1;
}


if (judge_sLP3_2<ring_sL3_2) {
    r=judge_sLP3_2/ring_sL3_2;
    dis=r;
    r=clamp(r,0,max_ring_sL3_2);
    r/=max_ring_sL3_2;
    if (type_s3_2=="clamp") {
             r=r-smoothstep(max_ring_sL3_2,1,dis);
     }
     if (type_s3_2=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL3_2,1,dis);
     }
     if (type_s3_2=="step") {
             r=1-smoothstep(max_ring_sL3_2,1,dis);
     }
     if (type_s3_2=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL3_2;
     ci_s3_2=mix(Ci,cs_s3_2,r);
     
     Ci=ci_s3_2;
}

if (judge_sLP4_1<ring_sL4_1) {
    r=judge_sLP4_1/ring_sL4_1;
    dis=r;
    r=clamp(r,0,max_ring_sL4_1);
    r/=max_ring_sL4_1;
    if (type_s4_1=="clamp") {
             r=r-smoothstep(max_ring_sL4_1,1,dis);
     }
     if (type_s4_1=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL4_1,1,dis);
     }
     if (type_s4_1=="step") {
             r=1-smoothstep(max_ring_sL4_1,1,dis);
     }
     if (type_s4_1=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL4_1;
     ci_s4_1=mix(Ci,cs_s4_1,r);
     
     Ci=ci_s4_1;
}


if (judge_sLP4_3<ring_sL4_3) {
    r=judge_sLP4_3/ring_sL4_3;
    dis=r;
    r=clamp(r,0,max_ring_sL4_3);
    r/=max_ring_sL4_3;
    if (type_s4_3=="clamp") {
             r=r-smoothstep(max_ring_sL4_3,1,dis);
     }
     if (type_s4_3=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL4_3,1,dis);
     }
     if (type_s4_3=="step") {
             r=1-smoothstep(max_ring_sL4_3,1,dis);
     }
     if (type_s4_3=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL4_3;
     ci_s4_3=mix(Ci,cs_s4_3,r);
     
     Ci=ci_s4_3;
}


if (judge_sLP4_2<ring_sL4_2) {
    r=judge_sLP4_2/ring_sL4_2;
    dis=r;
    r=clamp(r,0,max_ring_sL4_2);
    r/=max_ring_sL4_2;
    if (type_s4_2=="clamp") {
             r=r-smoothstep(max_ring_sL4_2,1,dis);
     }
     if (type_s4_2=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL4_2,1,dis);
     }
     if (type_s4_2=="step") {
             r=1-smoothstep(max_ring_sL4_2,1,dis);
     }
     if (type_s4_2=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL4_2;
     ci_s4_2=mix(Ci,cs_s4_2,r);
     
     Ci=ci_s4_2;
}


if (judge_sLP5_1<ring_sL5_1) {
    r=judge_sLP5_1/ring_sL5_1;
    dis=r;
    r=clamp(r,0,max_ring_sL5_1);
    r/=max_ring_sL5_1;
    if (type_s5_1=="clamp") {
             r=r-smoothstep(max_ring_sL5_1,1,dis);
     }
     if (type_s5_1=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL5_1,1,dis);
     }
     if (type_s5_1=="step") {
             r=1-smoothstep(max_ring_sL5_1,1,dis);
     }
     if (type_s5_1=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL5_1;
     ci_s5_1=mix(Ci,cs_s5_1,r);
     
     Ci=ci_s5_1;
}
     

if (judge_sLP5_2<ring_sL5_2) {
    r=judge_sLP5_2/ring_sL5_2;
    dis=r;
    r=clamp(r,0,max_ring_sL5_2);
    r/=max_ring_sL5_2;
    if (type_s5_2=="clamp") {
             r=r-smoothstep(max_ring_sL5_2,1,dis);
     }
     if (type_s5_2=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL5_2,1,dis);
     }
     if (type_s5_2=="step") {
             r=1-smoothstep(max_ring_sL5_2,1,dis);
     }
     if (type_s5_2=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL5_2;
     ci_s5_2=mix(Ci,cs_s5_2,r);
     
     Ci=ci_s5_2;
}


if (judge_sLP6_1<ring_sL6_1) {
    r=judge_sLP6_1/ring_sL6_1;
    dis=r;
    r=clamp(r,0,max_ring_sL6_1);
    r/=max_ring_sL6_1;
    if (type_s6_1=="clamp") {
             r=r-smoothstep(max_ring_sL6_1,1,dis);
     }
     if (type_s6_1=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL6_1,1,dis);
     }
     if (type_s6_1=="step") {
             r=1-smoothstep(max_ring_sL6_1,1,dis);
     }
     if (type_s6_1=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL6_1;
     ci_s6_1=mix(Ci,cs_s6_1,r);
     
     Ci=ci_s6_1;
}


if (judge_sLP6_3<ring_sL6_3) {
    r=judge_sLP6_3/ring_sL6_3;
    dis=r;
    r=clamp(r,0,max_ring_sL6_3);
    r/=max_ring_sL6_3;
    if (type_s6_3=="clamp") {
             r=r-smoothstep(max_ring_sL6_3,1,dis);
     }
     if (type_s6_3=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL6_3,1,dis);
     }
     if (type_s6_3=="step") {
             r=1-smoothstep(max_ring_sL6_3,1,dis);
     }
     if (type_s6_3=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL6_3;
     ci_s6_3=mix(Ci,cs_s6_3,r);
     
     Ci=ci_s6_3;
}


if (judge_sLP6_2<ring_sL6_2_1) {
    r=judge_sLP6_2/ring_sL6_2_1;
    dis=r;
    r=clamp(r,0,max_ring_sL6_2_1);
    r/=max_ring_sL6_2_1;
    if (type_s6_2_1=="clamp") {
             r=r-smoothstep(max_ring_sL6_2_1,1,dis);
     }
     if (type_s6_2_1=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL6_2_1,1,dis);
     }
     if (type_s6_2_1=="step") {
             r=1-smoothstep(max_ring_sL6_2_1,1,dis);
     }
     if (type_s6_2_1=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL6_2_1;
     ci_s6_2_1=mix(Ci,cs_s6_2_1,r);
     
     Ci=ci_s6_2_1;
}


if (judge_sLP6_2<ring_sL6_2_2) {
    r=judge_sLP6_2/ring_sL6_2_2;
    dis=r;
    r=clamp(r,0,max_ring_sL6_2_2);
    r/=max_ring_sL6_2_2;
    if (type_s6_2_2=="clamp") {
             r=r-smoothstep(max_ring_sL6_2_2,1,dis);
     }
     if (type_s6_2_2=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL6_2_2,1,dis);
     }
     if (type_s6_2_2=="step") {
             r=1-smoothstep(max_ring_sL6_2_2,1,dis);
     }
     if (type_s6_2_2=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL6_2_2;
     ci_s6_2_2=mix(Ci,cs_s6_2_2,r);
     
     Ci=ci_s6_2_2;
}


if (judge_sLP7<ring_sL7) {
    r=judge_sLP7/ring_sL7;
    dis=r;
    r=clamp(r,0,max_ring_sL7);
    r/=max_ring_sL7;
    if (type_s7=="clamp") {
             r=r-smoothstep(max_ring_sL7,1,dis);
     }
     if (type_s7=="sine") {
             r=sin(r*(PI/2));
             r=r-smoothstep(max_ring_sL7,1,dis);
     }
     if (type_s7=="step") {
             r=1-smoothstep(max_ring_sL7,1,dis);
     }
     if (type_s7=="specular") {
             r=pow(1+dis,s_roughness);
     }
     r*=max_mix_sL7;
     ci_s7=mix(Ci,cs_s7,r);

     Ci=ci_s7;
}

/*  small flare  END  */



/*  highlight   BEGIN  */
Highlight=pow(judgeS+1,-(Highlight_roughness));
Highlight*=Highlight_intensity;
Ci=mix(Ci,Highlightcolor,Highlight);

LightP=SWP-SWSunP;
lightdis=length(LightP);
norLightP=normalize(LightP);
PorM=sign(ycomp(norLightP));
Cosin=norLightP.(point(1,0,0));
lighttheta=PorM*acos(Cosin);
if (lighttheta<0) {
     lighttheta=PI-lighttheta;
}
lighttheta=noise(lightscale*lighttheta);
lighttheta=clamp(lighttheta,light_min,light_max);
lighttheta=(lighttheta-light_min)/(light_max-light_min);
lighttheta=lighttheta*pow(lightdis+1,-(light_roughness));
lighttheta*=light_intensity;

Ci=Oi*mix(Ci,lightcolor,lighttheta);
/*  highlight   END  */

}