Page 1 Source Listing EFILT 2014-12-17 20:47 EFILT.f 1 SUBROUTINE EFILT(EGRID) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: EFILT HEAVY HANDED E-GRID SMOOTHER 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-03-24 6 C 7 C ABSTRACT: 8 C THIS ROUTINE IS A HEAVY HANDED SMOOTHER THAT WORKS ON 9 C E-GRID MASS POINTS. IT IS BASED ON CODE IN THE OLD 10 C ETA POST PROCESSOR SUBROUTINE, OUTMAP. CURRENTLY, 11 C THIS SMOOTHER IS HARDCODED TO BE USED ONLY ON VORTICITY 12 C FIELDS IF THE USER ACTIVATES E-GRID SMOOTHING OF 13 C VORTICITY. IT DOES A VERY NICE JOB OF SMOOTHING THE 14 C VORTICITY FIELD. 15 C . 16 C 17 C PROGRAM HISTORY LOG: 18 C 93-03-24 RUSS TREADON 19 C 98-06-04 BLACK - CONVERSION TO 2-D 20 C 21 C USAGE: CALL EFILT(EGRID) 22 C INPUT ARGUMENT LIST: 23 C EGRID - ARRAY OF DATA ON E-GRID 24 C 25 C OUTPUT ARGUMENT LIST: 26 C EGRID - ARRAY OF SMOOTHED DATA ON E-GRID 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 COMMON - MASKS 36 C DYNAMD 37 C INDX 38 C 39 C ATTRIBUTES: 40 C LANGUAGE: FORTRAN 41 C MACHINE : CRAY Y-MP 42 C$$$ 43 C 44 C 45 C INCLUDE PARAMETERS. 46 INCLUDE "parmeta" 64 INCLUDE "params" 65 C 66 C DECLARE VARIABLES. 108 REAL EGRID(IM,JM),WORK1(IM,JM),WORK2(IM,JM),WORK3(IM,JM) 109 C 110 C INCLUDE COMMON BLOCKS. 111 INCLUDE "MASKS.comm" 119 INCLUDE "DYNAMD.comm" 131 INCLUDE "INDX.comm" 132 C 133 C******************************************************************* Page 2 Source Listing EFILT 2014-12-17 20:47 EFILT.f 134 C START EFILT HERE. 135 C 136 C LOAD PASSED E-GRID FIELD INTO WORK ARRAY. 137 138 !$omp parallel do 144 DO J=1,JM 145 DO I=1,IM 146 WORK1(I,J)=EGRID(I,J) 147 ENDDO 148 ENDDO 149 C 150 C LOAD WEIGHT ARRAY. 151 !$omp parallel do 152 DO J=3,JM-2 153 DO I=2,IM-1 154 WORK3(I,J)=1./(4.* HBM2(I,J)+ 155 1 2.*(HBM2(I+IHE(J),J+1)+HBM2(I+IHW(J),J+1)+ 156 2 HBM2(I+IHE(J),J-1)+HBM2(I+IHW(J),J-1))+ 157 3 HBM2(I+1,J)+HBM2(I-1,J)+ 158 4 HBM2(I,J+2)+HBM2(I,J-2)) 159 ENDDO 160 ENDDO 161 C 162 C SMOOTHING LOOP. SMTHA AND SMTHB ARE PARAMETERS SET IN 163 C INCLUDE FILE PARAMS. 164 MXFILT = NINT(SMTHA*DY+SMTHB) 165 C 166 DO 50 NFILT = 1,MXFILT 167 !$omp parallel do 168 DO J=3,JM-2 169 DO I=2,IM-1 170 WORK2(I,J) = 4.*WORK1(I,J)+ 171 1 2.*(WORK1(I+IHE(J),J+1)+WORK1(I+IHW(J),J+1)+ 172 2 WORK1(I+IHE(J),J-1)+WORK1(I+IHW(J),J-1))+ 173 3 WORK1(I+1,J)+WORK1(I-1,J)+ 174 4 WORK1(I,J+2)+WORK1(I,J-2) 175 ENDDO 176 ENDDO 177 C 178 !$omp parallel do 179 DO J=3,JM-2 180 DO I=2,IM-1 181 WORK1(I,J)=WORK3(I,J)*WORK2(I,J)*HBM2(I,J) 182 ENDDO 183 ENDDO 184 50 CONTINUE 185 C 186 C LOAD SMOOTHED EGRID INTO OUTPUT ARRAY. 187 !$omp parallel do 188 DO J=1,JM 189 DO I=1,IM 190 EGRID(I,J)=WORK1(I,J)*HBM2(I,J) 191 ENDDO 192 ENDDO 193 C 194 C END OF ROUTINE. 195 RETURN Page 3 Source Listing EFILT 2014-12-17 20:47 EFILT.f 196 END ENTRY POINTS Name efilt_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 50 Label 103 85 A2 Param 41 R(4) 4 scalar A3 Param 41 R(4) 4 scalar A4 Param 41 R(4) 4 scalar AD05 Param 33 R(4) 4 scalar AETAD Scalar 4 R(4) 4 1 45 COM BETA Param 37 R(4) 4 scalar 37 BTG Param 37 R(4) 4 scalar CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CM1 Param 23 R(4) 4 scalar CM2 Param 23 R(4) 4 scalar CM3 Param 23 R(4) 4 scalar CP Param 34 R(4) 4 scalar 36,40 CPDR Param 36 R(4) 4 scalar CPGFU Scalar 8 R(4) 4 2 91719 COM CPGFV Scalar 3 R(4) 4 scalar COM CURV Scalar 8 R(4) 4 2 91719 COM D00 Param 25 R(4) 4 scalar D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar D05 Param 29 R(4) 4 scalar 33 D1 Param 25 R(4) 4 scalar D115 Param 29 R(4) 4 scalar D125 Param 27 R(4) 4 scalar D25 Param 27 R(4) 4 scalar D50 Param 25 R(4) 4 scalar D608 Param 29 R(4) 4 scalar DDMPU Scalar 9 R(4) 4 2 91719 COM DDMPV Scalar 9 R(4) 4 2 91719 COM DETAD Scalar 4 R(4) 4 1 45 COM DFL Scalar 4 R(4) 4 1 46 COM DP10M Param 35 R(4) 4 scalar DTR Param 34 R(4) 4 scalar 34 DX Scalar 8 R(4) 4 2 91719 COM DY Scalar 3 R(4) 4 scalar COM 83 DYNAM Common 2 3675964 EF4T Scalar 3 R(4) 4 scalar COM EFILT Subr 1 EGRID Dummy 1 R(4) 4 2 91719 ARG,INOUT 65,109 ELDR Param 36 R(4) 4 scalar ELWV Param 36 R(4) 4 scalar 36 Page 4 Source Listing EFILT 2014-12-17 20:47 Symbol Table EFILT.f Name Object Declared Type Bytes Dimen Elements Attributes References EM Scalar 6 R(4) 4 1 760 COM EMT Scalar 6 R(4) 4 1 760 COM EN Scalar 3 R(4) 4 scalar COM ENT Scalar 3 R(4) 4 scalar COM EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar ETAD Scalar 4 R(4) 4 1 46 COM F Scalar 9 R(4) 4 2 91719 COM F4D Scalar 3 R(4) 4 scalar COM F4Q Scalar 3 R(4) 4 scalar COM F4Q2 Scalar 4 R(4) 4 1 45 COM FAD Scalar 9 R(4) 4 2 91719 COM FCP Scalar 8 R(4) 4 2 91719 COM FDIV Scalar 9 R(4) 4 2 91719 COM FMD Param 39 R(4) 4 scalar 39 FMW Param 39 R(4) 4 scalar 39 G Param 32 R(4) 4 scalar 32 GI Param 32 R(4) 4 scalar 36 H1 Param 26 R(4) 4 scalar 33 H100 Param 27 R(4) 4 scalar H1000 Param 27 R(4) 4 scalar H10000 Param 28 R(4) 4 scalar H10E5 Param 28 R(4) 4 scalar H1E2 Param 30 R(4) 4 scalar H1E4 Param 30 R(4) 4 scalar H1M12 Param 26 R(4) 4 scalar H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar H2 Param 26 R(4) 4 scalar H4 Param 26 R(4) 4 scalar H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar HBM2 Scalar 3 R(4) 4 2 91719 COM 73,74,75,76,77,100,109 HTM Scalar 5 R(4) 4 3 4127355 COM I Local 64 I(4) 4 scalar 64,65,72,73,74,75,76,77,88,89,90,9 1,92,93,99,100,108,109 IDIM1 Param 14 I(4) 4 scalar IDIM2 Param 14 I(4) 4 scalar 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 74,75,90,91 IHW Scalar 3 I(4) 4 1 387 COM 74,75,90,91 IM Param 3 I(4) 4 scalar 3,4,5,6,7,8,9,11,14,50,64,72,88,99 ,108 IM1 Param 6 I(4) 4 scalar IMT Param 4 I(4) 4 scalar INDX Common 2 6192 INPES Param 6 I(4) 4 scalar 11,14 ITAIL Param 11 I(4) 4 scalar IVE Scalar 3 I(4) 4 1 387 COM IVW Scalar 3 I(4) 4 1 387 COM J Local 63 I(4) 4 scalar 63,65,71,73,74,75,76,77,87,89,90,9 1,92,93,98,100,107,109 JAM Param 6 I(4) 4 scalar 6 Page 5 Source Listing EFILT 2014-12-17 20:47 Symbol Table EFILT.f Name Object Declared Type Bytes Dimen Elements Attributes References JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JGSTL Param 10 I(4) 4 scalar 15 JGSTR Param 10 I(4) 4 scalar 15 JM Param 3 I(4) 4 scalar 3,4,5,6,8,9,12,15,50,63,71,87,98,1 07 JM2 Param 6 I(4) 4 scalar JMT Param 4 I(4) 4 scalar JNPES Param 6 I(4) 4 scalar 12,15 JTAIL Param 12 I(4) 4 scalar KSAV Param 14 I(4) 4 scalar KSFI Param 13 I(4) 4 scalar KSLI Param 14 I(4) 4 scalar KSLP Param 14 I(4) 4 scalar KSLPD Param 13 I(4) 4 scalar KSMUD Param 13 I(4) 4 scalar KSRH Param 14 I(4) 4 scalar KSSLP Param 13 I(4) 4 scalar KST Param 13 I(4) 4 scalar KSUV Param 14 I(4) 4 scalar LB Param 9 I(4) 4 scalar LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LM Param 3 I(4) 4 scalar 4,5,8,19 LM1 Param 8 I(4) 4 scalar LP1 Param 8 I(4) 4 scalar 4 LSM Param 3 I(4) 4 scalar MASKS Common 2 34853220 MXFILT Local 83 I(4) 4 scalar 83,85 NFILT Local 85 I(4) 4 scalar 85 NINC Param 7 I(4) 4 scalar NINT Func 83 scalar 83 NSTAT Param 19 I(4) 4 scalar NSUMD Param 15 I(4) 4 scalar NWORD Param 19 I(4) 4 scalar 19 NWORD3 Param 19 I(4) 4 scalar NWORD6 Param 19 I(4) 4 scalar ONEPS Param 39 R(4) 4 scalar OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar PQ0 Param 41 R(4) 4 scalar PTDYN Scalar 3 R(4) 4 scalar COM RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDETA Scalar 4 R(4) 4 1 45 COM RDOCP Param 40 R(4) 4 scalar RDYN Scalar 3 R(4) 4 scalar COM RG Param 36 R(4) 4 scalar ROG Param 32 R(4) 4 scalar RTD Param 34 R(4) 4 scalar SICE Scalar 3 R(4) 4 2 91719 COM SM Scalar 3 R(4) 4 2 91719 COM SMTHA Param 15 R(4) 4 scalar 83 SMTHB Param 15 R(4) 4 scalar 83 TFRZ Param 39 R(4) 4 scalar Page 6 Source Listing EFILT 2014-12-17 20:47 Symbol Table EFILT.f Name Object Declared Type Bytes Dimen Elements Attributes References VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VTM Scalar 5 R(4) 4 3 4127355 COM WORK1 Local 50 R(4) 4 2 91719 65,89,90,91,92,93,100,109 WORK2 Local 50 R(4) 4 2 91719 89,100 WORK3 Local 50 R(4) 4 2 91719 73,100 WPDAR Scalar 8 R(4) 4 2 91719 COM Page 7 Source Listing EFILT 2014-12-17 20:47 EFILT.f 197 198 199 Page 8 Source Listing EFILT 2014-12-17 20:47 Subprograms/Common Blocks EFILT.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References DYNAM Common 2 3675964 EFILT Subr 1 INDX Common 2 6192 MASKS Common 2 34853220 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 9 Source Listing EFILT 2014-12-17 20:47 EFILT.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 : EFILT.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100