Page 1 Source Listing P2FLTVF 2025-03-12 18:24 P2FLTVF.f 1 SUBROUTINE P2FLTVF(MKOUNT,VBM2,FIELD) 2 CFPP$ NOCONCUR R 3 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 4 C . . . 5 C SUBPROGRAM: P2FLTVF 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. DOES FULL FIELD 13 C . 14 C 15 C PROGRAM HISTORY LOG: 16 C 93-03-02 RUSS TREADON - MODIFIED SUBROUTINE P2FILT. 17 C 98-06-04 BLACK - CONVERSION TO 2-D 18 C 19 C USAGE: CALL P2FLTV(MKOUNT,VBM2,FIELD) 20 C INPUT ARGUMENT LIST: 21 C MKOUNT - NUMBER OF APPLICATIONS OF THE SMOOTHER. 22 C HBM2 - VELOCITY POINT BOUNDARY MASK. 23 C FIELD - VELOCITY POINT FIELD TO BE SMOOTHED. 24 C 25 C OUTPUT ARGUMENT LIST: 26 C FIELD - SMOOTHED VELOCITY POINT FIELD. 27 C 28 C OUTPUT FILES: 29 C NONE 30 C 31 C SUBPROGRAMS CALLED: 32 C UTILITIES: 33 C NONE 34 C LIBRARY: 35 C NONE 36 C 37 C ATTRIBUTES: 38 C LANGUAGE: FORTRAN 39 C MACHINE : CRAY Y-MP 40 C$$$ 41 C 42 C 43 C INCLUDE GLOBAL PARAMETERS. SET DEPENDENT VARIABLES. 44 C 45 C---------------------------------------------------------------------- 46 INCLUDE "parmeta" 47 C 65 INCLUDE "INDX.comm" 66 C---------------------------------------------------------------------- 72 D I M E N S I O N 73 & VBM2(IM,JM),FIELD(IM,JM),WFIELD(IM,JM),FNE(IM,JM),FSE(IM,JM) 74 C 75 C*********************************************************************** 76 C START P2FLTV HERE. 77 C 78 C IF NUMBER OF PASSES IS LESS THAN 1, EXIT ROUTINE. 79 C Page 2 Source Listing P2FLTVF 2025-03-12 18:24 P2FLTVF.f 80 IF(MKOUNT.LT.1) RETURN 81 C 82 C OTHERWISE, APPLY SMOOTHER MKOUND TIMES. 83 C 84 C----------------------------------------------------------------------- 85 DO 400 KOUNT=1,MKOUNT 86 C----------------------------------------------------------------------- 87 C 88 C FILTERING STARTS 89 C 90 C----------------------------------------------------------------------- 91 !$omp parallel do 92 DO J=1,JM 93 DO I=1,IM 94 FNE (I,J)=0.0 95 FSE (I,J)=0.0 96 WFIELD(I,J)=0.0 97 ENDDO 98 ENDDO 99 C 100 C 2-ND ORDER DIAGONAL CONTRIBUTIONS 101 C 102 !$omp parallel do 103 DO J=1,JM-1 104 DO I=1,IM-1 105 FNE(I,J)=FIELD(I+IVE(J),J+1)-FIELD(I,J) 106 ENDDO 107 ENDDO 108 C 109 !$omp parallel do 110 DO J=2,JM 111 DO I=1,IM-1 112 FSE(I,J)=FIELD(I+IVE(J),J-1)-FIELD(I,J) 113 ENDDO 114 ENDDO 115 C 116 !$omp parallel do 117 DO J=3,JM-2 118 DO I=2,IM-1 119 WFIELD(I,J)=(FNE(I,J)-FNE(I+IVW(J),J-1) 120 1 +FSE(I,J)-FSE(I+IVW(J),J+1)) 121 2 *VBM2(I,J)*0.125 122 ENDDO 123 ENDDO 124 C 125 C 4-TH ORDER DIAGONAL CONTRIBUTIONS 126 C 127 !$omp parallel do 128 DO J=1,JM-1 129 DO I=1,IM-1 130 FNE(I,J)=WFIELD(I+IHE(J),J+1)-WFIELD(I,J) 131 ENDDO 132 ENDDO 133 C 134 !$omp parallel do 135 DO J=2,JM 136 DO I=1,IM-1 Page 3 Source Listing P2FLTVF 2025-03-12 18:24 P2FLTVF.f 137 FSE(I,J)=WFIELD(I+IVE(J),J-1)-WFIELD(I,J) 138 ENDDO 139 ENDDO 140 C 141 !$omp parallel do 142 DO J=3,JM-2 143 DO I=2,IM-1 144 FIELD(I,J)=FIELD(I,J)-(FNE(I,J)-FNE(I+IVW(J),J-1) 145 1 +FSE(I,J)-FSE(I+IVW(J),J+1)) 146 2 *VBM2(I,J)*0.125 147 ENDDO 148 ENDDO 149 C 150 400 CONTINUE 151 C 152 C END OF ROUTINE 153 RETURN 154 END ENTRY POINTS Name p2fltvf_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 400 Label 150 85 FIELD Dummy 1 R(4) 4 2 91719 ARG,INOUT 105,112,144 FNE Local 73 R(4) 4 2 91719 94,105,119,130,144 FSE Local 73 R(4) 4 2 91719 95,112,120,137,145 I Local 93 I(4) 4 scalar 93,94,95,96,104,105,111,112,118,11 9,120,121,129,130,136,137,143,144, 145,146 IDIM1 Param 61 I(4) 4 scalar IDIM2 Param 61 I(4) 4 scalar IGSTL Param 56 I(4) 4 scalar 61 IGSTR Param 56 I(4) 4 scalar 61 IM Param 50 I(4) 4 scalar 58,61,73,93,104,111,118,129,136,14 3 INDX Common 68 6192 SAVE INPES Param 53 I(4) 4 scalar 58,61 ITAIL Param 58 I(4) 4 scalar J Local 92 I(4) 4 scalar 92,94,95,96,103,105,110,112,117,11 9,120,121,128,130,135,137,142,144, 145,146 JDIM1 Param 62 I(4) 4 scalar JDIM2 Param 62 I(4) 4 scalar JGSTL Param 57 I(4) 4 scalar 62 JGSTR Param 57 I(4) 4 scalar 62 JM Param 50 I(4) 4 scalar 59,62,69,73,92,103,110,117,128,135 ,142 JNPES Param 53 I(4) 4 scalar 59,62 Page 4 Source Listing P2FLTVF 2025-03-12 18:24 Symbol Table P2FLTVF.f Name Object Declared Type Bytes Dimen Elements Attributes References JTAIL Param 59 I(4) 4 scalar KOUNT Local 85 I(4) 4 scalar 85 LM Param 50 I(4) 4 scalar LSM Param 50 I(4) 4 scalar MKOUNT Dummy 1 I(4) 4 scalar ARG,INOUT 80,85 P2FLTVF Subr 1 VBM2 Dummy 1 R(4) 4 2 91719 ARG,INOUT 121,146 WFIELD Local 73 R(4) 4 2 91719 96,119,130,137 TYPE COMPONENTS/COMMON VARIABLES Name Type Bytes Offset Dimen Elements Attributes References IHE I(4) 4 0 1 387 COM 130 IHW I(4) 4 1548 1 387 COM IVE I(4) 4 3096 1 387 COM 105,112,137 IVW I(4) 4 4644 1 387 COM 119,120,144,145 Page 5 Source Listing P2FLTVF 2025-03-12 18:24 Subprograms/Common Blocks P2FLTVF.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References INDX Common 68 6192 SAVE P2FLTVF Subr 1 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__ -D __x86_64 -D __x86_64__ -D __amd64 -D __amd64__ Page 6 Source Listing P2FLTVF 2025-03-12 18:24 P2FLTVF.f -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 -w usage no -wrap-margins Page 7 Source Listing P2FLTVF 2025-03-12 18:24 P2FLTVF.f -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 : P2FLTVF.lst no -o COMPILER: Intel(R) Fortran 19.1-1655