Page 1 Source Listing P2FILT 2014-12-17 20:48 P2FILT.f 1 SUBROUTINE P2FILT(MKOUNT,HBM2,FIELD) 2 CFPP$ NOCONCUR R 3 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 4 C . . . 5 C SUBPROGRAM: P2FILT 4-TH ORDER SMOOTHER FOR MASS 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 MASS POINTS ON THE STAGGERED 11 C E-GRID. MULTIPLE APPLICATIONS OF THIS FILTER MAY BE 12 C MADE. 13 C 14 C IN THIS VERSION, EACH MPI TASK WORKS ON ITS OWN DATA. 15 C . 16 C 17 C PROGRAM HISTORY LOG: 18 C ??-??-?? ??? - SUBROUTINE P2FILT IN ETA MODEL. 19 C 93-03-02 RUSS TREADON - ADDED DOCBLOC 20 C 98-06-04 BLACK - CONVERSION TO 2-D 21 C 00-01-04 JIM TUCCILLO - MPI VERSION 22 C 23 C USAGE: CALL P2FILT(MKOUNT,HBM2,FIELD) 24 C INPUT ARGUMENT LIST: 25 C MKOUNT - NUMBER OF APPLICATIONS OF THE SMOOTHER. 26 C HBM2 - MASS POINT BOUNDARY MASK. 27 C FIELD - MASS POINT FIELD TO BE SMOOTHED. 28 C 29 C OUTPUT ARGUMENT LIST: 30 C FIELD - SMOOTHED MASS POINT FIELD. 31 C 32 C OUTPUT FILES: 33 C NONE 34 C 35 C SUBPROGRAMS CALLED: 36 C UTILITIES: 37 C NONE 38 C LIBRARY: 39 C NONE 40 C 41 C ATTRIBUTES: 42 C LANGUAGE: FORTRAN 43 C MACHINE : CRAY C-90 44 C$$$ 45 C 46 C 47 C INCLUDE GLOBAL PARAMETERS. SET DEPENDENT VARIABLES. 48 C 49 C----------------------------------------------------------------------- 50 INCLUDE "parmeta" 51 C----------------------------------------------------------------------- 69 INCLUDE "INDX.comm" 75 INCLUDE "CTLBLK.comm" 76 C----------------------------------------------------------------------- 91 D I M E N S I O N 92 1 HBM2(IM,JM),FIELD(IM,JM),WFIELD(IM,JM),FNE(IM,JM),FSE(IM,JM) 93 C Page 2 Source Listing P2FILT 2014-12-17 20:48 P2FILT.f 94 C*********************************************************************** 95 C START P2FILT HERE. 96 C 97 C IF NUMBER OF PASSES IS LESS THAN 1, EXIT ROUTINE. 98 C 99 IF(MKOUNT.LT.1) RETURN 100 C 101 C OTHERWISE, APPLY SMOOTHER MKOUNT TIMES. 102 C 103 DO 400 KOUNT=1,MKOUNT 104 C 105 C FILTERING STARTS 106 C 107 !$omp parallel do 108 DO J=JSTA,JEND 109 DO I=1,IM 110 FNE (I,J)=0.0 111 FSE (I,J)=0.0 112 WFIELD(I,J)=0.0 113 ENDDO 114 ENDDO 115 C 116 C 2-ND ORDER DIAGONAL CONTRIBUTIONS 117 C 118 CALL EXCH(FIELD) 119 C 120 !$omp parallel do 121 DO J=JSTA,JEND_M 122 DO I=1,IM-1 123 FNE(I,J)=FIELD(I+IHE(J),J+1)-FIELD(I,J) 124 ENDDO 125 ENDDO 126 C 127 !$omp parallel do 128 DO J=JSTA_M,JEND 129 DO I=1,IM-1 130 FSE(I,J)=FIELD(I+IHE(J),J-1)-FIELD(I,J) 131 ENDDO 132 ENDDO 133 C 134 CALL EXCH(FNE) 135 CALL EXCH(FSE) 136 C 137 !$omp parallel do 138 DO J=JSTA_M2,JEND_M2 139 DO I=2,IM-1 140 WFIELD(I,J)=(FNE(I,J)-FNE(I+IHW(J),J-1) 141 1 +FSE(I,J)-FSE(I+IHW(J),J+1))*HBM2(I,J)*0.125 142 ENDDO 143 ENDDO 144 C 145 C 4-TH ORDER DIAGONAL CONTRIBUTIONS 146 C 147 CALL EXCH(WFIELD) 148 C 149 !$omp parallel do 150 DO J=JSTA,JEND_M Page 3 Source Listing P2FILT 2014-12-17 20:48 P2FILT.f 151 DO I=1,IM-1 152 FNE(I,J)=WFIELD(I+IHE(J),J+1)-WFIELD(I,J) 153 ENDDO 154 ENDDO 155 C 156 !$omp parallel do 157 DO J=JSTA_M,JEND 158 DO I=1,IM-1 159 FSE(I,J)=WFIELD(I+IHE(J),J-1)-WFIELD(I,J) 160 ENDDO 161 ENDDO 162 C 163 CALL EXCH(FNE) 164 CALL EXCH(FSE) 165 C 166 !$omp parallel do 167 DO J=JSTA_M2,JEND_M2 168 DO I=2,IM-1 169 FIELD(I,J)=FIELD(I,J)-(FNE(I,J)-FNE(I+IHW(J),J-1) 170 1 +FSE(I,J)-FSE(I+IHW(J),J+1)) 171 2 *HBM2(I,J)*0.125 172 ENDDO 173 ENDDO 174 C 175 400 CONTINUE 176 C 177 CALL EXCH(FIELD) 178 C 179 C END OF ROUTINE 180 RETURN 181 END Page 4 Source Listing P2FILT 2014-12-17 20:48 Entry Points P2FILT.f ENTRY POINTS Name p2filt_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 400 Label 139 67 CTLBLK Common 2 96 DT Scalar 4 R(4) 4 scalar COM EXCH Subr 82 82,98,99,111,127,128,141 FIELD Dummy 1 R(4) 4 2 91719 ARG,INOUT 82,87,94,133,141 FIRST Scalar 3 R(4) 4 scalar COM FNE Local 56 R(4) 4 2 91719 74,87,98,104,116,127,133 FSE Local 56 R(4) 4 2 91719 75,94,99,105,123,128,134 HBM2 Dummy 1 R(4) 4 2 91719 ARG,INOUT 105,135 I Local 73 I(4) 4 scalar 73,74,75,76,86,87,93,94,103,104,10 5,115,116,122,123,132,133,134,135 ICNT Scalar 7 I(4) 4 1 1024 COM IDAT Scalar 3 I(4) 4 1 3 COM IDIM1 Param 14 I(4) 4 scalar IDIM2 Param 14 I(4) 4 scalar IDN Scalar 7 I(4) 4 scalar COM IDSP Scalar 7 I(4) 4 1 1024 COM IDTAD Scalar 4 I(4) 4 scalar COM IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 IHE Scalar 3 I(4) 4 1 387 COM 87,94,116,123 IHRST Scalar 3 I(4) 4 scalar COM IHW Scalar 3 I(4) 4 1 387 COM 104,105,133,134 IM Param 3 I(4) 4 scalar 11,14,56,73,86,93,103,115,122,132 INDX Common 2 6192 INPES Param 6 I(4) 4 scalar 11,14 IOUT Scalar 3 I(4) 4 scalar COM ITAIL Param 11 I(4) 4 scalar IUP Scalar 7 I(4) 4 scalar COM IVE Scalar 3 I(4) 4 1 387 COM IVW Scalar 3 I(4) 4 1 387 COM J Local 72 I(4) 4 scalar 72,74,75,76,85,87,92,94,102,104,10 5,114,116,121,123,131,133,134,135 JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 72,92,121 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 85,114 JEND_M2 Scalar 7 I(4) 4 scalar COM 102,131 JGSTL Param 10 I(4) 4 scalar 15 JGSTR Param 10 I(4) 4 scalar 15 JM Param 3 I(4) 4 scalar 3,12,15,56 JNPES Param 6 I(4) 4 scalar 12,15 JSTA Scalar 6 I(4) 4 scalar COM 72,85,114 Page 5 Source Listing P2FILT 2014-12-17 20:48 Symbol Table P2FILT.f Name Object Declared Type Bytes Dimen Elements Attributes References JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 92,121 JSTA_M2 Scalar 7 I(4) 4 scalar COM 102,131 JTAIL Param 12 I(4) 4 scalar KOUNT Local 67 I(4) 4 scalar 67 LIST Scalar 3 I(4) 4 scalar COM LM Param 3 I(4) 4 scalar LSM Param 3 I(4) 4 scalar ME Scalar 6 I(4) 4 scalar COM MKOUNT Dummy 1 I(4) 4 scalar ARG,INOUT 63,67 NBC Scalar 3 I(4) 4 scalar COM NBOCO Scalar 4 I(4) 4 scalar COM NCP Scalar 4 I(4) 4 scalar COM NDDAMP Scalar 4 I(4) 4 scalar COM NFCST Scalar 3 I(4) 4 scalar COM NPHS Scalar 5 I(4) 4 scalar COM NPREC Scalar 4 I(4) 4 scalar COM NRADL Scalar 5 I(4) 4 scalar COM NRADS Scalar 5 I(4) 4 scalar COM NSHDE Scalar 4 I(4) 4 scalar COM NTSD Scalar 4 I(4) 4 scalar COM NTSTM Scalar 4 I(4) 4 scalar COM NUM_PROCS Scalar 6 I(4) 4 scalar COM P2FILT Subr 1 PARALLEL Common 9 8240 RESTRT Scalar 3 R(4) 4 scalar COM RUN Scalar 3 R(4) 4 scalar COM SIGMA Scalar 3 R(4) 4 scalar COM WFIELD Local 56 R(4) 4 2 91719 76,104,111,116,123 Page 6 Source Listing P2FILT 2014-12-17 20:48 Subprograms/Common Blocks P2FILT.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CTLBLK Common 2 96 INDX Common 2 6192 P2FILT Subr 1 PARALLEL Common 9 8240 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 Page 7 Source Listing P2FILT 2014-12-17 20:48 P2FILT.f -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 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/tp1/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/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/.f, /gpfs/tp1/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.7/include/.f, /usr/include/.f,/usr/include/.f -list filename : P2FILT.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100