Page 1 Source Listing INTERP 2013-05-14 19:54 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 2013-05-14 19:54 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 DATAI Scalar 10 R(4) 4 3 135000 COM 55,56,57 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 SINMI Scalar 10 R(4) 4 2 450 COM 48,49,50 SISTR Dummy 1 R(4) 4 1 0 ARG,INOUT 30,31,32,44,58 X0 Scalar 4 R(4) 4 scalar COM 45,48 X1 Scalar 4 R(4) 4 scalar COM 46,49 Page 3 Source Listing INTERP 2013-05-14 19:54 Symbol Table interp.f Name Object Declared Type Bytes Dimen Elements Attributes References X2 Scalar 4 R(4) 4 scalar COM 47,50 Y0 Scalar 4 R(4) 4 scalar COM 52,55 Y1 Scalar 4 R(4) 4 scalar COM 53,56 Y2 Scalar 4 R(4) 4 scalar COM 54,57 Page 4 Source Listing INTERP 2013-05-14 19:54 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 X0 Scalar 73 R(4) 4 scalar COM 75,76,77 X1 Scalar 73 R(4) 4 scalar COM 75,76,77 X2 Scalar 73 R(4) 4 scalar COM 75,76,77 Y0 Scalar 73 R(4) 4 scalar COM 75 Y1 Scalar 73 R(4) 4 scalar COM 76 Y2 Scalar 73 R(4) 4 scalar COM 77 Page 5 Source Listing FCHB 2013-05-14 19:54 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 nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores no -auto -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd Page 6 Source Listing FCHB 2013-05-14 19:54 interp.f -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest no -ftz -fp_model noprecise -fp_model fast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -heap_arrays 0 no -threadprivate_compat -g0 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /gpfs/gp2/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/,.f,./.f,/usrx/local/intel/composerxe/mkl/include/.f, /usrx/local/intel/composerxe/tbb/include/.f,/gpfs/gp2/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/.f, /gpfs/gp2/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/.f,/usr/local/include/.f,/usr/lib/gcc/x86_64-redhat-linux/4.4.6/include/.f, /usr/include/.f,/usr/include/.f -list filename : interp.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100