istp/gray
2
0
Fork 0
Code Issues Pull Requests Packages Projects Releases 1 Wiki Activity

solved bug in freq value assignment for ibeam>0

Browse Source
This commit is contained in:
Lorenzo Figini 2015-05-21 10:49:44 +00:00
parent 72b3c5f511
commit 81b03f9df3
1 changed files with 3 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
src/gray-externals.f
View File

@ -522,7 +522,7 @@ c
common/scal/iscal common/scal/iscal
common/fghz/fghz common/fghz/fghz
c c
open(602,file='gray.data',status= 'unknown') open(602,file='gray.data',status= 'old')
c c
c alpha0, beta0 (cartesian) launching angles c alpha0, beta0 (cartesian) launching angles
c fghz wave frequency (GHz) c fghz wave frequency (GHz)
@ -677,8 +677,6 @@ c
write(*,*) ' nrayr < 5 ! => OPTICAL CASE ONLY' write(*,*) ' nrayr < 5 ! => OPTICAL CASE ONLY'
write(*,*) write(*,*)
end if end if
c
fhz=fghz*1.0d9
c c
c read data for beam from file if ibeam>0 c read data for beam from file if ibeam>0
c c
@ -709,6 +707,7 @@ c
end if end if
end if end if
c c
fhz=fghz*1.0d9
akinv=1.0d0/ak0 akinv=1.0d0/ak0
c c
bresg=2.0d0*pi*fhz*me*vc/qe bresg=2.0d0*pi*fhz*me*vc/qe
@ -717,7 +716,7 @@ c
c xg=xgcn*dens19 c xg=xgcn*dens19
c c
xgcn=1.0d-5*qe**2/(pi*me*fghz**2) xgcn=1.0d-5*qe**2/(pi*me*fghz**2)
c c
sox=-1.0d0 sox=-1.0d0
if(iox.eq.2) sox=1.0d0 if(iox.eq.2) sox=1.0d0
c c
@ -996,7 +995,6 @@ c c============================================================================
c c
c beam of interest c beam of interest
read(nfbeam,*) beamname, iox, fghz, nalpha, nbeta read(nfbeam,*) beamname, iox, fghz, nalpha, nbeta
fhz=fghz*1.0d9
c c
c c==================================================================================== c c====================================================================================
c unused beam data grids c unused beam data grids