diff options
Diffstat (limited to 'Voigt.f')
| -rwxr-xr-x | Voigt.f | 14 |
1 files changed, 7 insertions, 7 deletions
@@ -9,9 +9,9 @@ c****************************************************************************** implicit real*8 (a-h,o-z) a2 = a*a v2 = v*v - if(a .eq. 0.) go to 62 - if(a .le. 0.2) go to 50 - if(a.le.1.4 .and. a+v.le.3.2) go to 30 + if(a == 0.) go to 62 + if(a <= 0.2) go to 50 + if(a<=1.4 .and. a+v<=3.2) go to 30 c case 1, a .gt. 1.4 or (a .gt. 0.2 and a+v .gt. 3.2) @@ -22,8 +22,8 @@ c case 1, a .gt. 1.4 or (a .gt. 0.2 and a+v .gt. 3.2) return 30 u=0.979895023-0.962846325*a+0.532770573*a2-0.122727278*a*a2 31 h0 = dexp(-v2) - if(v .ge. 2.4) go to 35 - if(v .ge. 1.3) go to 33 + if(v >= 2.4) go to 35 + if(v >= 1.3) go to 33 c case 2, 0.2 .lt. a .le. 1.4 and a+v .le. 3.2 @@ -36,7 +36,7 @@ c case 2, 0.2 .lt. a .le. 1.4 and a+v .le. 3.2 35 h1=(0.554153432+0.278711796*v-0.188325687*v2+0.042991293*v*v2 1 -0.003278278*v2*v2)/(v2 - 1.5) 36 h2 = (1.0 - 2.0*v2)*h0 - if(a .le. 0.2) go to 52 + if(a <= 0.2) go to 52 h1 = h1 + 1.1283790*h0 h2p = h2 h2 = h2 - h0 + 1.1283790*h1 @@ -48,7 +48,7 @@ c case 2, 0.2 .lt. a .le. 1.4 and a+v .le. 3.2 c case 3, a .le. 0.2 and v .lt. 5.0 -50 if(v .ge. 5.0) go to 60 +50 if(v >= 5.0) go to 60 go to 31 52 voigt = h0 + h1*a + h2*a2 voigt = voigt/1.772454 |