26 static long int nZoneEval=-1;
27 if(
nzone==nZoneEval )
34 for(
long nelem=ipISO; nelem<
LIMELM; ++nelem )
36 if( nelem >= 2 && !
dense.lgElmtOn[nelem] )
47 sp->
ex[ipHi][ipLo].pestrk = 0.;
48 sp->
ex[ipHi][ipLo].pestrk_up = 0.;
63 ah *= pow( (
double)(nelem+1), -4.5 );
67 stark = 0.264*pow(ah,0.4);
73 long ipHi =
iso_ctrl.nLyaLevel[ipISO];
80 sp->
ex[ipHi][0].pestrk_up = strkla/2.*
MAX2(1.,aa);
86 sp->
ex[ipHi][0].pestrk_up =
87 MIN2(.01,sp->
ex[ipHi][0].pestrk_up);
88 sp->
ex[ipHi][0].pestrk_up = 0.;
89 sp->
ex[ipHi][0].pestrk = sp->
ex[ipHi][0].pestrk_up *
98 sp->
ex[ipHi][0].pestrk_up = stark / 2. *
102 sp->
ex[ipHi][0].pestrk_up =
MIN2(.01,sp->
ex[ipHi][0].pestrk_up);
104 sp->
ex[ipHi][0].pestrk_up;
110 sp->
ex[ipHi][0].pestrk_up = 0.;
111 sp->
ex[ipHi][0].pestrk = 0.;
123 strkar( sp->
st[ipLo].n(), sp->
st[ipHi].n() ) *
125 sp->
ex[ipHi][ipLo].pestrk_up =
129 sp->
ex[ipHi][ipLo].pestrk_up;
138 sp->
ex[ipHi][ipLo].pestrk_up = 0.;
139 sp->
ex[ipHi][ipLo].pestrk = 0.;