Page 1 Source Listing CALMXW 2014-12-17 20:47 CALMXW.f 1 SUBROUTINE CALMXW(MXWP,MXWZ,MXWU,MXWV) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CALMXW COMPUTE MAX WIND LEVEL 5 C PRGRMMR: MANIKIN ORG: W/NP2 DATE: 97-03-04 6 C 7 C ABSTRACT: 8 C THIS ROUTINE COMPUTES MAX WIND LEVEL. AT EACH POINT, 9 C IT FINDS THE MAX WIND ABOVE 500 MB AND DETERMINES THE 10 C PRESSURE AND HEIGHT AT THAT LEVEL. 11 C 12 C 13 C PROGRAM HISTORY LOG: 14 C 97-03-04 GEOFF MANIKIN 15 C 98-06-15 T BLACK - CONVERSION FROM 1-D TO 2-D 16 C 00-01-02 JIM TUCCILLO - MPI VERSION 17 C 18 C USAGE: CALL CALMXW(MXWP,MXWZ,MXWU,MXWV) 19 C INPUT ARGUMENT LIST: 20 C NONE 21 C 22 C OUTPUT ARGUMENT LIST: 23 C MXWP - PRESSURE LEVEL OF THE MAX WIND 24 C MXWZ - HEIGHT OF THE MAX WIND 25 C MXWU - U COMPONENT OF THE ACTUAL MAX WIND 26 C MXWV - V COMPONENT OF THE ACTUAL MAX WIND 27 C 28 C OUTPUT FILES: 29 C NONE 30 C 31 C SUBPROGRAMS CALLED: 32 C UTILITIES: 33 C H2V 34 C 35 C LIBRARY: 36 C COMMON - VRBLS 37 C LOOPS 38 C EXTRA 39 C OPTIONS 40 C MASKS 41 C INDX 42 C 43 C ATTRIBUTES: 44 C LANGUAGE: FORTRAN 90 45 C MACHINE : CRAY C-90 46 C$$$ 47 C 48 C 49 C INCLUDE ETA GRID DIMENSIONS. SET/DERIVE PARAMETERS. 50 C 51 INCLUDE "parmeta" 69 INCLUDE "params" 70 C 112 INCLUDE "MASKS.comm" 120 INCLUDE "VRBLS.comm" 127 INCLUDE "OPTIONS.comm" 131 INCLUDE "LOOPS.comm" Page 2 Source Listing CALMXW 2014-12-17 20:47 CALMXW.f 139 INCLUDE "EXTRA.comm" 151 INCLUDE "INDX.comm" 157 INCLUDE "CTLBLK.comm" 158 C 159 C DECLARE VARIABLES. 160 C 175 REAL MXWP(IM,JM),MXWZ(IM,JM),MXWU(IM,JM),MXWV(IM,JM),MXWW 176 C 177 C 178 C***************************************************************************** 179 C START CALMXW HERE. 180 C 181 C LOOP OVER THE GRID. 182 C 183 CRITP=5.0E4 184 C 185 DO J=JSTA,JEND 186 DO I=1,IM 187 MXWU(I,J) = SPVAL 188 MXWV(I,J) = SPVAL 189 MXWP(I,J) = SPVAL 190 MXWZ(I,J) = SPVAL 191 ENDDO 192 ENDDO 193 C 194 C ASSUME THAT U AND V HAVE UPDATED HALOS 195 C 196 !$omp parallel do 197 !$omp& private(ie,iw,mxww,u0,v0,wind) 198 DO 20 J=JSTA_M,JEND_M 199 DO 20 I=2,IM-1 200 IE=I+IHE(J) 201 IW=I+IHW(J) 202 MXWW = -1000. 203 LLMH=LMH(I,J) 204 C 205 DO 10 L= LLMH-1,1,-1 206 U0 = D25*(U(I,J-1,L)+U(IW,J,L)+ 207 X U(IE,J,L)+U(I,J+1,L)) 208 V0 = D25*(V(I,J-1,L)+V(IW,J,L)+ 209 X V(IE,J,L)+V(I,J+1,L)) 210 WIND = SQRT(U0**2 + V0**2) 211 212 C MAX WIND LEVEL MUST BE ABOVE THE 500 MB 213 C ***NOTE, HOWEVER, THAT THE CHECK FOR 500 IS WITH THE 214 C INTERFACE PRESSURE. THEREFORE, WHEN THE PRESSURE OF 215 C THE ACTUAL LEVEL OF THE MAX WIND IT COMPUTED, IT MAY 216 C END UP WITH A PRESSURE SLIGHTLY GREATER THAN 500 217 218 IF (WIND .GT. MXWW .and. PINT(I,J,L) .LT. CRITP) THEN 219 MXWU(I,J) = U0 220 MXWV(I,J) = V0 221 MXWW = WIND 222 MXWP(I,J) = (PINT(I,J,L) + PINT(I,J,L+1)) * 0.5 223 MXWZ(I,J)=HTM(I,J,L+1)*T(I,J,L+1)*(Q(I,J,L+1)*D608+H1)*ROG* 224 X (LOG(PINT(I,J,L+1))-LOG(MXWP(I,J)))+ZINT(I,J,L+1) 225 ENDIF Page 3 Source Listing CALMXW 2014-12-17 20:47 CALMXW.f 226 10 CONTINUE 227 20 CONTINUE 228 229 C END OF ROUTINE. 230 C 231 RETURN 232 END ENTRY POINTS Name calmxw_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 115 94 20 Label 116 87,88 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 ALPINT Scalar 4 R(4) 4 3 4219074 COM BETA Param 37 R(4) 4 scalar 37 BTG Param 37 R(4) 4 scalar CALMXW Subr 1 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 CRITP Local 72 R(4) 4 scalar 72,107 CTLBLK Common 2 96 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 95,97 D50 Param 25 R(4) 4 scalar D608 Param 29 R(4) 4 scalar 112 DP10M Param 35 R(4) 4 scalar DT Scalar 4 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 ELDR Param 36 R(4) 4 scalar ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar Page 4 Source Listing CALMXW 2014-12-17 20:47 Symbol Table CALMXW.f Name Object Declared Type Bytes Dimen Elements Attributes References ERAD Param 34 R(4) 4 scalar EXTRA Common 3 71506268 FIRST Scalar 3 R(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM FMD Param 39 R(4) 4 scalar 39 FMW Param 39 R(4) 4 scalar 39 FSL Scalar 8 R(4) 4 2 91719 COM G Param 32 R(4) 4 scalar 32 GI Param 32 R(4) 4 scalar 36 H1 Param 26 R(4) 4 scalar 33,112 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 HTM Scalar 5 R(4) 4 3 4127355 COM 112 I Local 75 I(4) 4 scalar 75,76,77,78,79,88,89,90,92,95,96,9 7,98,107,108,109,111,112,113 IBESSL Scalar 2 I(4) 4 scalar COM 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 IE Local 89 I(4) 4 scalar 89,96,98 IFLAG Scalar 2 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 89 IHHA Scalar 3 I(4) 4 1 760 COM IHLA Scalar 3 I(4) 4 1 760 COM IHRST Scalar 3 I(4) 4 scalar COM IHW Scalar 3 I(4) 4 1 387 COM 90 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,14,64,75,88 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMT Param 4 I(4) 4 scalar INDX Common 2 6192 INPES Param 6 I(4) 4 scalar 11,14 IOFFS Scalar 2 I(4) 4 scalar COM 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 Page 5 Source Listing CALMXW 2014-12-17 20:47 Symbol Table CALMXW.f Name Object Declared Type Bytes Dimen Elements Attributes References IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IVW Scalar 3 I(4) 4 1 387 COM IW Local 90 I(4) 4 scalar 90,95,97 J Local 74 I(4) 4 scalar 74,76,77,78,79,87,89,90,92,95,96,9 7,98,107,108,109,111,112,113 JAM Param 6 I(4) 4 scalar 3 JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 74 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 87 JEND_M2 Scalar 7 I(4) 4 scalar COM JGSTL Param 10 I(4) 4 scalar 15 JGSTR Param 10 I(4) 4 scalar 15 JM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,12,15,64 JM2 Param 6 I(4) 4 scalar JMT Param 4 I(4) 4 scalar JNPES Param 6 I(4) 4 scalar 12,15 JRA Scalar 3 I(4) 4 1 760 COM JSTA Scalar 6 I(4) 4 scalar COM 74 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 87 JSTA_M2 Scalar 7 I(4) 4 scalar COM JTAIL Param 12 I(4) 4 scalar KMNT Scalar 5 I(4) 4 2 4118670 COM KMNTM Scalar 5 I(4) 4 1 45 COM KSAV Param 14 I(4) 4 scalar KSB Scalar 2 I(4) 4 scalar COM 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 L Local 94 I(4) 4 scalar 94,95,96,97,98,107,111,112,113 LB Param 9 I(4) 4 scalar LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LIST Scalar 3 I(4) 4 scalar COM LLMH Local 92 I(4) 4 scalar 92,94 LM Param 3 I(4) 4 scalar 5,8,19 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 92 LMV Scalar 5 I(4) 4 2 91719 COM LOG Func 113 scalar 113 LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4 LSM Param 3 I(4) 4 scalar MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM MXWP Dummy 1 R(4) 4 2 91719 ARG,INOUT 78,111,113 Page 6 Source Listing CALMXW 2014-12-17 20:47 Symbol Table CALMXW.f Name Object Declared Type Bytes Dimen Elements Attributes References MXWU Dummy 1 R(4) 4 2 91719 ARG,INOUT 76,108 MXWV Dummy 1 R(4) 4 2 91719 ARG,INOUT 77,109 MXWW Local 64 R(4) 4 scalar 91,107,110 MXWZ Dummy 1 R(4) 4 2 91719 ARG,INOUT 79,112 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 NINC Param 7 I(4) 4 scalar 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 NSTAT Param 19 I(4) 4 scalar NSUMD Param 15 I(4) 4 scalar 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 NWORD Param 19 I(4) 4 scalar 19 NWORD3 Param 19 I(4) 4 scalar NWORD6 Param 19 I(4) 4 scalar OLDRD Scalar 10 R(4) 4 scalar COM ONEPS Param 39 R(4) 4 scalar OPTIONS Common 2 24 OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar PARALLEL Common 9 8240 PD Scalar 3 R(4) 4 2 91719 COM PDSL Scalar 6 R(4) 4 2 91719 COM PDVP1 Scalar 6 R(4) 4 2 91719 COM PFSL Scalar 7 R(4) 4 2 91719 COM PINT Scalar 4 R(4) 4 3 4219074 COM 107,111,113 PQ0 Param 41 R(4) 4 scalar PSLP Scalar 6 R(4) 4 2 91719 COM PTSL Scalar 7 R(4) 4 2 91719 COM Q Scalar 5 R(4) 4 3 4127355 COM 112 QSL Scalar 9 R(4) 4 2 91719 COM RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDOCP Param 40 R(4) 4 scalar RES Scalar 3 R(4) 4 2 91719 COM RESTRT Scalar 3 R(4) 4 scalar COM RG Param 36 R(4) 4 scalar ROG Param 32 R(4) 4 scalar 112 RTD Param 34 R(4) 4 scalar RUN Scalar 3 R(4) 4 scalar COM SATDEL Scalar 2 R(4) 4 scalar COM SICE Scalar 3 R(4) 4 2 91719 COM SIGMA Scalar 3 R(4) 4 scalar COM SLP Scalar 5 R(4) 4 2 91719 COM SM Scalar 3 R(4) 4 2 91719 COM SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar Page 7 Source Listing CALMXW 2014-12-17 20:47 Symbol Table CALMXW.f Name Object Declared Type Bytes Dimen Elements Attributes References SPVAL Scalar 2 R(4) 4 scalar COM 76,77,78,79 SQRT Func 99 scalar 99 STDRD Scalar 10 R(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM 112 T500 Scalar 6 R(4) 4 2 91719 COM TFRZ Param 39 R(4) 4 scalar TSL Scalar 8 R(4) 4 2 91719 COM TTV Scalar 8 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM 95,96 U0 Local 95 R(4) 4 scalar 95,99,108 V Scalar 5 R(4) 4 3 4127355 COM 97,98 V0 Local 97 R(4) 4 scalar 97,99,109 VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VRBLS Common 2 67138308 VTM Scalar 5 R(4) 4 3 4127355 COM WIND Local 99 R(4) 4 scalar 99,107,110 Z1000 Scalar 9 R(4) 4 2 91719 COM ZINT Scalar 4 R(4) 4 3 4219074 COM 113 Page 8 Source Listing CALMXW 2014-12-17 20:47 Subprograms/Common Blocks CALMXW.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALMXW Subr 1 CTLBLK Common 2 96 EXTRA Common 3 71506268 INDX Common 2 6192 LOOPS Common 2 748952 MASKS Common 2 34853220 OPTIONS Common 2 24 PARALLEL Common 9 8240 VRBLS Common 2 67138308 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__ Page 9 Source Listing CALMXW 2014-12-17 20:47 CALMXW.f -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 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 : CALMXW.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100