Page 1 Source Listing INTERP
2022-05-02 15:46 interp.f
1 subroutine interp(datas,dataw,sistr,sigtop,jstr,idim,jdim)
2
3 dimension sistr(jstr),sigtop(jstr)
4 common/quadcb/ y0,y1,y2,x0,x1,x2
5
6 dimension datas(idim,jdim),sinm(jdim)
7 dimension dataw(idim,jstr)
8
9 c dimension datai(idim,jstr,3),sinmi(jstr,3)
10 dimension datai(300,150,3),sinmi(150,3)
11 common/itstcb/ datai,sinmi
12
13 c open(unit=60,file='intout.d',form='formatted')
14
15 pi=acos(-1.0)
16 dpsieq=1./36.
17 dsieqd=dpsieq*(180./pi)
18 dsinm=pi/72.
19 hdsi=dsinm/2.0
20 dsinmd=dsinm*(180./pi)
21 c write(60,*) 'idim,jdim,pi,jstr,dsieqd,dsinmd'
22 c write(60,*) idim,jdim,pi,jstr,dsieqd,dsinmd
23
24 do j=1,jdim
25 sinm(j)=(pi/2.)-(j-1)*dsinm
26 end do
27
28 do j=1,jstr
29 jeq=37-j
30 sistr(j)=asin(jeq*dpsieq)
31 sigtop(j)=(pi/2.)-sistr(j)
32 sid=sistr(j)*180./pi
33 sigd=sigtop(j)*180./pi
34 c write(60,*) j,jeq,sistr(j),sid,sigtop(j),sigd
35 end do
36
37 do i=1,idim
38 dataw(i,1)=datas(i,1)/jdim
39 dataw(i,jstr)=datas(i,jstr)/jdim
40 end do
41 do js=2,jstr-1
42 do j=2,jdim-1
43
44 if((sinm(j)-hdsi).lt.sistr(js))then
45 x0=sinm(j-1)
46 x1=sinm(j)
47 x2=sinm(j+1)
48 sinmi(js,1)=x0
49 sinmi(js,2)=x1
50 sinmi(js,3)=x2
51 do i=1,idim
52 y0=datas(i,j-1)
53 y1=datas(i,j)
54 y2=datas(i,j+1)
55 datai(i,js,1)=y0
56 datai(i,js,2)=y1
57 datai(i,js,3)=y2
Page 2 Source Listing INTERP
2022-05-02 15:46 interp.f
58 dataw(i,js)=fchb(sistr(js))
59 end do
60 go to 10
61 end if
62
63 end do
64 c end do loop over j
65 10 continue
66 end do
67 c end do over js
68 return
69 end
ENTRY POINTS
Name
interp_
SYMBOL CROSS REFERENCE
Name Object Declared Type Bytes Dimen Elements Attributes References
10 Label 65 60
ACOS Func 15 scalar 15
ASIN Func 30 scalar 30
DATAS Dummy 1 R(4) 4 2 0 ARG,INOUT 38,39,52,53,54
DATAW Dummy 1 R(4) 4 2 0 ARG,INOUT 38,39,58
DPSIEQ Local 16 R(4) 4 scalar 16,17,30
DSIEQD Local 17 R(4) 4 scalar 17
DSINM Local 18 R(4) 4 scalar 18,19,20,25
DSINMD Local 20 R(4) 4 scalar 20
FCHB Func 58 R(4) 4 scalar 58
HDSI Local 19 R(4) 4 scalar 19,44
I Local 37 I(4) 4 scalar 37,38,39,51,52,53,54,55,56,57,58
IDIM Dummy 1 I(4) 4 scalar ARG,INOUT 6,7,37,51
INTERP Subr 1
ITSTCB Common 11 541800
J Local 24 I(4) 4 scalar 24,25,28,29,30,31,32,33,42,44,45,4
6,47,52,53,54
JDIM Dummy 1 I(4) 4 scalar ARG,INOUT 6,24,38,39,42
JEQ Local 29 I(4) 4 scalar 29,30
JS Local 41 I(4) 4 scalar 41,44,48,49,50,55,56,57,58
JSTR Dummy 1 I(4) 4 scalar ARG,INOUT 3,7,28,39,41
PI Local 15 R(4) 4 scalar 15,17,18,20,25,31,32,33
QUADCB Common 4 24
SID Local 32 R(4) 4 scalar 32
SIGD Local 33 R(4) 4 scalar 33
SIGTOP Dummy 1 R(4) 4 1 0 ARG,INOUT 31,33
SINM Local 6 R(4) 4 1 0 25,44,45,46,47
SISTR Dummy 1 R(4) 4 1 0 ARG,INOUT 30,31,32,44,58
Page 3 Source Listing INTERP
2022-05-02 15:46 Symbol Table interp.f
TYPE COMPONENTS/COMMON VARIABLES
Name Type Bytes Offset Dimen Elements Attributes References
DATAI R(4) 4 0 3 135000 COM 55,56,57
SINMI R(4) 4 540000 2 450 COM 48,49,50
X0 R(4) 4 12 scalar COM 45,48
X1 R(4) 4 16 scalar COM 46,49
X2 R(4) 4 20 scalar COM 47,50
Y0 R(4) 4 0 scalar COM 52,55
Y1 R(4) 4 4 scalar COM 53,56
Y2 R(4) 4 8 scalar COM 54,57
Page 4 Source Listing INTERP
2022-05-02 15:46 interp.f
70
71 function fchb(si)
72
73 common/quadcb/ y0,y1,y2,x0,x1,x2
74
75 si0=y0*(si-x1)*(si-x2)/((x0-x1)*(x0-x2))
76 si1=y1*(si-x2)*(si-x0)/((x1-x0)*(x1-x2))
77 si2=y2*(si-x0)*(si-x1)/((x2-x0)*(x2-x1))
78
79 fchb=si0+si1+si2
80
81 return
82 end
ENTRY POINTS
Name
fchb_
SYMBOL CROSS REFERENCE
Name Object Declared Type Bytes Dimen Elements Attributes References
FCHB Func 71 R(4) 4 scalar 79
FCHB@0 Local 71 R(4) 4 scalar
QUADCB Common 73 24
SI Dummy 71 R(4) 4 scalar ARG,INOUT 75,76,77
SI0 Local 75 R(4) 4 scalar 75,79
SI1 Local 76 R(4) 4 scalar 76,79
SI2 Local 77 R(4) 4 scalar 77,79
TYPE COMPONENTS/COMMON VARIABLES
Name Type Bytes Offset Dimen Elements Attributes References
X0 R(4) 4 12 scalar COM 75,76,77
X1 R(4) 4 16 scalar COM 75,76,77
X2 R(4) 4 20 scalar COM 75,76,77
Y0 R(4) 4 0 scalar COM 75
Y1 R(4) 4 4 scalar COM 76
Y2 R(4) 4 8 scalar COM 77
Page 5 Source Listing FCHB
2022-05-02 15:46 Subprograms/Common Blocks interp.f
SUBPROGRAMS/COMMON BLOCKS
Name Object Declared Type Bytes Dimen Elements Attributes References
FCHB Func 71 R(4) 4 scalar 79
INTERP Subr 1
ITSTCB Common 11 541800
QUADCB Common 4 24
QUADCB Common 73 24
COMPILER OPTIONS BEING USED
-align noall -align nonone
-align nocommons -align nodcommons
-align noqcommons -align nozcommons
-align records -align sequence
-align norec1byte -align norec2byte
-align norec4byte -align norec8byte
-align norec16byte -align norec32byte
-align norec64byte -align noarray8byte
-align noarray16byte -align noarray32byte
-align noarray64byte -align noarray128byte
-align noarray256byte -altparam
-assume accuracy_sensitive -assume nobscc
-assume nobuffered_io -assume nobuffered_stdout
-assume byterecl -assume nocontiguous_assumed_shape
-assume nocontiguous_pointer -assume nocc_omp
-assume nocstring -assume nodummy_aliases
-assume nofpe_summary -assume noieee_fpe_flags
-assume nominus0 -assume noold_boz
-assume old_complex_align -assume old_unit_star
-assume old_inquire_recl -assume old_ldout_format
-assume old_ldout_zero -assume noold_logical_assign
-assume noold_logical_ldio -assume old_maxminloc
-assume old_xor -assume noprotect_allocates
-assume protect_constants -assume noprotect_parens
-assume split_common -assume source_include
-assume nostd_intent_in -assume std_minus0_rounding
-assume nostd_mod_proc_name -assume std_value
-assume realloc_lhs -assume underscore
-assume no2underscores -assume norecursion
no -auto -auto_scalar
no -bintext -ccdefault default
-check noarg_temp_created -check noassume
-check nobounds -check nocontiguous
-check noformat -check nooutput_conversion
-check nooverflow -check nopointers
-check noshape -check nostack
-check nouninitialized -check noudio_iostat
-coarray-num-procs 0 no -coarray-config-file
-convert big_endian -cross_reference
-D __INTEL_COMPILER=1910 -D __INTEL_COMPILER_UPDATE=3
-D __unix__ -D __unix
-D __linux__ -D __linux
-D __gnu_linux__ -D unix
Page 6 Source Listing FCHB
2022-05-02 15:46 interp.f
-D linux -D __ELF__
-D __x86_64 -D __x86_64__
-D __amd64 -D __amd64__
-D __INTEL_COMPILER_BUILD_DATE=20200925 -D __INTEL_OFFLOAD
-D __MMX__ -D __SSE__
-D __SSE_MATH__ -D __SSE2__
-D __SSE2_MATH__ -D __SSE3__
-D __SSSE3__ -D __SSE4_1__
-D __SSE4_2__ -D __POPCNT__
-D __PCLMUL__ -D __AES__
-D __AVX__ -D __F16C__
-D __AVX_I__ -D __RDRND__
-D __FMA__ -D __FP_FAST_FMA
-D __FP_FAST_FMAF -D __BMI__
-D __LZCNT__ -D __AVX2__
-D __haswell -D __haswell__
-D __tune_haswell__ -D __core_avx2
-D __core_avx2__ -D __tune_core_avx2__
-double_size 64 no -d_lines
no -Qdyncom -error_limit 30
no -f66 no -f77rtl
no -fast -fpscomp nofilesfromcmd
-fpscomp nogeneral -fpscomp noioformat
-fpscomp noldio_spacing -fpscomp nologicals
-fixed no -fpconstant
-fpe3 -fprm nearest
-ftz -fp_model precise
-fp_model nofast -fp_model nostrict
-fp_model nosource -fp_model nodouble
-fp_model noextended -fp_model novery_fast
-fp_model noexcept -fp_model nono_except
-fp_modbits nofp_contract -fp_modbits nono_fp_contract
-fp_modbits nofenv_access -fp_modbits nono_fenv_access
-fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range
-fp_modbits noprec_div -fp_modbits no_prec_div
-fp_modbits noprec_sqrt -fp_modbits no_prec_sqrt
-fp_modbits ftz -fp_modbits nono_ftz
-fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range
-fp_modbits notrunc_compares -fp_modbits nono_trunc_compares
-fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares
-fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion
-fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion
-fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy
-fp_modbits noexception_semantics -fp_modbits nono_exception_semantics
-fp_modbits noprecise_libm_functions -fp_modbits no_precise_libm_functions
-heap_arrays 0 no -threadprivate_compat
-g0 -iface nomixed_str_len_arg
-iface nono_mixed_str_len_arg -init noarrays
-init nohuge -init noinfinity
-init nominus_huge -init nominus_infinity
-init nominus_tiny -init nonan
-init nosnan -init notiny
-init nozero no -intconstant
-integer_size 32 no -mixed_str_len_arg
no -module -names lowercase
no -noinclude no -o
-offload-build=host -openmp-simd
Page 7 Source Listing FCHB
2022-05-02 15:46 interp.f
-O3 no -pad_source
-real_size 32 no -recursive
-reentrancy threaded -vec=simd
-show nofullpath -show noinclude
-show map -show options
no -syntax_only no -threadcom
no -U no -vms
-w noall -w nonone
-w alignments -w nodeclarations
-w noexternals -w general
-w noignore_bounds -w noignore_loc
-w nointerfaces -w noshape
-w notruncated_source -w uncalled
-w uninitialized -w nounused
-w usage no -wrap-margins
-includepath : /pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/,
/usr/include/,.FOR,./.FOR,/opt/cray/pe/mpich/8.1.9/ofi/intel/19.0/include/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/intel64/lp64/.FOR,
/pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/.FOR,
/pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/.FOR,
/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/.FOR,
/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/.FOR,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/.FOR,
/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/.FOR,/usr/lib64/gcc/x86_64-suse-linux/7/include/.FOR,
/usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/.FOR,/usr/include/.FOR,/usr/include/.FOR,/usr/include/.FOR
-list filename : interp.lst
no -o
COMPILER: Intel(R) Fortran 19.1-1655