Page 1 Source Listing VDIFQ 2014-12-17 20:47 /tmpnwprd/ifort8T41Y9.i 1 # 1 "VDIFQ.F" 2 C&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 3 SUBROUTINE VDIFQ(LMHK,KTM,DTQ2,Q2,EL,Z) 4 C ****************************************************************** 5 C * * 6 C * VERTICAL DIFFUSION * 7 C * * 8 C ****************************************************************** 9 C----------------------------------------------------------------------- 10 INCLUDE "parmeta" 11 # 1 "./sp.h" 1 12 # 4 13 14 # 11 "VDIFQ.F" 2 15 C----------------------------------------------------------------------- 38 P A R A M E T E R 39 &(LP1=LM+1,LM1=LM-1,LM2=LM-2) 40 P A R A M E T E R 41 &(ESQ=0.20,ELZ0=0.) 42 C----------------------------------------------------------------------- 43 D I M E N S I O N 44 & Q2 (LM) 45 D I M E N S I O N 46 & EL (LM1) 47 &,Z (LP1) 48 D I M E N S I O N 49 & CM (LM2),CR (LM2),RSQ2 (LM2),AKQ (LM2),DTOZ (LM2) 50 C----------------------------------------------------------------------- 51 C*********************************************************************** 52 DTDIF=DTQ2/FLOAT(KTM) 53 LMHM=LMHK-1 54 LMH2=LMHK-2 55 LMHP=LMHK+1 56 C----------------------------------------------------------------------- 57 DO 300 KT=1,KTM 58 C----------------------------------------------------------------------- 59 DO 100 L=1,LMH2 60 DTOZ(L)=(DTDIF+DTDIF)/(Z(L)-Z(L+2)) 61 AKQ(L)=SQRT((Q2(L)+Q2(L+1))*0.5)*(EL(L)+EL(L+1))*0.5*ESQ 62 & /(Z(L+1)-Z(L+2)) 63 CR(L)=-DTOZ(L)*AKQ(L) 64 100 CONTINUE 65 C 66 CM(1)=DTOZ(1)*AKQ(1)+1. 67 RSQ2(1)=Q2(1) 68 C----------------------------------------------------------------------- 69 DO 110 L=2,LMH2 70 CF=-DTOZ(L)*AKQ(L-1)/CM(L-1) 71 CM(L)=-CR(L-1)*CF+(AKQ(L-1)+AKQ(L))*DTOZ(L)+1. 72 RSQ2(L)=-RSQ2(L-1)*CF+Q2(L) 73 110 CONTINUE 74 C----------------------------------------------------------------------- 75 DTOZS=(DTDIF+DTDIF)/(Z(LMHM)-Z(LMHP)) 76 AKQS=SQRT((Q2(LMHM)+Q2(LMHK))*0.5)*(EL(LMHM)+ELZ0)*0.5*ESQ 77 & /(Z(LMHK)-Z(LMHP)) 78 C 79 CF=-DTOZS*AKQ(LMH2)/CM(LMH2) Page 2 Source Listing VDIFQ 2014-12-17 20:47 VDIFQ.F 80 C----------------------------------------------------------------------- 81 Q2(LMHM)=(DTOZS*AKQS*Q2(LMHK)-RSQ2(LMH2)*CF+Q2(LMHM)) 82 & /((AKQ(LMH2)+AKQS)*DTOZS-CR(LMH2)*CF+1.) 83 C 84 DO 120 IVI=1,LMH2 85 L=LMHM-IVI 86 Q2(L)=(-CR(L)*Q2(L+1)+RSQ2(L))/CM(L) 87 120 CONTINUE 88 C----------------------------------------------------------------------- 89 300 CONTINUE 90 C----------------------------------------------------------------------- 91 RETURN 92 END ENTRY POINTS Name vdifq_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 100 Label 38 33 110 Label 47 43 120 Label 61 58 300 Label 63 31 AKQ Local 23 R(4) 4 1 43 35,37,40,44,45,53,56 AKQS Local 50 R(4) 4 scalar 50,55,56 CF Local 44 R(4) 4 scalar 44,45,46,53,55,56 CM Local 23 R(4) 4 1 43 40,44,45,53,60 CR Local 23 R(4) 4 1 43 37,45,56,60 DTDIF Local 26 R(4) 4 scalar 26,34,49 DTOZ Local 23 R(4) 4 1 43 34,37,40,44,45 DTOZS Local 49 R(4) 4 scalar 49,53,55,56 DTQ2 Dummy 2 R(4) 4 scalar ARG,INOUT 26 EL Dummy 2 R(4) 4 1 44 ARG,INOUT 35,50 ELZ0 Param 15 R(4) 4 scalar 50 ESQ Param 15 R(4) 4 scalar 35,50 FLOAT Func 26 scalar 26 IDIM1 Param 19 I(4) 4 scalar IDIM2 Param 19 I(4) 4 scalar IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 IM Param 3 I(4) 4 scalar 16,19 INPES Param 11 I(4) 4 scalar 16,19 ITAIL Param 16 I(4) 4 scalar IVI Local 58 I(4) 4 scalar 58,59 JDIM1 Param 20 I(4) 4 scalar JDIM2 Param 20 I(4) 4 scalar JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 JM Param 3 I(4) 4 scalar 17,20 JNPES Param 11 I(4) 4 scalar 17,20 Page 3 Source Listing VDIFQ 2014-12-17 20:47 Symbol Table VDIFQ.F Name Object Declared Type Bytes Dimen Elements Attributes References JTAIL Param 17 I(4) 4 scalar KT Local 31 I(4) 4 scalar 31 KTM Dummy 2 I(4) 4 scalar ARG,INOUT 26,31 L Local 33 I(4) 4 scalar 33,34,35,36,37,43,44,45,46,59,60 LM Param 3 I(4) 4 scalar 13,18 LM1 Param 13 I(4) 4 scalar 20 LM2 Param 13 I(4) 4 scalar 23 LMH2 Local 28 I(4) 4 scalar 28,33,43,53,55,56,58 LMHK Dummy 2 I(4) 4 scalar ARG,INOUT 27,28,29,50,51,55 LMHM Local 27 I(4) 4 scalar 27,49,50,55,59 LMHP Local 29 I(4) 4 scalar 29,49,51 LP1 Param 13 I(4) 4 scalar 21 LSM Param 3 I(4) 4 scalar Q2 Dummy 2 R(4) 4 1 45 ARG,INOUT 35,41,46,50,55,60 RSQ2 Local 23 R(4) 4 1 43 41,46,55,60 SQRT Func 35 scalar 35,50 VDIFQ Subr 2 Z Dummy 2 R(4) 4 1 46 ARG,INOUT 34,36,49,51 Page 4 Source Listing VDIFQ 2014-12-17 20:47 VDIFQ.F 93 Page 5 Source Listing VDIFQ 2014-12-17 20:47 Subprograms/Common Blocks VDIFQ.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References VDIFQ Subr 2 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 -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest Page 6 Source Listing VDIFQ 2014-12-17 20:47 VDIFQ.F no -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 nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_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 nono_precise_libm_functions -heap_arrays 0 no -threadprivate_compat -g2 -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 : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : VDIFQ.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100