diff options
Diffstat (limited to 'Inlines.f')
| -rwxr-xr-x | Inlines.f | 78 |
1 files changed, 39 insertions, 39 deletions
@@ -17,38 +17,38 @@ c****************************************************************************** integer n2 - if (num .eq. 2) go to 4 - if (num .eq. 6) go to 340 + if (num == 2) go to 4 + if (num == 6) go to 340 n1marker = 1 n2 = 0 c*****decide if certain element abundances need to be modified. - if (numpecatom .gt. 0) then + if (numpecatom > 0) then do iatom=3,95 xabund(iatom) = 10.**pecabund(iatom,isynth)* . 10.**abfactor(isynth)*xabu(iatom) enddo endif - if (num .ne. 5) then + if (num /= 5) then write (nf1out,1004) xmetals = abscale + abfactor(isynth) - if (ninetynineflag .eq. 1) then + if (ninetynineflag == 1) then write (nf1out,1005) xmetals - if (nf2out .gt. 0) write (nf2out,1005) xmetals + if (nf2out > 0) write (nf2out,1005) xmetals else - if (nf2out .gt. 0) write (nf2out,1006) abscale + if (nf2out > 0) write (nf2out,1006) abscale endif do j=1,93 - if (pec(j) .gt. 0 ) then + if (pec(j) > 0 ) then dummy1(j) = dlog10(xabund(j)) + 12.0 - if (dummy1(j) .le. -10.) then + if (dummy1(j) <= -10.) then write (nf1out,1008) names(j),dummy1(j) - if (nf2out .gt. 0) + if (nf2out > 0) . write (nf2out,1008) names(j),dummy1(j) else write (nf1out,1007) names(j),dummy1(j) - if (nf2out .gt. 0) + if (nf2out > 0) . write (nf2out,1007) names(j),dummy1(j) endif endif @@ -57,12 +57,12 @@ c*****decide if certain element abundances need to be modified. c*****output information about the isotopic ratios - if (numiso .gt. 0) then + if (numiso > 0) then write (nf1out,1014) do i=1,numiso iiso = isotope(i) write (nf1out,1015) iiso, isotope(i), isoabund(i,isorun) - if (nf2out .gt. 0) write (nf2out,1015) + if (nf2out > 0) write (nf2out,1015) . iiso, isotope(i), isoabund(i,isorun) enddo endif @@ -74,7 +74,7 @@ c if 'printstrong' gt 0 then the strong lines have c been printed printstrong = -1 - if (num .ne. 4) then + if (num /= 4) then rewind nflines wave = start read (nflines,1001) linitle @@ -83,10 +83,10 @@ c been printed c*****read in the strong lines if needed 302 nstrong = 0 - if (dostrong .gt. 0 ) then + if (dostrong > 0 ) then rewind nfslines do j=1,41 - if (linfileopt .eq. 0) then + if (linfileopt == 0) then read (nfslines,1002,end=340) swave1(j),satom1(j),se(j), . sgf(j),sdampnum(j),sd0(j),swidth(j) else @@ -97,12 +97,12 @@ c*****read in the strong lines if needed iatom = satom1(j) scharge(j) = 1.0 + dble(int(10.0*(satom1(j) - iatom) . +0.0001)) - if (scharge(j) .gt. 3.) then + if (scharge(j) > 3.) then write (*,1003) swave1(i), satom1(i) stop endif enddo - if (nstrong .gt. 40) then + if (nstrong > 40) then write(*,*) 'STRONG LINE LIST HAS MORE THAN 40 LINES. THIS' write(*,*) 'IS NOT ALLOWED. I QUIT!' stop @@ -111,7 +111,7 @@ c*****read in the strong lines if needed 340 nlines = 2500 - nstrong j = 1 -333 if (linfileopt .eq. 0) then +333 if (linfileopt == 0) then read (nflines,1002,end=311) wave1(j),atom1(j),e(j,1),gf(j), . dampnum(j),d0(j),width(j) else @@ -120,26 +120,26 @@ c*****read in the strong lines if needed endif iatom = atom1(j) charge(j) = 1.0 + dble(int(10.0*(atom1(j) - iatom)+0.0001)) - if (charge(j) .gt. 3.) then + if (charge(j) > 3.) then write (*,1003) wave1(j), atom1(j) stop endif - if (width(j) .lt. 0.) then - if (control .eq. 'blends ') then + if (width(j) < 0.) then + if (control == 'blends ') then write (*,*) 'BLENDS cannot have negative EWs! I QUIT!' stop else go to 333 endif endif - if (iunits .eq. 1) wave1(j) = 1.d+4*wave1(j) + if (iunits == 1) wave1(j) = 1.d+4*wave1(j) j = j + 1 - if (j .le. nlines) go to 333 + if (j <= nlines) go to 333 311 nlines = j - 1 c*****append the strong lines here if necessary - if (dostrong .gt. 0) then + if (dostrong > 0) then do k=1,nstrong wave1(nlines+k) = swave1(k) atom1(nlines+k) = satom1(k) @@ -155,7 +155,7 @@ c*****append the strong lines here if necessary c*****here groups of lines for blended features are defined do j=1,nlines+nstrong - if (wave1(j) .lt. 0.) then + if (wave1(j) < 0.) then group(j) = 1 wave1(j) = dabs(wave1(j)) width(j) = width(j-1) @@ -167,7 +167,7 @@ c*****here groups of lines for blended features are defined c*****here excitation potentials are changed from cm^-1 to eV, if needed do j=1,nlines+nstrong - if (e(j,1) .gt. 50.) then + if (e(j,1) > 50.) then do jj=1,nlines+nstrong e(jj,1) = 1.2389e-4*e(jj,1) enddo @@ -178,7 +178,7 @@ c*****here excitation potentials are changed from cm^-1 to eV, if needed c*****here log(gf) values are turned into gf values, if needed do j=1,nlines+nstrong - if (gfstyle.eq.0 .or. gf(j) .lt. 0) then + if (gfstyle==0 .or. gf(j) < 0) then do jj=1,nlines+nstrong gf(jj) = 10.**gf(jj) enddo @@ -190,7 +190,7 @@ c*****here log(gf) values are turned into gf values, if needed c*****turn log(RW) values and EW values in mA into EW values in A. Stuff c duplicate EW values of the first line of a blend into all blend members. do j=1,nlines+nstrong - if (width(j) .lt. 0.) then + if (width(j) < 0.) then width(j) = 10.**width(j)*wave1(j) else width(j) = width(j)/1000. @@ -208,13 +208,13 @@ c and a different one for atomic lines c*****here are the calculations specific to molecular lines - if (iatom .ge. 100) then + if (iatom >= 100) then call sunder (atom1(j),ia,ib) - if (ia .gt. ib) then + if (ia > ib) then write (*,1010) ia,ib stop endif - if (atom10-int(atom10) .le. 0.0) then + if (atom10-int(atom10) <= 0.0) then amass(j) = xam(ia) + xam(ib) mas1 = xam(ia) + 0.0000001 mas2 = xam(ib) + 0.0000001 @@ -223,8 +223,8 @@ c*****here are the calculations specific to molecular lines mas1 = jat100 - 100*int(atom10) jat10000 = int(10000.*(atom10+0.00001)) mas2 = jat10000 - 100*jat100 - if (mas1.gt.mas2 .or. mas1.le.0.0 .or. - . mas2.le.0.0) then + if (mas1>mas2 .or. mas1<=0.0 .or. + . mas2<=0.0) then write (*,1011) mas1, mas2 stop endif @@ -232,9 +232,9 @@ c*****here are the calculations specific to molecular lines endif c*****use an internal dissociation energy for molecules if the user c does not read one in - if (d0(j) .eq. 0.) then + if (d0(j) == 0.) then do k=1,110 - if (int(datmol(1,k)+0.01) .eq. + if (int(datmol(1,k)+0.01) == . int(atom1(j)+0.01)) then d0(j) = datmol(2,k) go to 390 @@ -251,7 +251,7 @@ c does not read one in c*****here are the calculations specific to atomic lines else - if (atom10-int(atom10) .le. 0.0) then + if (atom10-int(atom10) <= 0.0) then amass(j) = xam(iatom) else atom10 = atom10 + 0.00001 @@ -266,8 +266,8 @@ c*****here are the calculations specific to atomic lines c*****quit the routine normally - if (nlines+nstrong .lt. 2500) then - if (sstop .gt. wave1(nlines)+10.) sstop = wave1(nlines)+10. + if (nlines+nstrong < 2500) then + if (sstop > wave1(nlines)+10.) sstop = wave1(nlines)+10. endif lim1line = 1 return |