41 static double zLimit = 35.5;
46 if(
nMatch(
"TRAC",chCard) )
65 fprintf(
ioQQQ,
" Abundances relative to H, Z\n" );
68 fprintf(
ioQQQ,
" Z ");
75 fprintf(
ioQQQ,
" \n" );
79 else if( lgEOL && !lgDebug )
82 fprintf(
ioQQQ,
" The metallicity must appear on this line.\n" );
92 if(
nMatch(
" LOG",chCard) )
98 else if(
nMatch(
"LINE",chCard) )
102 fprintf(
ioQQQ,
" Z .le.0 not allowed, Z=%10.2e\n",
122 while( zed <= zHigh )
124 if( zed < 1e-3 || zed > zLimit )
126 fprintf(
ioQQQ,
" The metallicity must be between 1E-3 and%10.2e\n",
136 zh = 1.081646723 - 0.04600492*zed + 8.6569e-6*zed2 + 1.922e-5*
137 zed3 - 2.0087e-7*zed4;
140 zhe = 0.864675891 + 0.044423807*zed + 7.10886e-5*zed2 - 5.3242e-5*
141 zed3 + 5.70194e-7*zed4;
150 zc = 0.347489799 + 0.016054107*zed - 0.00185855*zed2 + 5.43015e-5*
151 zed3 - 5.3123e-7*zed4;
155 zn = -0.06549567 + 0.332073984*zed + 0.009146066*zed2 - 0.00054441*
156 zed3 + 6.16363e-6*zed4;
163 zn = -0.00044731103 + 0.00026453554*zed + 0.52354843*zed2 -
164 0.41156655*zed3 + 0.1290967*zed4;
167 zn = 0.000344828*zed;
173 zo = 1.450212747 - 0.05379041*zed + 0.000498919*zed2 + 1.09646e-5*
174 zed3 - 1.8147e-7*zed4;
178 zne = 1.110139023 + 0.002551998*zed - 2.09516e-7*zed3 - 0.00798157*
186 zna = 1.072721387 - 0.02391599*
POW2(zedlog) + .068644972*
187 zedlog + 0.017030935/sqrzed;
189 zna =
MAX2(1e-12,zna);
193 zmg = 1.147209646 - 7.9491e-7*
POW3(zed) - .00264458*
POW2(zedlog) -
198 zal = 1.068116358 - 0.00520227*sqrzed*zedlog - 0.01403851*
199 POW2(zedlog) + 0.066186787*zedlog;
201 zal =
MAX2(1e-12,zal);
205 zsi = 1.067372578 + 0.011818743*zed - 0.00107725*zed2 + 3.66056e-5*
206 zed3 - 3.556e-7*zed4;
221 zar = 1.091067724 + 2.51124e-6*zed3 - 0.0039589*sqrzed*zedlog +
229 zca = 1.077553875 - 0.00888806*zed + 0.001479866*zed2 - 6.5689e-5*
230 zed3 + 1.16935e-6*zed4;
234 zfe = 0.223713045 + 0.001400746*zed + 0.000624734*zed2 - 3.5629e-5*
235 zed3 + 8.13184e-7*zed4;
272 for( i=2; i <
LIMELM; i++ )
280 fprintf(
ioQQQ,
"%10.2e", zed );
281 for( i=0; i <
LIMELM; i++ )
285 fprintf(
ioQQQ,
"\n" );
296 zed =
MIN2(zed*1.5,zLimit);