Page 1 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.f 1 SUBROUTINE CALCAPE(ITYPE,P1D,T1D,Q1D,L1D,CAPE,CINS) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CALCAPE COMPUTES CAPE AND CINS 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-02-10 6 C 7 C ABSTRACT: 8 C 9 C THIS ROUTINE COMPUTES CAPE AND CINS GIVEN TEMPERATURE, 10 C PRESSURE, AND SPECIFIC HUMIDTY. IN "STORM AND CLOUD 11 C DYNAMICS" (1989, ACADEMIC PRESS) COTTON AND ANTHES DEFINE 12 C CAPE (EQUATION 9.16, P501) AS 13 C 14 C EL 15 C CAPE = SUM G * LN(THETAP/THETAA) DZ 16 C LCL 17 C 18 C WHERE, 19 C EL = EQUILIBRIUM LEVEL, 20 C LCL = LIFTING CONDENSTATION LEVEL, 21 C G = GRAVITATIONAL ACCELERATION, 22 C THETAP = LIFTED PARCEL POTENTIAL TEMPERATURE, 23 C THETAA = AMBIENT POTENTIAL TEMPERATURE. 24 C 25 C NOTE THAT THE INTEGRAND LN(THETAP/THETAA) APPROXIMATELY 26 C EQUALS (THETAP-THETAA)/THETAA. THIS RATIO IS OFTEN USED 27 C IN THE DEFINITION OF CAPE/CINS. 28 C 29 C TWO TYPES OF CAPE/CINS CAN BE COMPUTED BY THIS ROUTINE. THE 30 C SUMMATION PROCESS IS THE SAME FOR BOTH CASES. WHAT DIFFERS 31 C IS THE DEFINITION OF THE PARCEL TO LIFT. FOR ITYPE=1 THE 32 C PARCEL WITH THE WARMEST THETA-E IN A DPBND PASCAL LAYER ABOVE 33 C THE MODEL SURFACE IS LIFTED. THE ARRAYS P1D, T1D, Q1D, AND 34 C L1D ARE NOT USED. FOR ITYPE=2 THE ARRAYS P1D, T1D, Q1D, AND 35 C L1D DEFINE THE PARCEL TO LIFT IN EACH COLUMN. BOTH TYPES OF 36 C CAPE/CINS MAY BE COMPUTED IN A SINGLE EXECUTION OF THE POST 37 C PROCESSOR. 38 C 39 C THIS ALGORITHM PROCEEDS AS FOLLOWS. 40 C FOR EACH COLUMN, 41 C (1) INITIALIZE RUNNING CAPE AND CINS SUM TO 0.0 42 C (2) COMPUTE TEMPERATURE AND PRESSURE AT THE LCL USING 43 C LOOK UP TABLE (PTBL). USE EITHER PARCEL THAT GIVES 44 C MAX THETAE IN LOWEST DPBND ABOVE GROUND (ITYPE=1) 45 C OR GIVEN PARCEL FROM T1D,Q1D,...(ITYPE=2). 46 C (3) COMPUTE THE TEMP OF A PARCEL LIFTED FROM THE LCL. 47 C WE KNOW THAT THE PARCEL'S 48 C EQUIVALENT POTENTIAL TEMPERATURE (THESP) REMAINS 49 C CONSTANT THROUGH THIS PROCESS. WE CAN 50 C COMPUTE TPAR USING THIS KNOWLEDGE USING LOOK 51 C UP TABLE (SUBROUTINE TTBLEX). 52 C (4) FIND THE EQUILIBRIUM LEVEL. THIS IS DEFINED AS THE 53 C HIGHEST POSITIVELY BUOYANT LAYER. 54 C (IF THERE IS NO POSITIVELY BUOYANT LAYER, CAPE/CINS 55 C WILL BE ZERO) 56 C (5) COMPUTE CAPE/CINS. 57 C (A) COMPUTE THETAP. WE KNOW TPAR AND P. Page 2 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.f 58 C (B) COMPUTE THETAA. WE KNOW T AND P. 59 C (6) ADD G*(THETAP-THETAA)*DZ TO THE RUNNING CAPE OR CINS SUM. 60 C (A) IF THETAP > THETAA, ADD TO THE CAPE SUM. 61 C (B) IF THETAP < THETAA, ADD TO THE CINS SUM. 62 C (7) ARE WE AT EQUILIBRIUM LEVEL? 63 C (A) IF YES, STOP THE SUMMATION. 64 C (B) IF NO, CONTIUNUE THE SUMMATION. 65 C (8) ENFORCE LIMITS ON CAPE AND CINS (I.E. NO NEGATIVE CAPE) 66 C 67 C PROGRAM HISTORY LOG: 68 C 93-02-10 RUSS TREADON 69 C 93-06-19 RUSS TREADON - GENERALIZED ROUTINE TO ALLOW FOR 70 C TYPE 2 CAPE/CINS CALCULATIONS. 71 C 94-09-23 MIKE BALDWIN - MODIFIED TO USE LOOK UP TABLES 72 C INSTEAD OF COMPLICATED EQUATIONS. 73 C 94-10-13 MIKE BALDWIN - MODIFIED TO CONTINUE CAPE/CINS CALC 74 C UP TO AT HIGHEST BUOYANT LAYER. 75 C 98-06-12 T BLACK - CONVERSION FROM 1-D TO 2-D 76 C 98-08-18 T BLACK - COMPUTE APE INTERNALLY 77 C 00-01-04 JIM TUCCILLO - MPI VERSION 78 C 79 C USAGE: CALL CALCAPE(ITYPE,P1D,T1D,Q1D,L1D,CAPE,CINS) 80 C INPUT ARGUMENT LIST: 81 C ITYPE - INTEGER FLAG SPECIFYING HOW PARCEL TO LIFT IS 82 C IDENTIFIED. SEE COMMENTS ABOVE. 83 C P1D - ARRAY OF PRESSURE OF PARCELS TO LIFT. 84 C T1D - ARRAY OF TEMPERATURE OF PARCELS TO LIFT. 85 C Q1D - ARRAY OF SPECIFIC HUMIDITY OF PARCELS TO LIFT. 86 C L1D - ARRAY OF MODEL LEVEL OF PARCELS TO LIFT. 87 C 88 C OUTPUT ARGUMENT LIST: 89 C CAPE - CONVECTIVE AVAILABLE POTENTIAL ENERGY (J/KG) 90 C CINS - CONVECTIVE INHIBITION (J/KG) 91 C 92 C OUTPUT FILES: 93 C STDOUT - RUN TIME STANDARD OUT. 94 C 95 C SUBPROGRAMS CALLED: 96 C UTILITIES: 97 C P2FILT - SMOOTH DATA ON STAGGERED E-GRID. 98 C BOUND - BOUND (CLIP) DATA BETWEEN UPPER AND LOWER LIMTS. 99 C TTBLEX - LOOKUP TABLE ROUTINE TO GET T FROM THETAE AND P 100 C 101 C LIBRARY: 102 C COMMON - VRBLS 103 C PHYS 104 C EXTRA 105 C IOUNIT 106 C LOOPS 107 C MASKS 108 C 109 C ATTRIBUTES: 110 C LANGUAGE: FORTRAN 90 111 C MACHINE : CRAY C-90 112 C$$$ 113 C 114 C Page 3 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.f 115 C 116 C INCLUDE/SET PARAMETERS. CONSTANTS ARE FROM BOLTON (MWR, 1980). 117 INCLUDE "parmeta" 135 INCLUDE "parm.tbl" 137 INCLUDE "params" 179 PARAMETER (DPBND=70.E2, 180 & ELIVW=2.72E6,ELOCP=ELIVW/CP) 181 PARAMETER (ISMTHP=2,ISMTHT=2,ISMTHQ=2) 182 C 183 C DECLARE VARIABLES. 184 C 185 INTEGER L1D(IM,JM),IEQL(IM,JM),IPTB(IM,JM),ITHTB(IM,JM) 186 INTEGER KLRES(IM,JM),KHRES(IM,JM),LCL(IM,JM) 187 C INTEGER KNUML(LM),KNUMH(LM) 188 C 189 REAL P1D(IM,JM),T1D(IM,JM),Q1D(IM,JM) 190 REAL CAPE(IM,JM),CINS(IM,JM) 191 REAL THESP(IM,JM),TPAR(IM,JM,LM),PSP(IM,JM) 192 REAL QQ(IM,JM),PP(IM,JM) 193 INTEGER IDX(IM,JM) 194 C 195 C INCLUDE COMMON BLOCKS. 196 INCLUDE "VRBLS.comm" 203 INCLUDE "EXTRA.comm" 215 INCLUDE "IOUNIT.comm" 223 INCLUDE "LOOPS.comm" 231 INCLUDE "PHYS.comm" 252 INCLUDE "MASKS.comm" 260 INCLUDE "CTLBLK.comm" 261 C 262 C************************************************************** 263 C START CALCAPE HERE. 264 C 265 C 266 C COMPUTE CAPE/CINS 267 C 268 C WHICH IS: THE SUM FROM THE LCL TO THE EQ LEVEL OF 269 C G * (LN(THETAP) - LN(THETAA)) * DZ 270 C 271 C (POSITIVE AREA FOR CAPE, NEGATIVE FOR CINS) 272 C 273 C WHERE: 274 C THETAP IS THE PARCEL THETA 275 C THETAA IS THE AMBIENT THETA 276 C DZ IS THE THICKNESS OF THE LAYER 277 C 278 C USING LCL AS LEVEL DIRECTLY BELOW SATURATION POINT 279 C AND EQ LEVEL IS THE HIGHEST POSITIVELY BUOYANT LEVEL. 280 C 281 C IEQL = EQ LEVEL 282 C 283 C INITIALIZE CAPE AND CINS ARRAYS 284 C 285 !$omp parallel do 300 DO J=JSTA,JEND 301 DO I=1,IM 302 CAPE(I,J) = D00 Page 4 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.f 303 CINS(I,J) = D00 304 LCL(I,J) = D00 305 THESP(I,J)= D00 306 IEQL(I,J) = LM+1 307 ENDDO 308 ENDDO 309 C 310 !$omp parallel do 311 DO L=1,LM 312 DO J=JSTA,JEND 313 DO I=1,IM 314 TPAR(I,J,L)= D00 315 ENDDO 316 ENDDO 317 ENDDO 318 C 319 C FOR TYPE 2 CAPE/CINS SMOOTH DATA PRIOR TO CALCULATION. 320 C NOTE THAT FOR TYPE 1 CAPE/CINS ARRAYS P1D, T1D, Q1D 321 C ARE DUMMY ARRAYS. 322 C 323 IF (ITYPE.EQ.2) THEN 324 CALL P2FILT(ISMTHP,HBM2,P1D) 325 CALL P2FILT(ISMTHT,HBM2,T1D) 326 CALL P2FILT(ISMTHQ,HBM2,Q1D) 327 CALL BOUNDL(Q1D,H1M12,H99999,IM,JM) 328 ENDIF 329 C-------FOR ITYPE=1----------------------------------------------------- 330 C---------FIND MAXIMUM THETA E LAYER IN LOWEST DPBND ABOVE GROUND------- 331 C-------FOR ITYPE=2----------------------------------------------------- 332 C---------FIND THETA E LAYER OF GIVEN T1D, Q1D, P1D--------------------- 333 C--------------TRIAL MAXIMUM BUOYANCY LEVEL VARIABLES------------------- 334 DO 20 KB=1,LM 335 IF (ITYPE.EQ.2.AND.KB.GT.1) GOTO 20 336 !$omp parallel do 337 !$omp& private(apebtk,apespk,bqk,bqs00k,bqs10k,iq,it,ittbk,lmhk, 338 !$omp& p00k,p01k,p10k,p11k,pkl,psfck,qbtk,sqk,sqs00k, 339 !$omp& sqs10k,tbtk,tpspk,tqk,tthbtk,tthesk,tthk) 340 DO 10 J=JSTA_M,JEND_M 341 DO 10 I=2,IM-1 342 LMHK =LMH(I,J) 343 PSFCK =AETA(LMHK)*PDSL(I,J)+PT 344 PKL = AETA(KB)*PDSL(I,J)+PT 345 IF (ITYPE.EQ.1.AND.(PKL.LT.PSFCK-DPBND.OR.PKL.GT.PSFCK)) 346 & GOTO 10 347 IF (ITYPE.EQ.1) THEN 348 TBTK =T(I,J,KB) 349 QBTK =Q(I,J,KB) 350 APEBTK =(H10E5/(PDSL(I,J)*AETA(KB)+PT))**CAPA 351 ELSE 352 PKL =P1D(I,J) 353 TBTK =T1D(I,J) 354 QBTK =Q1D(I,J) 355 APEBTK =(H10E5/PKL)**CAPA 356 ENDIF 357 C--------------SCALING POTENTIAL TEMPERATURE & TABLE INDEX-------------- 358 TTHBTK =TBTK*APEBTK 359 TTHK =(TTHBTK-THL)*RDTH Page 5 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.f 360 QQ(I,J)=TTHK-AINT(TTHK) 361 ITTBK =INT(TTHK)+1 362 C--------------KEEPING INDICES WITHIN THE TABLE------------------------- 363 IF(ITTBK.LT.1) THEN 364 ITTBK =1 365 QQ(I,J)=D00 366 ENDIF 367 IF(ITTBK.GE.JTB) THEN 368 ITTBK =JTB-1 369 QQ(I,J)=D00 370 ENDIF 371 C--------------BASE AND SCALING FACTOR FOR SPEC. HUMIDITY--------------- 372 BQS00K=QS0(ITTBK) 373 SQS00K=SQS(ITTBK) 374 BQS10K=QS0(ITTBK+1) 375 SQS10K=SQS(ITTBK+1) 376 C--------------SCALING SPEC. HUMIDITY & TABLE INDEX--------------------- 377 BQK =(BQS10K-BQS00K)*QQ(I,J)+BQS00K 378 SQK =(SQS10K-SQS00K)*QQ(I,J)+SQS00K 379 TQK =(QBTK-BQK)/SQK*RDQ 380 PP(I,J)=TQK-AINT(TQK) 381 IQ =INT(TQK)+1 382 C--------------KEEPING INDICES WITHIN THE TABLE------------------------- 383 IF(IQ.LT.1) THEN 384 IQ =1 385 PP(I,J)=D00 386 ENDIF 387 IF(IQ.GE.ITB) THEN 388 IQ =ITB-1 389 PP(I,J)=D00 390 ENDIF 391 C--------------SATURATION PRESSURE AT FOUR SURROUNDING TABLE PTS.------- 392 IT=ITTBK 393 P00K=PTBL(IQ ,IT ) 394 P10K=PTBL(IQ+1,IT ) 395 P01K=PTBL(IQ ,IT+1) 396 P11K=PTBL(IQ+1,IT+1) 397 C--------------SATURATION POINT VARIABLES AT THE BOTTOM----------------- 398 TPSPK=P00K+(P10K-P00K)*PP(I,J)+(P01K-P00K)*QQ(I,J) 399 2 +(P00K-P10K-P01K+P11K)*PP(I,J)*QQ(I,J) 400 APESPK=(H10E5/TPSPK)**CAPA 401 TTHESK=TTHBTK*EXP(ELOCP*QBTK*APESPK/TTHBTK) 402 C--------------CHECK FOR MAXIMUM THETA E-------------------------------- 403 IF(TTHESK.GT.THESP(I,J)) THEN 404 PSP (I,J)=TPSPK 405 THESP(I,J)=TTHESK 406 ENDIF 407 10 CONTINUE 408 20 CONTINUE 409 C-----CHOOSE LAYER DIRECTLY BELOW PSP AS LCL AND------------------------ 410 C-----ENSURE THAT THE LCL IS ABOVE GROUND.------------------------------ 411 C-------(IN SOME RARE CASES FOR ITYPE=2, IT IS NOT)--------------------- 412 DO L=1,LM 413 !$omp parallel do 414 !$omp& private(pkl) 415 DO J=JSTA_M,JEND_M 416 DO I=2,IM-1 Page 6 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.f 417 PKL = AETA(L)*PDSL(I,J)+PT 418 IF (PKL.LT.PSP(I,J)) LCL(I,J)=L+1 419 ENDDO 420 ENDDO 421 ENDDO 422 !$omp parallel do 423 DO J=JSTA_M,JEND_M 424 DO I=2,IM-1 425 IF (LCL(I,J).GT.LMH(I,J)) LCL(I,J)=LMH(I,J) 426 ENDDO 427 ENDDO 428 C----------------------------------------------------------------------- 429 C---------FIND TEMP OF PARCEL LIFTED ALONG MOIST ADIABAT (TPAR)--------- 430 C----------------------------------------------------------------------- 431 !!$omp parallel do 432 DO 30 L=LM,1,-1 433 C--------------SCALING PRESSURE & TT TABLE INDEX------------------------ 434 KNUML=0 435 KNUMH=0 436 DO J=JSTA_M,JEND_M 437 DO I=2,IM-1 438 KLRES(I,J)=0 439 KHRES(I,J)=0 440 PKL=AETA(L)*PDSL(I,J)+PT 441 IF(L.LE.LCL(I,J)) THEN 442 IF(PKL.LT.PLQ)THEN 443 KNUML=KNUML+1 444 KLRES(I,J)=1 445 ELSE 446 KNUMH=KNUMH+1 447 KHRES(I,J)=1 448 ENDIF 449 ENDIF 450 ENDDO 451 ENDDO 452 C*** 453 C*** COMPUTE PARCEL TEMPERATURE ALONG MOIST ADIABAT FOR PRESSUREPLQ 462 C** 463 IF(KNUMH.GT.0)THEN 464 CALL TTBLEX(TPAR(1,1,L),TTBLQ,ITBQ,JTBQ,KHRES,PDSL,AETA(L) 465 1, HTM(1,1,L),PT,PLQ,QQ,PP,RDPQ,THE0Q,STHEQ,RDTHEQ 466 2, THESP,IPTB,ITHTB) 467 ENDIF 468 469 C------------SEARCH FOR EQ LEVEL---------------------------------------- 470 DO J=JSTA_M,JEND_M 471 DO I=2,IM-1 472 IF(KHRES(I,J).GT.0)THEN 473 IF(TPAR(I,J,L).GT.T(I,J,L)) IEQL(I,J)=L Page 7 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.f 474 ENDIF 475 ENDDO 476 ENDDO 477 C 478 DO J=JSTA_M,JEND_M 479 DO I=2,IM-1 480 IF(KLRES(I,J).GT.0)THEN 481 IF(TPAR(I,J,L).GT.T(I,J,L)) IEQL(I,J)=L 482 ENDIF 483 ENDDO 484 ENDDO 485 C----------------------------------------------------------------------- 486 30 CONTINUE 487 C------------COMPUTE CAPE AND CINS-------------------------------------- 488 c DO K=KHL01,KHH01 489 c LCLK=LCL(K) 490 c IEQK=IEQL(K) 491 cCDIR$ SHORTLOOP 492 c DO L=IEQK,LCLK 493 c PRESK=AETA(L)*PDSL(K)+PT 494 c DZKL=T(K,L)*(Q(K,L)*D608+H1)*ROG*PDSL(K)*DETA(L)/PRESK 495 c THETAP=TPAR(K,L)*(H10E5/PRESK)**CAPA 496 c THETAA=T(K,L)*(H10E5/PRESK)**CAPA 497 c IF (THETAP.LT.THETAA) 498 c & CINS(K)=CINS(K)+G*(ALOG(THETAP)-ALOG(THETAA))*DZKL 499 c IF (THETAP.GT.THETAA) 500 c & CAPE(K)=CAPE(K)+G*(ALOG(THETAP)-ALOG(THETAA))*DZKL 501 c ENDDO 502 c ENDDO 503 LBEG=100 504 LEND=0 505 C 506 !$omp parallel do 507 !$omp& private(lbeg,lend) 508 DO J=JSTA_M,JEND_M 509 DO I=2,IM-1 510 LBEG=MIN(IEQL(I,J),LBEG) 511 LEND=MAX(LCL(I,J),LEND) 512 ENDDO 513 ENDDO 514 C 515 DO L=LBEG,LEND 516 DO J=JSTA_M,JEND_M 517 DO I=2,IM-1 518 IDX(I,J)=0 519 IF(L.GE.IEQL(I,J).AND.L.LE.LCL(I,J))THEN 520 IDX(I,J)=1 521 ENDIF 522 ENDDO 523 C 524 ENDDO 525 C 526 !$omp parallel do 527 !$omp& private(dzkl,presk,thetaa,thetap) 528 DO J=JSTA_M,JEND_M 529 DO I=2,IM-1 530 IF(IDX(I,J).GT.0)THEN Page 8 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.f 531 PRESK=AETA(L)*PDSL(I,J)+PT 532 DZKL=T(I,J,L)*(Q(I,J,L)*D608+H1) 533 * *ROG*PDSL(I,J)*DETA(L)/PRESK 534 THETAP=TPAR(I,J,L)*(H10E5/PRESK)**CAPA 535 THETAA=T(I,J,L)*(H10E5/PRESK)**CAPA 536 IF(THETAP.LT.THETAA)THEN 537 CINS(I,J)=CINS(I,J) 538 cjjt * +G*(ALOG(THETAP/THETAA))*DZKL 539 * +G*(ALOG(THETAP)-ALOG(THETAA))*DZKL 540 ELSEIF(THETAP.GT.THETAA)THEN 541 CAPE(I,J)=CAPE(I,J) 542 cjjt * +G*(ALOG(THETAP/THETAA))*DZKL 543 * +G*(ALOG(THETAP)-ALOG(THETAA))*DZKL 544 ENDIF 545 ENDIF 546 ENDDO 547 ENDDO 548 ENDDO 549 550 C 551 C ENFORCE LOWER LIMIT OF 0.0 ON CAPE AND UPPER 552 C LIMIT OF 0.0 ON CINS. 553 C 554 !$omp parallel do 555 DO 40 J=JSTA,JEND 556 DO 40 I=1,IM 557 c IF (CAPE(K).LT.D00) 558 c X WRITE(STDOUT,*) 559 c X 'CALCAPE: NEGATIVE CAPE AT K=',K,CAPE(K) 560 c IF (CINS(K).GT.D00) 561 c X WRITE(STDOUT,*) 562 c X 'CALCAPE: POSTIVE CINS AT K = ',K,CINS(K) 563 CAPE(I,J) = AMAX1(D00,CAPE(I,J)) 564 CINS(I,J) = AMIN1(CINS(I,J),D00) 565 40 CONTINUE 566 C 567 C END OF ROUTINE. 568 C 569 RETURN 570 END Page 9 Source Listing CALCAPE 2014-12-17 20:47 Entry Points CALCAPE.f ENTRY POINTS Name calcape_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 276 209,210,215 20 Label 277 203,204 30 Label 355 301 40 Label 434 424,425 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 212,213,219,286,309,325,333,400 AINT Func 229 scalar 229,249 ALBASE Scalar 14 R(4) 4 2 91719 COM ALBEDO Scalar 15 R(4) 4 2 91719 COM ALOG Func 408 scalar 408,412 ALPINT Scalar 4 R(4) 4 3 4219074 COM AMAX1 Func 432 scalar 432 AMIN1 Func 433 scalar 433 APEBTK Local 219 R(4) 4 scalar 219,224,227 APESPK Local 269 R(4) 4 scalar 269,270 BETA Param 37 R(4) 4 scalar 37 BOUNDL Subr 196 196 BQK Local 246 R(4) 4 scalar 246,248 BQS00K Local 241 R(4) 4 scalar 241,246 BQS10K Local 243 R(4) 4 scalar 243,246 BTG Param 37 R(4) 4 scalar CALCAPE Subr 1 CAPA Param 28 R(4) 4 scalar 35,219,224,269,403,404 CAPE Dummy 1 R(4) 4 2 91719 ARG,INOUT 171,410,432 CFT0 Param 33 R(4) 4 scalar CI Scalar 4 R(4) 4 scalar COM CINS Dummy 1 R(4) 4 2 91719 ARG,INOUT 172,406,433 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,121 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 171,172,173,174,183,234,238,254,25 8,432,433 D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar Page 10 Source Listing CALCAPE 2014-12-17 20:47 Symbol Table CALCAPE.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 401 DETA Scalar 7 R(4) 4 1 45 COM 402 DFRLG Scalar 7 R(4) 4 1 46 COM DI Scalar 4 R(4) 4 scalar COM DP10M Param 35 R(4) 4 scalar DPBND Param 120 R(4) 4 scalar 214 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 DZKL Local 401 R(4) 4 scalar 401,408,412 ELDR Param 36 R(4) 4 scalar ELIVW Param 121 R(4) 4 scalar 121 ELOCP Param 121 R(4) 4 scalar 270 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 EXP Func 270 scalar 270 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,408,412 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,401 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 219,224,269,403,404 H1E2 Param 30 R(4) 4 scalar H1E4 Param 30 R(4) 4 scalar H1M12 Param 26 R(4) 4 scalar 196 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 196 HBM2 Scalar 3 R(4) 4 2 91719 COM 193,194,195 HBOT Scalar 13 R(4) 4 2 91719 COM Page 11 Source Listing CALCAPE 2014-12-17 20:47 Symbol Table CALCAPE.f Name Object Declared Type Bytes Dimen Elements Attributes References HDAC Scalar 16 R(4) 4 2 91719 COM HDACV Scalar 16 R(4) 4 2 91719 COM HTM Scalar 5 R(4) 4 3 4127355 COM 326,334 HTOP Scalar 13 R(4) 4 2 91719 COM I Local 170 I(4) 4 scalar 170,171,172,173,174,175,182,183,21 0,211,212,213,217,218,219,221,222, 223,229,234,238,246,247,249,254,25 8,267,268,272,273,274,285,286,287, 293,294,306,307,308,309,310,313,31 6,340,341,342,348,349,350,378,379, 380,386,387,388,389,398,399,400,40 1,402,403,404,406,410,425,432,433 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 IDX Local 134 I(4) 4 2 91719 387,389,399 IEQL Local 126 I(4) 4 2 91719 175,342,350,379,388 IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 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 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,12,13,14,15,16, 126,127,130,131,132,133,134,170,18 2,196,210,285,293,306,340,348,378, 386,398,425 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMT Param 4 I(4) 4 scalar INPES Param 6 I(4) 4 scalar 11,14 INT Func 230 scalar 230,250 IOUNIT Common 4 20 IOUT Scalar 3 I(4) 4 scalar COM IPTB Local 126 I(4) 4 2 91719 327,335 IQ Local 250 I(4) 4 scalar 250,252,253,256,257,262,263,264,26 5 ISMTHP Param 122 I(4) 4 scalar 193 ISMTHQ Param 122 I(4) 4 scalar 195 ISMTHT Param 122 I(4) 4 scalar 194 IT Local 261 I(4) 4 scalar 261,262,263,264,265 ITAIL Param 11 I(4) 4 scalar ITB Param 1 I(4) 4 scalar 8,18,256,257,325 ITBQ Param 1 I(4) 4 scalar 9,18,333 ITHTB Local 126 I(4) 4 2 91719 327,335 ITTBK Local 230 I(4) 4 scalar 230,232,233,236,237,241,242,243,24 4,261 ITYPE Dummy 1 I(4) 4 scalar ARG,INOUT 192,204,214,216 IUP Scalar 7 I(4) 4 scalar COM IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM J Local 169 I(4) 4 scalar 169,171,172,173,174,175,181,183,20 Page 12 Source Listing CALCAPE 2014-12-17 20:47 Symbol Table CALCAPE.f Name Object Declared Type Bytes Dimen Elements Attributes References 9,211,212,213,217,218,219,221,222, 223,229,234,238,246,247,249,254,25 8,267,268,272,273,274,284,286,287, 292,294,305,307,308,309,310,313,31 6,339,341,342,347,349,350,377,379, 380,385,387,388,389,397,399,400,40 1,402,403,404,406,410,424,432,433 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 169,181,424 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 209,284,292,305,339,347,377,385,39 7 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, 126,127,130,131,132,133,134,196 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 169,181,424 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 209,284,292,305,339,347,377,385,39 7 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,236,237,325 JTBQ Param 1 I(4) 4 scalar 18,333 KB Local 203 I(4) 4 scalar 203,204,213,217,218,219 KHRES Local 127 I(4) 4 2 91719 308,316,333,341 KLRES Local 127 I(4) 4 2 91719 307,313,325,349 KMNT Scalar 5 I(4) 4 2 4118670 COM KMNTM Scalar 5 I(4) 4 1 45 COM KNUMH Local 304 I(4) 4 scalar 304,315,332 KNUML Local 303 I(4) 4 scalar 303,312,324 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 KTM Scalar 3 I(4) 4 scalar COM L Local 180 I(4) 4 scalar 180,183,281,286,287,301,309,310,32 5,326,333,334,342,350,384,388,400, 401,402,403,404 L1D Dummy 1 I(4) 4 2 91719 ARG,INOUT LB Param 9 I(4) 4 scalar LBEG Local 372 I(4) 4 scalar 372,379,384 Page 13 Source Listing CALCAPE 2014-12-17 20:47 Symbol Table CALCAPE.f Name Object Declared Type Bytes Dimen Elements Attributes References LCL Local 127 I(4) 4 2 91719 173,287,294,310,380,388 LCNTRL Scalar 5 I(4) 4 scalar COM LEND Local 373 I(4) 4 scalar 373,380,384 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,132,175,180,203,281,301 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 211,294 LMHK Local 211 I(4) 4 scalar 211,212 LMV Scalar 5 I(4) 4 2 91719 COM LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4,7 LSM Param 3 I(4) 4 scalar LUNCO Scalar 5 I(4) 4 scalar COM LUNLL Scalar 5 I(4) 4 scalar COM LUNOUT Scalar 5 I(4) 4 scalar COM MASKS Common 2 34853220 MAX Func 380 scalar 380 ME Scalar 6 I(4) 4 scalar COM MIN Func 379 scalar 379 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 OVERRC Param 33 R(4) 4 scalar 33 P00K Local 262 R(4) 4 scalar 262,267,268 P01K Local 264 R(4) 4 scalar 264,267,268 P1000 Param 35 R(4) 4 scalar P10K Local 263 R(4) 4 scalar 263,267,268 P11K Local 265 R(4) 4 scalar 265,268 P1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 193,221 P2FILT Subr 193 193,194,195 PARALLEL Common 9 8240 PD Scalar 3 R(4) 4 2 91719 COM PDSL Scalar 6 R(4) 4 2 91719 COM 212,213,219,286,309,325,333,400,40 2 PDVP1 Scalar 6 R(4) 4 2 91719 COM Page 14 Source Listing CALCAPE 2014-12-17 20:47 Symbol Table CALCAPE.f Name Object Declared Type Bytes Dimen Elements Attributes References PFSL Scalar 7 R(4) 4 2 91719 COM PHYS Common 2 7690036 PINT Scalar 4 R(4) 4 3 4219074 COM PKL Local 213 R(4) 4 scalar 213,214,221,224,286,287,309,311 PL Scalar 5 R(4) 4 scalar COM 326 PLQ Scalar 6 R(4) 4 scalar COM 311,334 PP Local 133 R(4) 4 2 91719 249,254,258,267,268,326,334 PQ0 Param 41 R(4) 4 scalar PRESK Local 400 R(4) 4 scalar 400,402,403,404 PSFCK Local 212 R(4) 4 scalar 212,214 PSLP Scalar 6 R(4) 4 2 91719 COM PSP Local 132 R(4) 4 2 91719 273,287 PT Scalar 3 R(4) 4 scalar COM 212,213,219,286,309,326,334,400 PTBL Scalar 18 R(4) 4 2 10184 COM 262,263,264,265 PTSL Scalar 7 R(4) 4 2 91719 COM Q Scalar 5 R(4) 4 3 4127355 COM 218,401 Q1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 195,196,223 QBTK Local 218 R(4) 4 scalar 218,223,248,270 QQ Local 133 R(4) 4 2 91719 229,234,238,246,247,267,268,326,33 4 QS0 Scalar 8 R(4) 4 1 134 COM 241,243 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 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 326 RDPQ Scalar 6 R(4) 4 scalar COM 334 RDQ Scalar 5 R(4) 4 scalar COM 248 RDTH Scalar 5 R(4) 4 scalar COM 228 RDTHE Scalar 5 R(4) 4 scalar COM 326 RDTHEQ Scalar 6 R(4) 4 scalar COM 334 RES Scalar 3 R(4) 4 2 91719 COM 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 402 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 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 SQK Local 247 R(4) 4 scalar 247,248 SQS Scalar 8 R(4) 4 1 134 COM 242,244 SQS00K Local 242 R(4) 4 scalar 242,247 SQS10K Local 244 R(4) 4 scalar 244,247 SST Scalar 14 R(4) 4 2 91719 COM STDOUT Scalar 2 I(4) 4 scalar COM Page 15 Source Listing CALCAPE 2014-12-17 20:47 Symbol Table CALCAPE.f Name Object Declared Type Bytes Dimen Elements Attributes References STDRD Scalar 10 R(4) 4 scalar COM STHE Scalar 8 R(4) 4 1 76 COM 326 STHEQ Scalar 9 R(4) 4 1 152 COM 334 SWNET Scalar 12 R(4) 4 2 91719 COM T Scalar 5 R(4) 4 3 4127355 COM 217,342,350,401,404 T1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 194,222 T500 Scalar 6 R(4) 4 2 91719 COM TBTK Local 217 R(4) 4 scalar 217,222,227 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 326 THE0Q Scalar 9 R(4) 4 1 152 COM 334 THESP Local 132 R(4) 4 2 91719 174,272,274,326,335 THETAA Local 404 R(4) 4 scalar 404,405,408,409,412 THETAP Local 403 R(4) 4 scalar 403,405,408,409,412 THL Scalar 5 R(4) 4 scalar COM 228 TPAR Local 132 R(4) 4 3 4127355 183,325,333,342,350,403 TPSPK Local 267 R(4) 4 scalar 267,269,273 TQK Local 248 R(4) 4 scalar 248,249,250 TSL Scalar 8 R(4) 4 2 91719 COM TTBL Scalar 18 R(4) 4 2 10184 COM 325 TTBLEX Subr 325 325,333 TTBLQ Scalar 18 R(4) 4 2 66880 COM 333 TTHBTK Local 227 R(4) 4 scalar 227,228,270 TTHESK Local 270 R(4) 4 scalar 270,272,274 TTHK Local 228 R(4) 4 scalar 228,229,230 TTV Scalar 8 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM V Scalar 5 R(4) 4 3 4127355 COM 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 Z1000 Scalar 9 R(4) 4 2 91719 COM ZINT Scalar 4 R(4) 4 3 4219074 COM Page 16 Source Listing CALCAPE 2014-12-17 20:47 Subprograms/Common Blocks CALCAPE.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALCAPE Subr 1 CTLBLK Common 2 96 EXTRA Common 3 71506268 IOUNIT Common 4 20 LOOPS Common 2 748952 MASKS Common 2 34853220 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__ -D __SSE__ -D __MMX__ Page 17 Source Listing CALCAPE 2014-12-17 20:47 CALCAPE.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 : CALCAPE.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100