Page 1 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 1 SUBROUTINE ETA2P(IMOUT,JMOUT) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: ETA2P VERT INTRP OF ETA TO PRESSURE 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 92-12-21 6 C 7 C ABSTRACT: 8 C FOR MOST APPLICATIONS THIS ROUTINE IS THE WORKHORSE 9 C OF THE POST PROCESSOR. IN A NUTSHELL IT INTERPOLATES 10 C DATA FROM ETA TO PRESSURE SURFACES. IT ORIGINATED 11 C FROM THE VERTICAL INTERPOLATION CODE IN THE OLD ETA 12 C POST PROCESSOR SUBROUTINE OUTMAP. 13 C . 14 C 15 C PROGRAM HISTORY LOG: 16 C 92-12-21 RUSS TREADON 17 C 96-03-21 GEOFF MANIKIN - ADDED CLOUD ICE ON P 18 C 98-06-16 T BLACK - CONVERSION FROM 1-D TO 2-D 19 C 98-07-17 MIKE BALDWIN - REMOVED LABL84 20 C 98-08-18 T BLACK - REMOVED MOST 3-D ARRAYS FROM 21 C COMMON BLOCK JIMA 22 C 98-12-22 MIKE BALDWIN - BACK OUT RH OVER ICE 23 C 00-01-04 JIM TUCCILLO - MPI VERSION 24 C 25 C USAGE: CALL ETA2P(IMOUT,JMOUT) 26 C INPUT ARGUMENT LIST: 27 C IMOUT - FIRST DIMENSION OF OUTPUT GRID 28 C JMOUT - SECOND DIMENSION OF OUTPUT GRID 29 C 30 C OUTPUT ARGUMENT LIST: 31 C NONE 32 C 33 C OUTPUT FILES: 34 C NONE 35 C 36 C SUBPROGRAMS CALLED: 37 C UTILITIES: 38 C SCLFLD - SCALE ARRAY ELEMENTS BY CONSTANT. 39 C E2OUT - E-GRID TO OUTPUT GRID INTERPOLATION/SMOOTHING. 40 C OUTPUT - POST ARRAY TO OUTPUT FILE. 41 C CALPOT2 - COMPUTE POTENTIAL TEMPERATURE. 42 C CALRH2 - COMPUTE RELATIVE HUMIDITY. 43 C CALDWP2 - COMPUTE DEWPOINT TEMPERATURE. 44 C BOUND - BOUND ARRAY ELEMENTS BETWEEN LOWER AND UPPER LIMITS. 45 C CALMCVG - COMPUTE MOISTURE CONVERGENCE. 46 C CALVOR - COMPUTE ABSOLUTE VORTICITY. 47 C CALSTRM - COMPUTE GEOSTROPHIC STREAMFUNCTION. 48 C 49 C LIBRARY: 50 C COMMON - OMGAOT 51 C LOOPS 52 C MASKS 53 C MAPOT 54 C VRBLS 55 C PVRBLS 56 C RQSTFLD 57 C EXTRA Page 2 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 58 C CLDWTR 59 C 60 C ATTRIBUTES: 61 C LANGUAGE: FORTRAN 90 62 C MACHINE : IBM SP 63 C$$$ 64 C 65 C 66 C 67 C INCLUDE ETA MODEL DIMENSIONS. SET/DERIVE OTHER PARAMETERS. 68 C GAMMA AND RGAMOG ARE USED IN THE EXTRAPOLATION OF VIRTUAL 69 C TEMPERATURES BEYOND THE UPPER OF LOWER LIMITS OF ETA DATA. 70 C 71 INCLUDE "parmeta" 89 INCLUDE "parmout" 93 INCLUDE "params" 94 C 136 PARAMETER (GAMMA=6.5E-3,RGAMOG=RD*GAMMA/G) 137 C 138 C DECLARE VARIABLES. 139 C 140 LOGICAL RUN,FIRST,RESTRT,SIGMA,OLDRD,STDRD 141 LOGICAL IOOMG,IOALL 142 REAL OSL(IM,JM),USL(IM,JM),VSL(IM,JM) 143 REAL PRSI(IM,JM),QCSL(IM,JM),Q2SL(IM,JM) 144 REAL ICE(IM,JM),IW(IM,JM,LM),IWU,IWL 145 146 REAL EGRID1(IM,JM),EGRID2(IM,JM) 147 REAL GRID1(IMOUT,JMOUT),GRID2(IMOUT,JMOUT) 148 C 149 C INCLUDE COMMON BLOCKS. 150 INCLUDE "CTLBLK.comm" 165 INCLUDE "OMGAOT.comm" 169 INCLUDE "LOOPS.comm" 177 INCLUDE "MASKS.comm" 185 INCLUDE "MAPOT.comm" 198 INCLUDE "VRBLS.comm" 205 INCLUDE "PVRBLS.comm" 217 INCLUDE "RQSTFLD.comm" 230 INCLUDE "EXTRA.comm" 242 INCLUDE "CLDWTR.comm" 248 INCLUDE "E2PFLG.comm" 249 C 255 COMMON/JIMA/NL1X(IM,JM),ALPETUX(IM,JM),ALPET2X(IM,JM) 256 257 C 258 C****************************************************************************** 259 C 260 C START ETA2P. 261 C 262 C SET TOTAL NUMBER OF POINTS ON OUTPUT GRID. 263 C 264 C--------------------------------------------------------------- 265 C 266 C *** PART I *** 267 C 268 C VERTICAL INTERPOLATION OF EVERYTHING ELSE. EXECUTE ONLY Page 3 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 269 C IF THERE'S SOMETHING WE WANT. 270 C 271 IF((IGET(012).GT.0).OR.(IGET(013).GT.0).OR. 272 X (IGET(014).GT.0).OR.(IGET(015).GT.0).OR. 273 X (IGET(016).GT.0).OR.(IGET(017).GT.0).OR. 274 X (IGET(018).GT.0).OR.(IGET(019).GT.0).OR. 275 X (IGET(020).GT.0).OR.(IGET(030).GT.0).OR. 276 X (IGET(021).GT.0).OR.(IGET(022).GT.0).OR. 277 X (IGET(153).GT.0).OR.(IGET(166).GT.0))THEN 278 C 279 C SET UP UTIM FOR THIS TIME STEP 280 C 281 UTIM=1. 282 CLIMIT =1.0E-20 283 C 284 DO 975 L=1,LM 285 IF(L.EQ.1)THEN 286 !$omp parallel do 287 DO J=JSTA,JEND 288 DO I=1,IM 289 IW(I,J,L)=0. 290 ENDDO 291 ENDDO 292 GO TO 975 293 ENDIF 294 C 295 !$omp parallel do 296 !$omp& private(cwmkl,fiq,hh,lml,pp,qi,qkl,qw,tkl,tmt0,tmt15,u00kl) 297 DO 970 J=JSTA,JEND 298 DO 970 I=1,IM 299 LML=LM-LMH(I,J) 300 HH=HTM(I,J,L)*HBM2(I,J) 301 TKL=T(I,J,L) 302 QKL=Q(I,J,L) 303 CWMKL=CWM(I,J,L) 304 TMT0=(TKL-273.16)*HH 305 TMT15=AMIN1(TMT0,-15.)*HH 306 PP=PDSL(I,J)*AETA(L)+PT 307 QW=HH*PQ0/PP*EXP(HH*A2*(TKL-A3)/(TKL-A4)) 308 QI=QW*(1.+0.01*AMIN1(TMT0,0.)) 309 U00KL=U00(I,J)+UL(L+LML)*(0.95-U00(I,J))*UTIM 310 C 311 IF(TMT0.LT.-15.0)THEN 312 FIQ=QKL-U00KL*QI 313 IF(FIQ.GT.D00.OR.CWMKL.GT.CLIMIT) THEN 314 IW(I,J,L)=1. 315 ELSE 316 IW(I,J,L)=0. 317 ENDIF 318 ENDIF 319 C 320 IF(TMT0.GE.0.0)IW(I,J,L)=0. 321 IF(TMT0.LT.0.0.AND.TMT0.GE.-15.0)THEN 322 IW(I,J,L)=0. 323 IF(IW(I,J,L-1).EQ.1.0.AND.CWMKL.GT.CLIMIT)IW(I,J,L)=1. 324 ENDIF 325 970 CONTINUE Page 4 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 326 975 CONTINUE 327 C 328 C 329 C VERTICAL INTERPOLATION OF GEOPOTENTIAL, SPECIFIC HUMIDITY, 330 C AND TEMPERATURE. START AT THE UPPERMOST TARGET PRESSURE LEVEL. 331 C 332 ALPTH = ALOG(1.E5) 333 DO 310 L=1,LSL 334 C 335 C 336 c IF(IOALL)GO TO 225 337 C 338 TRF=H2*ALSL(L) 339 C 340 C LOOP OVER HORIZONTAL GRID. 341 C 342 !$omp parallel do 343 !$omp& private(lma,lmap1,ppdsl) 344 DO 180 J=JSTA,JEND 345 DO 180 I=1,IM 346 LMA=LM 347 CX IF(OLDRD)LMA=LMH(I,J) 348 LMAP1=LMA+1 349 C 350 C SET PRESSURE DEPTH IN THIS COLUMN. 351 C 352 PPDSL=PDSL(I,J) 353 C 354 C LOCATE VERTICAL INDEX OF ETA INTERFACE PRESSURES BOUNDING 355 C THE STANDARD PRESSURE LEVEL TO WHICH WE'RE INTERPOLATING. 356 C 357 DO 170 IL=2,LMAP1 358 IF((ALSL(L)-ALPINT(I,J,IL)).GT.D00)GO TO 170 359 NL1X(I,J)=IL 360 GO TO 180 361 170 CONTINUE 362 NL1X(I,J)=LMAP1 363 180 CONTINUE 364 365 !$omp parallel do 366 !$omp& private(ahf,ahfo,ahfq,ahfq2,ahfqc,ahfqi,ai,b,bi,bom, 367 !$omp& bq,bq2,bqc,bqc_2,bqi_2,bqi,fac,gmiw,gmiw_2,iwl,iwu, 368 !$omp& lma,lmap1,pl,pnl1,pu,q2a,q2b,qabv,qi,qint,ql, 369 !$omp& qsat,qu,qw,rhu,tabv,tblo,tl,tmt0,tmt15,tu, 370 !$omp& tvrabv,tvrblo,tvrl,tvru,zl,zu) 371 DO 220 J=JSTA,JEND 372 DO 220 I=1,IM 373 LMA =LM 374 LMAP1=LMA+1 375 IF((TRF-ALPINT(I,J,NL1X(I,J)) 376 1 -ALPINT(I,J,NL1X(I,J)-1)).LE.D00) 377 2 NL1X(I,J)=NL1X(I,J)-1 378 PNL1=PINT(I,J,NL1X(I,J)) 379 C 380 C BRANCH TO APPROPRIATE BLOCK TO COMPUTE COEFFICIENTS. 381 C 382 IF(NL1X(I,J).EQ.1)THEN Page 5 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 383 PU=PINT(I,J,2) 384 ZU=ZINT(I,J,2) 385 TU=D50*(T(I,J,1)+T(I,J,2)) 386 QU=D50*(Q(I,J,1)+Q(I,J,2)) 387 C 388 IWU=D50*(IW(I,J,1)+IW(I,J,2)) 389 TMT0=TU-273.16 390 TMT15=AMIN1(TMT0,-15.) 391 AI=0.008855 392 BI=1. 393 IF(TMT0.LT.-20.)THEN 394 AI=0.007225 395 BI=0.9674 396 ENDIF 397 QW=PQ0/PU 398 1 *EXP(A2*(TU-A3)/(TU-A4)) 399 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 400 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 401 IF(TMT0.LT.-15.)THEN 402 QSAT=QI 403 ELSEIF(TMT0.GE.0.)THEN 404 QSAT=QINT 405 ELSE 406 IF(IWU.GT.0.0) THEN 407 QSAT=QI 408 ELSE 409 QSAT=QINT 410 ENDIF 411 ENDIF 412 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 413 C DELETE THIS LINE TO SWITCH BACK TO RH VS ICE 414 QSAT=QW 415 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 416 RHU =QU/QSAT 417 C 418 IF(RHU.GT.H1)THEN 419 RHU=H1 420 QU =RHU*QSAT 421 ENDIF 422 C 423 IF(RHU.LT.D01)THEN 424 RHU=D01 425 QU =RHU*QSAT 426 ENDIF 427 C 428 TVRU=TU*(H1+D608*QU) 429 TVRABV=TVRU*(SPL(L)/PU)**RGAMOG 430 TABV=TVRABV/(H1+D608*QU) 431 C 432 C 433 TMT0=TABV-273.16 434 TMT15=AMIN1(TMT0,-15.) 435 AI=0.008855 436 BI=1. 437 IF(TMT0.LT.-20.)THEN 438 AI=0.007225 439 BI=0.9674 Page 6 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 440 ENDIF 441 QW=PQ0/SPL(L) 442 1 *EXP(A2*(TABV-A3)/(TABV-A4)) 443 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 444 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 445 IF(TMT0.LT.-15.)THEN 446 QSAT=QI 447 ELSEIF(TMT0.GE.0.)THEN 448 QSAT=QINT 449 ELSE 450 IF(IWU.GT.0.0) THEN 451 QSAT=QI 452 ELSE 453 QSAT=QINT 454 ENDIF 455 ENDIF 456 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 457 C DELETE THIS LINE TO SWITCH BACK TO RH VS ICE 458 QSAT=QW 459 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 460 QABV =RHU*QSAT 461 QABV =AMAX1(H1M12,QABV) 462 B =TABV 463 BQ =QABV 464 BOM =OMGA(I,J,1) 465 GMIW =IW(I,J,1) 466 BQC =(1.-GMIW)*CWM(I,J,1) 467 BQI =GMIW*CWM(I,J,1) 468 Q2A =D50*(Q2(I,J,1)+Q2(I,J,2)) 469 BQ2 =Q2A 470 AHF =D00 471 AHFQ =D00 472 AHFO =D00 473 AHFQC=D00 474 AHFQI=D00 475 AHFQ2=D00 476 FAC =D00 477 C 478 ELSEIF(NL1X(I,J).EQ.LMAP1)THEN 479 C 480 C EXTRAPOLATION AT LOWER BOUND. THE LOWER BOUND IS 481 C LM IF OLDRD=.FALSE. IF OLDRD=.TRUE. THE LOWER 482 C BOUND IS THE FIRST ATMOSPHERIC ETA LAYER. 483 C 484 PL=PINT(I,J,LMA-1) 485 ZL=ZINT(I,J,LMA-1) 486 TL=D50*(T(I,J,LMA-2)+T(I,J,LMA-1)) 487 QL=D50*(Q(I,J,LMA-2)+Q(I,J,LMA-1)) 488 C 489 IWL=D50*(IW(I,J,LMA-2)+IW(I,J,LMA-1)) 490 491 TMT0=TL-273.16 492 TMT15=AMIN1(TMT0,-15.) 493 AI=0.008855 494 BI=1. 495 IF(TMT0.LT.-20.)THEN 496 AI=0.007225 Page 7 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 497 BI=0.9674 498 ENDIF 499 QW=PQ0/PL 500 1 *EXP(A2*(TL-A3)/(TL-A4)) 501 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 502 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 503 IF(TMT0.LT.-15.)THEN 504 QSAT=QI 505 ELSEIF(TMT0.GE.0.)THEN 506 QSAT=QINT 507 ELSE 508 IF(IWL.GT.0.0) THEN 509 QSAT=QI 510 ELSE 511 QSAT=QINT 512 ENDIF 513 ENDIF 514 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 515 C DELETE THIS LINE TO SWITCH BACK TO RH VS ICE 516 QSAT=QW 517 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 518 RHL=QL/QSAT 519 C 520 IF(RHL.GT.H1)THEN 521 RHL=H1 522 QL =RHL*QSAT 523 ENDIF 524 C 525 IF(RHL.LT.D01)THEN 526 RHL=D01 527 QL =RHL*QSAT 528 ENDIF 529 C 530 TVRL =TL*(H1+D608*QL) 531 TVRBLO=TVRL*(SPL(L)/PL)**RGAMOG 532 TBLO =TVRBLO/(H1+D608*QL) 533 C 534 TMT0=TBLO-273.16 535 TMT15=AMIN1(TMT0,-15.) 536 AI=0.008855 537 BI=1. 538 IF(TMT0.LT.-20.)THEN 539 AI=0.007225 540 BI=0.9674 541 ENDIF 542 QW=PQ0/SPL(L) 543 1 *EXP(A2*(TBLO-A3)/(TBLO-A4)) 544 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 545 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 546 IF(TMT0.LT.-15.)THEN 547 QSAT=QI 548 ELSEIF(TMT0.GE.0.)THEN 549 QSAT=QINT 550 ELSE 551 IF(IWL.GT.0.0) THEN 552 QSAT=QI 553 ELSE Page 8 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 554 QSAT=QINT 555 ENDIF 556 ENDIF 557 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 558 C DELETE THIS LINE TO SWITCH BACK TO RH VS ICE 559 QSAT=QW 560 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 561 QBLO =RHL*QSAT 562 QBLO =AMAX1(H1M12,QBLO) 563 B =TBLO 564 BQ =QBLO 565 BOM =OMGA(I,J,LMA) 566 GMIW =IW(I,J,LMH(I,J)) 567 BQC =(1.-GMIW)*CWM(I,J,LMH(I,J)) 568 BQI =GMIW*CWM(I,J,LMH(I,J)) 569 Q2A =D50*(Q2(I,J,LMH(I,J)-1)+Q2(I,J,LMH(I,J))) 570 BQ2 =Q2A 571 AHF =D00 572 AHFQ =D00 573 AHFO =D00 574 AHFQC=D00 575 AHFQI=D00 576 AHFQ2=D00 577 FAC =D00 578 C 579 ELSE 580 C 581 C INTERPOLATION BETWEEN LOWER AND UPPER BOUNDS. 582 C 583 B =T(I,J,NL1X(I,J)) 584 BQ =Q(I,J,NL1X(I,J)) 585 BOM =OMGA(I,J,NL1X(I,J)) 586 GMIW =IW(I,J,NL1X(I,J)) 587 BQC =(1.-GMIW)*CWM(I,J,NL1X(I,J)) 588 BQI =GMIW*CWM(I,J,NL1X(I,J)) 589 GMIW_2=IW(I,J,NL1X(I,J)-1) 590 BQC_2 =(1.-GMIW)*CWM(I,J,NL1X(I,J)-1) 591 BQI_2 =GMIW*CWM(I,J,NL1X(I,J)-1) 592 Q2B =D50*(Q2(I,J,NL1X(I,J)-1)+Q2(I,J,NL1X(I,J))) 593 C 594 IF(NL1X(I,J).GT.2)THEN 595 Q2A=D50*(Q2(I,J,NL1X(I,J)-2)+Q2(I,J,NL1X(I,J)-1)) 596 ELSE 597 Q2A=Q2B 598 ENDIF 599 C 600 BQ2=Q2B*HTM(I,J,NL1X(I,J)) 601 FAC =H2*ALOG(PT+PDSL(I,J)*AETA(NL1X(I,J))) 602 AHF =(B-T(I,J,NL1X(I,J)-1))/ 603 1 (ALPINT(I,J,NL1X(I,J)+1)-ALPINT(I,J,NL1X(I,J)-1)) 604 AHFQ =(BQ-Q(I,J,NL1X(I,J)-1))/ 605 1 (ALPINT(I,J,NL1X(I,J)+1)-ALPINT(I,J,NL1X(I,J)-1)) 606 AHFO =(BOM-OMGA(I,J,NL1X(I,J)-1))/ 607 1 (ALPINT(I,J,NL1X(I,J)+1)-ALPINT(I,J,NL1X(I,J)-1)) 608 AHFQC=(BQC-BQC_2)/ 609 1 (ALPINT(I,J,NL1X(I,J)+1)-ALPINT(I,J,NL1X(I,J)-1)) 610 AHFQI=(BQI-BQI_2)/ Page 9 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 611 1 (ALPINT(I,J,NL1X(I,J)+1)-ALPINT(I,J,NL1X(I,J)-1)) 612 AHFQ2=(BQ2-Q2A*HTM(I,J,NL1X(I,J)-1))/ 613 1 (ALPINT(I,J,NL1X(I,J)+1)-ALPINT(I,J,NL1X(I,J)-1)) 614 ENDIF 615 C 616 TSL(I,J)=B+AHF*(TRF-FAC) 617 QSL(I,J)=BQ+AHFQ*(TRF-FAC) 618 QSL(I,J)=AMAX1(QSL(I,J),H1M12) 619 IF (QSL(I,J).LT.1.5E-12) QSL(I,J)=0.0 ! NARR 620 OSL(I,J)=BOM+AHFO*(TRF-FAC) 621 QCSL(I,J)=BQC+AHFQC*(TRF-FAC) 622 QCSL(I,J)=AMAX1(QCSL(I,J),H1M12) 623 IF (QCSL(I,J).LT.1.5E-12) QCSL(I,J)=0.0 ! NARR 624 ICE(I,J)=BQI+AHFQI*(TRF-FAC) 625 ICE(I,J)=AMAX1(ICE(I,J),H1M12) 626 IF (ICE(I,J).LT.1.5E-12) ICE(I,J)=0.0 ! NARR 627 Q2SL(I,J)=BQ2+AHFQ2*(TRF-FAC) 628 Q2SL(I,J)=AMAX1(Q2SL(I,J),D00) 629 FSL(I,J)=(PNL1-SPL(L))/(SPL(L)+PNL1) 630 1 *((ALSL(L)+ALPINT(I,J,NL1X(I,J))-FAC)*AHF+B)*R*H2 631 2 +ZINT(I,J,NL1X(I,J))*G 632 220 CONTINUE 633 C 634 C LOAD GEOPOTENTIAL AND TEMPERATURE INTO STANDARD LEVEL 635 C ARRAYS FOR THE NEXT PASS. 636 C 637 225 CONTINUE 638 C 639 C SAVE 500MB TEMPERATURE FOR LIFTED INDEX. 640 C 641 IF((NINT(SPL(L)).EQ.50000).AND. 642 1 ((IGET(030).GT.0).OR.(IGET(031).GT.0).OR. 643 2 (IGET(075).GT.0)))THEN 644 !$omp parallel do 645 DO J=JSTA,JEND 646 DO I=1,IM 647 T500(I,J)=TSL(I,J) 648 ENDDO 649 ENDDO 650 ENDIF 651 C 652 C CALCULATE 1000MB GEOPOTENTIALS CONSISTENT WITH SLP OBTAINED 653 C FROM THE MESINGER OR NWS SHUELL SLP REDUCTION. 654 C 655 IF(NINT(SPL(L)).EQ.100000)THEN 656 C 657 C MESINGER SLP 658 C 659 IF(IGET(023).GT.0)THEN 660 !$omp parallel do 661 DO J=JSTA,JEND 662 DO I=1,IM 663 ALPSL=ALOG(PSLP(I,J)) 664 IF(FIS(I,J).GT.H1)THEN 665 FSL(I,J)=FIS(I,J)/(ALPSL-ALOG(PD(I,J)+PT))* 666 1 (ALPSL-ALPTH) 667 ELSE Page 10 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 668 FSL(I,J)=R*T(I,J,LM)*(ALPSL-ALPTH) 669 ENDIF 670 Z1000(I,J)=FSL(I,J)*GI 671 ENDDO 672 ENDDO 673 C 674 C NWS SHUELL SLP. NGMSLP2 COMPUTES 1000MB GEOPOTENTIAL. 675 C 676 ELSE 677 !$omp parallel do 678 DO J=JSTA,JEND 679 DO I=1,IM 680 FSL(I,J)=Z1000(I,J)*G 681 ENDDO 682 ENDDO 683 ENDIF 684 ENDIF 685 C 686 C INTERPOLATE WIND COMPONENTS FROM ETA TO PRESSURE. 687 C 688 IF((IGET(018).GT.0).OR.(IGET(019).GT.0).OR. 689 1 (IGET(021).GT.0).OR.(IGET(085).GT.0))THEN 690 C 691 !$omp parallel do 692 DO J=JSTA,JEND 693 DO I=1,IM 694 USL(I,J)=D00 695 VSL(I,J)=D00 696 ENDDO 697 ENDDO 698 C 699 !!$omp parallel do 700 !!$omp& private(alpet2,alpetl,alpetu,lmb,petal,petau) 701 !!$omp& shared (alpet2x,alpetux,nl1x) 702 DO 281 J=JSTA,JEND 703 DO 281 I=1,IM 704 CNOTE 705 CNOTE 29 JANUARY 1993, RUSS TREADON. 706 CNOTE - AS FOR THE OTHER FIELDS WE INTERPOLATE ONLY 707 CNOTE BETWEEN THE FAL AND THE MODEL TOP. BELOW 708 CNOTE SURFACE VALUES ARE FAL VALUES. 709 C 710 LMB = LMV(I,J) 711 C 712 PETAU=PT+PDVP1(I,J)*ETA(1) 713 ALPETU=ALOG(PETAU) 714 DO 280 IL=2,LMB 715 PETAL=PT+PDVP1(I,J)*ETA(IL) 716 c PETAU=PT+PDVP1(I,J)*ETA(IL-1) 717 ALPETL=ALOG(PETAL) 718 c ALPETU=ALOG(PETAU) 719 ALPET2=SQRT(0.5E0*(ALPETL*ALPETL+ALPETU*ALPETU)) 720 C 721 C SEARCH FOR HIGHEST MID-LAYER ETA SURFACE (NOT SUBMERGED) 722 C THAT IS BELOW THE GIVEN STANDARD PRESSURE LEVEL. 723 IF(ALSL(L).LT.ALPET2)THEN 724 NL1X(I,J)=IL-1 Page 11 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 725 ALPETUX(I,J)=ALPETU 726 ALPET2X(I,J)=ALPET2 727 GO TO 281 728 ENDIF 729 C If we arent on the last iterate of the 280 loop, reset PETAU and ALPETU 730 if ( il .eq. lmb ) goto 280 731 PETAU=PETAL 732 ALPETU=ALPETL 733 280 CONTINUE 734 NL1X(I,J)=LMB+1 735 ALPETUX(I,J)=ALPETU 736 ALPET2X(I,J)=ALPET2 737 281 CONTINUE 738 C 739 C BELOW GROUND USE FAL WINDS. 740 C 741 !$omp parallel do 742 !$omp& private(alpet1,alpetl,alpetu,fact,petau) 743 DO 290 J=JSTA,JEND 744 DO 290 I=1,IM 745 IF(NL1X(I,J).GT.LMV(I,J))THEN 746 USL(I,J)=U(I,J,LMV(I,J)) 747 VSL(I,J)=V(I,J,LMV(I,J)) 748 C 749 C IF REQUESTED PRESSURE LEVEL IS NOT BELOW THE LOCAL GROUND 750 C THEN WE HAVE TWO POSSIBILITIES. IF THE REQUESTED PRESSURE 751 C LEVEL IS BETWEEN THE LOCAL SURFACE PRESSURE AND TOP OF 752 C MODEL PRESSURE, VERTICALLY INTERPOLATE BETWEEN NEAREST 753 C BOUNDING ETA LEVELS TO GET THE WIND COMPONENTS. IF THE 754 C REQUESTED PRESSURE LEVEL IS ABOVE THE MODEL TOP, USE 755 C CONSTANT EXTRAPOLATION OF TOP ETA LAYER (L=1) WINDS. 756 C 757 ELSE 758 IF(NL1X(I,J).GT.1)THEN 759 ALPETL=ALPETUX(I,J) 760 PETAU=PT+PDVP1(I,J)*ETA(NL1X(I,J)-1) 761 ALPETU=ALOG(PETAU) 762 ALPET1=SQRT(0.5*(ALPETL*ALPETL+ALPETU*ALPETU)) 763 FACT=(ALPET2X(I,J)-ALSL(L))/(ALPET2X(I,J)-ALPET1) 764 USL(I,J)=U(I,J,NL1X(I,J)) 765 1 +(U(I,J,NL1X(I,J)-1)-U(I,J,NL1X(I,J)))*FACT 766 VSL(I,J)=V(I,J,NL1X(I,J))+(V(I,J,NL1X(I,J)-1) 767 1 -V(I,J,NL1X(I,J)))*FACT 768 ELSE 769 USL(I,J)=U(I,J,NL1X(I,J)) 770 VSL(I,J)=V(I,J,NL1X(I,J)) 771 ENDIF 772 C 773 C ALPET2 IS MID-LAYER ETA SURFACE JUST BELOW STANDARD PRESSURE 774 C LEVEL AND ALPET1 IS DASHED ETA SURFACE JUST ABOVE. 775 C NOTE THAT IF THE STANDARD PRESSURE SURFACE IS SUBMERGED, THEN 776 C ALPET2 AND ALPET1 ARE THE LOWEST AND 2ND LOWEST MID-LAYER 777 C ETA SURFACES ABOVE THE TOPOGRAPHY (WITH OLDRD=.TRUE., ZJ). 778 C 779 ENDIF 780 290 CONTINUE 781 ENDIF Page 12 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 782 C 783 C 784 C *** PART II *** 785 C 786 C INTERPOLATE/OUTPUT SELECTED FIELDS. 787 C 788 C GEOPOTENTIAL (SCALE BY GI) 789 IF (IGET(012).GT.0) THEN 790 IF (LVLS(L,IGET(012)).GT.0) THEN 791 !$omp parallel do 792 DO J=JSTA,JEND 793 DO I=1,IM 794 EGRID1(I,J)=FSL(I,J)*GI 795 ENDDO 796 ENDDO 797 C 798 CALL E2OUT(012,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 799 ID(1:25)=0 800 CALL OUTPUT(IOUTYP,IGET(012),L,GRID1,IMOUT,JMOUT) 801 ENDIF 802 ENDIF 803 C 804 C TEMPERATURE. 805 IF (IGET(013).GT.0) THEN 806 IF (LVLS(L,IGET(013)).GT.0) THEN 807 CALL E2OUT(013,000,TSL,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 808 ID(1:25)=0 809 CALL OUTPUT(IOUTYP,IGET(013),L,GRID1,IMOUT,JMOUT) 810 ENDIF 811 ENDIF 812 C 813 C POTENTIAL TEMPERATURE. 814 IF (IGET(014).GT.0) THEN 815 IF (LVLS(L,IGET(014)).GT.0) THEN 816 !$omp parallel do 817 DO J=JSTA,JEND 818 DO I=1,IM 819 EGRID2(I,J)=SPL(L) 820 ENDDO 821 ENDDO 822 CALL CALPOT2(EGRID2,TSL,EGRID1,IM,JM) 823 CALL E2OUT(014,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 824 ID(1:25)=0 825 CALL OUTPUT(IOUTYP,IGET(014),L,GRID1,IMOUT,JMOUT) 826 ENDIF 827 ENDIF 828 C 829 C RELATIVE HUMIDITY. 830 IF (IGET(017).GT.0) THEN 831 IF (LVLS(L,IGET(017)).GT.0) THEN 832 !$omp parallel do 833 DO J=JSTA,JEND 834 DO I=1,IM 835 EGRID2(I,J)=SPL(L) 836 ENDDO 837 ENDDO 838 CALL CALRH2(EGRID2,TSL,QSL,ICE,EGRID1,IM,JM) Page 13 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 839 CALL E2OUT(017,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 840 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 841 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 842 ID(1:25)=0 843 CALL OUTPUT(IOUTYP,IGET(017),L,GRID1,IMOUT,JMOUT) 844 ENDIF 845 ENDIF 846 C 847 C DEWPOINT TEMPERATURE. 848 IF (IGET(015).GT.0) THEN 849 IF (LVLS(L,IGET(015)).GT.0) THEN 850 !$omp parallel do 851 DO J=JSTA,JEND 852 DO I=1,IM 853 EGRID2(I,J)=SPL(L) 854 ENDDO 855 ENDDO 856 CALL CALDWP2(EGRID2,QSL,EGRID1,TSL) 857 CALL E2OUT(015,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 858 ID(1:25)=0 859 CALL OUTPUT(IOUTYP,IGET(015),L,GRID1,IMOUT,JMOUT) 860 ENDIF 861 ENDIF 862 C 863 C SPECIFIC HUMIDITY. 864 IF (IGET(016).GT.0) THEN 865 IF (LVLS(L,IGET(016)).GT.0) THEN 866 CALL E2OUT(016,000,QSL,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 867 CALL BOUND(GRID1,H1M12,H99999,IMOUT,JMOUT) 868 ID(1:25)=0 869 CALL OUTPUT(IOUTYP,IGET(016),L,GRID1,IMOUT,JMOUT) 870 ENDIF 871 ENDIF 872 C 873 C OMEGA 874 IF (IGET(020).GT.0) THEN 875 IF (LVLS(L,IGET(020)).GT.0) THEN 876 CALL E2OUT(020,000,OSL,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 877 ID(1:25)=0 878 CALL OUTPUT(IOUTYP,IGET(020),L,GRID1,IMOUT,JMOUT) 879 ENDIF 880 ENDIF 881 C 882 C MOISTURE CONVERGENCE. 883 IF (IGET(085).GT.0) THEN 884 IF (LVLS(L,IGET(085)).GT.0) THEN 885 CALL CALMCVG(QSL,USL,VSL,-1,EGRID1) 886 CALL E2OUT(085,000,EGRID1,EGRID2, 887 X GRID1,GRID2,IMOUT,JMOUT) 888 C CONVERT TO DIVERGENCE FOR GRIB UNITS 889 CALL SCLFLD(GRID1,-1.0,IMOUT,JMOUT) 890 ID(1:25)=0 891 CALL OUTPUT(IOUTYP,IGET(085),L,GRID1,IMOUT,JMOUT) 892 ENDIF 893 ENDIF 894 C 895 C U AND/OR V WIND. Page 14 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 896 IF (IGET(018).GT.0.OR.IGET(019).GT.0) THEN 897 IF (LVLS(L,IGET(018)).GT.0.OR.LVLS(L,IGET(019)).GT.0) THEN 898 CALL E2OUT(018,019,USL,VSL,GRID1,GRID2,IMOUT,JMOUT) 899 ID(1:25)=0 900 IF (IGET(018).GT.0)then 901 CALL OUTPUT(IOUTYP,IGET(018),L,GRID1,IMOUT,JMOUT) 902 endif 903 ID(1:25)=0 904 IF (IGET(019).GT.0) 905 X CALL OUTPUT(IOUTYP,IGET(019),L,GRID2,IMOUT,JMOUT) 906 ENDIF 907 ENDIF 908 C 909 C ABSOLUTE VORTICITY. 910 IF (IGET(021).GT.0) THEN 911 IF (LVLS(L,IGET(021)).GT.0) THEN 912 CALL CALVOR(USL,VSL,EGRID1) 913 CALL E2OUT(021,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 914 ID(1:25)=0 915 CALL OUTPUT(IOUTYP,IGET(021),L,GRID1,IMOUT,JMOUT) 916 ENDIF 917 ENDIF 918 C 919 C GEOSTROPHIC STREAMFUNCTION. 920 IF (IGET(086).GT.0) THEN 921 IF (LVLS(L,IGET(086)).GT.0) THEN 922 !$omp parallel do 923 DO J=JSTA,JEND 924 DO I=1,IM 925 EGRID2(I,J)=FSL(I,J)*GI 926 ENDDO 927 ENDDO 928 CALL CALSTRM(EGRID2,EGRID1) 929 CALL E2OUT(086,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 930 ID(1:25)=0 931 CALL OUTPUT(IOUTYP,IGET(086),L,GRID1,IMOUT,JMOUT) 932 ENDIF 933 ENDIF 934 C 935 C TURBULENT KINETIC ENERGY. 936 IF (IGET(022).GT.0) THEN 937 IF (LVLS(L,IGET(022)).GT.0) THEN 938 CALL E2OUT(022,000,Q2SL,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 939 ID(1:25)=0 940 CALL OUTPUT(IOUTYP,IGET(022),L,GRID1,IMOUT,JMOUT) 941 ENDIF 942 ENDIF 943 C 944 C TOTAL CLOUD WATER. 945 IF (IGET(153).GT.0) THEN 946 IF (LVLS(L,IGET(153)).GT.0) THEN 947 CALL E2OUT(153,000,QCSL,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 948 CALL BOUND(GRID1,H1M12,H99999,IMOUT,JMOUT) 949 ID(1:25)=0 950 CALL OUTPUT(IOUTYP,IGET(153),L,GRID1,IMOUT,JMOUT) 951 ENDIF 952 ENDIF Page 15 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f 953 C 954 C TOTAL CLOUD ICE 955 IF (IGET(166).GT.0) THEN 956 IF (LVLS(L,IGET(166)).GT.0) THEN 957 CALL E2OUT(166,000,ICE,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 958 CALL BOUND(GRID1,H1M12,H99999,IMOUT,JMOUT) 959 ID(1:25)=0 960 CALL OUTPUT(IOUTYP,IGET(166),L,GRID1,IMOUT,JMOUT) 961 ENDIF 962 ENDIF 963 964 C 965 C END OF MAIN VERTICAL LOOP. 966 C 967 310 CONTINUE 968 IOALL=.TRUE. 969 C 970 C ENDIF FOR IF TEST SEEING IF WE WANT ANY OTHER VARIABLES 971 ENDIF 972 C 973 C END OF ROUTINE. 974 C 975 RETURN 976 END ENTRY POINTS Name eta2p_ Page 16 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 170 Label 207 203,204 180 Label 209 190,191,206 220 Label 478 217,218 225 Label 483 280 Label 579 560,576 281 Label 583 548,549,573 290 Label 626 589,590 310 Label 813 179 970 Label 171 143,144 975 Label 172 130,138 A2 Param 41 R(4) 4 scalar 153,244,288,346,389 A3 Param 41 R(4) 4 scalar 153,244,288,346,389 A4 Param 41 R(4) 4 scalar 153,244,288,346,389 ACCLIQ Scalar 7 R(4) 4 2 91719 COM ACPREC Scalar 7 R(4) 4 2 91719 COM AD05 Param 33 R(4) 4 scalar AETA Scalar 8 R(4) 4 1 45 COM 152,447 AHF Local 316 R(4) 4 scalar 316,417,448,462,476 AHFO Local 318 R(4) 4 scalar 318,419,452,466 AHFQ Local 317 R(4) 4 scalar 317,418,450,463 AHFQ2 Local 321 R(4) 4 scalar 321,422,458,473 AHFQC Local 319 R(4) 4 scalar 319,420,454,467 AHFQI Local 320 R(4) 4 scalar 320,421,456,470 AI Local 237 R(4) 4 scalar 237,240,245,281,284,289,339,342,34 7,382,385,390 AKHS Scalar 3 R(4) 4 2 91719 COM AKMS Scalar 3 R(4) 4 2 91719 COM ALOG Func 178 scalar 178,447,509,511,559,563,607 ALPET1 Local 608 R(4) 4 scalar 608,609 ALPET2 Local 565 R(4) 4 scalar 565,569,572,582 ALPET2X Scalar 101 R(4) 4 2 91719 COM 572,582,609 ALPETL Local 563 R(4) 4 scalar 563,565,578,605,608 ALPETU Local 559 R(4) 4 scalar 559,565,571,578,581,607,608 ALPETUX Scalar 101 R(4) 4 2 91719 COM 571,581,605 ALPINT Scalar 4 R(4) 4 3 4219074 COM 204,221,222,449,451,453,455,457,45 9,476 ALPSL Local 509 R(4) 4 scalar 509,511,512,514 ALPTH Local 178 R(4) 4 scalar 178,512,514 ALSL Scalar 7 R(4) 4 1 39 COM 184,204,476,569,609 AMAX1 Func 307 scalar 307,408,464,468,471,474 AMIN1 Func 151 scalar 151,154,236,245,280,289,338,347,38 1,390 AVBL Scalar 2 CHAR 20 1 250 COM B Local 308 R(4) 4 scalar 308,409,429,448,462,476 BETA Param 37 R(4) 4 scalar 37 BI Local 238 R(4) 4 scalar 238,241,245,282,285,289,340,343,34 7,383,386,390 BOM Local 310 R(4) 4 scalar 310,411,431,452,466 BOUND Subr 687 687,713,794,804 BQ Local 309 R(4) 4 scalar 309,410,430,450,463 BQ2 Local 315 R(4) 4 scalar 315,416,446,458,473 Page 17 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References BQC Local 312 R(4) 4 scalar 312,413,433,454,467 BQC_2 Local 436 R(4) 4 scalar 436,454 BQI Local 313 R(4) 4 scalar 313,414,434,456,470 BQI_2 Local 437 R(4) 4 scalar 437,456 BTG Param 37 R(4) 4 scalar CALDWP2 Subr 702 702 CALMCVG Subr 731 731 CALPOT2 Subr 668 668 CALRH2 Subr 684 684 CALSTRM Subr 774 774 CALVOR Subr 758 758 CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CLDEFI Scalar 6 R(4) 4 2 91719 COM CLDWTR Common 2 17610408 CLIMIT Local 128 R(4) 4 scalar 128,159,169 CM1 Param 23 R(4) 4 scalar CM2 Param 23 R(4) 4 scalar CM3 Param 23 R(4) 4 scalar CMLD Scalar 6 R(4) 4 scalar COM CP Param 34 R(4) 4 scalar 36,40 CPDR Param 36 R(4) 4 scalar CPHI0 Scalar 9 R(4) 4 scalar COM CTLBLK Common 2 96 CUPREC Scalar 7 R(4) 4 2 91719 COM CWM Scalar 3 R(4) 4 3 4127355 COM 149,312,313,413,414,433,434,436,43 7 CWMKL Local 149 R(4) 4 scalar 149,159,169 D00 Param 25 R(4) 4 scalar 159,204,222,316,317,318,319,320,32 1,322,417,418,419,420,421,422,423, 474,540,541 D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar 269,270,371,372 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 231,232,234,314,332,333,335,415,43 8,441 D608 Param 29 R(4) 4 scalar 274,276,376,378 DATSET Scalar 3 CHAR 6 scalar COM DEC Scalar 8 R(4) 4 1 250 COM DETA Scalar 8 R(4) 4 1 45 COM DISLP Scalar 7 R(4) 4 scalar COM DLMD Scalar 6 R(4) 4 scalar COM DP10M Param 35 R(4) 4 scalar DP30 Scalar 6 R(4) 4 scalar COM DPHD Scalar 6 R(4) 4 scalar COM DT Scalar 4 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 E2OUT Subr 644 644,653,669,685,703,712,722,732,74 4,759,775,784,793,803 E2PFL Common 2 8 Page 18 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References EGRID1 Local 85 R(4) 4 2 91719 640,644,668,669,684,685,702,703,73 1,732,758,759,774,775 EGRID2 Local 85 R(4) 4 2 91719 644,653,665,668,669,681,684,685,69 9,702,703,712,722,732,759,771,774, 775,784,793,803 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 ERAD Param 34 R(4) 4 scalar ERLAM0 Scalar 9 R(4) 4 scalar COM ETA Scalar 8 R(4) 4 1 46 COM 558,561,606 ETA2P Subr 1 EXP Func 153 scalar 153,244,288,346,389 EXTRA Common 3 71506268 FAC Local 322 R(4) 4 scalar 322,423,447,462,463,466,467,470,47 3,476 FACT Local 609 R(4) 4 scalar 609,611,613 FIELD Scalar 2 CHAR 20 1 250 COM FIQ Local 158 R(4) 4 scalar 158,159 FIRST Scalar 79 L(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM 510,511 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 475,511,514,516,526,640,771 G Param 32 R(4) 4 scalar 32,75,477,526 GAMMA Param 75 R(4) 4 scalar 75 GI Param 32 R(4) 4 scalar 36,516,640,771 GMIW Local 311 R(4) 4 scalar 311,312,313,412,413,414,432,433,43 4,436,437 GMIW_2 Local 435 R(4) 4 scalar 435 GRID1 Local 86 R(4) 4 2 0 644,646,653,655,669,671,685,686,68 7,689,703,705,712,713,715,722,724, 733,735,737,744,747,759,761,775,77 7,784,786,793,794,796,803,804,806 GRID2 Local 86 R(4) 4 2 0 644,653,669,685,703,712,722,733,74 4,751,759,775,784,793,803 H1 Param 26 R(4) 4 scalar 33,264,265,274,276,366,367,376,378 ,510,687 H100 Param 27 R(4) 4 scalar 686,687 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 307,408,464,468,471,713,794,804 H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar H2 Param 26 R(4) 4 scalar 184,447,476 H4 Param 26 R(4) 4 scalar H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar 713,794,804 HBM2 Scalar 3 R(4) 4 2 91719 COM 146 HH Local 146 R(4) 4 scalar 146,150,151,153 HTM Scalar 5 R(4) 4 3 4127355 COM 146,446,458 Page 19 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References I Local 134 I(4) 4 scalar 134,135,144,145,146,147,148,149,15 2,155,160,162,166,168,169,191,198, 204,205,208,218,221,222,223,224,22 8,229,230,231,232,234,310,311,312, 313,314,324,330,331,332,333,335,41 1,412,413,414,415,429,430,431,432, 433,434,435,436,437,438,440,441,44 6,447,448,449,450,451,452,453,455, 457,458,459,462,463,464,465,466,46 7,468,469,470,471,472,473,474,475, 476,477,492,493,508,509,510,511,51 4,516,525,526,539,540,541,549,556, 558,561,570,571,572,580,581,582,59 0,591,592,593,604,605,606,609,610, 611,612,613,615,616,639,640,664,66 5,680,681,698,699,770,771 ICE Local 83 R(4) 4 2 91719 470,471,472,684,803 ICNT Scalar 7 I(4) 4 1 1024 COM ID Scalar 11 I(4) 4 1 25 COM 645,654,670,688,704,714,723,736,74 5,749,760,776,785,795,805 IDAT Scalar 3 I(4) 4 1 3 COM IDENT Scalar 10 I(4) 4 1 250 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 IFILV Scalar 11 I(4) 4 1 250 COM IGET Scalar 8 I(4) 4 1 250 COM 117,118,119,120,121,122,123,488,48 9,505,534,535,635,636,646,651,652, 655,660,661,671,676,677,689,694,69 5,705,710,711,715,720,721,724,729, 730,737,742,743,746,747,750,751,75 6,757,761,766,767,777,782,783,786, 791,792,796,801,802,806 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 IL Local 203 I(4) 4 scalar 203,204,205,560,561,570,576 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,14,81,82,83,85, 101,134,144,191,218,492,508,525,53 9,549,590,639,664,668,680,684,698, 770 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 86,644,646,653,655,669,671,685,686 ,687,689,703,705,712,713,715,722,7 24,733,735,737,744,747,751,759,761 ,775,777,784,786,793,794,796,803,8 04,806 IMT Param 4 I(4) 4 scalar IMX Param 3 I(4) 4 scalar INPES Param 6 I(4) 4 scalar 11,14 Page 20 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References IOALL Scalar 80 L(4) 4 scalar COM 814 IOOMG Scalar 80 L(4) 4 scalar COM IOUT Scalar 3 I(4) 4 scalar COM IOUTYP Scalar 8 I(4) 4 scalar COM 646,655,671,689,705,715,724,737,74 7,751,761,777,786,796,806 IQ Scalar 9 I(4) 4 1 250 COM IS Scalar 9 I(4) 4 1 250 COM ISHDE Scalar 4 I(4) 4 1 61 COM ISMFUL Scalar 9 I(4) 4 1 250 COM ISMOUT Scalar 10 I(4) 4 1 250 COM ISMSTG Scalar 9 I(4) 4 1 250 COM ITAIL Param 11 I(4) 4 scalar IUNIT Scalar 10 I(4) 4 scalar COM 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 IW Local 83 R(4) 4 3 4127355 135,160,162,166,168,169,234,311,33 5,412,432,435 IWL Local 83 R(4) 4 scalar 335,354,397 IWU Local 83 R(4) 4 scalar 234,252,296 IXM Scalar 3 I(4) 4 scalar COM IYM Scalar 3 I(4) 4 scalar COM J Local 133 I(4) 4 scalar 133,135,143,145,146,147,148,149,15 2,155,160,162,166,168,169,190,198, 204,205,208,217,221,222,223,224,22 8,229,230,231,232,234,310,311,312, 313,314,324,330,331,332,333,335,41 1,412,413,414,415,429,430,431,432, 433,434,435,436,437,438,440,441,44 6,447,448,449,450,451,452,453,455, 457,458,459,462,463,464,465,466,46 7,468,469,470,471,472,473,474,475, 476,477,491,493,507,509,510,511,51 4,516,524,526,538,540,541,548,556, 558,561,570,571,572,580,581,582,58 9,591,592,593,604,605,606,609,610, 611,612,613,615,616,638,640,663,66 5,679,681,697,699,769,771 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 133,143,190,217,491,507,524,538,54 8,589,638,663,679,697,769 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 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 JIMA Common 101 1100628 JM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,12,15,81,82,83,85, 101,668,684 JM2 Param 6 I(4) 4 scalar JMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 86,644,646,653,655,669,671,685,686 ,687,689,703,705,712,713,715,722,7 24,733,735,737,744,747,751,759,761 Page 21 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References ,775,777,784,786,793,794,796,803,8 04,806 JMT Param 4 I(4) 4 scalar JMX Param 3 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 133,143,190,217,491,507,524,538,54 8,589,638,663,679,697,769 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM JSTA_M2 Scalar 7 I(4) 4 scalar COM JTAIL Param 12 I(4) 4 scalar KGTYPE Scalar 8 I(4) 4 scalar COM 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 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 KUNIT Scalar 10 I(4) 4 scalar COM L Local 130 I(4) 4 scalar 130,131,135,146,147,148,149,152,15 5,160,162,166,168,169,179,184,204, 275,287,377,388,475,476,487,501,56 9,609,636,646,652,655,661,665,671, 677,681,689,695,699,705,711,715,72 1,724,730,737,743,747,751,757,761, 767,777,783,786,792,796,802,806 LB Param 9 I(4) 4 scalar LC Scalar 3 I(4) 4 2 91719 COM 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 2,3,5,8,9,19,83,130,145,192,219,51 4 LM1 Param 8 I(4) 4 scalar LMA Local 192 I(4) 4 scalar 192,194,219,220,330,331,332,333,33 5,411 LMAP1 Local 194 I(4) 4 scalar 194,203,208,220,324 LMB Local 556 I(4) 4 scalar 556,560,576,580 LMH Scalar 5 I(4) 4 2 91719 COM 145,412,413,414,415 LML Local 145 I(4) 4 scalar 145,155 LMV Scalar 5 I(4) 4 2 91719 COM 556,591,592,593 LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4,8 LSL Scalar 3 I(4) 4 scalar COM 179 LSM Param 3 I(4) 4 scalar 7 LUNIT Scalar 10 I(4) 4 scalar COM LVLS Scalar 10 I(4) 4 2 11250 COM 636,652,661,677,695,711,721,730,74 3,757,767,783,792,802 Page 22 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References MAPOT Common 2 1448 MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM MUNIT Scalar 10 I(4) 4 scalar COM MXFLD Param 3 I(4) 4 scalar 2,8,9,10,11 MXLVL Param 3 I(4) 4 scalar 10 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 NFLD Scalar 8 I(4) 4 scalar COM NINC Param 7 I(4) 4 scalar NINT Func 487 scalar 487,501 NL1X Scalar 101 I(4) 4 2 91719 COM 205,208,221,222,223,224,228,324,42 9,430,431,432,433,434,435,436,437, 438,440,441,446,447,448,449,450,45 1,452,453,455,457,458,459,476,477, 570,580,591,604,606,610,611,612,61 3,615,616 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 NUNIT Scalar 10 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 79 L(4) 4 scalar COM OMGA Scalar 2 R(4) 4 3 4127355 COM 310,411,431,452 OMGAOT Common 2 16509420 ONEPS Param 39 R(4) 4 scalar OSL Local 81 R(4) 4 2 91719 466,722 OUTPUT Subr 646 646,655,671,689,705,715,724,737,74 7,751,761,777,786,796,806 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 511 PDSL Scalar 6 R(4) 4 2 91719 COM 152,198,447 PDVP1 Scalar 6 R(4) 4 2 91719 COM 558,561,606 PETAL Local 561 R(4) 4 scalar 561,563,577 PETAU Local 558 R(4) 4 scalar 558,559,577,606,607 PFSL Scalar 7 R(4) 4 2 91719 COM PINT Scalar 4 R(4) 4 3 4219074 COM 224,229,330 PL Local 330 R(4) 4 scalar 330,345,377 PNL1 Local 224 R(4) 4 scalar 224,475 PP Local 152 R(4) 4 scalar 152,153 PPDSL Local 198 R(4) 4 scalar 198 Page 23 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References PQ0 Param 41 R(4) 4 scalar 153,243,287,345,388 PREC Scalar 7 R(4) 4 2 91719 COM PRSI Local 82 R(4) 4 2 91719 PSLP Scalar 6 R(4) 4 2 91719 COM 509 PT Scalar 5 R(4) 4 scalar COM 152,447,511,558,561,606 PTSL Scalar 7 R(4) 4 2 91719 COM PU Local 229 R(4) 4 scalar 229,243,275 PVRBLS Common 2 23480064 Q Scalar 5 R(4) 4 3 4127355 COM 148,232,333,430,450 Q2 Scalar 9 R(4) 4 3 4127355 COM 314,415,438,441 Q2A Local 314 R(4) 4 scalar 314,315,415,416,441,443,458 Q2B Local 438 R(4) 4 scalar 438,443,446 Q2SL Local 82 R(4) 4 2 91719 473,474,784 QABV Local 306 R(4) 4 scalar 306,307,309 QBLO Local 407 R(4) 4 scalar 407,408,410 QCSL Local 82 R(4) 4 2 91719 467,468,469,793 QI Local 154 R(4) 4 scalar 154,158,245,248,253,289,292,297,34 7,350,355,390,393,398 QINT Local 246 R(4) 4 scalar 246,250,255,290,294,299,348,352,35 7,391,395,400 QKL Local 148 R(4) 4 scalar 148,158 QL Local 333 R(4) 4 scalar 333,364,368,373,376,378 QS Scalar 3 R(4) 4 2 91719 COM QSAT Local 248 R(4) 4 scalar 248,250,253,255,260,262,266,271,29 2,294,297,299,304,306,350,352,355, 357,362,364,368,373,393,395,398,40 0,405,407 QSL Scalar 9 R(4) 4 2 91719 COM 463,464,465,684,702,712,731 QU Local 232 R(4) 4 scalar 232,262,266,271,274,276 QW Local 153 R(4) 4 scalar 153,154,243,245,246,260,287,289,29 0,304,345,347,348,362,388,390,391, 405 QWBS Scalar 5 R(4) 4 2 91719 COM QZ0 Scalar 4 R(4) 4 2 91719 COM R Scalar 5 R(4) 4 scalar COM 476,514 RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40,75 RDOCP Param 40 R(4) 4 scalar RES Scalar 3 R(4) 4 2 91719 COM RESTRT Scalar 79 L(4) 4 scalar COM RF Scalar 5 R(4) 4 2 91719 COM RG Param 36 R(4) 4 scalar RGAMOG Param 75 R(4) 4 scalar 275,377 RHL Local 364 R(4) 4 scalar 364,366,367,368,371,372,373,407 RHU Local 262 R(4) 4 scalar 262,264,265,266,269,270,271,306 RITE2 Scalar 5 L(4) 4 scalar COM RITEHD Scalar 5 L(4) 4 scalar COM ROG Param 32 R(4) 4 scalar RQSTFLD Common 7 64130 RTD Param 34 R(4) 4 scalar RUN Scalar 79 L(4) 4 scalar COM SBD Scalar 6 R(4) 4 scalar COM SCLFLD Subr 686 686,735 SI Scalar 6 R(4) 4 2 91719 COM SICE Scalar 3 R(4) 4 2 91719 COM Page 24 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References SIGMA Scalar 79 L(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 SNO Scalar 6 R(4) 4 2 91719 COM SPHI0 Scalar 9 R(4) 4 scalar COM SPL Scalar 7 R(4) 4 1 39 COM 275,287,377,388,475,487,501,665,68 1,699 SQRT Func 565 scalar 565,608 SR Scalar 3 R(4) 4 2 91719 COM STDRD Scalar 79 L(4) 4 scalar COM SVALUE Scalar 8 R(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM 147,231,332,429,448,514 T500 Scalar 6 R(4) 4 2 91719 COM 493 TABV Local 276 R(4) 4 scalar 276,279,288,308 TBLO Local 378 R(4) 4 scalar 378,380,389,409 TFRZ Param 39 R(4) 4 scalar THS Scalar 3 R(4) 4 2 91719 COM THZ0 Scalar 4 R(4) 4 2 91719 COM TKL Local 147 R(4) 4 scalar 147,150,153 TL Local 332 R(4) 4 scalar 332,337,346,376 TLM0D Scalar 6 R(4) 4 scalar COM TMT0 Local 150 R(4) 4 scalar 150,151,154,157,166,167,235,236,23 9,245,247,249,279,280,283,289,291, 293,337,338,341,347,349,351,380,38 1,384,390,392,394 TMT15 Local 151 R(4) 4 scalar 151,236,246,280,290,338,348,381,39 1 TPH0D Scalar 6 R(4) 4 scalar COM TRF Local 184 R(4) 4 scalar 184,221,462,463,466,467,470,473 TSHDE Scalar 8 R(4) 4 1 61 COM TSL Scalar 8 R(4) 4 2 91719 COM 462,493,653,668,684,702 TSPH Scalar 5 R(4) 4 scalar COM TTV Scalar 8 R(4) 4 2 91719 COM TU Local 231 R(4) 4 scalar 231,235,244,274 TVRABV Local 275 R(4) 4 scalar 275,276 TVRBLO Local 377 R(4) 4 scalar 377,378 TVRL Local 376 R(4) 4 scalar 376,377 TVRU Local 274 R(4) 4 scalar 274,275 TWBS Scalar 5 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM 592,610,611,615 U00 Scalar 3 R(4) 4 2 91719 COM 155 U00KL Local 155 R(4) 4 scalar 155,158 UL Scalar 3 R(4) 4 1 90 COM 155 USL Local 81 R(4) 4 2 91719 540,592,610,615,731,744,758 UTIM Local 127 R(4) 4 scalar 127,155 UZ0 Scalar 4 R(4) 4 2 91719 COM V Scalar 5 R(4) 4 3 4127355 COM 593,612,613,616 VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VRBLS Common 2 67138308 VSL Local 81 R(4) 4 2 91719 541,593,612,616,731,744,758 VTM Scalar 5 R(4) 4 3 4127355 COM VZ0 Scalar 4 R(4) 4 2 91719 COM Page 25 Source Listing ETA2P 2014-12-17 20:47 Symbol Table ETA2P.f Name Object Declared Type Bytes Dimen Elements Attributes References WBD Scalar 6 R(4) 4 scalar COM X1P Scalar 6 R(4) 4 scalar COM Y1P Scalar 6 R(4) 4 scalar COM Z0 Scalar 3 R(4) 4 2 91719 COM Z0SLP Scalar 7 R(4) 4 scalar COM Z1000 Scalar 9 R(4) 4 2 91719 COM 516,526 ZINT Scalar 4 R(4) 4 3 4219074 COM 230,331,477 ZL Local 331 R(4) 4 scalar 331 ZU Local 230 R(4) 4 scalar 230 Page 26 Source Listing ETA2P 2014-12-17 20:47 Subprograms/Common Blocks ETA2P.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CLDWTR Common 2 17610408 CTLBLK Common 2 96 E2PFL Common 2 8 ETA2P Subr 1 EXTRA Common 3 71506268 JIMA Common 101 1100628 LOOPS Common 2 748952 MAPOT Common 2 1448 MASKS Common 2 34853220 OMGAOT Common 2 16509420 PARALLEL Common 9 8240 PVRBLS Common 2 23480064 RQSTFLD Common 7 64130 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 Page 27 Source Listing ETA2P 2014-12-17 20:47 ETA2P.f -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 -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 : ETA2P.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100