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/mkl/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/, /usr/include/,./,/opt/cray/pe/mpich/8.1.9/ofi/intel/19.0/include/,/opt/cray/pe/libsci/21.08.1.2/INTEL/19.0/x86_64/include/, /pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/, /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/, /pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/, /usr/lib64/gcc/x86_64-suse-linux/7/include/,/usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/,/usr/include/, /usr/include/,/usr/include/ -list filename : interp.lst no -o COMPILER: Intel(R) Fortran 19.1-1655