Page 1 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 1 SUBROUTINE ETAFLD2(IMOUT,JMOUT) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: ETAFLD2 SLP AND ETA LEVEL POSTING 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 92-12-21 6 C 7 C ABSTRACT: 8 C THIS ROUTINE DOES SEVERAL THINGS. IT IS THE FIRST 9 C ROUTINE CALLED BY POST PROCESSOR SUBROUTINE PROCESS 10 C WHICH SETS THE ORDER IN WHICH FIELDS ARE POSTED. THE 11 C INTERNAL POST VERSION OF ETAFLD2 FIRST LOADS FUNDAMENTAL 12 C MODEL VARIABLES (T,Q,U,V) INTO POST-ONLY ARRAYS. 13 C NEGATIVE SPECIFIC HUMIDITY IS CLIPPED. NEXT PRESSURE AND 14 C LOG PRESSURE ARE COMPUTED AT ALL MASS GRID POINTS. THE 15 C MODEL ITSELF COMPUTES THE MESINGER SEA LEVEL PRESSURE 16 C PRIOR TO POSTING THE PROFILE DATA. THIS ROUTINE WILL 17 C COMPUTE THE STANDARD NMC SEA LEVEL PRESSURE IF THIS OPTION 18 C IS ACTIVATED. BY EITHER METHOD WE COMPUTE BELOW SURFACE 19 C TEMPERATURES. AFTER COMPUTING OMEGA ON THE ETA LEVELS, 20 C WE SET BELOW SURFACE FIELDS ON ETA LEVELS. USING THE 21 C HYDROSTATIC EQUATION WE COMPUTE THE HEIGHT AT ETA LAYER 22 C INTERFACES. FINALLY WE COMPUTE/POST REQUESTED FIELDS ON 23 C ETA LAYERS. 24 C 25 C . 26 C 27 C PROGRAM HISTORY LOG: 28 C 92-12-21 RUSS TREADON 29 C 93-09-01 RUSS TREADON - ADDED ADDITIONAL OUTPUT FIELDS. 30 C 96-03-20 MIKE BALDWIN - ADDED CLOUD TOP TEMPS, CHANGE CLOUD WATER 31 C TO CONTAIN WATER ONLY 32 C 97-04-29 GEOFF MANIKIN - MOVED CLOUD TOP TEMPS TO CLDRAD 33 C 98-06-01 T BLACK - CONVERSION FROM 1-D TO 2-D 34 C 98-07-20 MIKE BALDWIN - REMOVED LABL84 35 C 98-08-18 T BLACK - REMOVED EXCESS SPACE IN EXTRA.com 36 C 00-01-04 JIM TUCCILLO - MPI VERSION 37 C 03-01-22 H CHUANG & M EK - ADDED PBL HEIGHT CALCULATION 38 C 39 C USAGE: CALL ETAFLD2(IMOUT,JMOUT) 40 C INPUT ARGUMENT LIST: 41 C IMOUT - FIRST DIMENSION OF OUTPUT GRID. 42 C JMOUT - SECOND DIMENSION OF OUTPUT GRID. 43 C 44 C OUTPUT ARGUMENT LIST: 45 C NONE 46 C 47 C OUTPUT FILES: 48 C NONE 49 C 50 C SUBPROGRAMS CALLED: 51 C UTILITIES: 52 C BOUND - BOUND ARRAY ELEMENTS BETWEEN LOWER AND UPPER LIMITS. 53 C E2OUT - INTERPOLATE/SMOOTH E-GRID TO OUTPUT GRID. 54 C SCLFLD - SCALE ARRAY ELEMENTS BY SCALAR CONSTANT. 55 C OUTPUT - POST DATA TO OUTPUT GRID IN SPECIFIED FORMAT. 56 C NGMSLP2 - COMPUTE SLP USING STANDARD NMC REDUCTION METHOD. 57 C BLOSFC2 - SET BELOW SURFACE ETA LEVEL DATA. Page 2 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 58 C NETAL - EXTRACT DATA ON CONSTANT ETA LAYER OR CONSTANT 59 C ATMOSPHERIC ETA LAYER. 60 C CALPOT2 - COMPUTE POTENTIAL TEMPERATURE. 61 C CALRH2 - COMPUTE RELATIVE HUMIDITY. 62 C CALDWP2 - COMPUTE DEWPOINT TEMPERATURE. 63 C CALMCVG - COMPUTE MOISTURE CONVERGENCE. 64 C CALVOR - COMPUTE ABSOLUTE VORTICITY. 65 C CALSTRM - COMPUTE GEOSTROPHIC STREAMFUNCTION. 66 C LIBRARY: 67 C COMMON - VRBLS 68 C PVRBLS 69 C EXTRA 70 C MASKS 71 C MAPOT 72 C DYNAMD 73 C OMGAOT 74 C RQSTFLD 75 C CTLBLK 76 C LOOPS 77 C ACMCLH 78 C ACMRDL 79 C ACMRDS 80 C CLDWTR 81 C IOUNIT 82 C 83 C ATTRIBUTES: 84 C LANGUAGE: FORTRAN 85 C MACHINE : CRAY C-90 86 C$$$ 87 C 88 C 89 C INCLUDE ETA MODEL DIMENSIONS. SET/DERIVE OTHER PARAMETERS. 90 C 91 INCLUDE "parmeta" 109 INCLUDE "parmout" 113 INCLUDE "params" 155 INCLUDE "parm.tbl" 157 parameter(CPBLT=10.) 158 C 159 PARAMETER (RAINCON=1.1787E4) 160 PARAMETER (SNOCON=1.4594E5) 161 PARAMETER (VCON1=1.66476,VCON2=0.55683) 162 C 163 C DECLARE VARIABLES. 164 C 165 LOGICAL RUN,FIRST,RESTRT,SIGMA,OLDRD,STDRD 166 LOGICAL NORTH,NEED(IM,JM) 167 REAL EGRID1(IM,JM),EGRID2(IM,JM),EGRID3(IM,JM),EGRID4(IM,JM) 168 REAL HGT(IM,JM),EL0(IM,JM),FI(IM,JM,2) 169 REAL P1D(IM,JM),T1D(IM,JM),Q1D(IM,JM),EGRID5(IM,JM) 170 REAL GRID1(IMOUT,JMOUT), GRID2(IMOUT,JMOUT) 171 REAL PMID(IM,JM,LM),ZMID(IM,JM,LM),IW(IM,JM,LM) 172 REAL EL(IM,JM,LM),RICHNO(IM,JM,LM) 173 REAL QI(IM,JM),QINT(IM,JM) 174 REAL TT(IM,JM),PPP(IM,JM),QV(IM,JM),QCD(IM,JM),QICE(IM,JM) 175 REAL QRAIN(IM,JM),QSNO(IM,JM),VIS(IM,JM) 176 REAL HPBL(IM,JM) Page 3 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 177 C 178 C INCLUDE REQUIRED COMMONS. 179 INCLUDE "VRBLS.comm" 186 INCLUDE "PVRBLS.comm" 198 INCLUDE "EXTRA.comm" 210 INCLUDE "MASKS.comm" 218 INCLUDE "MAPOT.comm" 231 INCLUDE "DYNAMD.comm" 243 INCLUDE "OMGAOT.comm" 247 INCLUDE "RQSTFLD.comm" 260 INCLUDE "CTLBLK.comm" 275 INCLUDE "ACMCLH.comm" 280 INCLUDE "ACMRDL.comm" 286 INCLUDE "ACMRDS.comm" 292 INCLUDE "CLDWTR.comm" 298 INCLUDE "IOUNIT.comm" 306 INCLUDE "OUTFIL.comm" 312 INCLUDE "LOOPS.comm" 313 C 321 REAL TTND(IM,JM),TRAIN(IM,JM),TCUCN(IM,JM) 322 C 323 C DECLARE EQUIVALENCES. 324 EQUIVALENCE (EGRID1(1,1),P1D(1,1),EL0(1,1)) 325 EQUIVALENCE (EGRID2(1,1),T1D(1,1),HGT(1,1)) 326 EQUIVALENCE (PMID(1,1,1),EL(1,1,1)) 327 EQUIVALENCE (ZMID(1,1,1),RICHNO(1,1,1)) 328 C 329 C***************************************************************************** 330 C START SUBROUTINE ETAFLD. 331 C 332 print *, " CWM start of ETAFLD" 333 do l=1,lm 334 do j=jm,1,-1 335 ! write(*,"(10E12.7)") (CWM(I,J,L),I=1,IM) 336 end do 337 end do 338 print *, " CWM end of ETAFLD" 339 C 340 C SET UP UTIM FOR THIS TIME STEP 341 C 342 CTIM1=0. 343 CTIM2=24.*3600. 344 CTIM =NTSD*DT 345 IF(CTIM.LT.CTIM1)THEN 346 UTIM=0. 347 ELSE 348 IF(CTIM.LE.CTIM2)THEN 349 UTIM=(CTIM-CTIM1)/(CTIM2-CTIM1) 350 ELSE 351 UTIM=1. 352 ENDIF 353 ENDIF 354 C SET TOTAL NUMBER OF OUTPUT GRID POINTS. 355 C 356 C FROM THE ETA MODEL WE COMPUTE SEA LEVEL PRESSURE USING 357 C MESINGER'S ALGORITHM. THIS CREATES AN UNDERGROUND 358 C TEMPERATURE FIELD. SUBROUTINE ETA2P MAKES USE OF THIS Page 4 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 359 C UNDERGROUND TEMPERATURE FIELD IN COMPUTING GEOPOTENTIAL 360 C BELOW THE GROUND. IF THE USER WANTS SEA LEVEL PRESSURE 361 C VIA THE SHUELL SCHEME, MAKE IT SO. THIS SCHEME CREATES 362 C ITS OWN UNDERGROUND TEMPERATURES WHICH ETA2P WILL USE. 363 C THE ROUTINE ALSO COMPUTES ITS OWN 1000MB GEOPOTENTIAL. 364 C 365 C OUTPUT SEA LEVEL PRESSURE IF REQUESTED. 366 C FIRST, MESINGER'S SEA LEVEL PRESSURE. 367 IF (IGET(023).GT.0) THEN 368 CALL E2OUT(023,000,PSLP,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 369 ID(1:25) = 0 370 CALL OUTPUT(IOUTYP,IGET(023),LVLS(1,IGET(023)), 371 X GRID1,IMOUT,JMOUT) 372 ENDIF 373 C 374 C SECOND, STANDARD NGM SEA LEVEL PRESSURE. 375 IF (IGET(105).GT.0) THEN 376 CALL NGMSLP2 377 CALL E2OUT(105,000,SLP,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 378 ID(1:25) = 0 379 CALL OUTPUT(IOUTYP,IGET(105),LVLS(1,IGET(105)), 380 X GRID1,IMOUT,JMOUT) 381 ENDIF 382 C 383 C SET BELOW GROUND Q, U, V, AND OMEGA 384 C 385 CALL BLOSFC2 386 C 387 C COMPUTE HEIGHT AT INTERFACES. 388 C SET SURFACE VALUES. 389 !$omp parallel do 390 DO J=JSTA,JEND 391 DO I=1,IM 392 ZINT(I,J,LP1)=FIS(I,J)*GI 393 FI(I,J,1)=FIS(I,J) 394 ENDDO 395 ENDDO 396 C 397 C COMPUTE VALUES FROM THE SURFACE UP. 398 C 399 DO 80 L=LM,1,-1 400 !$omp parallel do 401 DO J=JSTA,JEND 402 DO I=1,IM 403 FI(I,J,2)=HTM(I,J,L)*T(I,J,L)*(Q(I,J,L)*D608+H1)*R* 404 1 (ALPINT(I,J,L+1)-ALPINT(I,J,L))+FI(I,J,1) 405 ZINT(I,J,L)=FI(I,J,2)*GI 406 FI(I,J,1)=FI(I,J,2) 407 ENDDO 408 ENDDO 409 80 CONTINUE 410 C 411 C COMPUTE VALUES FROM THE SURFACE BELOW. 412 KMM=KMNTM(LM) 413 !$omp parallel do 414 !$omp& private(i,iadd,j,k,lftov1,lmap1,ndrow) 415 DO 100 KM=1,KMM Page 5 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 416 K=KMNT(KM,LM) 417 NDROW=K/IMT 418 LFTOV1=MOD(K,IMT) 419 IF(LFTOV1-IM.GT.0)THEN 420 I=K-NDROW*IMT-IM 421 IADD=2 422 ELSEIF(LFTOV1.GT.0)THEN 423 I=K-NDROW*IMT 424 IADD=1 425 ELSEIF(LFTOV1.EQ.0)THEN 426 I=IM-1 427 IADD=0 428 ENDIF 429 J=2*NDROW+IADD 430 C IF ( J .GE. JSTA .AND. J .LE. JEND ) THEN 431 C 432 LMAP1=LMH(I,J)+1 433 DO L=LMAP1,LM 434 ZINT(I,J,L+1)=DFL(L+1)*GI 435 END DO 436 C END IF 437 100 CONTINUE 438 C 439 C COMPUTE MIDLAYER HEIGHTS AND PRESSURES. 440 !$omp parallel do 441 DO 110 L = 1,LM 442 DO J=JSTA,JEND 443 DO I=1,IM 444 PMID(I,J,L)=D50*(PINT(I,J,L+1)+PINT(I,J,L)) 445 ZMID(I,J,L)=D50*(ZINT(I,J,L+1)+ZINT(I,J,L)) 446 ENDDO 447 ENDDO 448 110 CONTINUE 449 C 450 C OUTPUT/CALCULATE PRESSURE, OMEGA, POTENTIAL TEMPERATURE, 451 C DEWPOINT TEMPERATURE, RELATIVE HUMIDITY, AND 452 C ABSOLUTE VORTICITY ON ETA SURFACES. 453 C 454 NREC0=1 455 C 456 IF ( (IGET(001).GT.0).OR.(IGET(077).GT.0).OR. 457 X (IGET(002).GT.0).OR.(IGET(003).GT.0).OR. 458 X (IGET(004).GT.0).OR.(IGET(005).GT.0).OR. 459 X (IGET(006).GT.0).OR.(IGET(083).GT.0).OR. 460 X (IGET(007).GT.0).OR.(IGET(008).GT.0).OR. 461 X (IGET(009).GT.0).OR.(IGET(010).GT.0).OR. 462 X (IGET(084).GT.0).OR.(IGET(011).GT.0).OR. 463 X (IGET(041).GT.0).OR.(IGET(124).GT.0).OR. 464 X (IGET(125).GT.0).OR.(IGET(145).GT.0).OR. 465 X (IGET(078).GT.0).OR.(IGET(079).GT.0).OR. 466 X (IGET(140).GT.0).OR.(IGET(040).GT.0).OR. 467 X (IGET(180).GT.0) ) THEN 468 C 469 470 C IF ANY OF THE CLOUD ARRAYS ARE REQUESTED, COMPUTE ICE/WATER 471 C ALSO NEED THIS IF VISIBILITY OR RH IS REQUESTED 472 Page 6 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 473 IF(IGET(124).GT.0.OR.IGET(125).GT.0.OR.IGET(145).GT.0 474 X .OR.IGET(006).GT.0.OR.IGET(180).GT.0)THEN 475 CLIMIT =1.0E-20 476 477 DO J=JSTA,JEND 478 DO I=1,IM 479 IW(I,J,1)=0. 480 ENDDO 481 ENDDO 482 C 483 DO 125 L=2,LM 484 DO J=JSTA,JEND 485 DO I=1,IM 486 LML=LM-LMH(I,J) 487 HH=HTM(I,J,L)*HBM2(I,J) 488 TKL=T(I,J,L) 489 QKL=Q(I,J,L) 490 CWMKL=CWM(I,J,L) 491 TMT0=(TKL-273.15)*HH 492 TMT15=AMIN1(TMT0,-15.)*HH 493 PP=PDSL(I,J)*AETA(L)+PT 494 QW=HH*PQ0/PP*EXP(HH*A2*(TKL-A3)/(TKL-A4)) 495 QI(I,J)=QW*(1.+0.01*AMIN1(TMT0,0.)) 496 C 497 U00KL=U00(I,J)+UL(L+LML)*(0.95-U00(I,J))*UTIM 498 IF(TMT0.LT.-15.0)THEN 499 FIQ=QKL-U00KL*QI(I,J) 500 IF(FIQ.GT.D00.OR.CWMKL.GT.CLIMIT) THEN 501 IW(I,J,L)=1. 502 ELSE 503 IW(I,J,L)=0. 504 ENDIF 505 ENDIF 506 IF(TMT0.GE.0.0)IW(I,J,L)=0. 507 IF(TMT0.LT.0.0.AND.TMT0.GE.-15.0)THEN 508 IW(I,J,L)=0. 509 IF(IW(I,J,L-1).EQ.1.0.AND.CWMKL.GT.CLIMIT)IW(I,J,L)=1. 510 ENDIF 511 ENDDO 512 ENDDO 513 125 CONTINUE 514 ENDIF 515 C 516 DO 190 L=1,LM 517 C 518 C PRESSURE ON ETA SURFACES. 519 IF (IGET(001).GT.0) THEN 520 IF (LVLS(L,IGET(001)).GT.0) THEN 521 ITYPE = LVLS(L,IGET(001)) 522 CALL NETAL(PMID,ITYPE,L,LMH,EGRID1) 523 CALL E2OUT(001,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 524 ID(1:25) = 0 525 LYR = L 526 IF (ITYPE.GT.1) ID(9) = 109 527 CALL OUTPUT(IOUTYP,IGET(001),LYR,GRID1,IMOUT,JMOUT) 528 ENDIF 529 ENDIF Page 7 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 530 C 531 C HEIGHTS ON ETA SURFACES. 532 IF (IGET(077).GT.0) THEN 533 IF (LVLS(L,IGET(077)).GT.0) THEN 534 ITYPE = LVLS(L,IGET(077)) 535 CALL NETAL(ZMID,ITYPE,L,LMH,EGRID1) 536 CALL E2OUT(077,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 537 ID(1:25) = 0 538 LYR = L 539 IF (ITYPE.GT.1) ID(9) = 109 540 CALL OUTPUT(IOUTYP,IGET(077),LYR,GRID1,IMOUT,JMOUT) 541 ENDIF 542 ENDIF 543 C 544 C TEMPERATURE ON ETA SURFACES. 545 IF (IGET(002).GT.0) THEN 546 IF (LVLS(L,IGET(002)).GT.0) THEN 547 ITYPE = LVLS(L,IGET(002)) 548 CALL NETAL(T,ITYPE,L,LMH,EGRID1) 549 CALL E2OUT(002,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 550 ID(1:25) = 0 551 LYR = L 552 IF (ITYPE.GT.1) ID(9) = 109 553 CALL OUTPUT(IOUTYP,IGET(002),LYR,GRID1,IMOUT,JMOUT) 554 ENDIF 555 ENDIF 556 C 557 C POTENTIAL TEMPERATURE ON ETA SURFACES. 558 IF (IGET(003).GT.0) THEN 559 IF (LVLS(L,IGET(003)).GT.0) THEN 560 ITYPE = LVLS(L,IGET(003)) 561 CALL NETAL(PMID,ITYPE,L,LMH,P1D) 562 CALL NETAL(T,ITYPE,L,LMH,T1D) 563 CALL CALPOT2(P1D,T1D,EGRID3,IM,JM) 564 CALL E2OUT(003,000,EGRID3,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 565 ID(1:25) = 0 566 LYR = L 567 IF (ITYPE.GT.1) ID(9) = 109 568 CALL OUTPUT(IOUTYP,IGET(003),LYR,GRID1,IMOUT,JMOUT) 569 ENDIF 570 ENDIF 571 C 572 C RELATIVE HUMIDITY ON ETA SURFACES. 573 IF (IGET(006).GT.0) THEN 574 IF (LVLS(L,IGET(006)).GT.0) THEN 575 ITYPE = LVLS(L,IGET(006)) 576 CALL NETAL(PMID,ITYPE,L,LMH,P1D) 577 CALL NETAL(T,ITYPE,L,LMH,T1D) 578 CALL NETAL(Q,ITYPE,L,LMH,Q1D) 579 CALL NETAL(IW,ITYPE,L,LMH,EGRID3) 580 CALL CALRH2(P1D,T1D,Q1D,EGRID3,EGRID4,IM,JM) 581 CALL E2OUT(006,000,EGRID4,EGRID2, 582 X GRID1,GRID2,IMOUT,JMOUT) 583 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 584 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 585 ID(1:25) = 0 586 LYR = L Page 8 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 587 IF (ITYPE.GT.1) ID(9) = 109 588 CALL OUTPUT(IOUTYP,IGET(006),LYR,GRID1,IMOUT,JMOUT) 589 ENDIF 590 ENDIF 591 C 592 C DEWPOINT ON ETA SURFACES. 593 IF (IGET(004).GT.0) THEN 594 IF (LVLS(L,IGET(004)).GT.0) THEN 595 ITYPE = LVLS(L,IGET(004)) 596 CALL NETAL(PMID,ITYPE,L,LMH,P1D) 597 CALL NETAL(Q,ITYPE,L,LMH,Q1D) 598 CALL NETAL(T,ITYPE,L,LMH,T1D) 599 CALL CALDWP2(P1D,Q1D,EGRID3,T1D) 600 CALL E2OUT(004,000,EGRID3,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 601 ID(1:25) = 0 602 LYR = L 603 IF (ITYPE.GT.1) ID(9) = 109 604 CALL OUTPUT(IOUTYP,IGET(004),LYR,GRID1,IMOUT,JMOUT) 605 ENDIF 606 ENDIF 607 C 608 C SPECIFIC HUMIDITY ON ETA SURFACES. 609 IF (IGET(005).GT.0) THEN 610 IF (LVLS(L,IGET(005)).GT.0) THEN 611 ITYPE = LVLS(L,IGET(005)) 612 CALL NETAL(Q,ITYPE,L,LMH,EGRID1) 613 CALL E2OUT(005,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 614 CALL BOUND(GRID1,H1M12,H99999,IMOUT,JMOUT) 615 ID(1:25) = 0 616 LYR = L 617 IF (ITYPE.GT.1) ID(9) = 109 618 CALL OUTPUT(IOUTYP,IGET(005),LYR,GRID1,IMOUT,JMOUT) 619 ENDIF 620 ENDIF 621 C 622 C MOISTURE CONVERGENCE ON ETA SURFACES. 623 IF (IGET(083).GT.0) THEN 624 IF (LVLS(L,IGET(083)).GT.0) THEN 625 ITYPE = LVLS(L,IGET(083)) 626 CALL NETAL(Q,ITYPE,L,LMH,Q1D) 627 CALL NETAL(U,ITYPE,L,LMV,EGRID1) 628 CALL NETAL(V,ITYPE,L,LMV,EGRID2) 629 CALL CALMCVG(Q1D,EGRID1,EGRID2,-1,EGRID3) 630 CALL E2OUT(083,000,EGRID3,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 631 ID(1:25) = 0 632 LYR = L 633 IF (ITYPE.GT.1) ID(9) = 109 634 CALL OUTPUT(IOUTYP,IGET(083),LYR,GRID1,IMOUT,JMOUT) 635 ENDIF 636 ENDIF 637 C 638 C U AND/OR V WIND ON ETA SURFACES. 639 640 IF (IGET(007).GT.0.OR.IGET(008).GT.0) THEN 641 IF (LVLS(L,IGET(007)).GT.0.OR.LVLS(L,IGET(008)).GT.0) THEN 642 ITYPE = LVLS(L,IGET(007)) 643 CALL NETAL(U,ITYPE,L,LMV,EGRID1) Page 9 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 644 ITYPE = LVLS(L,IGET(008)) 645 CALL NETAL(V,ITYPE,L,LMV,EGRID2) 646 CALL E2OUT(007,008,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 647 ID(1:25) = 0 648 LYR = L 649 IF (ITYPE.GT.1) ID(9) = 109 650 IF (IGET(007).GT.0) 651 X CALL OUTPUT(IOUTYP,IGET(007),LYR,GRID1,IMOUT,JMOUT) 652 ID(1:25) = 0 653 LYR = L 654 IF (ITYPE.GT.1) ID(9) = 109 655 IF (IGET(008).GT.0) 656 X CALL OUTPUT(IOUTYP,IGET(008),LYR,GRID2,IMOUT,JMOUT) 657 ENDIF 658 ENDIF 659 C 660 C OMEGA ON ETA SURFACES. 661 IF (IGET(009).GT.0) THEN 662 IF (LVLS(L,IGET(009)).GT.0) THEN 663 ITYPE = LVLS(L,IGET(009)) 664 CALL NETAL(OMGA,ITYPE,L,LMH,EGRID1) 665 CALL E2OUT(009,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 666 ID(1:25) = 0 667 LYR = L 668 IF (ITYPE.GT.1) ID(9) = 109 669 CALL OUTPUT(IOUTYP,IGET(009),LYR,GRID1,IMOUT,JMOUT) 670 ENDIF 671 ENDIF 672 C 673 C ABSOLUTE VORTICITY ON ETA SURFACES. 674 IF (IGET(010).GT.0) THEN 675 IF (LVLS(L,IGET(010)).GT.0) THEN 676 ITYPE = LVLS(L,IGET(010)) 677 CALL NETAL(U,ITYPE,L,LMV,EGRID1) 678 CALL NETAL(V,ITYPE,L,LMV,EGRID2) 679 CALL CALVOR(EGRID1,EGRID2,EGRID3) 680 CALL E2OUT(010,000,EGRID3,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 681 ID(1:25) = 0 682 LYR = L 683 IF (ITYPE.GT.1) ID(9) = 109 684 CALL OUTPUT(IOUTYP,IGET(010),LYR,GRID1,IMOUT,JMOUT) 685 ENDIF 686 ENDIF 687 C 688 C GEOSTROPHIC STREAMFUNCTION ON ETA SURFACES. 689 IF (IGET(084).GT.0) THEN 690 IF (LVLS(L,IGET(084)).GT.0) THEN 691 ITYPE = LVLS(L,IGET(084)) 692 CALL NETAL(ZMID,ITYPE,L,LMH,EGRID1) 693 CALL CALSTRM(EGRID1,EGRID2) 694 CALL E2OUT(084,000,EGRID2,EGRID3,GRID1,GRID2,IMOUT,JMOUT) 695 ID(1:25) = 0 696 LYR = L 697 IF (ITYPE.GT.1) ID(9) = 109 698 CALL OUTPUT(IOUTYP,IGET(084),LYR,GRID1,IMOUT,JMOUT) 699 ENDIF 700 ENDIF Page 10 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 701 C 702 C TURBULENT KINETIC ENERGY ON ETA SURFACES. 703 IF (IGET(011).GT.0) THEN 704 IF (LVLS(L,IGET(011)).GT.0) THEN 705 ITYPE = LVLS(L,IGET(011)) 706 CALL NETAL(Q2,ITYPE,L,LMH,EGRID1) 707 CALL E2OUT(011,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 708 ID(1:25) = 0 709 LYR = L 710 IF (ITYPE.GT.1) ID(9) = 109 711 CALL OUTPUT(IOUTYP,IGET(011),LYR,GRID1,IMOUT,JMOUT) 712 ENDIF 713 ENDIF 714 C 715 C CLOUD WATER CONTENT 716 IF (IGET(124).GT.0) THEN 717 IF (LVLS(L,IGET(124)).GT.0) THEN 718 DO J=JSTA,JEND 719 DO I=1,IM 720 IF(CWM(I,J,L).LT.0..AND.CWM(I,J,L).GT.-1.E-10) 721 1 CWM(I,J,L)=0. 722 ENDDO 723 ENDDO 724 C 725 ITYPE = LVLS(L,IGET(124)) 726 CALL NETAL(CWM,ITYPE,L,LMH,EGRID1) 727 CALL NETAL(IW,ITYPE,L,LMH,EGRID2) 728 C 729 DO J=JSTA,JEND 730 DO I=1,IM 731 EGRID1(I,J)=EGRID1(I,J)*(1.-EGRID2(I,J)) 732 ENDDO 733 ENDDO 734 C 735 CALL E2OUT(124,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 736 ID(1:25) = 0 737 LYR = L 738 IF (ITYPE.GT.1) ID(9) = 109 739 CALL OUTPUT(IOUTYP,IGET(124),LYR,GRID1,IMOUT,JMOUT) 740 ENDIF 741 ENDIF 742 C 743 C CLOUD ICE CONTENT. 744 IF (IGET(125).GT.0) THEN 745 IF (LVLS(L,IGET(125)).GT.0) THEN 746 ITYPE = LVLS(L,IGET(125)) 747 CALL NETAL(CWM,ITYPE,L,LMH,EGRID1) 748 CALL NETAL(IW,ITYPE,L,LMH,EGRID2) 749 C 750 DO J=JSTA,JEND 751 DO I=1,IM 752 EGRID1(I,J)=EGRID1(I,J)*EGRID2(I,J) 753 ENDDO 754 ENDDO 755 C 756 CALL E2OUT(125,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 757 ID(1:25) = 0 Page 11 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 758 LYR = L 759 IF (ITYPE.GT.1) ID(9) = 109 760 CALL OUTPUT(IOUTYP,IGET(125),LYR,GRID1,IMOUT,JMOUT) 761 ENDIF 762 ENDIF 763 C 764 C CLOUD FRACTION 765 C 766 IF (IGET(145).GT.0) THEN 767 IF (LVLS(L,IGET(145)).GT.0) THEN 768 ITYPE = LVLS(L,IGET(145)) 769 US=H1 770 CCLIMIT=1.0E-3 771 CLIMIT =1.0E-20 772 C 773 IF(L.EQ.1)THEN 774 !$omp parallel do 775 DO J=JSTA,JEND 776 DO I=1,IM 777 EGRID1(I,J)=0. 778 ENDDO 779 ENDDO 780 GO TO 180 781 ENDIF 782 C 783 !$omp parallel do 784 !$omp& private(cwmkl,fiw,hh,lml,pp,qc,qkl,qw,rqkl,tkl,tmt0, 785 !$omp& tmt15,u00kl) 786 DO J=JSTA,JEND 787 DO I=1,IM 788 LML=LM-LMH(I,J) 789 HH=HTM(I,J,L)*HBM2(I,J) 790 TKL=T(I,J,L) 791 QKL=Q(I,J,L) 792 CWMKL=CWM(I,J,L) 793 TMT0=(TKL-273.15)*HH 794 TMT15=AMIN1(TMT0,-15.)*HH 795 PP=PDSL(I,J)*AETA(L)+PT 796 QW =HH*PQ0/PP*EXP(HH*A2*(TKL-A3)/(TKL-A4)) 797 QI(I,J)=QW *(1.+0.01*AMIN1(TMT0,0.)) 798 QINT(I,J)=QW *(1.-0.0004*TMT15*(TMT15+15.)) 799 IF(TMT0.LE.-40.)QINT(I,J)=QI(I,J) 800 C 801 U00KL=U00(I,J)+UL(L+LML)*(0.95-U00(I,J))*UTIM 802 C -----------THE SATUATION SPECIFIC HUMIDITY--------- 803 FIW=IW(I,J,L) 804 QC =(H1-FIW)*QINT(I,J)+FIW*QI(I,J) 805 C -----------THE RELATIVE HUMIDITY------------------- 806 IF(QC.LE.D00) THEN 807 RQKL=D00 808 ELSE 809 RQKL=QKL/QC 810 ENDIF 811 C -----------CLOUD COVER RATIO (EGRID1)-------------- 812 IF(RQKL.GE.0.9999) THEN 813 EGRID1(I,J)=AMIN1(1.0,RQKL) 814 ELSE Page 12 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 815 ARG=-1000.0*CWMKL/(1.-RQKL) 816 ARG=AMAX1(ARG,-12.) 817 EGRID1(I,J)=RQKL*(1.-EXP(ARG)) 818 ENDIF 819 C 820 ENDDO 821 ENDDO 822 180 CONTINUE 823 C 824 CALL E2OUT(145,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 825 ID(1:25) = 0 826 LYR = L 827 IF (ITYPE.GT.1) ID(9) = 109 828 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 829 CALL OUTPUT(IOUTYP,IGET(145),LYR,GRID1,IMOUT,JMOUT) 830 ENDIF 831 ENDIF 832 C 833 C TEMPERATURE TENDENCY DUE TO RADIATIVE FLUX CONVERGENCE 834 IF (IGET(140).GT.0) THEN 835 IF (LVLS(L,IGET(140)).GT.0) THEN 836 NREC1=LM+11+9*(L-1) 837 CALL RDRST2D(TTND,IM,JM,LRSTRT,NREC0,NREC1,1,.FALSE.) 838 NREC0=NREC1+1 839 ITYPE = LVLS(L,IGET(140)) 840 CALL E2OUT(140,000,TTND,EGRID2, 841 X GRID1,GRID2,IMOUT,JMOUT) 842 ID(1:25) = 0 843 LYR = L 844 IF (ITYPE.GT.1) ID(9) = 109 845 CALL OUTPUT(IOUTYP,IGET(140),LYR,GRID1,IMOUT,JMOUT) 846 ENDIF 847 ENDIF 848 C 849 C TEMPERATURE TENDENCY DUE TO SHORT WAVE RADIATION. 850 IF (IGET(040).GT.0) THEN 851 IF (LVLS(L,IGET(040)).GT.0) THEN 852 c NREC1=LM+11+9*(L-1) 853 c CALL RDRST2D(TTND,IM,JM,LRSTRT,NREC0,NREC1,1,.FALSE.) 854 c NREC0=NREC1+1 855 ITYPE = LVLS(L,IGET(040)) 856 c CALL NETAL(TTND,ITYPE,L,LMH,EGRID1) 857 C FILLED WILL -H999 FOR NOW 858 !$omp parallel do 859 DO J=JSTA,JEND 860 DO I=1,IM 861 EGRID1(I,J) = -H999 862 ENDDO 863 ENDDO 864 C 865 CALL E2OUT(040,000,EGRID1,EGRID2, 866 X GRID1,GRID2,IMOUT,JMOUT) 867 ID(1:25) = 0 868 LYR = L 869 IF (ITYPE.GT.1) ID(9) = 109 870 CALL OUTPUT(IOUTYP,IGET(040),LYR,GRID1,IMOUT,JMOUT) 871 ENDIF Page 13 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 872 ENDIF 873 C 874 C TEMPERATURE TENDENCY DUE TO LONG WAVE RADIATION. 875 IF (IGET(041).GT.0) THEN 876 IF (LVLS(L,IGET(041)).GT.0) THEN 877 c NREC1=LM+11+9*(L-1) 878 c CALL RDRST2D(TTND,IM,JM,LRSTRT,NREC0,NREC1,1,.FALSE.) 879 c NREC0=NREC1+1 880 ITYPE = LVLS(L,IGET(041)) 881 c CALL NETAL(TTND,ITYPE,L,LMH,EGRID1) 882 C FILLED WILL -H999 FOR NOW 883 !$omp parallel do 884 DO J=JSTA,JEND 885 DO I=1,IM 886 EGRID1(I,J)=-H999 887 ENDDO 888 ENDDO 889 C 890 CALL E2OUT(041,000,EGRID1,EGRID2, 891 1 GRID1,GRID2,IMOUT,JMOUT) 892 ID(1:25) = 0 893 LYR = L 894 IF (ITYPE.GT.1) ID(9) = 109 895 CALL OUTPUT(IOUTYP,IGET(041),LYR,GRID1,IMOUT,JMOUT) 896 ENDIF 897 ENDIF 898 C 899 C LATENT HEATING FROM GRID SCALE RAIN/EVAP. (TIME AVE) 900 IF (IGET(078).GT.0) THEN 901 IF (LVLS(L,IGET(078)).GT.0) THEN 902 NREC1=LM+13+9*(L-1) 903 CALL RDRST2D(TRAIN,IM,JM,LRSTRT,NREC0,NREC1,1,.FALSE.) 904 NREC0=NREC1+1 905 ITYPE = LVLS(L,IGET(078)) 906 IF(AVRAIN.GT.0.)THEN 907 RRNUM=1./AVRAIN 908 ELSE 909 RRNUM=0. 910 ENDIF 911 !$omp parallel do 912 DO J=JSTA,JEND 913 DO I=1,IM 914 EGRID1(I,J)=TRAIN(I,J)*RRNUM 915 ENDDO 916 ENDDO 917 C 918 CALL E2OUT(078,000,EGRID1,EGRID2, 919 1 GRID1,GRID2,IMOUT,JMOUT) 920 ID(1:25) = 0 921 IFHR = NTSD/TSPH+0.5 922 ITHEAT = INT(THEAT) 923 IFINCR = MOD(IFHR,ITHEAT) 924 ID(19) = IFHR 925 ID(20) = 3 926 IF (IFINCR.EQ.0) THEN 927 ID(18) = IFHR-ITHEAT 928 ELSE Page 14 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 929 ID(18) = IFHR-IFINCR 930 ENDIF 931 IF (ID(18).LT.0) ID(18) = 0 932 LYR = L 933 IF (ITYPE.GT.1) ID(9) = 109 934 CALL OUTPUT(IOUTYP,IGET(078),LYR,GRID1,IMOUT,JMOUT) 935 ENDIF 936 ENDIF 937 C 938 C LATENT HEATING FROM CONVECTION. (TIME AVE) 939 IF (IGET(079).GT.0) THEN 940 IF (LVLS(L,IGET(079)).GT.0) THEN 941 NREC1=LM+14+9*(L-1) 942 CALL RDRST2D(TCUCN,IM,JM,LRSTRT,NREC0,NREC1,1,.FALSE.) 943 NREC0=NREC1+1 944 ITYPE = LVLS(L,IGET(079)) 945 IF(AVCNVC.GT.0.)THEN 946 RRNUM=1./AVCNVC 947 ELSE 948 RRNUM=0. 949 ENDIF 950 !$omp parallel do 951 DO J=JSTA,JEND 952 DO I=1,IM 953 EGRID1(I,J) = TCUCN(I,J)*RRNUM 954 ENDDO 955 ENDDO 956 C 957 CALL E2OUT(079,000,EGRID1,EGRID2, 958 1 GRID1,GRID2,IMOUT,JMOUT) 959 ID(1:25) = 0 960 IFHR = NTSD/TSPH+0.5 961 ITHEAT = INT(THEAT) 962 IFINCR = MOD(IFHR,ITHEAT) 963 ID(19) = IFHR 964 ID(20) = 3 965 IF (IFINCR.EQ.0) THEN 966 ID(18) = IFHR-ITHEAT 967 ELSE 968 ID(18) = IFHR-IFINCR 969 ENDIF 970 IF (ID(18).LT.0) ID(18) = 0 971 LYR = L 972 IF (ITYPE.GT.1) ID(9) = 109 973 CALL OUTPUT(IOUTYP,IGET(079),LYR,GRID1,IMOUT,JMOUT) 974 ENDIF 975 ENDIF 976 C 977 C PROCESS NEXT ETA LEVEL. 978 C 979 190 CONTINUE 980 C 981 C END OF ETA SURFACE OUTPUT BLOCK. 982 C 983 ENDIF 984 C VISIBILITY 985 IF (IGET(180).GT.0) THEN Page 15 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 986 RDTPHS= 1./(NPHS*DT) 987 988 C NEED TO CALCULATE RAIN WATER AND SNOW MIXING RATIOS 989 DO J=JSTA,JEND 990 DO I=1,IM 991 IF (PREC(I,J).EQ.0) THEN 992 QSNO(I,J)=0. 993 QRAIN(I,J)=0. 994 ELSE 995 LLMH=LMH(I,J) 996 SNORATE=SR(I,J)*PREC(I,J)*RDTPHS 997 RAINRATE=(1-SR(I,J))*PREC(I,J)*RDTPHS 998 TERM1=(T(I,J,LM)/PSLP(I,J))**0.4167 999 TERM2=(T(I,J,LLMH)/(PDSL(I,J)*AETA(LMH(I,J))+PT))**0.5833 1000 TERM3=RAINRATE**0.8333 1001 QRAIN(I,J)=RAINCON*TERM1*TERM2*TERM3 1002 TERM4=(T(I,J,LM)/PSLP(I,J))**0.47 1003 TERM5=(T(I,J,LLMH)/(PDSL(I,J)*AETA(LMH(I,J))+PT))**0.53 1004 TERM6=SNORATE**0.94 1005 QSNO(I,J)=SNOCON*TERM4*TERM5*TERM6 1006 ENDIF 1007 LLMH=LMH(I,J) 1008 TT(I,J)=T(I,J,LLMH) 1009 QV(I,J)=Q(I,J,LLMH) 1010 QCD(I,J)=(1-IW(I,J,LLMH))*CWM(I,J,LLMH) 1011 QICE(I,J)=IW(I,J,LLMH)*CWM(I,J,LLMH) 1012 PPP(I,J)=PDSL(I,J)*AETA(LMH(I,J))+PT 1013 ENDDO 1014 ENDDO 1015 c CALL CALVIS(QV,QCD,QR,QICE,QS,TT,PPP,PRSNOW,METH,IICE,VIS) 1016 CALL CALVIS(QV,QCD,QRAIN,QICE,QSNO,TT,PPP,VIS) 1017 CALL E2OUT(180,000,VIS,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1018 ID(1:25) = 0 1019 CALL OUTPUT(IOUTYP,IGET(180),LVLS(1,IGET(180)), 1020 X GRID1,IMOUT,JMOUT) 1021 ENDIF 1022 C 1023 C 1024 C ASYMPTOTIC AND FREE ATMOSPHERE MASTER LENGTH SCALE (EL), PLUS 1025 C GRADIENT RICHARDSON NUMBER. THESE FIELDS ARE PLACED OUTSIDE 1026 C THE ABOVE ETA LAYER LOOP SO THAT WE CAN EQUIVALENCE THE 3-D 1027 C EL ARRAY WITH 3-D ARRAY PMID. THIS IS DONE TO MAKE THE POST 1028 C A BIT SMALLER. 1029 C 1030 IF ( (IGET(111).GT.0) .OR. (IGET(146).GT.0) .OR. 1031 X (IGET(147).GT.0) .OR. (IGET(221).GT.0) ) THEN 1032 C 1033 C COMPUTE ASYMPTOTIC MASTER LENGTH SCALE. 1034 CALL CLMAX(DETA,PDSL,HTM,Q2,ZINT,SM,ZINT(1,1,LP1), 1035 X LMH,IM,JM,LM,LP1, 1036 X EL0,EGRID2,EGRID3,EGRID4,EGRID5) 1037 C 1038 C IF REQUESTED, POST ASYMPTOTIC MASTER LENGTH SCALE. 1039 IF (IGET(147).GT.0) THEN 1040 CALL E2OUT(147,000,EL0,EGRID2, 1041 X GRID1,GRID2,IMOUT,JMOUT) 1042 ID(1:25) = 0 Page 16 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 1043 CALL OUTPUT(IOUTYP,IGET(147),LM,GRID1,IMOUT,JMOUT) 1044 ENDIF 1045 ! H CHUANG AND M EK LIFTED HPBL COMPUTATION FROM ETA MODEL 1046 IF (IGET(221).GT.0) THEN 1047 DO J=JSTA,JEND 1048 DO I=1,IM 1049 LMHK=LMH(I,J) 1050 CALL HPBLCAL(I,J,LM,LMHK,LPBL,HPBL(I,J) 1051 + ,Q2(I,J,1:LM),ZINT(I,J,1:LM+1)) 1052 if(i.eq.10.and.j.eq.10)print*,'Debug:sample HPBLCAL inp' 1053 +, I,J,LM,LMHK,LPBL,HPBL(I,J),(Q2(I,J,L),l=1,lm) 1054 + ,(ZINT(I,J,l),l=1,lm) 1055 if(i.eq.550.and.j.eq.144)print* 1056 +, 'Debug:sample HPBLCAL in ETAFLD2' 1057 +, I,J,LM,LMHK,LPBL,HPBL(I,J),SM(I,J),(Q2(I,J,L),l=1,lm) 1058 + ,(ZINT(I,J,l),l=1,lm+1) 1059 C----------------------------------------------------------------------- 1060 END DO 1061 END DO 1062 CALL E2OUT(221,000,HPBL,EGRID2, 1063 X GRID1,GRID2,IMOUT,JMOUT) 1064 ID(1:25) = 0 1065 CALL OUTPUT(IOUTYP,IGET(221),LVLS(1,IGET(221)), 1066 X GRID1,IMOUT,JMOUT) 1067 ENDIF 1068 1069 C 1070 C IF REQUESTED, POST FREE ATMOSPHERE MASTER LENGTH SCALE 1071 C AND/OR THE GRADIENT RICHARDSON NUMBER. 1072 C 1073 IF ( (IGET(111).GT.0) .OR. (IGET(146).GT.0) ) THEN 1074 C 1075 C COMPUTE FREE ATMOSPHERE MASTER LENGTH SCALE. 1076 !$omp parallel do 1077 DO J=JSTA,JEND 1078 DO I=1,IM 1079 HGT(I,J)=ZINT(I,J,LP1) 1080 ENDDO 1081 ENDDO 1082 !$omp parallel do 1083 DO L=1,LM 1084 DO J=JSTA,JEND 1085 DO I=1,IM 1086 EL(I,J,L)=D00 1087 ENDDO 1088 ENDDO 1089 ENDDO 1090 CALL MIXLEN(ZINT,T,PDSL,AETA,PT,Q2,HGT,HTM,EL0, 1091 X LM,LM1,LP1,IM,JM, 1092 X EL) 1093 C 1094 C COMPUTE GRADIENT RICHARDSON NUMBER IF REQUESTED. 1095 C 1096 IF ( (IGET(111).GT.0) ) CALL CALRCH(EL,RICHNO) 1097 C 1098 C LOOP OVER ETA LAYERS. 1099 DO 200 L = 1,LM Page 17 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f 1100 C 1101 C POST MIXING LENGTH. 1102 C 1103 IF (IGET(146).GT.0) THEN 1104 IF (LVLS(L,IGET(146)).GT.0) THEN 1105 ITYPE = LVLS(L,IGET(146)) 1106 CALL NETAL(EL,ITYPE,L,LMH,EGRID1) 1107 CALL E2OUT(146,000,EGRID1,EGRID2, 1108 X GRID1,GRID2,IMOUT,JMOUT) 1109 ID(1:25) = 0 1110 LYR = L 1111 IF (ITYPE.GT.1) ID(9) = 109 1112 CALL OUTPUT(IOUTYP,IGET(146),LYR,GRID1,IMOUT,JMOUT) 1113 ENDIF 1114 ENDIF 1115 C 1116 C POST GRADIENT RICHARDSON NUMBER. 1117 C 1118 IF (IGET(111).GT.0) THEN 1119 IF (LVLS(L,IGET(111)).GT.0) THEN 1120 ITYPE = LVLS(L,IGET(111)) 1121 CALL NETAL(RICHNO,ITYPE,L,LMH,EGRID1) 1122 CALL E2OUT(111,000,EGRID1,EGRID2, 1123 X GRID1,GRID2,IMOUT,JMOUT) 1124 ID(1:25) = 0 1125 LYR = L 1126 IF (ITYPE.GT.1) ID(9) = 109 1127 CALL OUTPUT(IOUTYP,IGET(111),LYR,GRID1,IMOUT,JMOUT) 1128 ENDIF 1129 ENDIF 1130 200 CONTINUE 1131 C 1132 ENDIF 1133 ENDIF 1134 C 1135 C 1136 C END OF ROUTINE. 1137 C 1138 print *, " CWM start of ETAFLD" 1139 do l=1,lm 1140 do j=jm,1,-1 1141 ! write(*,"(10E12.7)") (CWM(I,J,L),I=1,IM) 1142 end do 1143 end do 1144 print *, " CWM end of ETAFLD" 1145 RETURN 1146 END Page 18 Source Listing ETAFLD2 2014-12-17 20:47 Entry Points ETAFLD2.f ENTRY POINTS Name etafld2_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 100 Label 250 228 110 Label 261 254 125 Label 326 296 180 Label 635 593 190 Label 792 329 200 Label 943 912 80 Label 222 212 A2 Param 41 R(4) 4 scalar 307,609 A3 Param 41 R(4) 4 scalar 307,609 A4 Param 41 R(4) 4 scalar 307,609 ACCLIQ Scalar 7 R(4) 4 2 91719 COM ACMCLH Common 2 24 ACMRDL Common 2 2201268 ACMRDS Common 2 2201268 ACPREC Scalar 7 R(4) 4 2 91719 COM ACUTIM Scalar 3 R(4) 4 scalar COM AD05 Param 33 R(4) 4 scalar AETA Scalar 8 R(4) 4 1 45 COM 306,608,812,816,825,903 AETAD Scalar 4 R(4) 4 1 45 COM AKHS Scalar 3 R(4) 4 2 91719 COM AKMS Scalar 3 R(4) 4 2 91719 COM ALPINT Scalar 4 R(4) 4 3 4219074 COM 217 ALSL Scalar 7 R(4) 4 1 39 COM ALWIN Scalar 4 R(4) 4 2 91719 COM ALWOUT Scalar 4 R(4) 4 2 91719 COM ALWTOA Scalar 4 R(4) 4 2 91719 COM AMAX1 Func 629 scalar 629 AMIN1 Func 305 scalar 305,308,607,610,626 ARATIM Scalar 3 R(4) 4 scalar COM ARDLW Scalar 3 R(4) 4 scalar COM ARDSW Scalar 3 R(4) 4 scalar COM ARG Local 628 R(4) 4 scalar 628,629,630 ASWIN Scalar 4 R(4) 4 2 91719 COM ASWOUT Scalar 4 R(4) 4 2 91719 COM ASWTOA Scalar 4 R(4) 4 2 91719 COM AVBL Scalar 2 CHAR 20 1 250 COM AVCNVC Scalar 3 R(4) 4 scalar COM 758,759 AVRAIN Scalar 3 R(4) 4 scalar COM 719,720 BETA Param 37 R(4) 4 scalar 37 BLOSFC2 Subr 198 198 BOUND Subr 397 397,427 BTG Param 37 R(4) 4 scalar CALDWP2 Subr 412 412 CALMCVG Subr 442 442 Page 19 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References CALPOT2 Subr 376 376 CALRCH Subr 909 909 CALRH2 Subr 393 393 CALSTRM Subr 506 506 CALVIS Subr 829 829 CALVOR Subr 492 492 CAPA Param 28 R(4) 4 scalar 35 CCLIMIT Local 583 R(4) 4 scalar 583 CFT0 Param 33 R(4) 4 scalar CLDEFI Scalar 6 R(4) 4 2 91719 COM CLDWTR Common 2 17610408 CLIMIT Local 288 R(4) 4 scalar 288,313,322,584 CLMAX Subr 847 847 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 CPBLT Param 95 R(4) 4 scalar CPDR Param 36 R(4) 4 scalar CPGFU Scalar 8 R(4) 4 2 91719 COM CPGFV Scalar 3 R(4) 4 scalar COM CPHI0 Scalar 9 R(4) 4 scalar COM CTIM Local 157 R(4) 4 scalar 157,158,161,162 CTIM1 Local 155 R(4) 4 scalar 155,158,162 CTIM2 Local 156 R(4) 4 scalar 156,161,162 CTLBLK Common 2 96 CUPREC Scalar 7 R(4) 4 2 91719 COM CURV Scalar 8 R(4) 4 2 91719 COM CWM Scalar 3 R(4) 4 3 4127355 COM 303,533,534,539,560,605,823,824 CWMKL Local 303 R(4) 4 scalar 303,313,322,605,628 D00 Param 25 R(4) 4 scalar 313,619,620,899 D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar D05 Param 29 R(4) 4 scalar 33 D1 Param 25 R(4) 4 scalar D115 Param 29 R(4) 4 scalar D125 Param 27 R(4) 4 scalar D25 Param 27 R(4) 4 scalar D50 Param 25 R(4) 4 scalar 257,258 D608 Param 29 R(4) 4 scalar 216 DATSET Scalar 3 CHAR 6 scalar COM DDMPU Scalar 9 R(4) 4 2 91719 COM DDMPV Scalar 9 R(4) 4 2 91719 COM DEC Scalar 8 R(4) 4 1 250 COM DETA Scalar 8 R(4) 4 1 45 COM 847 DETAD Scalar 4 R(4) 4 1 45 COM DFL Scalar 4 R(4) 4 1 46 COM 247 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 157,799 Page 20 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References DTR Param 34 R(4) 4 scalar 34 DX Scalar 8 R(4) 4 2 91719 COM DY Scalar 3 R(4) 4 scalar COM DYNAM Common 2 3675964 E2OUT Subr 181 181,190,336,349,362,377,394,413,42 6,443,459,478,493,507,520,548,569, 637,653,678,703,731,770,830,853,87 5,920,935 EF4T Scalar 3 R(4) 4 scalar COM EGRID1 Local 105 R(4) 4 2 91719 335,336,348,349,361,362,425,426,44 0,442,456,459,477,478,490,492,505, 506,519,520,539,544,548,560,565,56 9,590,626,630,637,674,678,699,703, 727,731,766,770,919,920,934,935 EGRID2 Local 105 R(4) 4 2 91719 181,190,336,349,362,377,394,413,42 6,441,442,443,458,459,478,491,492, 493,506,507,520,540,544,548,561,56 5,569,637,653,678,703,731,770,830, 849,853,875,920,935 EGRID3 Local 105 R(4) 4 2 91719 376,377,392,393,412,413,442,443,49 2,493,507,849 EGRID4 Local 105 R(4) 4 2 91719 393,394,849 EGRID5 Local 107 R(4) 4 2 91719 849 EL Local 110 R(4) 4 3 4127355 899,905,909,919 EL0 Local 106 R(4) 4 2 91719 849,853,903 ELDR Param 36 R(4) 4 scalar ELWV Param 36 R(4) 4 scalar 36 EM Scalar 6 R(4) 4 1 760 COM EMT Scalar 6 R(4) 4 1 760 COM EN Scalar 3 R(4) 4 scalar COM ENT Scalar 3 R(4) 4 scalar COM EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar ERLAM0 Scalar 9 R(4) 4 scalar COM ETA Scalar 8 R(4) 4 1 46 COM ETAD Scalar 4 R(4) 4 1 46 COM ETAFLD2 Subr 1 EXP Func 307 scalar 307,609,630 EXTRA Common 3 71506268 F Scalar 9 R(4) 4 2 91719 COM F4D Scalar 3 R(4) 4 scalar COM F4Q Scalar 3 R(4) 4 scalar COM F4Q2 Scalar 4 R(4) 4 1 45 COM FAD Scalar 9 R(4) 4 2 91719 COM FCP Scalar 8 R(4) 4 2 91719 COM FDIV Scalar 9 R(4) 4 2 91719 COM FI Local 106 R(4) 4 3 183438 206,216,217,218,219 FIELD Scalar 2 CHAR 20 1 250 COM FIQ Local 312 R(4) 4 scalar 312,313 FIRST Scalar 103 L(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM 205,206 FIW Local 616 R(4) 4 scalar 616,617 FMD Param 39 R(4) 4 scalar 39 FMW Param 39 R(4) 4 scalar 39 Page 21 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References FSL Scalar 8 R(4) 4 2 91719 COM G Param 32 R(4) 4 scalar 32 GI Param 32 R(4) 4 scalar 36,205,218,247 GRID1 Local 108 R(4) 4 2 0 181,184,190,193,336,340,349,353,36 2,366,377,381,395,396,397,401,413, 417,426,427,431,443,447,459,464,47 8,482,493,497,507,511,520,524,548, 552,569,573,637,641,642,654,658,67 9,683,704,708,732,747,771,786,830, 833,854,856,876,879,921,925,936,94 0 GRID2 Local 108 R(4) 4 2 0 181,190,336,349,362,377,395,413,42 6,443,459,469,478,493,507,520,548, 569,637,654,679,704,732,771,830,85 4,876,921,936 H1 Param 26 R(4) 4 scalar 33,216,397,582,617 H100 Param 27 R(4) 4 scalar 396,397 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 427 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 674,699 H99999 Param 26 R(4) 4 scalar 427 HBM2 Scalar 3 R(4) 4 2 91719 COM 300,602 HGT Local 106 R(4) 4 2 91719 892,903 HH Local 300 R(4) 4 scalar 300,304,305,307,602,606,607,609 HPBL Local 114 R(4) 4 2 91719 863,866,870,875 HPBLCAL Subr 863 863 HTM Scalar 5 R(4) 4 3 4127355 COM 216,300,602,847,903 I Local 204 I(4) 4 scalar 204,205,206,215,216,217,218,219,23 3,236,239,245,247,256,257,258,291, 292,298,299,300,301,302,303,306,30 8,310,312,314,316,319,321,322,532, 533,534,543,544,564,565,589,590,60 0,601,602,603,604,605,608,610,611, 612,614,616,617,626,630,673,674,69 8,699,726,727,765,766,803,804,805, 806,808,809,810,811,812,814,815,81 6,818,820,821,822,823,824,825,861, 862,863,864,865,866,867,868,870,87 1,891,892,898,899 IADD Local 234 I(4) 4 scalar 234,237,240,242 ICNT Scalar 7 I(4) 4 1 1024 COM ID Scalar 11 I(4) 4 1 25 COM 182,191,337,339,350,352,363,365,37 8,380,398,400,414,416,428,430,444, 446,460,462,465,467,479,481,494,49 6,508,510,521,523,549,551,570,572, 638,640,655,657,680,682,705,707,73 3,737,738,740,742,744,746,772,776, Page 22 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References 777,779,781,783,785,831,855,877,92 2,924,937,939 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 IFHR Local 734 I(4) 4 scalar 734,736,737,740,742,773,775,776,77 9,781 IFILV Scalar 11 I(4) 4 1 250 COM IFINCR Local 736 I(4) 4 scalar 736,739,742,775,778,781 IGET Scalar 8 I(4) 4 1 250 COM 180,183,188,192,269,270,271,272,27 3,274,275,276,277,278,279,280,286, 287,332,333,334,340,345,346,347,35 3,358,359,360,366,371,372,373,381, 386,387,388,401,406,407,408,417,42 2,423,424,431,436,437,438,447,453, 454,455,457,463,464,468,469,474,47 5,476,482,487,488,489,497,502,503, 504,511,516,517,518,524,529,530,53 8,552,557,558,559,573,579,580,581, 642,647,648,652,658,663,664,668,68 3,688,689,693,708,713,714,718,747, 752,753,757,786,798,832,843,844,85 2,856,859,878,886,909,916,917,918, 925,931,932,933,940 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,14,104,105,106, 107,109,110,111,112,113,114,134,20 4,215,232,233,239,256,291,298,376, 393,532,543,564,589,600,650,673,69 8,716,726,755,765,803,848,861,891, 898,904 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 108,181,184,190,193,336,340,349,35 3,362,366,377,381,395,396,397,401, 413,417,426,427,431,443,447,459,46 4,469,478,482,493,497,507,511,520, 524,548,552,569,573,637,641,642,65 4,658,679,683,704,708,732,747,771, 786,830,833,854,856,876,879,921,92 5,936,940 IMT Param 4 I(4) 4 scalar 230,231,233,236 IMX Param 3 I(4) 4 scalar INPES Param 6 I(4) 4 scalar 11,14 INT Func 735 scalar 735,774 IOUNIT Common 4 20 IOUT Scalar 3 I(4) 4 scalar COM Page 23 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References IOUTYP Scalar 8 I(4) 4 scalar COM 183,192,340,353,366,381,401,417,43 1,447,464,469,482,497,511,524,552, 573,642,658,683,708,747,786,832,85 6,878,925,940 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 ITAG Scalar 3 I(4) 4 scalar COM ITAIL Param 11 I(4) 4 scalar ITB Param 1 I(4) 4 scalar ITBQ Param 1 I(4) 4 scalar ITHEAT Local 735 I(4) 4 scalar 735,736,740,774,775,779 ITYPE Local 334 I(4) 4 scalar 334,335,339,347,348,352,360,361,36 5,373,374,375,380,388,389,390,391, 392,400,408,409,410,411,416,424,42 5,430,438,439,440,441,446,455,456, 457,458,462,467,476,477,481,489,49 0,491,496,504,505,510,518,519,523, 538,539,540,551,559,560,561,572,58 1,640,652,657,668,682,693,707,718, 746,757,785,918,919,924,933,934,93 9 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 109 R(4) 4 3 4127355 292,314,316,319,321,322,392,540,56 1,616,823,824 IXM Scalar 3 I(4) 4 scalar COM IYM Scalar 3 I(4) 4 scalar COM J Local 147 I(4) 4 scalar 147,203,205,206,214,216,217,218,21 9,242,245,247,255,257,258,290,292, 297,299,300,301,302,303,306,308,31 0,312,314,316,319,321,322,531,533, 534,542,544,563,565,588,590,599,60 1,602,603,604,605,608,610,611,612, 614,616,617,626,630,672,674,697,69 9,725,727,764,766,802,804,805,806, 808,809,810,811,812,814,815,816,81 8,820,821,822,823,824,825,860,862, 863,864,865,866,867,868,870,871,89 0,892,897,899,953 JAM Param 6 I(4) 4 scalar 3,6 JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 203,214,255,290,297,531,542,563,58 8,599,672,697,725,764,802,860,890, 897 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 Page 24 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References JGSTR Param 10 I(4) 4 scalar 15 JM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,12,15,104,105,106, 107,109,110,111,112,113,114,134,14 7,376,393,650,716,755,848,904,953 JM2 Param 6 I(4) 4 scalar JMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 108,181,184,190,193,336,340,349,35 3,362,366,377,381,395,396,397,401, 413,417,426,427,431,443,447,459,46 4,469,478,482,493,497,507,511,520, 524,548,552,569,573,637,641,642,65 4,658,679,683,704,708,732,747,771, 786,830,833,854,856,876,879,921,92 5,936,940 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 203,214,255,290,297,531,542,563,58 8,599,672,697,725,764,802,860,890, 897 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 JTB Param 1 I(4) 4 scalar JTBQ Param 1 I(4) 4 scalar K Local 229 I(4) 4 scalar 229,230,231,233,236 KGTYPE Scalar 8 I(4) 4 scalar COM KM Local 228 I(4) 4 scalar 228,229 KMM Local 225 I(4) 4 scalar 225,228 KMNT Scalar 5 I(4) 4 2 4118670 COM 229 KMNTM Scalar 5 I(4) 4 1 45 COM 225 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 146 I(4) 4 scalar 146,212,216,217,218,246,247,254,25 7,258,296,300,301,302,303,306,310, 314,316,319,321,322,329,333,334,33 5,338,346,347,348,351,359,360,361, 364,372,373,374,375,379,387,388,38 9,390,391,392,399,407,408,409,410, 411,415,423,424,425,429,437,438,43 9,440,441,445,454,455,456,457,458, 461,466,475,476,477,480,488,489,49 0,491,495,503,504,505,509,517,518, 519,522,530,533,534,538,539,540,55 0,558,559,560,561,571,580,581,586, Page 25 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References 602,603,604,605,608,614,616,639,64 8,649,652,656,664,668,681,689,693, 706,714,715,718,745,753,754,757,78 4,866,867,870,871,896,899,912,917, 918,919,923,932,933,934,938,952 LB Param 9 I(4) 4 scalar LC Scalar 3 I(4) 4 2 91719 COM LCNTRL Scalar 5 I(4) 4 scalar COM LFTOV1 Local 231 I(4) 4 scalar 231,232,235,238 LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LIST Scalar 3 I(4) 4 scalar COM LLMH Local 808 I(4) 4 scalar 808,812,816,820,821,822,823,824 LM Param 3 I(4) 4 scalar 2,3,4,5,8,9,19,109,110,146,212,225 ,229,246,254,296,299,329,601,649,7 15,754,811,815,848,856,863,864,866 ,867,870,871,896,904,912,952 LM1 Param 8 I(4) 4 scalar 904 LMAP1 Local 245 I(4) 4 scalar 245,246 LMH Scalar 5 I(4) 4 2 91719 COM 245,299,335,348,361,374,375,389,39 0,391,392,409,410,411,425,439,477, 505,519,539,540,560,561,601,808,81 2,816,820,825,848,862,919,934 LMHK Local 862 I(4) 4 scalar 862,863,866,870 LML Local 299 I(4) 4 scalar 299,310,601,614 LMV Scalar 5 I(4) 4 2 91719 COM 440,441,456,458,490,491 LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4,8,205,847,848,892,904 LPBL Local 863 I(4) 4 scalar 863,866,870 LRSTRT Scalar 3 I(4) 4 scalar COM 650,716,755 LSL Scalar 3 I(4) 4 scalar COM LSM Param 3 I(4) 4 scalar 7 LUNCO Scalar 5 I(4) 4 scalar COM LUNIT Scalar 10 I(4) 4 scalar COM LUNLL Scalar 5 I(4) 4 scalar COM LUNOUT Scalar 5 I(4) 4 scalar COM LVLS Scalar 10 I(4) 4 2 11250 COM 183,192,333,334,346,347,359,360,37 2,373,387,388,407,408,423,424,437, 438,454,455,457,475,476,488,489,50 3,504,517,518,530,538,558,559,580, 581,648,652,664,668,689,693,714,71 8,753,757,832,878,917,918,932,933 LYR Local 338 I(4) 4 scalar 338,340,351,353,364,366,379,381,39 9,401,415,417,429,431,445,447,461, 464,466,469,480,482,495,497,509,51 1,522,524,550,552,571,573,639,642, 656,658,681,683,706,708,745,747,78 4,786,923,925,938,940 MAPOT Common 2 1448 MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM MIXLEN Subr 903 903 MOD Func 231 scalar 231,736,775 MUNIT Scalar 10 I(4) 4 scalar COM MXFLD Param 3 I(4) 4 scalar 2,8,9,10,11 Page 26 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 NDROW Local 230 I(4) 4 scalar 230,233,236,242 NEED Local 104 L(4) 4 2 91719 NETAL Subr 335 335,348,361,374,375,389,390,391,39 2,409,410,411,425,439,440,441,456, 458,477,490,491,505,519,539,540,56 0,561,919,934 NFCST Scalar 3 I(4) 4 scalar COM NFLD Scalar 8 I(4) 4 scalar COM NGMSLP2 Subr 189 189 NHEAT Scalar 3 I(4) 4 scalar COM NINC Param 7 I(4) 4 scalar NORTH Local 104 L(4) 4 scalar NPHS Scalar 5 I(4) 4 scalar COM 799 NPREC Scalar 4 I(4) 4 scalar COM NRADL Scalar 5 I(4) 4 scalar COM NRADS Scalar 5 I(4) 4 scalar COM NRDLW Scalar 3 I(4) 4 scalar COM NRDSW Scalar 3 I(4) 4 scalar COM NREC0 Local 267 I(4) 4 scalar 267,650,651,716,717,755,756 NREC1 Local 649 I(4) 4 scalar 649,650,651,715,716,717,754,755,75 6 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 157,734,773 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 103 L(4) 4 scalar COM OMGA Scalar 2 R(4) 4 3 4127355 COM 477 OMGAOT Common 2 16509420 ONEPS Param 39 R(4) 4 scalar OUTFIL Common 3 21 OUTPUT Subr 183 183,192,340,353,366,381,401,417,43 1,447,464,469,482,497,511,524,552, 573,642,658,683,708,747,786,832,85 6,878,925,940 OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar P1D Local 107 R(4) 4 2 91719 374,376,389,393,409,412 PARALLEL Common 9 8240 PD Scalar 3 R(4) 4 2 91719 COM PDSL Scalar 6 R(4) 4 2 91719 COM 306,608,812,816,825,847,903 PDVP1 Scalar 6 R(4) 4 2 91719 COM PFSL Scalar 7 R(4) 4 2 91719 COM PINT Scalar 4 R(4) 4 3 4219074 COM 257 PMID Local 109 R(4) 4 3 4127355 257,335,374,389,409 Page 27 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References PP Local 306 R(4) 4 scalar 306,307,608,609 PPP Local 112 R(4) 4 2 91719 825,829 PQ0 Param 41 R(4) 4 scalar 307,609 PREC Scalar 7 R(4) 4 2 91719 COM 804,809,810 PSLP Scalar 6 R(4) 4 2 91719 COM 181,811,815 PT Scalar 5 R(4) 4 scalar COM 306,608,812,816,825,903 PTDYN Scalar 3 R(4) 4 scalar COM PTSL Scalar 7 R(4) 4 2 91719 COM PVRBLS Common 2 23480064 Q Scalar 5 R(4) 4 3 4127355 COM 216,302,391,410,425,439,604,822 Q1D Local 107 R(4) 4 2 91719 391,393,410,412,439,442 Q2 Scalar 9 R(4) 4 3 4127355 COM 519,847,864,866,870,903 QC Local 617 R(4) 4 scalar 617,619,622 QCD Local 112 R(4) 4 2 91719 823,829 QI Local 111 R(4) 4 2 91719 308,312,610,612,617 QICE Local 112 R(4) 4 2 91719 824,829 QINT Local 111 R(4) 4 2 91719 611,612,617 QKL Local 302 R(4) 4 scalar 302,312,604,622 QRAIN Local 113 R(4) 4 2 91719 806,814,829 QS Scalar 3 R(4) 4 2 91719 COM QSL Scalar 9 R(4) 4 2 91719 COM QSNO Local 113 R(4) 4 2 91719 805,818,829 QV Local 112 R(4) 4 2 91719 822,829 QW Local 307 R(4) 4 scalar 307,308,609,610,611 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 216 RAINCON Param 97 R(4) 4 scalar 814 RAINRATE Local 810 R(4) 4 scalar 810,813 RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDETA Scalar 4 R(4) 4 1 45 COM RDOCP Param 40 R(4) 4 scalar RDRST2D Subr 650 650,716,755 RDTPHS Local 799 R(4) 4 scalar 799,809,810 RDYN Scalar 3 R(4) 4 scalar COM RES Scalar 3 R(4) 4 2 91719 COM RESTRT Scalar 103 L(4) 4 scalar COM RF Scalar 5 R(4) 4 2 91719 COM RG Param 36 R(4) 4 scalar RICHNO Local 110 R(4) 4 3 4127355 909,934 RITE2 Scalar 5 L(4) 4 scalar COM RITEHD Scalar 5 L(4) 4 scalar COM RLWIN Scalar 3 R(4) 4 2 91719 COM RLWOUT Scalar 3 R(4) 4 2 91719 COM RLWTOA Scalar 3 R(4) 4 2 91719 COM ROG Param 32 R(4) 4 scalar RQKL Local 620 R(4) 4 scalar 620,622,625,626,628,630 RQSTFLD Common 7 64130 RRNUM Local 720 R(4) 4 scalar 720,722,727,759,761,766 RSTFIL Scalar 2 CHAR 13 scalar COM RSWIN Scalar 3 R(4) 4 2 91719 COM RSWOUT Scalar 3 R(4) 4 2 91719 COM RSWTOA Scalar 3 R(4) 4 2 91719 COM RTD Param 34 R(4) 4 scalar Page 28 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References RUN Scalar 103 L(4) 4 scalar COM SBD Scalar 6 R(4) 4 scalar COM SCLFLD Subr 396 396,641 SI Scalar 6 R(4) 4 2 91719 COM SICE Scalar 3 R(4) 4 2 91719 COM SIGMA Scalar 103 L(4) 4 scalar COM SLP Scalar 5 R(4) 4 2 91719 COM 190 SM Scalar 3 R(4) 4 2 91719 COM 847,870 SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar SNO Scalar 6 R(4) 4 2 91719 COM SNOCON Param 98 R(4) 4 scalar 818 SNORATE Local 809 R(4) 4 scalar 809,817 SPHI0 Scalar 9 R(4) 4 scalar COM SPL Scalar 7 R(4) 4 1 39 COM SR Scalar 3 R(4) 4 2 91719 COM 809,810 STDOUT Scalar 2 I(4) 4 scalar COM STDRD Scalar 103 L(4) 4 scalar COM SVALUE Scalar 8 R(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM 216,301,361,375,390,411,603,811,81 2,815,816,821,903 T1D Local 107 R(4) 4 2 91719 375,376,390,393,411,412 T500 Scalar 6 R(4) 4 2 91719 COM TCUCN Local 134 R(4) 4 2 91719 755,766 TERM1 Local 811 R(4) 4 scalar 811,814 TERM2 Local 812 R(4) 4 scalar 812,814 TERM3 Local 813 R(4) 4 scalar 813,814 TERM4 Local 815 R(4) 4 scalar 815,818 TERM5 Local 816 R(4) 4 scalar 816,818 TERM6 Local 817 R(4) 4 scalar 817,818 TFRZ Param 39 R(4) 4 scalar THEAT Scalar 3 R(4) 4 scalar COM 735,774 THS Scalar 3 R(4) 4 2 91719 COM THZ0 Scalar 4 R(4) 4 2 91719 COM TKL Local 301 R(4) 4 scalar 301,304,307,603,606,609 TLM0D Scalar 6 R(4) 4 scalar COM TMT0 Local 304 R(4) 4 scalar 304,305,308,311,319,320,606,607,61 0,612 TMT15 Local 305 R(4) 4 scalar 305,607,611 TPH0D Scalar 6 R(4) 4 scalar COM TRAIN Local 134 R(4) 4 2 91719 716,727 TRDLW Scalar 3 R(4) 4 scalar COM TRDSW Scalar 3 R(4) 4 scalar COM TSHDE Scalar 8 R(4) 4 1 61 COM TSL Scalar 8 R(4) 4 2 91719 COM TSPH Scalar 5 R(4) 4 scalar COM 734,773 TT Local 112 R(4) 4 2 91719 821,829 TTND Local 134 R(4) 4 2 91719 650,653 TTV Scalar 8 R(4) 4 2 91719 COM TWBS Scalar 5 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM 440,456,490 U00 Scalar 3 R(4) 4 2 91719 COM 310,614 U00KL Local 310 R(4) 4 scalar 310,312,614 UL Scalar 3 R(4) 4 1 90 COM 310,614 US Local 582 R(4) 4 scalar 582 Page 29 Source Listing ETAFLD2 2014-12-17 20:47 Symbol Table ETAFLD2.f Name Object Declared Type Bytes Dimen Elements Attributes References UTIM Local 159 R(4) 4 scalar 159,162,164,310,614 UZ0 Scalar 4 R(4) 4 2 91719 COM V Scalar 5 R(4) 4 3 4127355 COM 441,458,491 VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VCON1 Param 99 R(4) 4 scalar VCON2 Param 99 R(4) 4 scalar VIS Local 113 R(4) 4 2 91719 829,830 VRBLS Common 2 67138308 VTM Scalar 5 R(4) 4 3 4127355 COM VZ0 Scalar 4 R(4) 4 2 91719 COM WBD Scalar 6 R(4) 4 scalar COM WPDAR Scalar 8 R(4) 4 2 91719 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 ZINT Scalar 4 R(4) 4 3 4219074 COM 205,218,247,258,847,864,867,871,89 2,903 ZMID Local 109 R(4) 4 3 4127355 258,348,505 Page 30 Source Listing ETAFLD2 2014-12-17 20:47 Subprograms/Common Blocks ETAFLD2.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References ACMCLH Common 2 24 ACMRDL Common 2 2201268 ACMRDS Common 2 2201268 CLDWTR Common 2 17610408 CTLBLK Common 2 96 DYNAM Common 2 3675964 ETAFLD2 Subr 1 EXTRA Common 3 71506268 IOUNIT Common 4 20 LOOPS Common 2 748952 MAPOT Common 2 1448 MASKS Common 2 34853220 OMGAOT Common 2 16509420 OUTFIL Common 3 21 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 Page 31 Source Listing ETAFLD2 2014-12-17 20:47 ETAFLD2.f -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -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 : ETAFLD2.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100