From 81b03f9df3a290fa0a7a81456a1d7f303fa11380 Mon Sep 17 00:00:00 2001 From: Lorenzo Figini Date: Thu, 21 May 2015 10:49:44 +0000 Subject: [PATCH] solved bug in freq value assignment for ibeam>0 --- src/gray-externals.f | 8 +++----- 1 file changed, 3 insertions(+), 5 deletions(-) diff --git a/src/gray-externals.f b/src/gray-externals.f index c0699bb..b38cdbb 100644 --- a/src/gray-externals.f +++ b/src/gray-externals.f @@ -522,7 +522,7 @@ c common/scal/iscal common/fghz/fghz c - open(602,file='gray.data',status= 'unknown') + open(602,file='gray.data',status= 'old') c c alpha0, beta0 (cartesian) launching angles c fghz wave frequency (GHz) @@ -677,8 +677,6 @@ c write(*,*) ' nrayr < 5 ! => OPTICAL CASE ONLY' write(*,*) end if -c - fhz=fghz*1.0d9 c c read data for beam from file if ibeam>0 c @@ -709,6 +707,7 @@ c end if end if c + fhz=fghz*1.0d9 akinv=1.0d0/ak0 c bresg=2.0d0*pi*fhz*me*vc/qe @@ -717,7 +716,7 @@ c c xg=xgcn*dens19 c xgcn=1.0d-5*qe**2/(pi*me*fghz**2) -c +c sox=-1.0d0 if(iox.eq.2) sox=1.0d0 c @@ -996,7 +995,6 @@ c c============================================================================ c c beam of interest read(nfbeam,*) beamname, iox, fghz, nalpha, nbeta - fhz=fghz*1.0d9 c c c==================================================================================== c unused beam data grids