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

trunk/beams: fix to alpha0 value if outside table range in read_beam1

Browse Source
This commit is contained in:
Daniele Micheletti 2018-09-26 12:58:54 +00:00
parent 9fcc54804f
commit 2c46d3327e
2 changed files with 4 additions and 3 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
Makefile.sum
View File

@ -17,7 +17,7 @@ FC=gfortran
FFLAGS=-O3 FFLAGS=-O3
#FFLAGS=-Wall -g -finit-real=nan -ffpe-trap=invalid -fcheck=all -fbounds-check #FFLAGS=-Wall -g -finit-real=nan -ffpe-trap=invalid -fcheck=all -fbounds-check
DIRECTIVES = -DREVISION="'$(shell svnversion src)'" DIRECTIVES = -DREVISION="178M"
all: $(EXE) all: $(EXE)
@ -65,7 +65,7 @@ utils.o: const_and_precisions.o
.PHONY: clean install .PHONY: clean install
# Remove output files # Remove output files
clean: clean:
rm -rf *.o *.mod $(EXE) rm -rf *.o *.mod
install: install:
@if [ -f $(EXE) ]; then \ @if [ -f $(EXE) ]; then \

3
src/beams.f90
View File

@ -57,7 +57,7 @@ contains
implicit none implicit none
! arguments ! arguments
character(len=*), intent(in) :: file_beam character(len=*), intent(in) :: file_beam
real(wp_), intent(in) :: alpha0 real(wp_), intent(inout) :: alpha0
real(wp_), intent(out) :: fghz,x00,y00,z00,beta0 real(wp_), intent(out) :: fghz,x00,y00,z00,beta0
real(wp_), intent(out) :: wcsi,weta,rcicsi,rcieta,phir,phiw real(wp_), intent(out) :: wcsi,weta,rcicsi,rcieta,phir,phiw
integer, intent(in), optional :: unit integer, intent(in), optional :: unit
@ -131,6 +131,7 @@ contains
! alpha0 outside table range ! alpha0 outside table range
if(alpha0 >= alphastv(nisteer)) ii=nisteer if(alpha0 >= alphastv(nisteer)) ii=nisteer
if(alpha0 <= alphastv(1)) ii=1 if(alpha0 <= alphastv(1)) ii=1
alpha0=alphastv(ii)
beta0=betastv(ii) beta0=betastv(ii)
x00=x00v(ii) x00=x00v(ii)
y00=y00v(ii) y00=y00v(ii)