Properties {_MainTex ("Base (RGB) Gloss (A)", 2D) = "white" {}//刮风的方向（世界坐标系下）_Wind("Wind params",Vector) = (1,1,1,1)//风的频率_WindEdgeFlutter("Wind edge fultter factor", float) = 0.5//风的频率的缩放_WindEdgeFlutterFreqScale("Wind edge fultter freq scale",float) = 0.5
}

//计算风的一些参数
float4  windParams  = float4(0,_WindEdgeFlutter,bendingFact.xx);
float2  windTime = _Time.y * float2(_WindEdgeFlutterFreqScale,1);
float4  mdlPos  = AnimateVertex2(v.vertex,v.normal,windParams,wind,windTime);
//mvp矩阵变换
o.pos   = mul(UNITY_MATRIX_MVP,mdlPos);

inline float4 AnimateVertex2(float4 pos, float3 normal, float4 animParams,float4 wind,float2 time)
{   // animParams stored in color// animParams.x = branch phase// animParams.y = edge flutter factor// animParams.z = primary factor// animParams.w = secondary factorfloat fDetailAmp = 0.1f;float fBranchAmp = 0.3f;// Phases (object, vertex, branch)float fObjPhase = dot(_Object2World[3].xyz, 1);float fBranchPhase = fObjPhase + animParams.x;float fVtxPhase = dot(pos.xyz, animParams.y + fBranchPhase);// x is used for edges; y is used for branchesfloat2 vWavesIn = time.yy  + float2(fVtxPhase, fBranchPhase );// 1.975, 0.793, 0.375, 0.193 are good frequenciesfloat4 vWaves = (frac( vWavesIn.xxyy * float4(1.975, 0.793, 0.375, 0.193) ) * 2.0 - 1.0);vWaves = SmoothTriangleWave( vWaves );float2 vWavesSum = vWaves.xz + vWaves.yw;// Edge (xz) and branch bending (y)float3 bend = animParams.y * fDetailAmp * normal.xyz;bend.y = animParams.w * fBranchAmp;pos.xyz += ((vWavesSum.xyx * bend) + (wind.xyz * vWavesSum.y * animParams.w)) * wind.w; // Primary bending// Displace positionpos.xyz += animParams.z * wind.xyz;return pos;
}

o.uv = TRANSFORM_TEX(v.texcoord.xy,_MainTex) + frac(float2(_ScrollX, _ScrollY) * _Time);
o.uv2 = TRANSFORM_TEX(v.texcoord.xy,_DetailTex) + frac(float2(_Scroll2X, _Scroll2Y) * _Time);

fixed4 frag (v2f i) : COLOR
{fixed4 o;fixed4 tex = tex2D (_MainTex, i.uv);fixed4 tex2 = tex2D (_DetailTex, i.uv2);o = (tex * tex2) * i.color;return o;
}

float3   viewPos = mul(UNITY_MATRIX_MV,v.vertex);
float       dist        = length(viewPos);
float       nfadeout    = saturate(dist / _FadeOutDistNear);
float       ffadeout    = 1 - saturate(max(dist - _FadeOutDistFar,0) * 0.2);

float4 vpos = v.vertex;
vpos.xyz -=   v.normal * saturate(1 - nfadeout) * v.color.a * _ContractionAmount;

float    fracTime    = fmod(time,_TimeOnDuration + _TimeOffDuration);
float   wave        = smoothstep(0,_TimeOnDuration * 0.25,fracTime)  * (1 - smoothstep(_TimeOnDuration * 0.75,_TimeOnDuration,fracTime));
float   noiseTime   = time *  (6.2831853f / _TimeOnDuration);
float   noise       = sin(noiseTime) * (0.5f * cos(noiseTime * 0.6366f + 56.7272f) + 0.5f);
float   noiseWave   = _NoiseAmount * noise + (1 - _NoiseAmount);wave = _NoiseAmount < 0.01f ? wave : noiseWave;
o.color = nfadeout * _Color * _Multiplier * wave;

float3 centerOffs      = float3(float(0.5).xx - v.color.rg,0) * v.texcoord1.xyy;
float3  centerLocal = v.vertex.xyz + centerOffs.xyz;
float3  viewerLocal = mul(_World2Object,float4(_WorldSpaceCameraPos,1));
float3  localDir            = viewerLocal - centerLocal;localDir[1] = lerp(0,localDir[1],_VerticalBillboarding);float     localDirLength=length(localDir);
float3  rightLocal;
float3  upLocal;CalcOrthonormalBasis(localDir / localDirLength,rightLocal,upLocal);float        distScale       = CalcDistScale(localDirLength) * v.color.a;
float3  BBNormal        = rightLocal * v.normal.x + upLocal * v.normal.y;
float3  BBLocalPos  = centerLocal - (rightLocal * centerOffs.x + upLocal * centerOffs.y) + BBNormal * distScale;BBLocalPos += _ViewerOffset * localDir;

#if 1// quite suprisinly this looks better (probably there is some error in AO calculation)ao = 1 - saturate(SphereAO(_Sphere0,wrldPos,wrldNormal) + SphereAO(_Sphere1,wrldPos,wrldNormal) + SphereAO(_Sphere2,wrldPos,wrldNormal));
#elseao = 1 - max(max(SphereAO(_Sphere0,wrldPos,wrldNormal),SphereAO(_Sphere1,wrldPos,wrldNormal)),SphereAO(_Sphere2,wrldPos,wrldNormal));
#endifao = max(ao,1 - _Intensity) + (1 - v.color.r);o.color = fixed4(ao,ao,ao,ao);#endif

float SphereAO(float4 sphere,float3 pos,float3 normal)
{float3 dir = sphere.xyz - pos;float   d   = length(dir);float    v;dir /= d;v = (sphere.w / d);return dot(normal,dir) * v * v;
}