40 static double zLimit = 35.5;
52 fprintf(
ioQQQ,
" Abundances relative to H, Z\n" );
55 fprintf(
ioQQQ,
" Z ");
62 fprintf(
ioQQQ,
" \n" );
75 fprintf(
ioQQQ,
" Z .le.0 not allowed, Z=%10.2e\n",
90 if( zed < 1e-3 || zed > zLimit )
92 fprintf(
ioQQQ,
" The metallicity must be between 1E-3 and%10.2e\n",
102 zh = 1.081646723 - 0.04600492*zed + 8.6569e-6*zed2 + 1.922e-5*
103 zed3 - 2.0087e-7*zed4;
106 zhe = 0.864675891 + 0.044423807*zed + 7.10886e-5*zed2 - 5.3242e-5*
107 zed3 + 5.70194e-7*zed4;
116 zc = 0.347489799 + 0.016054107*zed - 0.00185855*zed2 + 5.43015e-5*
117 zed3 - 5.3123e-7*zed4;
121 zn = -0.06549567 + 0.332073984*zed + 0.009146066*zed2 - 0.00054441*
122 zed3 + 6.16363e-6*zed4;
129 zn = -0.00044731103 + 0.00026453554*zed + 0.52354843*zed2 -
130 0.41156655*zed3 + 0.1290967*zed4;
133 zn = 0.000344828*zed;
139 zo = 1.450212747 - 0.05379041*zed + 0.000498919*zed2 + 1.09646e-5*
140 zed3 - 1.8147e-7*zed4;
144 zne = 1.110139023 + 0.002551998*zed - 2.09516e-7*zed3 - 0.00798157*
152 zna = 1.072721387 - 0.02391599*
POW2(zedlog) + .068644972*
153 zedlog + 0.017030935/sqrzed;
155 zna =
MAX2(1e-12,zna);
159 zmg = 1.147209646 - 7.9491e-7*
POW3(zed) - .00264458*
POW2(zedlog) -
164 zal = 1.068116358 - 0.00520227*sqrzed*zedlog - 0.01403851*
165 POW2(zedlog) + 0.066186787*zedlog;
167 zal =
MAX2(1e-12,zal);
171 zsi = 1.067372578 + 0.011818743*zed - 0.00107725*zed2 + 3.66056e-5*
172 zed3 - 3.556e-7*zed4;
187 zar = 1.091067724 + 2.51124e-6*zed3 - 0.0039589*sqrzed*zedlog +
195 zca = 1.077553875 - 0.00888806*zed + 0.001479866*zed2 - 6.5689e-5*
196 zed3 + 1.16935e-6*zed4;
200 zfe = 0.223713045 + 0.001400746*zed + 0.000624734*zed2 - 3.5629e-5*
201 zed3 + 8.13184e-7*zed4;
238 for(
long i=2; i <
LIMELM; i++ )
246 fprintf(
ioQQQ,
"%10.2e", zed );
247 for(
long i=0; i <
LIMELM; i++ )
249 fprintf(
ioQQQ,
"%6.2f", log10(
abund.solar[i]) );
251 fprintf(
ioQQQ,
"\n" );
262 zed =
MIN2(zed*1.5,zLimit);