Page 1 Source Listing P2FLTV 2025-03-12 18:23 P2FLTV.f 1 SUBROUTINE P2FLTV(MKOUNT,VBM2,FIELD) 2 CFPP$ NOCONCUR R 3 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 4 C . . . 5 C SUBPROGRAM: P2FLTV 4-TH ORDER SMOOTHER for vel pts 6 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-03-02 7 C 8 C ABSTRACT: 9 C THIS ROUTINE WILL APPLY A FOURTH OTHER FILTER (SMOOTHER) 10 C OF A CONTINOUS FIELD AT VELOCITY POINTS ON THE STAGGERED 11 C E-GRID. MULTIPLE APPLICATIONS OF THIS FILTER MAY BE 12 C MADE. 13 C IN THIS VERSION, EACH TASK WORKS ON THE DATA IT OWNS. 14 C . 15 C 16 C PROGRAM HISTORY LOG: 17 C 93-03-02 RUSS TREADON - MODIFIED SUBROUTINE P2FILT. 18 C 98-06-04 BLACK - CONVERSION TO 2-D 19 C 00-01-04 JIM TUCCILLO - MPI VERSION 20 C 21 C USAGE: CALL P2FLTV(MKOUNT,VBM2,FIELD) 22 C INPUT ARGUMENT LIST: 23 C MKOUNT - NUMBER OF APPLICATIONS OF THE SMOOTHER. 24 C HBM2 - VELOCITY POINT BOUNDARY MASK. 25 C FIELD - VELOCITY POINT FIELD TO BE SMOOTHED. 26 C 27 C OUTPUT ARGUMENT LIST: 28 C FIELD - SMOOTHED VELOCITY POINT FIELD. 29 C 30 C OUTPUT FILES: 31 C NONE 32 C 33 C SUBPROGRAMS CALLED: 34 C UTILITIES: 35 C NONE 36 C LIBRARY: 37 C NONE 38 C 39 C ATTRIBUTES: 40 C LANGUAGE: FORTRAN 41 C MACHINE : CRAY Y-MP 42 C$$$ 43 C 44 C 45 C INCLUDE GLOBAL PARAMETERS. SET DEPENDENT VARIABLES. 46 C 47 C---------------------------------------------------------------------- 48 INCLUDE "parmeta" 49 C 67 INCLUDE "INDX.comm" 73 INCLUDE "CTLBLK.comm" 74 C---------------------------------------------------------------------- 89 D I M E N S I O N 90 & VBM2(IM,JM),FIELD(IM,JM),WFIELD(IM,JM),FNE(IM,JM),FSE(IM,JM) 91 C 92 C*********************************************************************** 93 C START P2FLTV HERE. Page 2 Source Listing P2FLTV 2025-03-12 18:23 P2FLTV.f 94 C 95 C IF NUMBER OF PASSES IS LESS THAN 1, EXIT ROUTINE. 96 C 97 IF(MKOUNT.LT.1) RETURN 98 C 99 C OTHERWISE, APPLY SMOOTHER MKOUND TIMES. 100 C 101 C----------------------------------------------------------------------- 102 DO 400 KOUNT=1,MKOUNT 103 C----------------------------------------------------------------------- 104 C 105 C FILTERING STARTS 106 C 107 C----------------------------------------------------------------------- 108 !$omp parallel do 109 DO J=JSTA,JEND 110 DO I=1,IM 111 FNE (I,J)=0.0 112 FSE (I,J)=0.0 113 WFIELD(I,J)=0.0 114 ENDDO 115 ENDDO 116 C 117 C 118 CALL EXCH(FIELD) 119 C 120 C 2-ND ORDER DIAGONAL CONTRIBUTIONS 121 C 122 !$omp parallel do 123 DO J=JSTA,JEND_M 124 DO I=1,IM-1 125 FNE(I,J)=FIELD(I+IVE(J),J+1)-FIELD(I,J) 126 ENDDO 127 ENDDO 128 C 129 !$omp parallel do 130 DO J=JSTA_M,JEND 131 DO I=1,IM-1 132 FSE(I,J)=FIELD(I+IVE(J),J-1)-FIELD(I,J) 133 ENDDO 134 ENDDO 135 C 136 CALL EXCH(FNE) 137 CALL EXCH(FSE) 138 C 139 !$omp parallel do 140 DO J=JSTA_M2,JEND_M2 141 DO I=2,IM-1 142 WFIELD(I,J)=(FNE(I,J)-FNE(I+IVW(J),J-1) 143 1 +FSE(I,J)-FSE(I+IVW(J),J+1)) 144 2 *VBM2(I,J)*0.125 145 ENDDO 146 ENDDO 147 C 148 C 4-TH ORDER DIAGONAL CONTRIBUTIONS 149 C 150 CALL EXCH(WFIELD) Page 3 Source Listing P2FLTV 2025-03-12 18:23 P2FLTV.f 151 C 152 !$omp parallel do 153 DO J=JSTA,JEND_M 154 DO I=1,IM-1 155 FNE(I,J)=WFIELD(I+IHE(J),J+1)-WFIELD(I,J) 156 ENDDO 157 ENDDO 158 C 159 !$omp parallel do 160 DO J=JSTA_M,JEND 161 DO I=1,IM-1 162 FSE(I,J)=WFIELD(I+IVE(J),J-1)-WFIELD(I,J) 163 ENDDO 164 ENDDO 165 C 166 CALL EXCH(FNE) 167 CALL EXCH(FSE) 168 C 169 !$omp parallel do 170 DO J=JSTA_M2,JEND_M2 171 DO I=2,IM-1 172 FIELD(I,J)=FIELD(I,J)-(FNE(I,J)-FNE(I+IVW(J),J-1) 173 1 +FSE(I,J)-FSE(I+IVW(J),J+1)) 174 2 *VBM2(I,J)*0.125 175 ENDDO 176 ENDDO 177 C 178 400 CONTINUE 179 C 180 CALL EXCH(FIELD) 181 C 182 C END OF ROUTINE 183 RETURN 184 END Page 4 Source Listing P2FLTV 2025-03-12 18:23 Entry Points P2FLTV.f ENTRY POINTS Name p2fltv_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 400 Label 178 102 CTLBLK Common 76 96 SAVE EXCH Subr 118 118,136,137,150,166,167,180 FIELD Dummy 1 R(4) 4 2 91719 ARG,INOUT 118,125,132,172,180 FNE Local 90 R(4) 4 2 91719 111,125,136,142,155,166,172 FSE Local 90 R(4) 4 2 91719 112,132,137,143,162,167,173 I Local 110 I(4) 4 scalar 110,111,112,113,124,125,131,132,14 1,142,143,144,154,155,161,162,171, 172,173,174 IDIM1 Param 63 I(4) 4 scalar IDIM2 Param 63 I(4) 4 scalar IGSTL Param 58 I(4) 4 scalar 63 IGSTR Param 58 I(4) 4 scalar 63 IM Param 52 I(4) 4 scalar 60,63,90,110,124,131,141,154,161,1 71 INDX Common 69 6192 SAVE INPES Param 55 I(4) 4 scalar 60,63 ITAIL Param 60 I(4) 4 scalar J Local 109 I(4) 4 scalar 109,111,112,113,123,125,130,132,14 0,142,143,144,153,155,160,162,170, 172,173,174 JDIM1 Param 64 I(4) 4 scalar JDIM2 Param 64 I(4) 4 scalar JGSTL Param 59 I(4) 4 scalar 64 JGSTR Param 59 I(4) 4 scalar 64 JM Param 52 I(4) 4 scalar 61,64,70,90 JNPES Param 55 I(4) 4 scalar 61,64 JTAIL Param 61 I(4) 4 scalar KOUNT Local 102 I(4) 4 scalar 102 LM Param 52 I(4) 4 scalar LSM Param 52 I(4) 4 scalar MKOUNT Dummy 1 I(4) 4 scalar ARG,INOUT 97,102 P2FLTV Subr 1 PARALLEL Common 83 8240 VBM2 Dummy 1 R(4) 4 2 91719 ARG,INOUT 144,174 WFIELD Local 90 R(4) 4 2 91719 113,142,150,155,162 Page 5 Source Listing P2FLTV 2025-03-12 18:23 Symbol Table P2FLTV.f TYPE COMPONENTS/COMMON VARIABLES Name Type Bytes Offset Dimen Elements Attributes References DT R(4) 4 48 scalar COM FIRST R(4) 4 24 scalar COM ICNT I(4) 4 40 1 1024 COM IDAT I(4) 4 8 1 3 COM IDN I(4) 4 36 scalar COM IDSP I(4) 4 4136 1 1024 COM IDTAD I(4) 4 68 scalar COM IHE I(4) 4 0 1 387 COM 155 IHRST I(4) 4 20 scalar COM IHW I(4) 4 1548 1 387 COM IOUT I(4) 4 44 scalar COM IUP I(4) 4 32 scalar COM IVE I(4) 4 3096 1 387 COM 125,132,162 IVW I(4) 4 4644 1 387 COM 142,143,172,173 JEND I(4) 4 12 scalar COM 109,130,160 JEND_2U I(4) 4 8236 scalar COM JEND_M I(4) 4 20 scalar COM 123,153 JEND_M2 I(4) 4 28 scalar COM 140,170 JSTA I(4) 4 8 scalar COM 109,123,153 JSTA_2L I(4) 4 8232 scalar COM JSTA_M I(4) 4 16 scalar COM 130,160 JSTA_M2 I(4) 4 24 scalar COM 140,170 LIST I(4) 4 40 scalar COM ME I(4) 4 4 scalar COM NBC I(4) 4 36 scalar COM NBOCO I(4) 4 72 scalar COM NCP I(4) 4 80 scalar COM NDDAMP I(4) 4 60 scalar COM NFCST I(4) 4 32 scalar COM NPHS I(4) 4 84 scalar COM NPREC I(4) 4 64 scalar COM NRADL I(4) 4 92 scalar COM NRADS I(4) 4 88 scalar COM NSHDE I(4) 4 76 scalar COM NTSD I(4) 4 52 scalar COM NTSTM I(4) 4 56 scalar COM NUM_PROCS I(4) 4 0 scalar COM RESTRT R(4) 4 28 scalar COM RUN R(4) 4 4 scalar COM SIGMA R(4) 4 0 scalar COM Page 6 Source Listing P2FLTV 2025-03-12 18:23 Subprograms/Common Blocks P2FLTV.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CTLBLK Common 76 96 SAVE INDX Common 69 6192 SAVE P2FLTV Subr 1 PARALLEL Common 83 8240 COMPILER OPTIONS BEING USED -align noall -align nonone -align nocommons -align nodcommons -align noqcommons -align nozcommons -align records -align nosequence -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 -D linux -D __ELF__ Page 7 Source Listing P2FLTV 2025-03-12 18:23 P2FLTV.f -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__ -D __CRAY_X86_ROME -D __CRAYXT_COMPUTE_LINUX_TARGET -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 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 -g2 -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 -O2 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 Page 8 Source Listing P2FLTV 2025-03-12 18:23 P2FLTV.f -w usage no -wrap-margins -includepath : /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/, .f90,./.f90,/opt/cray/pe/mpich/8.1.12/ofi/intel/19.0/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/.f90,/usr/lib64/gcc/x86_64-suse-linux/7/include/.f90, /usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/.f90,/usr/include/.f90,/usr/include/.f90,/usr/include/.f90 -list filename : P2FLTV.lst no -o COMPILER: Intel(R) Fortran 19.1-1655