Page 1 Source Listing CALHEL 2014-12-17 20:47 CALHEL.f 1 SUBROUTINE CALHEL(UST,VST,HELI) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CALHEL COMPUTES STORM RELATIVE HELICITY 5 C PRGRMMR: BALDWIN ORG: W/NP2 DATE: 94-08-22 6 C 7 C ABSTRACT: 8 C THIS ROUTINE COMPUTES ESTIMATED STORM MOTION AND 9 C STORM-RELATIVE ENVIRONMENTAL HELICITY. 10 C (DAVIES-JONES ET AL 1990) THE ALGORITHM PROCEEDS AS 11 C FOLLOWS. 12 C 13 C THE STORM MOTION COMPUTATION NO LONGER EMPLOYS THE DAVIES AND 14 C JOHNS (1993) METHOD WHICH DEFINED STORM MOTION AS 30 DEGREES TO 15 C THE RIGHT OF THE 0-6 KM MEAN WIND AT 75% OF THE SPEED FOR MEAN 16 C SPEEDS LESS THAN 15 M/S AND 20 DEGREES TO THE RIGHT FOR SPEEDS 17 C GREATER THAN 15 M/S. INSTEAD, WE NOW USE THE DYNAMIC METHOD 18 C (BUNKERS ET AL. 1998) WHICH HAS BEEN FOUND TO DO BETTER IN 19 C CASES WITH 'NON-CLASSIC' HODOGRAPHS (SUCH AS NORTHWEST-FLOW 20 C EVENTS) AND DO AS WELL OR BETTER THAN THE OLD METHOD IN MORE 21 C CLASSIC SITUATIONS. 22 C 23 C PROGRAM HISTORY LOG: 24 C 94-08-22 MICHAEL BALDWIN 25 C 97-03-27 MICHAEL BALDWIN - SPEED UP CODE 26 C 98-06-15 T BLACK - CONVERSION FROM 1-D TO 2-D 27 C 00-01-10 G MANIKIN - CHANGED CODE TO USE BUNKERS METHOD 28 C AS DESCRIBED ABOVE 29 C 30 C USAGE: CALHEL(UST,VST,HELI) 31 C INPUT ARGUMENT LIST: 32 C NONE 33 C 34 C OUTPUT ARGUMENT LIST: 35 C UST - ESTIMATED U COMPONENT (M/S) OF STORM MOTION. 36 C VST - ESTIMATED V COMPONENT (M/S) OF STORM MOTION. 37 C HELI - STORM-RELATIVE HELICITY (M**2/S**2) 38 C 39 C OUTPUT FILES: 40 C NONE 41 C 42 C SUBPROGRAMS CALLED: 43 C UTILITIES: 44 C 45 C LIBRARY: 46 C COMMON - VRBLS 47 C LOOPS 48 C PHYS 49 C EXTRA 50 C MASKS 51 C OPTIONS 52 C INDX 53 C 54 C ATTRIBUTES: 55 C LANGUAGE: FORTRAN 90 56 C MACHINE : CRAY C-90 57 C$$$ Page 2 Source Listing CALHEL 2014-12-17 20:47 CALHEL.f 58 C 59 C 60 C INCLUDE PARAMETERS. 61 INCLUDE "parmeta" 79 INCLUDE "params" 121 INCLUDE "parm.tbl" 123 PARAMETER (PI=3.141592654) 124 PARAMETER (P150=15000.0,P300=30000.0,S15=15.0) 125 PARAMETER (D3000=3000.0,PI6=0.5235987756,PI9=0.34906585) 126 PARAMETER (D5500=5500.0,D6000=6000.0,D7000=7000.0) 127 PARAMETER (D500=500.0) 128 C 129 C DECLARE VARIABLES 130 C 131 REAL UST(IM,JM),VST(IM,JM),HELI(IM,JM),HTSFC(IM,JM) 132 REAL UST6(IM,JM),VST6(IM,JM),ETOT6(IM,JM) 133 REAL UST5(IM,JM),VST5(IM,JM),ETOT5(IM,JM) 134 REAL UST1(IM,JM),VST1(IM,JM),ETOT1(IM,JM) 135 C 136 C INCLUDE COMMON BLOCKS. 137 INCLUDE "VRBLS.comm" 144 INCLUDE "LOOPS.comm" 152 INCLUDE "EXTRA.comm" 164 INCLUDE "MASKS.comm" 172 INCLUDE "PHYS.comm" 193 INCLUDE "OPTIONS.comm" 197 INCLUDE "INDX.comm" 203 INCLUDE "CTLBLK.comm" 204 C 205 C**************************************************************** 206 C START CALHEL HERE 207 C 208 C INITIALIZE ARRAYS. 209 C 210 !$omp parallel do 225 DO J=JSTA,JEND 226 DO I=1,IM 227 UST(I,J) = 0.0 228 VST(I,J) = 0.0 229 HELI(I,J) = 0.0 230 UST1(I,J) = 0.0 231 VST1(I,J) = 0.0 232 ETOT1(I,J) = 0.0 233 UST5(I,J) = 0.0 234 VST5(I,J) = 0.0 235 ETOT5(I,J) = 0.0 236 UST6(I,J) = 0.0 237 VST6(I,J) = 0.0 238 ETOT6(I,J) = 0.0 239 240 ENDDO 241 ENDDO 242 C 243 C LOOP OVER HORIZONTAL GRID. 244 C 245 CALL EXCH(RES) 246 CALL EXCH(PD) Page 3 Source Listing CALHEL 2014-12-17 20:47 CALHEL.f 247 DO L = 1,LP1 248 CALL EXCH(ZINT(1,1,L)) 249 END DO 250 C 251 !$omp parallel do 252 !$omp& private(htsfc,ie,iw,pdslvk,pkl,psfck) 253 DO L = 1,LM 254 DO J=JSTA_M,JEND_M 255 DO I=2,IM-1 256 IE=I+IVE(J) 257 IW=I+IVW(J) 258 PDSLVK=(PD(IW,J)*RES(IW,J)+PD(IE,J)*RES(IE,J)+ 259 1 PD(I,J+1)*RES(I,J+1)+PD(I,J-1)*RES(I,J-1))*0.25 260 PSFCK=AETA(LMV(I,J))*PDSLVK+PT 261 HTSFC(I,J)=0.25*(ZINT(IW,J,LMV(I,J)+1)+ZINT(IE,J,LMV(I,J)+1)+ 262 1 ZINT(I,J+1,LMV(I,J)+1)+ZINT(I,J-1,LMV(I,J)+1)) 263 C 264 C COMPUTE MASS WEIGHTED MEAN WIND IN THE 0-6 KM LAYER, THE 265 C 0-0.5 KM LAYER, AND THE 5.5-6 KM LAYER 266 C 267 Z2=0.125*(ZINT(IW,J,L)+ZINT(IW,J,L+1)+ 268 1 ZINT(IE,J,L)+ZINT(IE,J,L+1)+ 269 1 ZINT(I,J+1,L)+ZINT(I,J+1,L+1)+ 270 1 ZINT(I,J-1,L)+ZINT(I,J-1,L+1)) 271 DZABV=Z2-HTSFC(I,J) 272 273 IF (DZABV.LE.D6000 .AND. L.LE.LMV(I,J)) THEN 274 UST6(I,J) = UST6(I,J) + U(I,J,L) * DETA(L) 275 VST6(I,J) = VST6(I,J) + V(I,J,L) * DETA(L) 276 ETOT6(I,J) = ETOT6(I,J) + DETA(L) 277 ENDIF 278 279 IF (DZABV.LT.D6000 .AND. DZABV.GE.D5500 .AND. 280 & L.LE.LMV(I,J)) THEN 281 UST5(I,J) = UST5(I,J) + U(I,J,L) * DETA(L) 282 VST5(I,J) = VST5(I,J) + V(I,J,L) * DETA(L) 283 ETOT5(I,J) = ETOT5(I,J) + DETA(L) 284 ENDIF 285 286 IF (DZABV.LT.D500 .AND. L.LE.LMV(I,J)) THEN 287 UST1(I,J) = UST1(I,J) + U(I,J,L) * DETA(L) 288 VST1(I,J) = VST1(I,J) + V(I,J,L) * DETA(L) 289 ETOT1(I,J) = ETOT1(I,J) + DETA(L) 290 ENDIF 291 292 ENDDO 293 ENDDO 294 ENDDO 295 296 C CASE WHERE THERE IS NO LEVEL WITH HEIGHT BETWEEN 5500 AND 6000 297 DO J=JSTA_M,JEND_M 298 DO I=2,IM-1 299 IF (ETOT5(I,J) .EQ. 0) THEN 300 DO L=LM,1,-1 301 IE=I+IVE(J) 302 IW=I+IVW(J) 303 Z2=0.125*(ZINT(IW,J,L)+ZINT(IW,J,L+1)+ Page 4 Source Listing CALHEL 2014-12-17 20:47 CALHEL.f 304 1 ZINT(IE,J,L)+ZINT(IE,J,L+1)+ 305 1 ZINT(I,J+1,L)+ZINT(I,J+1,L+1)+ 306 1 ZINT(I,J-1,L)+ZINT(I,J-1,L+1)) 307 DZABV=Z2-HTSFC(I,J) 308 IF (DZABV.LT.D7000 .AND. DZABV.GE.D6000) THEN 309 UST5(I,J) = U(I,J,L) * DETA(L) 310 VST5(I,J) = V(I,J,L) * DETA(L) 311 ETOT5(I,J) = DETA(L) 312 GOTO 30 313 ENDIF 314 ENDDO 315 ENDIF 316 30 CONTINUE 317 ENDDO 318 ENDDO 319 320 321 !$omp parallel do 322 !$omp& private(umean6,vmean6,umean5,vmean5,umean1,vmean1,ushr,vshr) 323 DO J=JSTA_M,JEND_M 324 DO I=2,IM-1 325 IF (ETOT6(I,J).GT.0.0 .AND. ETOT1(I,J) .GT. 0.0 326 1 .AND. ETOT5(I,J) .GT. 0.0) THEN 327 UMEAN6 = UST6(I,J) / ETOT6(I,J) 328 VMEAN6 = VST6(I,J) / ETOT6(I,J) 329 UMEAN5 = UST5(I,J) / ETOT5(I,J) 330 VMEAN5 = VST5(I,J) / ETOT5(I,J) 331 UMEAN1 = UST1(I,J) / ETOT1(I,J) 332 VMEAN1 = VST1(I,J) / ETOT1(I,J) 333 334 C 335 C COMPUTE STORM MOTION VECTOR 336 C IT IS DEFINED AS 7.5 M/S TO THE RIGHT OF THE 0-6 KM MEAN 337 C WIND CONSTRAINED ALONG A LINE WHICH IS BOTH PERPENDICULAR 338 C TO THE 0-6 KM MEAN VERTICAL WIND SHEAR VECTOR AND PASSES 339 C THROUGH THE 0-6 KM MEAN WIND. THE WIND SHEAR VECTOR IS 340 C SET AS THE DIFFERENCE BETWEEN THE 5.5-6 KM WIND (THE HEAD 341 C OF THE SHEAR VECTOR) AND THE 0-0.5 KM WIND (THE TAIL). 342 C THIS IS FOR THE RIGHT-MOVING CASE; WE IGNORE THE LEFT MOVER. 343 344 USHR = UMEAN5 - UMEAN1 345 VSHR = VMEAN5 - VMEAN1 346 347 UST(I,J) = UMEAN6 + (7.5*VSHR/SQRT(USHR*USHR+VSHR*VSHR)) 348 VST(I,J) = VMEAN6 - (7.5*USHR/SQRT(USHR*USHR+VSHR*VSHR)) 349 ELSE 350 UST(I,J) = 0.0 351 VST(I,J) = 0.0 352 ENDIF 353 ENDDO 354 ENDDO 355 C 356 C 357 C COMPUTE STORM-RELATIVE HELICITY 358 C 359 !$omp parallel do 360 !$omp& private(du1,du2,dv1,dv2,dz,dz1,dz2,dzabv,ie,iw,z1,z2,z3) Page 5 Source Listing CALHEL 2014-12-17 20:47 CALHEL.f 361 DO L = 2,LM-1 362 DO J=JSTA_M,JEND_M 363 DO I=2,IM-1 364 IW=I+IVW(J) 365 IE=I+IVE(J) 366 Z2=0.125*(ZINT(IW,J,L)+ZINT(IW,J,L+1)+ 367 & ZINT(IE,J,L)+ZINT(IE,J,L+1)+ 368 & ZINT(I,J+1,L)+ZINT(I,J+1,L+1)+ 369 & ZINT(I,J-1,L)+ZINT(I,J-1,L+1)) 370 DZABV=Z2-HTSFC(I,J) 371 C 372 IF(DZABV.LT.D3000.AND.L.LE.LMV(I,J))THEN 373 Z1=0.125*(ZINT(IW,J,L+1)+ZINT(IW,J,L+2)+ 374 & ZINT(IE,J,L+1)+ZINT(IE,J,L+2)+ 375 & ZINT(I,J+1,L+1)+ZINT(I,J+1,L+2)+ 376 & ZINT(I,J-1,L+1)+ZINT(I,J-1,L+2)) 377 Z3=0.125*(ZINT(IW,J,L-1)+ZINT(IW,J,L)+ 378 & ZINT(IE,J,L-1)+ZINT(IE,J,L)+ 379 & ZINT(I,J+1,L-1)+ZINT(I,J+1,L)+ 380 & ZINT(I,J-1,L-1)+ZINT(I,J-1,L)) 381 DZ=0.25*((ZINT(IW,J,L)+ZINT(IE,J,L)+ 382 & ZINT(I,J-1,L)+ZINT(I,J+1,L))- 383 & (ZINT(IW,J,L+1)+ZINT(IE,J,L+1)+ 384 & ZINT(I,J-1,L+1)+ZINT(I,J+1,L+1))) 385 DZ1=Z1-Z2 386 DZ2=Z2-Z3 387 DU1=U(I,J,L+1)-U(I,J,L) 388 DU2=U(I,J,L)-U(I,J,L-1) 389 DV1=V(I,J,L+1)-V(I,J,L) 390 DV2=V(I,J,L)-V(I,J,L-1) 391 HELI(I,J)=((V(I,J,L)-VST(I,J))* 392 1 (DZ2*(DU1/DZ1)+DZ1*(DU2/DZ2)) 393 2 -(U(I,J,L)-UST(I,J))* 394 3 (DZ2*(DV1/DZ1)+DZ1*(DV2/DZ2))) 395 4 *DZ/(DZ1+DZ2)+HELI(I,J) 396 ENDIF 397 ENDDO 398 ENDDO 399 ENDDO 400 C 401 C END OF ROUTINE. 402 C 403 RETURN 404 END Page 6 Source Listing CALHEL 2014-12-17 20:47 Entry Points CALHEL.f ENTRY POINTS Name calhel_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 30 Label 184 180 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 AETA Scalar 7 R(4) 4 1 45 COM 128 ALBASE Scalar 14 R(4) 4 2 91719 COM ALBEDO Scalar 15 R(4) 4 2 91719 COM ALPINT Scalar 4 R(4) 4 3 4219074 COM BETA Param 37 R(4) 4 scalar 37 BTG Param 37 R(4) 4 scalar CALHEL Subr 1 CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CI Scalar 4 R(4) 4 scalar COM 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 CS Scalar 4 R(4) 4 scalar COM CTLBLK Common 2 96 CZEN Scalar 12 R(4) 4 2 91719 COM CZMEAN Scalar 16 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 D3000 Param 66 R(4) 4 scalar 240 D50 Param 25 R(4) 4 scalar D500 Param 68 R(4) 4 scalar 154 D5500 Param 67 R(4) 4 scalar 147 D6000 Param 67 R(4) 4 scalar 141,147,176 D608 Param 29 R(4) 4 scalar D7000 Param 67 R(4) 4 scalar 176 DETA Scalar 7 R(4) 4 1 45 COM 142,143,144,149,150,151,155,156,15 7,177,178,179 DFRLG Scalar 7 R(4) 4 1 46 COM DI Scalar 4 R(4) 4 scalar COM Page 7 Source Listing CALHEL 2014-12-17 20:47 Symbol Table CALHEL.f Name Object Declared Type Bytes Dimen Elements Attributes References DP10M Param 35 R(4) 4 scalar DS Scalar 4 R(4) 4 scalar COM DT Scalar 4 R(4) 4 scalar COM DTD Scalar 3 R(4) 4 scalar COM DTQ2 Scalar 3 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 DU1 Local 255 R(4) 4 scalar 255,260 DU2 Local 256 R(4) 4 scalar 256,260 DV1 Local 257 R(4) 4 scalar 257,262 DV2 Local 258 R(4) 4 scalar 258,262 DZ Local 249 R(4) 4 scalar 249,263 DZ1 Local 253 R(4) 4 scalar 253,260,262,263 DZ2 Local 254 R(4) 4 scalar 254,260,262,263 DZABV Local 139 R(4) 4 scalar 139,141,147,154,175,176,238,240 ELDR Param 36 R(4) 4 scalar ELWDWN Scalar 13 R(4) 4 2 91719 COM ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar EPSR Scalar 11 R(4) 4 2 91719 COM ERAD Param 34 R(4) 4 scalar ETOT1 Local 75 R(4) 4 2 91719 100,157,193,199,200 ETOT5 Local 74 R(4) 4 2 91719 103,151,167,179,194,197,198 ETOT6 Local 73 R(4) 4 2 91719 106,144,193,195,196 EXCH Subr 113 113,114,116 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 GFFC Scalar 14 R(4) 4 2 91719 COM GI Param 32 R(4) 4 scalar 36 GLAT Scalar 12 R(4) 4 2 91719 COM GLON Scalar 12 R(4) 4 2 91719 COM 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 HBOT Scalar 13 R(4) 4 2 91719 COM HDAC Scalar 16 R(4) 4 2 91719 COM HDACV Scalar 16 R(4) 4 2 91719 COM HELI Dummy 1 R(4) 4 2 91719 ARG,INOUT 97,259,263 Page 8 Source Listing CALHEL 2014-12-17 20:47 Symbol Table CALHEL.f Name Object Declared Type Bytes Dimen Elements Attributes References HTM Scalar 5 R(4) 4 3 4127355 COM HTOP Scalar 13 R(4) 4 2 91719 COM HTSFC Local 72 R(4) 4 2 91719 129,139,175,238 I Local 94 I(4) 4 scalar 94,95,96,97,98,99,100,101,102,103, 104,105,106,123,124,125,127,128,12 9,130,137,138,139,141,142,143,144, 148,149,150,151,154,155,156,157,16 6,167,169,170,173,174,175,177,178, 179,192,193,194,195,196,197,198,19 9,200,215,216,218,219,231,232,233, 236,237,238,240,243,244,247,248,25 0,252,255,256,257,258,259,261,263 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 124 I(4) 4 scalar 124,126,129,136,169,172,233,235,24 2,246,249,251 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 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 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,12,13,14,15,16, 72,73,74,75,94,123,166,192,231 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 ITB Param 1 I(4) 4 scalar 8,18 ITBQ Param 1 I(4) 4 scalar 9,18 IUP Scalar 7 I(4) 4 scalar COM IVE Scalar 3 I(4) 4 1 387 COM 124,169,233 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 125,170,232 IW Local 125 I(4) 4 scalar 125,126,129,135,170,171,232,234,24 1,245,249,251 J Local 93 I(4) 4 scalar 93,95,96,97,98,99,100,101,102,103, 104,105,106,122,124,125,126,127,12 8,129,130,135,136,137,138,139,141, 142,143,144,148,149,150,151,154,15 5,156,157,165,167,169,170,171,172, 173,174,175,177,178,179,191,193,19 Page 9 Source Listing CALHEL 2014-12-17 20:47 Symbol Table CALHEL.f Name Object Declared Type Bytes Dimen Elements Attributes References 4,195,196,197,198,199,200,215,216, 218,219,230,232,233,234,235,236,23 7,238,240,241,242,243,244,245,246, 247,248,249,250,251,252,255,256,25 7,258,259,261,263 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 93 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 122,165,191,230 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,11,12,13,14,15,16, 72,73,74,75 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 93 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 122,165,191,230 JSTA_M2 Scalar 7 I(4) 4 scalar COM JTAIL Param 12 I(4) 4 scalar JTB Param 1 I(4) 4 scalar 8,18 JTBQ Param 1 I(4) 4 scalar 18 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 KTM Scalar 3 I(4) 4 scalar COM L Local 115 I(4) 4 scalar 115,116,121,135,136,137,138,141,14 2,143,144,148,149,150,151,154,155, 156,157,168,171,172,173,174,177,17 8,179,229,234,235,236,237,240,241, 242,243,244,245,246,247,248,249,25 0,251,252,255,256,257,258,259,261 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 LM Param 3 I(4) 4 scalar 5,7,8,19,121,168,229 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM LMV Scalar 5 I(4) 4 2 91719 COM 128,129,130,141,148,154,240 Page 10 Source Listing CALHEL 2014-12-17 20:47 Symbol Table CALHEL.f Name Object Declared Type Bytes Dimen Elements Attributes References LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4,7,115 LSM Param 3 I(4) 4 scalar MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM 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 P150 Param 65 R(4) 4 scalar P300 Param 65 R(4) 4 scalar PARALLEL Common 9 8240 PD Scalar 3 R(4) 4 2 91719 COM 114,126,127 PDSL Scalar 6 R(4) 4 2 91719 COM PDSLVK Local 126 R(4) 4 scalar 126,128 PDVP1 Scalar 6 R(4) 4 2 91719 COM PFSL Scalar 7 R(4) 4 2 91719 COM PHYS Common 2 7690036 PI Param 64 R(4) 4 scalar PI6 Param 66 R(4) 4 scalar PI9 Param 66 R(4) 4 scalar PINT Scalar 4 R(4) 4 3 4219074 COM PL Scalar 5 R(4) 4 scalar COM PLQ Scalar 6 R(4) 4 scalar COM PQ0 Param 41 R(4) 4 scalar PSFCK Local 128 R(4) 4 scalar 128 PSLP Scalar 6 R(4) 4 2 91719 COM PT Scalar 3 R(4) 4 scalar COM 128 PTBL Scalar 18 R(4) 4 2 10184 COM PTSL Scalar 7 R(4) 4 2 91719 COM Q Scalar 5 R(4) 4 3 4127355 COM QS0 Scalar 8 R(4) 4 1 134 COM QSL Scalar 9 R(4) 4 2 91719 COM RADIN Scalar 11 R(4) 4 2 91719 COM RADOT Scalar 11 R(4) 4 2 91719 COM Page 11 Source Listing CALHEL 2014-12-17 20:47 Symbol Table CALHEL.f Name Object Declared Type Bytes Dimen Elements Attributes References RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDOCP Param 40 R(4) 4 scalar RDP Scalar 5 R(4) 4 scalar COM RDPQ Scalar 6 R(4) 4 scalar COM RDQ Scalar 5 R(4) 4 scalar COM RDTH Scalar 5 R(4) 4 scalar COM RDTHE Scalar 5 R(4) 4 scalar COM RDTHEQ Scalar 6 R(4) 4 scalar COM RES Scalar 3 R(4) 4 2 91719 COM 113,126,127 RESTRT Scalar 3 R(4) 4 scalar COM RG Param 36 R(4) 4 scalar RMXSNAL Scalar 11 R(4) 4 2 91719 COM ROG Param 32 R(4) 4 scalar ROI Scalar 4 R(4) 4 scalar COM ROS Scalar 4 R(4) 4 scalar COM RTD Param 34 R(4) 4 scalar RUN Scalar 3 R(4) 4 scalar COM S15 Param 65 R(4) 4 scalar 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 SIGT4 Scalar 16 R(4) 4 2 91719 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 SPVAL Scalar 2 R(4) 4 scalar COM SQRT Func 215 scalar 215,216 SQS Scalar 8 R(4) 4 1 134 COM SST Scalar 14 R(4) 4 2 91719 COM STDRD Scalar 10 R(4) 4 scalar COM STHE Scalar 8 R(4) 4 1 76 COM STHEQ Scalar 9 R(4) 4 1 152 COM SWNET Scalar 12 R(4) 4 2 91719 COM T Scalar 5 R(4) 4 3 4127355 COM T500 Scalar 6 R(4) 4 2 91719 COM TDTD Scalar 3 R(4) 4 scalar COM TDTQ2 Scalar 3 R(4) 4 scalar COM TFRZ Param 39 R(4) 4 scalar TG Scalar 14 R(4) 4 2 91719 COM THE0 Scalar 8 R(4) 4 1 76 COM THE0Q Scalar 9 R(4) 4 1 152 COM THL Scalar 5 R(4) 4 scalar COM TSL Scalar 8 R(4) 4 2 91719 COM TTBL Scalar 18 R(4) 4 2 10184 COM TTBLQ Scalar 18 R(4) 4 2 66880 COM TTV Scalar 8 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM 142,149,155,177,255,256,261 UMEAN1 Local 199 R(4) 4 scalar 199,212 UMEAN5 Local 197 R(4) 4 scalar 197,212 UMEAN6 Local 195 R(4) 4 scalar 195,215 USHR Local 212 R(4) 4 scalar 212,215,216 UST Dummy 1 R(4) 4 2 91719 ARG,INOUT 95,215,218,261 UST1 Local 75 R(4) 4 2 91719 98,155,199 Page 12 Source Listing CALHEL 2014-12-17 20:47 Symbol Table CALHEL.f Name Object Declared Type Bytes Dimen Elements Attributes References UST5 Local 74 R(4) 4 2 91719 101,149,177,197 UST6 Local 73 R(4) 4 2 91719 104,142,195 V Scalar 5 R(4) 4 3 4127355 COM 143,150,156,178,257,258,259 VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VMEAN1 Local 200 R(4) 4 scalar 200,213 VMEAN5 Local 198 R(4) 4 scalar 198,213 VMEAN6 Local 196 R(4) 4 scalar 196,216 VRBLS Common 2 67138308 VSHR Local 213 R(4) 4 scalar 213,215,216 VST Dummy 1 R(4) 4 2 91719 ARG,INOUT 96,216,219,259 VST1 Local 75 R(4) 4 2 91719 99,156,200 VST5 Local 74 R(4) 4 2 91719 102,150,178,198 VST6 Local 73 R(4) 4 2 91719 105,143,196 VTM Scalar 5 R(4) 4 3 4127355 COM Z1 Local 241 R(4) 4 scalar 241,253 Z1000 Scalar 9 R(4) 4 2 91719 COM Z2 Local 135 R(4) 4 scalar 135,139,171,175,234,238,253,254 Z3 Local 245 R(4) 4 scalar 245,254 ZINT Scalar 4 R(4) 4 3 4219074 COM 116,129,130,135,136,137,138,171,17 2,173,174,234,235,236,237,241,242, 243,244,245,246,247,248,249,250,25 1,252 Page 13 Source Listing CALHEL 2014-12-17 20:47 Subprograms/Common Blocks CALHEL.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALHEL 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 PHYS Common 2 7690036 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__ Page 14 Source Listing CALHEL 2014-12-17 20:47 CALHEL.f -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 -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 : CALHEL.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100