36 #define DEBUGSTATE false
56 char *chToken,*chAtNo,*chIonStg;
61 static int nCalled = 0;
78 fprintf(
ioQQQ,
" Nemala_Start opening species.ini:");
80 atmolDATA =
open_data(
"species.ini",
"r" );
82 if(
read_whole_line( chLine , (
int)
sizeof(chLine) , atmolDATA ) == NULL )
84 fprintf(
ioQQQ,
" Nemala_Start could not read first line of species.ini.\n");
91 while(
read_whole_line( chLine , (
int)
sizeof(chLine) , atmolDATA ) != NULL )
95 if( (chLine[0] !=
'#')&&(chLine[0] !=
'\n')&&(chLine[0] !=
' '))
99 printf(
"The number of species is %li \n",
nSpecies);
102 if( fseek( atmolDATA , 0 , SEEK_SET ) != 0 )
104 fprintf(
ioQQQ,
" Nemala_Start could not rewind species.ini.\n");
142 while(
read_whole_line( chLine , (
int)
sizeof(chLine) , atmolDATA ) != NULL )
145 if ((chLine[0]!=
'#') && (chLine[0]!=
'\n')&&(chLine[0]!=
'\t')&&(chLine[0]!=
'\r'))
150 strcpy(chDLine, chLine);
151 chToken = strtok(chDLine,
" ");
152 los = (int)strlen(chToken);
154 printf(
"The length of the string is %d\n",los);
156 strcpy(
Species[i].chptrSpName,chToken);
160 printf(
"The name of the species is%s \n",
Species[i].chptrSpName);
161 if(
nMatch(
"LEID",chLine) )
165 else if(
nMatch(
"CHIA",chLine))
168 chAtNo = strtok(chToken,
"_");
170 chIonStg = strtok(NULL,
"_");
191 for( intNoSp=0; intNoSp<
nSpecies; intNoSp++ )
249 printf(
"The species is %s \n",
Species[i].chptrSpName);
250 printf(
"The data output is in the following format \n");
251 printf(
"Label Energy St.wt Pop Lifetime\n");
255 printf(
"This is the %ld state \n",j);
279 for( ipLo=0; ipLo<ipHi; ipLo++ )
281 if(
atmolTrans[intSpIndex][ipHi][ipLo].Emis!=NULL)
301 for( intNS = 0; intNS <
nSpecies; intNS++ )
322 else if(strcmp(p,
"HE")== 0)
326 else if(strcmp(p,
"C")== 0)
330 else if(strcmp(p,
"N")== 0)
334 else if(strcmp(p,
"O")== 0)
338 else if(strcmp(p,
"NE")== 0)
342 else if(strcmp(p,
"NA")== 0)
346 else if(strcmp(p,
"MG")== 0)
350 else if(strcmp(p,
"AL")== 0)
354 else if(strcmp(p,
"SI")== 0)
358 else if(strcmp(p,
"P")== 0)
362 else if(strcmp(p,
"S")== 0)
366 else if(strcmp(p,
"CL")== 0)
370 else if(strcmp(p,
"AR")== 0)
374 else if(strcmp(p,
"K")== 0)
378 else if(strcmp(p,
"CA")== 0)
382 else if(strcmp(p,
"SC")== 0)
386 else if(strcmp(p,
"TI")== 0)
390 else if(strcmp(p,
"V")== 0)
394 else if(strcmp(p,
"CR")== 0)
398 else if(strcmp(p,
"MN")== 0)
402 else if(strcmp(p,
"FE")== 0)
406 else if(strcmp(p,
"CO")== 0)
410 else if(strcmp(p,
"NI")== 0)
414 else if(strcmp(p,
"CU")== 0)
418 else if(strcmp(p,
"ZN")== 0)