Page 1 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1 SUBROUTINE SURFCE2(IMOUT,JMOUT) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: SURFCE2 POST SURFACE BASED FIELDS 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 92-12-21 6 C 7 C ABSTRACT: 8 C THIS ROUTINE POSTS SURFACE BASED FIELDS. 9 C . 10 C 11 C PROGRAM HISTORY LOG: 12 C 92-12-21 RUSS TREADON 13 C 94-08-04 MICHAEL BALDWIN - ADDED OUTPUT OF SFC FLUXES OF 14 C SENS AND LATENT HEAT AND THETA AT Z0 15 C 94-11-04 MICHAEL BALDWIN - ADDED INSTANTANEOUS PRECIP TYPE 16 C 96-03-19 MICHAEL BALDWIN - CHANGE SOIL PARAMETERS 17 C 96-09-25 MICHAEL BALDWIN - ADDED SNOW RATIO FROM EXPLICIT SCHEME 18 C 96-10-17 MICHAEL BALDWIN - CHANGED SFCEVP,POTEVP TO ACCUM. TOOK 19 C OUT -PTRACE FOR ACSNOW,SSROFF,BGROFF. 20 C 97-04-23 MICHAEL BALDWIN - TOOK OUT -PTRACE FOR ALL PRECIP FIELDS 21 C 98-06-12 T BLACK - CONVERSION FROM 1-D TO 2-D 22 C 98-07-17 MIKE BALDWIN - REMOVED LABL84 23 C 98-08-18 MIKE BALDWIN - COMPUTE RH OVER ICE 24 C 98-12-22 MIKE BALDWIN - BACK OUT RH OVER ICE 25 C 00-01-04 JIM TUCCILLO - MPI VERSION 26 C 03-01-15 H CHUANG - ADDED VEG AND SOIL TYPES 27 C 03-01-17 CHUANG AND EK- ADDED CANOPY CONDUCTANCE, SNOW COVER AND DEPTH, 28 C AND OTHER LAND SURFACE VARIABLES 29 C 03-02-20 DUSAN JOVIC - ADDED 30m T,Q,WIND AND INTEGRATED MOISTURE FLUXES 30 C 31 C USAGE: CALL SURFCE2(IMOUT,JMOUT) 32 C INPUT ARGUMENT LIST: 33 C IMOUT - FIRST DIMENSION OF OUTPUT GRID. 34 C JMOUT - SECOND DIMENSION OF OUTPUT GRID. 35 C 36 C OUTPUT ARGUMENT LIST: 37 C 38 C OUTPUT FILES: 39 C NONE 40 C 41 C SUBPROGRAMS CALLED: 42 C UTILITIES: 43 C E2OUT - INTERPOLATION/SMOOTHING ROUTINE. 44 C OUTPUT - DRIVER FOR OUTPUT ROUTINES. 45 C BOUND - ENFORCE LOWER AND UPPER LIMITS ON ARRAY ELEMENTS. 46 C SCLFLD - SCALE ARRAY ELEMENTS BY CONSTANT. 47 C SHELTR2 - COMPUTE 2M TEMPERATURE AND SPECIFIC HUMIDITY. 48 C ANEMLV6 - COMPUTE 10M U AND V WINDS. 49 C DEWPOINT - COMPUTE DEWPOINT TEMPERATURE. 50 C CALDRG - COMPUTE SURFACE LAYER DRAG COEFFICENT 51 C CALTAU - COMPUTE SURFACE LAYER U AND V WIND STRESSES. 52 C 53 C LIBRARY: 54 C COMMON - CTLBLK 55 C RQSTFLD 56 C EXTRA 57 C VRBLS Page 2 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 58 C MAPOT 59 C MASKS 60 C PVRBLS 61 C CLDWTR 62 C LOOPS 63 C PHYS2 64 C SRFDSP 65 C CNVCLD 66 C LLGRDS 67 C SOIL 68 C ACMSFC 69 C ACMPRE 70 C IOUNIT 71 C 72 C ATTRIBUTES: 73 C LANGUAGE: FORTRAN 74 C MACHINE : CRAY C-90 75 C$$$ 76 C 77 C 78 C INCLUDE GRID DIMENSIONS. SET/DERIVE OTHER PARAMETERS. 79 C 80 INCLUDE "parmeta" 98 INCLUDE "parmout" 102 INCLUDE "params" 144 INCLUDE "parm.tbl" 146 INCLUDE "parmsoil" 147 C 148 C IN NGM SUBROUTINE OUTPUT WE FIND THE FOLLOWING COMMENT. 149 C "IF THE FOLLOWING THRESHOLD VALUES ARE CHANGED, CONTACT 150 C TDL/SYNOPTIC-SCALE TECHNIQUES BRANCH (PAUL DALLAVALLE 151 C AND JOHN JENSENIUS). THEY MAY BE USING IT IN ONE OF 152 C THEIR PACKING CODES." THE THRESHOLD VALUE IS 0.01 INCH 153 C OR 2.54E-4 METER. PRECIPITATION VALUES LESS THAN THIS 154 C THRESHOLD ARE SET TO MINUS ONE TIMES THIS THRESHOLD. 156 PARAMETER (PTRACE = 0.000254E0) 157 C 158 C SET CELCIUS TO KELVIN AND SECOND TO HOUR CONVERSION. 159 PARAMETER (C2K = 273.15) 160 PARAMETER (SEC2HR = 1./3600.) 161 C 162 C DECLARE VARIABLES. 163 C 164 LOGICAL RUN,FIRST,RESTRT,SIGMA,OLDRD,STDRD 165 INTEGER IWX1(IM,JM),NROOTS(IM,JM) 166 REAL PSFC(IM,JM),TSFC(IM,JM),QSFC(IM,JM),RHSFC(IM,JM) 167 REAL ZSFC(IM,JM),THSFC(IM,JM),DWPSFC(IM,JM),EVP(IM,JM) 168 REAL ANCPRC(IM,JM),P1D(IM,JM),T1D(IM,JM),Q1D(IM,JM) 169 REAL EGRID1(IM,JM),EGRID2(IM,JM),UA(IM,JM),VA(IM,JM) 170 REAL GC(IM,JM) 171 REAL GRID1(IMOUT,JMOUT),GRID2(IMOUT,JMOUT),IW(IM,JM),IWM1 172 REAL SLEET(IM,JM),RAIN(IM,JM),FREEZR(IM,JM),SNOW(IM,JM) 173 REAL ECAN(IM,JM),EDIR(IM,JM),ETRANS(IM,JM),ESNOW(IM,JM) 174 &,SMCDRY(IM,JM),SMCMAX(IM,JM) 175 REAL RSMIN(IM,JM),SMCWLT(IM,JM),SMCREF(IM,JM) 176 & ,RCS(IM,JM),RCQ(IM,JM),RCT(IM,JM),RCSOIL(IM,JM) 177 REAL P10(IM,JM),P30(IM,JM),T10(IM,JM),T30(IM,JM) Page 3 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 178 C 179 C INCLUDE COMMON BLOCKS. 180 INCLUDE "CTLBLK.comm" 195 INCLUDE "RQSTFLD.comm" 208 INCLUDE "EXTRA.comm" 220 INCLUDE "VRBLS.comm" 227 INCLUDE "MAPOT.comm" 240 INCLUDE "MASKS.comm" 248 INCLUDE "PVRBLS.comm" 260 INCLUDE "CLDWTR.comm" 266 INCLUDE "LOOPS.comm" 274 INCLUDE "PHYS2.comm" 295 INCLUDE "SRFDSP.comm" 301 INCLUDE "CNVCLD.comm" 307 INCLUDE "LLGRDS.comm" 321 INCLUDE "SOIL.comm" 330 INCLUDE "ACMSFC.comm" 337 INCLUDE "ACMPRE.comm" 343 INCLUDE "IOUNIT.comm" 351 INCLUDE "ACMRDS.comm" 357 INCLUDE "QFLX.comm" 380 INCLUDE "PPTASM.comm" 381 C 382 C**************************************************************************** 383 C 384 C START SURFCE. 385 C 386 C COMPUTE IW AT SFC FOR SFC AND 2M RH 387 C 410 IF ( (IGET(076).GT.0).OR.(IGET(114).GT.0) ) THEN 411 CLIMIT =1.0E-20 412 IW=0. 413 C 414 DO L=2,LM 415 DO J=JSTA,JEND 416 DO I=1,IM 417 IF (L.LE.LMH(I,J)) THEN 418 IWM1=IW(I,J) 419 IF(CWM(I,J,L).GT.CLIMIT) THEN 420 IF(T(I,J,L).LT.258.15)THEN 421 IW(I,J)=1. 422 ELSEIF(T(I,J,L).GE.273.15)THEN 423 IW(I,J)=0. 424 ELSE 425 IF(IWM1.EQ.1.0)IW(I,J)=1. 426 ENDIF 427 ELSE 428 IW(I,J)=0. 429 ENDIF 430 ENDIF 431 ENDDO 432 ENDDO 433 ENDDO 434 C 435 ENDIF 436 C 437 C*** BLOCK 1. SURFACE BASED FIELDS. Page 4 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 438 C 439 C IF ANY OF THE FOLLOWING "SURFACE" FIELDS ARE REQUESTED, 440 C WE NEED TO COMPUTE THE FIELDS FIRST. 441 C 442 IF ( (IGET(024).GT.0).OR.(IGET(025).GT.0).OR. 443 X (IGET(026).GT.0).OR.(IGET(027).GT.0).OR. 444 X (IGET(028).GT.0).OR.(IGET(029).GT.0).OR. 445 X (IGET(154).GT.0).OR. 446 X (IGET(034).GT.0).OR.(IGET(076).GT.0) ) THEN 447 C 448 DO 40 J=JSTA,JEND 449 DO 40 I=1,IM 450 C 451 C PRESSURE AND TEMPERATURE AT 10 AND 30 M. 452 LMHK=LMH(I,J) 453 P10(I,J)=(PD(I,J)+PT)*EXP(-10.0*G/(287.04*T(I,J,LMHK))) 454 P30(I,J)=(PD(I,J)+PT)*EXP(-30.0*G/(287.04*T(I,J,LMHK))) 455 T10(I,J)=TH10(I,J)*(P10(I,J)/P1000)**CAPA 456 T30(I,J)=TH30(I,J)*(P30(I,J)/P1000)**CAPA 457 C 458 C SCALE ARRAY FIS BY GI TO GET SURFACE HEIGHT. 459 ZSFC(I,J)=FIS(I,J)*GI 460 C 461 C SURFACE PRESSURE. 462 PSFC(I,J)=PD(I,J)+PT 463 C 464 C SURFACE (SKIN) POTENTIAL TEMPERATURE AND TEMPERATURE. 465 THSFC(I,J)=THS(I,J) 466 TSFC(I,J) =THSFC(I,J)*(PSFC(I,J)/P1000)**CAPA 467 C 468 C SURFACE SPECIFIC HUMIDITY, RELATIVE HUMIDITY, 469 C AND DEWPOINT. ADJUST SPECIFIC HUMIDITY IF 470 C RELATIVE HUMIDITY EXCEEDS 0.1 OR 1.0. 471 C 472 QSFC(I,J)=QS(I,J) 473 QSFC(I,J)=AMAX1(H1M12,QSFC(I,J)) 474 TSFCK =TSFC(I,J) 475 C 476 TMT0=TSFCK-273.16 477 TMT15=AMIN1(TMT0,-15.) 478 AI=0.008855 479 BI=1. 480 IF(TMT0.LT.-20.)THEN 481 AI=0.007225 482 BI=0.9674 483 ENDIF 484 QW=PQ0/PSFC(I,J) 485 1 *EXP(A2*(TSFCK-A3)/(TSFCK-A4)) 486 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 487 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 488 IF(TMT0.LT.-15.)THEN 489 QSAT=QI 490 ELSEIF(TMT0.GE.0.)THEN 491 QSAT=QINT 492 ELSE 493 IF(IW(I,J).GT.0.0) THEN 494 QSAT=QI Page 5 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 495 ELSE 496 QSAT=QINT 497 ENDIF 498 ENDIF 499 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 500 C DELETE THIS LINE TO SWITCH BACK TO RH VS ICE 501 QSAT=QW 502 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 503 C 504 RHSFC(I,J)=QSFC(I,J)/QSAT 505 506 IF (RHSFC(I,J).GT.H1 ) RHSFC(I,J) = H1 507 IF (RHSFC(I,J).LT.D00) RHSFC(I,J) = D01 508 QSFC(I,J) = RHSFC(I,J)*QSAT 509 EVP(I,J) = PSFC(I,J)*QSFC(I,J)/(EPS+ONEPS*QSFC(I,J)) 510 EVP(I,J) = EVP(I,J)*D001 511 C 512 C ACCUMULATED NON-CONVECTIVE PRECIP. 513 IF(IGET(034).GT.0)THEN 514 IF(LVLS(1,IGET(034)).GT.0)THEN 515 ANCPRC(I,J)=ACPREC(I,J)-CUPREC(I,J) 516 ENDIF 517 ENDIF 518 519 40 CONTINUE 520 C 521 C INTERPOLATE/OUTPUT REQUESTED SURFACE FIELDS. 522 C 523 C SURFACE PRESSURE. 524 IF (IGET(024).GT.0) THEN 525 CALL E2OUT(024,000,PSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 526 ID(1:25) = 0 527 CALL OUTPUT(IOUTYP,IGET(024),LVLS(1,IGET(024)), 528 X GRID1,IMOUT,JMOUT) 529 ENDIF 530 C 531 C SURFACE PRESSURE NN. 532 IF (IGET(247).GT.0) THEN 533 CALL E2OUT(247,000,PSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 534 ID(1:25) = 0 535 CALL OUTPUT(IOUTYP,IGET(247),LVLS(1,IGET(247)), 536 X GRID1,IMOUT,JMOUT) 537 ENDIF 538 C 539 C SURFACE HEIGHT. 540 IF (IGET(025).GT.0) THEN 541 CALL E2OUT(025,000,ZSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 542 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) 543 ID(1:25) = 0 544 CALL OUTPUT(IOUTYP,IGET(025),LVLS(1,IGET(025)), 545 X GRID1,IMOUT,JMOUT) 546 ENDIF 547 C 548 C SURFACE HEIGHT NN. 549 IF (IGET(245).GT.0) THEN 550 CALL E2OUT(245,000,ZSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 551 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) Page 6 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 552 ID(1:25) = 0 553 CALL OUTPUT(IOUTYP,IGET(245),LVLS(1,IGET(245)), 554 X GRID1,IMOUT,JMOUT) 555 ENDIF 556 C 557 C SURFACE (SKIN) TEMPERATURE. 558 IF (IGET(026).GT.0) THEN 559 CALL E2OUT(026,000,TSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 560 ID(1:25) = 0 561 CALL OUTPUT(IOUTYP,IGET(026),LVLS(1,IGET(026)), 562 X GRID1,IMOUT,JMOUT) 563 ENDIF 564 C 565 C SURFACE (SKIN) POTENTIAL TEMPERATURE. 566 IF (IGET(027).GT.0) THEN 567 CALL E2OUT(027,000,THSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 568 ID(1:25) = 0 569 CALL OUTPUT(IOUTYP,IGET(027),LVLS(1,IGET(027)), 570 X GRID1,IMOUT,JMOUT) 571 ENDIF 572 C 573 C SURFACE SPECIFIC HUMIDITY. 574 IF (IGET(028).GT.0) THEN 575 CALL E2OUT(028,000,QSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 576 CALL BOUND(GRID1,H1M12,H99999,IMOUT,JMOUT) 577 ID(1:25) = 0 578 CALL OUTPUT(IOUTYP,IGET(028),LVLS(1,IGET(028)), 579 X GRID1,IMOUT,JMOUT) 580 ENDIF 581 C 582 C SURFACE DEWPOINT TEMPERATURE. 583 IF (IGET(029).GT.0) THEN 584 CALL DEWPOINT(EVP,DWPSFC,IM,JM) 585 CALL E2OUT(029,000,DWPSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 586 ID(1:25) = 0 587 CALL OUTPUT(IOUTYP,IGET(029),LVLS(1,IGET(029)), 588 X GRID1,IMOUT,JMOUT) 589 ENDIF 590 C 591 C SURFACE RELATIVE HUMIDITY. 592 IF (IGET(076).GT.0) THEN 593 CALL E2OUT(076,000,RHSFC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 594 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 595 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 596 ID(1:25) = 0 597 CALL OUTPUT(IOUTYP,IGET(076),LVLS(1,IGET(076)), 598 X GRID1,IMOUT,JMOUT) 599 ENDIF 600 C 601 ENDIF 602 C 603 C ADDITIONAL SURFACE-SOIL LEVEL FIELDS. 604 C 605 606 DO L=1,NSOIL 607 C SOIL TEMPERATURE. 608 IF (IGET(116).GT.0) THEN Page 7 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 609 IF (LVLS(L,IGET(116)).GT.0) THEN 610 DO J=JSTA,JEND 611 DO I=1,IM 612 EGRID1(I,J)=STC(I,J,L) 613 ENDDO 614 ENDDO 615 CALL E2OUT(116,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 616 ID(1:25) = 0 617 DTOP=0. 618 DO LS=1,L-1 619 DTOP=DTOP+SLDPTH(LS) 620 ENDDO 621 DBOT=DTOP+SLDPTH(L) 622 ID(10) = NINT(DTOP*100.) 623 ID(11) = NINT(DBOT*100.) 624 CALL OUTPUT(IOUTYP,IGET(116),L,GRID1,IMOUT,JMOUT) 625 ENDIF 626 ENDIF 627 C 628 C SOIL MOISTURE. 629 IF (IGET(117).GT.0) THEN 630 IF (LVLS(L,IGET(117)).GT.0) THEN 631 DO J=JSTA,JEND 632 DO I=1,IM 633 EGRID1(I,J)=SMC(I,J,L) 634 ENDDO 635 ENDDO 636 CALL E2OUT(117,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 637 ID(1:25) = 0 638 DTOP=0. 639 DO LS=1,L-1 640 DTOP=DTOP+SLDPTH(LS) 641 ENDDO 642 DBOT=DTOP+SLDPTH(L) 643 ID(10) = NINT(DTOP*100.) 644 ID(11) = NINT(DBOT*100.) 645 CALL OUTPUT(IOUTYP,IGET(117),L,GRID1,IMOUT,JMOUT) 646 ENDIF 647 ENDIF 648 ! ADD LIQUID SOIL MOISTURE 649 IF (IGET(225).GT.0) THEN 650 IF (LVLS(L,IGET(225)).GT.0) THEN 651 DO J=JSTA,JEND 652 DO I=1,IM 653 EGRID1(I,J)=SH2O(I,J,L) 654 ENDDO 655 ENDDO 656 CALL E2OUT(225,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 657 ID(1:25) = 0 658 DTOP=0. 659 DO LS=1,L-1 660 DTOP=DTOP+SLDPTH(LS) 661 ENDDO 662 DBOT=DTOP+SLDPTH(L) 663 ID(10) = NINT(DTOP*100.) 664 ID(11) = NINT(DBOT*100.) 665 ID(02) = 130 Page 8 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 666 CALL OUTPUT(IOUTYP,IGET(225),L,GRID1,IMOUT,JMOUT) 667 ENDIF 668 ENDIF 669 ! END OF NSOIL LOOP 670 ENDDO 671 C 672 C BOTTOM SOIL TEMPERATURE. 673 IF (IGET(115).GT.0) THEN 674 CALL E2OUT(115,000,SOILTB,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 675 ID(1:25) = 0 676 ISVALUE = 800 677 ID(11) = ISVALUE 678 CALL OUTPUT(IOUTYP,IGET(115),LVLS(1,IGET(115)), 679 X GRID1,IMOUT,JMOUT) 680 ENDIF 681 C 682 C SOIL MOISTURE AVAILABILITY 683 IF (IGET(171).GT.0) THEN 684 CALL E2OUT(171,000,SMSTAV,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 685 ID(1:25) = 0 686 ID(10) = 0 687 ID(11) = 100 688 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 689 CALL OUTPUT(IOUTYP,IGET(171),LVLS(1,IGET(171)), 690 X GRID1,IMOUT,JMOUT) 691 ENDIF 692 C 693 C TOTAL SOIL MOISTURE 694 IF (IGET(036).GT.0) THEN 695 CALL E2OUT(036,000,SMSTOT,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 696 ID(1:25) = 0 697 ID(10) = 0 698 ID(11) = 200 699 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 700 CALL OUTPUT(IOUTYP,IGET(036),LVLS(1,IGET(036)), 701 X GRID1,IMOUT,JMOUT) 702 ENDIF 703 C 704 C PLANT CANOPY SURFACE WATER. 705 IF ( IGET(118).GT.0 ) THEN 706 CALL E2OUT(118,000,CMC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 707 ID(1:25) = 0 708 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 709 CALL OUTPUT(IOUTYP,IGET(118),LVLS(1,IGET(118)), 710 X GRID1,IMOUT,JMOUT) 711 ENDIF 712 C 713 C SNOW WATER EQUIVALENT. 714 IF ( IGET(119).GT.0 ) THEN 715 CALL E2OUT(119,000,SNO,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 716 ID(1:25) = 0 717 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 718 CALL OUTPUT(IOUTYP,IGET(119),LVLS(1,IGET(119)), 719 X GRID1,IMOUT,JMOUT) 720 ENDIF 721 C 722 C PERCENT SNOW COVER. Page 9 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 723 IF ( IGET(120).GT.0 ) THEN 724 DO J=JSTA,JEND 725 DO I=1,IM 726 SNEQV=SNO(I,J) 727 IVEG=IVGTYP(I,J) 728 CALL SNFRAC (SNEQV,IVEG,SNCOVR) 729 PCTSNO(I,J)=SNCOVR 730 END DO 731 END DO 732 CALL E2OUT(120,000,PCTSNO,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 733 CALL BOUND(GRID1,D00,H1,IMOUT,JMOUT) 734 ID(1:25) = 0 735 CALL OUTPUT(IOUTYP,IGET(120),LVLS(1,IGET(120)), 736 X GRID1,IMOUT,JMOUT) 737 ENDIF 738 ! ADD SNOW DEPTH 739 IF ( IGET(224).GT.0 ) THEN 740 CALL E2OUT(224,000,SI,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 741 ID(1:25) = 0 742 CALL OUTPUT(IOUTYP,IGET(224),LVLS(1,IGET(224)), 743 X GRID1,IMOUT,JMOUT) 744 ENDIF 745 ! ADD POTENTIAL EVAPORATION 746 C IF ( IGET(242).GT.0 ) THEN 747 C CALL E2OUT(242,000,POTEVP,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 748 C ID(1:25) = 0 749 C CALL OUTPUT(IOUTYP,IGET(242),LVLS(1,IGET(242)), 750 C X GRID1,IMOUT,JMOUT) 751 C END IF 752 ! ADD EC,EDIR,ETRANS,ESNOW,SMCDRY,SMCMAX 753 IF ( IGET(228).GT.0 .OR. IGET(229).GT.0 754 + .OR.IGET(230).GT.0 .OR. IGET(231).GT.0 755 + .OR.IGET(232).GT.0 .OR. IGET(233).GT.0) THEN 756 DO J=JSTA,JEND 757 DO I=1,IM 758 C ---------------------------------------------------------------------- 759 ! IF(QWBS(I,J).gt.0.001)print*,'NONZERO QWBS',i,j,QWBS(I,J) 760 ! IF(abs(SM(I,J)-0.).lt.1.0E-5)THEN 761 IF( (abs(SM(I,J)-0.) .lt. 1.0E-5) .AND. 762 & (abs(SICE(I,J)-0.) .lt. 1.0E-5) ) THEN 763 CALL ETCALC(QWBS(I,J),POTEVP(I,J),SNO(I,J),VEGFRC(I,J) 764 & , ISLTYP(I,J),SH2O(I,J,1:1),CMC(I,J) 765 & , ECAN(I,J),EDIR(I,J),ETRANS(I,J),ESNOW(I,J),SMCDRY(I,J) 766 & , SMCMAX(I,J) ) 767 ELSE 768 ECAN(I,J)=0. 769 EDIR(I,J)=0. 770 ETRANS(I,J)=0. 771 ESNOW(I,J)=0. 772 SMCDRY(I,J)=0. 773 SMCMAX(I,J)=0. 774 END IF 775 END DO 776 END DO 777 IF ( IGET(228).GT.0 )THEN 778 CALL E2OUT(228,000,ECAN,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 779 ID(1:25) = 0 Page 10 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 780 CALL OUTPUT(IOUTYP,IGET(228),LVLS(1,IGET(228)), 781 X GRID1,IMOUT,JMOUT) 782 END IF 783 ! 784 IF ( IGET(229).GT.0 )THEN 785 CALL E2OUT(229,000,EDIR,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 786 ID(1:25) = 0 787 CALL OUTPUT(IOUTYP,IGET(229),LVLS(1,IGET(229)), 788 X GRID1,IMOUT,JMOUT) 789 END IF 790 ! 791 IF ( IGET(230).GT.0 )THEN 792 CALL E2OUT(230,000,ETRANS,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 793 ID(1:25) = 0 794 CALL OUTPUT(IOUTYP,IGET(230),LVLS(1,IGET(230)), 795 X GRID1,IMOUT,JMOUT) 796 END IF 797 ! 798 IF ( IGET(231).GT.0 )THEN 799 CALL E2OUT(231,000,ESNOW,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 800 ID(1:25) = 0 801 ID(02)= 130 802 CALL OUTPUT(IOUTYP,IGET(231),LVLS(1,IGET(231)), 803 X GRID1,IMOUT,JMOUT) 804 END IF 805 ! 806 IF ( IGET(232).GT.0 )THEN 807 CALL E2OUT(232,000,SMCDRY,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 808 ID(1:25) = 0 809 ID(02)= 130 810 CALL OUTPUT(IOUTYP,IGET(232),LVLS(1,IGET(232)), 811 X GRID1,IMOUT,JMOUT) 812 END IF 813 ! 814 IF ( IGET(233).GT.0 )THEN 815 CALL E2OUT(233,000,SMCMAX,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 816 ID(1:25) = 0 817 ID(02)= 130 818 CALL OUTPUT(IOUTYP,IGET(233),LVLS(1,IGET(233)), 819 X GRID1,IMOUT,JMOUT) 820 END IF 821 ! 822 ENDIF 823 C 824 C 825 C 826 C*** BLOCK 2. SHELTER (2M) LEVEL FIELDS. 827 C 828 C COMPUTE/POST SHELTER LEVEL FIELDS. 829 C 830 IF ( (IGET(106).GT.0).OR.(IGET(112).GT.0).OR. 831 X (IGET(113).GT.0).OR.(IGET(114).GT.0).OR. 832 X (IGET(138).GT.0) ) THEN 833 C 834 C CALL SHELTR2(PSHLTR,QSHLTR,TSHLTR) 835 C 836 C SHELTER LEVEL TEMPERATURE Page 11 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 837 IF (IGET(106).GT.0) THEN 838 CALL E2OUT(106,000,TSHLTR,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 839 ID(1:25) = 0 840 ISVALUE = 2 841 ID(10) = MOD(ISVALUE/256,256) 842 ID(11) = MOD(ISVALUE,256) 843 CALL OUTPUT(IOUTYP,IGET(106),LVLS(1,IGET(106)), 844 X GRID1,IMOUT,JMOUT) 845 ENDIF 846 C 847 C SHELTER LEVEL SPECIFIC HUMIDITY. 848 IF (IGET(112).GT.0) THEN 849 CALL E2OUT(112,000,QSHLTR,EGRID1,GRID1,GRID2,IMOUT,JMOUT) 850 CALL BOUND (GRID1,H1M12,H99999,IMOUT,JMOUT) 851 ID(1:25) = 0 852 ISVALUE = 2 853 ID(10) = MOD(ISVALUE/256,256) 854 ID(11) = MOD(ISVALUE,256) 855 CALL OUTPUT(IOUTYP,IGET(112),LVLS(1,IGET(112)), 856 X GRID1,IMOUT,JMOUT) 857 ENDIF 858 C 859 C SHELTER LEVEL DEWPOINT. 860 IF (IGET(113).GT.0) THEN 861 DO J=JSTA,JEND 862 DO I=1,IM 863 EVP(I,J)=PSHLTR(I,J)*QSHLTR(I,J)/(EPS+ONEPS*QSHLTR(I,J)) 864 EVP(I,J)=EVP(I,J)*D001 865 ENDDO 866 ENDDO 867 CALL DEWPOINT(EVP,EGRID1,IM,JM) 868 CALL E2OUT(113,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 869 ID(1:25) = 0 870 ISVALUE = 2 871 ID(10) = MOD(ISVALUE/256,256) 872 ID(11) = MOD(ISVALUE,256) 873 CALL OUTPUT(IOUTYP,IGET(113),LVLS(1,IGET(113)), 874 X GRID1,IMOUT,JMOUT) 875 ENDIF 876 C 877 C SHELTER LEVEL RELATIVE HUMIDITY. 878 IF (IGET(114).GT.0) THEN 879 CALL CALRH2(PSHLTR,TSHLTR,QSHLTR,IW,EGRID1,IM,JM) 880 CALL E2OUT(114,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 881 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 882 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 883 ID(1:25) = 0 884 ISVALUE = 2 885 ID(10) = MOD(ISVALUE/256,256) 886 ID(11) = MOD(ISVALUE,256) 887 CALL OUTPUT(IOUTYP,IGET(114),LVLS(1,IGET(114)), 888 X GRID1,IMOUT,JMOUT) 889 ENDIF 890 C 891 C SHELTER LEVEL PRESSURE. 892 IF (IGET(138).GT.0) THEN 893 CALL E2OUT(138,000,PSHLTR,EGRID1,GRID1,GRID2,IMOUT,JMOUT) Page 12 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 894 ID(1:25) = 0 895 ISVALUE = 2 896 ID(10) = MOD(ISVALUE/256,256) 897 ID(11) = MOD(ISVALUE,256) 898 CALL OUTPUT(IOUTYP,IGET(138),LVLS(1,IGET(138)), 899 X GRID1,IMOUT,JMOUT) 900 ENDIF 901 C 902 ENDIF 903 C 904 C 905 C BLOCK 3. ANEMOMETER LEVEL (10M) WINDS, THETA, AND Q. 906 C 907 IF ( (IGET(064).GT.0).OR.(IGET(065).GT.0) ) THEN 908 C CALL ANEMLV6(UA,VA) 909 ID(1:25) = 0 910 ISVALUE = 10 911 ID(10) = MOD(ISVALUE/256,256) 912 ID(11) = MOD(ISVALUE,256) 913 C 914 C ANEMOMETER LEVEL U WIND AND/OR V WIND. 915 IF ((IGET(064).GT.0).OR.(IGET(065).GT.0)) THEN 916 CALL E2OUT(064,065,U10,V10,GRID1,GRID2,IMOUT,JMOUT) 917 IF (IGET(064).GT.0) CALL OUTPUT(IOUTYP,IGET(064), 918 X LVLS(1,IGET(064)),GRID1,IMOUT,JMOUT) 919 IF (IGET(065).GT.0) CALL OUTPUT(IOUTYP,IGET(065), 920 X LVLS(1,IGET(065)),GRID2,IMOUT,JMOUT) 921 ENDIF 922 ENDIF 923 C 924 C ANEMOMETER LEVEL (10 M) POTENTIAL TEMPERATURE. 925 C 926 IF (IGET(158).GT.0) THEN 927 ID(1:25) = 0 928 ISVALUE = 10 929 ID(10) = MOD(ISVALUE/256,256) 930 ID(11) = MOD(ISVALUE,256) 931 CALL E2OUT(158,000,TH10,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 932 CALL OUTPUT(IOUTYP,IGET(158), 933 X LVLS(1,IGET(158)),GRID1,IMOUT,JMOUT) 934 ENDIF 935 C 936 C ANEMOMETER LEVEL (10 M) TEMPERATURE. 937 C 938 IF (IGET(185).GT.0) THEN 939 ID(1:25) = 0 940 ISVALUE = 10 941 ID(10) = MOD(ISVALUE/256,256) 942 ID(11) = MOD(ISVALUE,256) 943 CALL E2OUT(185,000,T10,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 944 CALL OUTPUT(IOUTYP,IGET(185), 945 X LVLS(1,IGET(185)),GRID1,IMOUT,JMOUT) 946 ENDIF 947 C 948 C ANEMOMETER LEVEL (10 M) PRESSURE. 949 C 950 IF (IGET(186).GT.0) THEN Page 13 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 951 ID(1:25) = 0 952 ISVALUE = 10 953 ID(10) = MOD(ISVALUE/256,256) 954 ID(11) = MOD(ISVALUE,256) 955 CALL E2OUT(186,000,P10,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 956 CALL OUTPUT(IOUTYP,IGET(186), 957 X LVLS(1,IGET(186)),GRID1,IMOUT,JMOUT) 958 ENDIF 959 C 960 C ANEMOMETER LEVEL (10 M) SPECIFIC HUMIDITY. 961 C 962 IF (IGET(159).GT.0) THEN 963 ID(1:25) = 0 964 ISVALUE = 10 965 ID(10) = MOD(ISVALUE/256,256) 966 ID(11) = MOD(ISVALUE,256) 967 CALL E2OUT(159,000,Q10,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 968 CALL OUTPUT(IOUTYP,IGET(159), 969 X LVLS(1,IGET(159)),GRID1,IMOUT,JMOUT) 970 ENDIF 971 C 972 C 973 C BLOCK 3.1 30M WINDS, THETA, AND Q. 974 C 975 IF ( (IGET(193).GT.0).OR.(IGET(194).GT.0) ) THEN 976 C CALL ANEMLV6(UA,VA) 977 ID(1:25) = 0 978 ISVALUE = 30 979 ID(10) = MOD(ISVALUE/256,256) 980 ID(11) = MOD(ISVALUE,256) 981 C 982 C ANEMOMETER LEVEL U WIND AND/OR V WIND. 983 IF ((IGET(193).GT.0).OR.(IGET(194).GT.0)) THEN 984 CALL E2OUT(193,194,U30,V30,GRID1,GRID2,IMOUT,JMOUT) 985 IF (IGET(193).GT.0) CALL OUTPUT(IOUTYP,IGET(193), 986 X LVLS(1,IGET(193)),GRID1,IMOUT,JMOUT) 987 IF (IGET(194).GT.0) CALL OUTPUT(IOUTYP,IGET(194), 988 X LVLS(1,IGET(194)),GRID2,IMOUT,JMOUT) 989 ENDIF 990 ENDIF 991 C 992 C 30 M POTENTIAL TEMPERATURE. 993 C 994 IF (IGET(191).GT.0) THEN 995 ID(1:25) = 0 996 ISVALUE = 30 997 ID(10) = MOD(ISVALUE/256,256) 998 ID(11) = MOD(ISVALUE,256) 999 CALL E2OUT(191,000,TH30,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1000 CALL OUTPUT(IOUTYP,IGET(191), 1001 X LVLS(1,IGET(191)),GRID1,IMOUT,JMOUT) 1002 ENDIF 1003 C 1004 C 30 M TEMPERATURE. 1005 C 1006 IF (IGET(189).GT.0) THEN 1007 ID(1:25) = 0 Page 14 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1008 ISVALUE = 30 1009 ID(10) = MOD(ISVALUE/256,256) 1010 ID(11) = MOD(ISVALUE,256) 1011 CALL E2OUT(189,000,T30,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1012 CALL OUTPUT(IOUTYP,IGET(189), 1013 X LVLS(1,IGET(189)),GRID1,IMOUT,JMOUT) 1014 ENDIF 1015 C 1016 C 30 M PRESSURE. 1017 C 1018 IF (IGET(190).GT.0) THEN 1019 ID(1:25) = 0 1020 ISVALUE = 30 1021 ID(10) = MOD(ISVALUE/256,256) 1022 ID(11) = MOD(ISVALUE,256) 1023 CALL E2OUT(190,000,P30,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1024 CALL OUTPUT(IOUTYP,IGET(190), 1025 X LVLS(1,IGET(190)),GRID1,IMOUT,JMOUT) 1026 ENDIF 1027 C 1028 C 30 M SPECIFIC HUMIDITY. 1029 C 1030 IF (IGET(192).GT.0) THEN 1031 ID(1:25) = 0 1032 ISVALUE = 30 1033 ID(10) = MOD(ISVALUE/256,256) 1034 ID(11) = MOD(ISVALUE,256) 1035 CALL E2OUT(192,000,Q30,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1036 CALL OUTPUT(IOUTYP,IGET(192), 1037 X LVLS(1,IGET(192)),GRID1,IMOUT,JMOUT) 1038 ENDIF 1039 C 1040 C 1041 C 1042 C*** BLOCK 4. PRECIPITATION RELATED FIELDS. 1043 C 1044 C SNOW FRACTION FROM EXPLICIT CLOUD SCHEME. LABELLED AS 1045 C 'PROB OF FROZEN PRECIP' IN GRIB, 1046 C DIDN'T KNOW WHAT ELSE TO CALL IT 1047 IF (IGET(172).GT.0) THEN 1048 DO J=JSTA,JEND 1049 DO I=1,IM 1050 EGRID1(I,J)=SR(I,J)*100. 1051 ENDDO 1052 ENDDO 1053 CALL E2OUT(172,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1054 ID(1:25) = 0 1055 CALL OUTPUT(IOUTYP,IGET(172),LVLS(1,IGET(172)), 1056 X GRID1,IMOUT,JMOUT) 1057 ENDIF 1058 C INSTANTANEOUS PRECIPITATION RATE. 1059 IF (IGET(167).GT.0) THEN 1060 RDTPHS=1./DTQ2 1061 write(*,*) 'prdgen surface2 prate' 1062 DO J=JSTA,JEND 1063 DO I=1,IM 1064 EGRID1(I,J)=PREC(I,J)*RDTPHS Page 15 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1065 if (egrid1(i,j).lt.0.0) write(*,*) 'neg prate ij=',i,j,egrid1(i,j) 1066 ENDDO 1067 ENDDO 1068 write(*,*) '>>e2out' 1069 CALL E2OUT(167,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1070 ID(1:25) = 0 1071 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1072 CALL OUTPUT(IOUTYP,IGET(167),LVLS(1,IGET(167)), 1073 X GRID1,IMOUT,JMOUT) 1074 ENDIF 1075 C 1076 C ACCUMULATED TOTAL PRECIPITATION. 1077 IF (IGET(087).GT.0) THEN 1078 CALL E2OUT(087,000,ACPREC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1079 ID(1:25) = 0 1080 IFHR = NTSD/TSPH + 0.5 1081 ITPREC = INT(TPREC) 1082 IFINCR = MOD(IFHR,ITPREC) 1083 ID(19) = IFHR 1084 ID(20) = 4 1085 IF (IFINCR.EQ.0) THEN 1086 ID(18) = IFHR-ITPREC 1087 ELSE 1088 ID(18) = IFHR-IFINCR 1089 ENDIF 1090 IF (ID(18).LT.0) ID(18) = 0 1091 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1092 CALL OUTPUT(IOUTYP,IGET(087),LVLS(1,IGET(087)), 1093 X GRID1,IMOUT,JMOUT) 1094 ENDIF 1095 C 1096 C ACCUMULATED TOTAL PRECIPITATION NN - NEAREST NEIGHBOR. 1097 IF (IGET(246).GT.0) THEN 1098 CALL E2OUT(246,000,ACPREC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1099 ID(1:25) = 0 1100 IFHR = NTSD/TSPH + 0.5 1101 ITPREC = INT(TPREC) 1102 IFINCR = MOD(IFHR,ITPREC) 1103 ID(19) = IFHR 1104 ID(20) = 4 1105 IF (IFINCR.EQ.0) THEN 1106 ID(18) = IFHR-ITPREC 1107 ELSE 1108 ID(18) = IFHR-IFINCR 1109 ENDIF 1110 IF (ID(18).LT.0) ID(18) = 0 1111 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1112 CALL OUTPUT(IOUTYP,IGET(246),LVLS(1,IGET(246)), 1113 X GRID1,IMOUT,JMOUT) 1114 ENDIF 1115 C 1116 C ACCUMULATED CONVECTIVE PRECIPITATION. 1117 IF (IGET(033).GT.0) THEN 1118 CALL E2OUT(033,000,CUPREC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1119 ID(1:25) = 0 1120 IFHR = NTSD/TSPH + 0.5 1121 ITPREC = INT(TPREC) Page 16 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1122 IFINCR = MOD(IFHR,ITPREC) 1123 ID(19) = IFHR 1124 ID(20) = 4 1125 IF (IFINCR.EQ.0) THEN 1126 ID(18) = IFHR-ITPREC 1127 ELSE 1128 ID(18) = IFHR-IFINCR 1129 ENDIF 1130 IF (ID(18).LT.0) ID(18) = 0 1131 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1132 CALL OUTPUT(IOUTYP,IGET(033),LVLS(1,IGET(033)), 1133 X GRID1,IMOUT,JMOUT) 1134 ENDIF 1135 C 1136 C ACCUMULATED GRID-SCALE PRECIPITATION. 1137 IF (IGET(034).GT.0) THEN 1138 CALL E2OUT(034,000,ANCPRC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1139 ID(1:25) = 0 1140 IFHR = NTSD/TSPH + 0.5 1141 ITPREC = INT(TPREC) 1142 IFINCR = MOD(IFHR,ITPREC) 1143 ID(19) = IFHR 1144 ID(20) = 4 1145 IF (IFINCR.EQ.0) THEN 1146 ID(18) = IFHR-ITPREC 1147 ELSE 1148 ID(18) = IFHR-IFINCR 1149 ENDIF 1150 IF (ID(18).LT.0) ID(18) = 0 1151 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1152 CALL OUTPUT(IOUTYP,IGET(034),LVLS(1,IGET(034)), 1153 X GRID1,IMOUT,JMOUT) 1154 ENDIF 1155 C 1156 C ACCUMULATED SNOWFALL. 1157 IF (IGET(035).GT.0) THEN 1158 CALL E2OUT(035,000,ACSNOW,EGRID2, 1159 x GRID1,GRID2,IMOUT,JMOUT) 1160 ID(1:25) = 0 1161 IFHR = NTSD/TSPH + 0.5 1162 ITPREC = INT(TPREC) 1163 IFINCR = MOD(IFHR,ITPREC) 1164 ID(19) = IFHR 1165 ID(20) = 4 1166 IF (IFINCR.EQ.0) THEN 1167 ID(18) = IFHR-ITPREC 1168 ELSE 1169 ID(18) = IFHR-IFINCR 1170 ENDIF 1171 IF (ID(18).LT.0) ID(18) = 0 1172 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1173 CALL OUTPUT(IOUTYP,IGET(035),LVLS(1,IGET(035)), 1174 X GRID1,IMOUT,JMOUT) 1175 ENDIF 1176 C 1177 C ACCUMULATED SNOW MELT. 1178 IF (IGET(121).GT.0) THEN Page 17 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1179 CALL E2OUT(121,000,ACSNOM,EGRID2, 1180 x GRID1,GRID2,IMOUT,JMOUT) 1181 ID(1:25) = 0 1182 IFHR = NTSD/TSPH + 0.5 1183 ITPREC = INT(TPREC) 1184 IFINCR = MOD(IFHR,ITPREC) 1185 ID(19) = IFHR 1186 ID(20) = 4 1187 IF (IFINCR.EQ.0) THEN 1188 ID(18) = IFHR-ITPREC 1189 ELSE 1190 ID(18) = IFHR-IFINCR 1191 ENDIF 1192 IF (ID(18).LT.0) ID(18) = 0 1193 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1194 CALL OUTPUT(IOUTYP,IGET(121),LVLS(1,IGET(121)), 1195 X GRID1,IMOUT,JMOUT) 1196 ENDIF 1197 C 1198 C ACCUMULATED STORM SURFACE RUNOFF. 1199 IF (IGET(122).GT.0) THEN 1200 CALL E2OUT(122,000,SSROFF,EGRID2, 1201 x GRID1,GRID2,IMOUT,JMOUT) 1202 ID(1:25) = 0 1203 IFHR = NTSD/TSPH + 0.5 1204 ITPREC = INT(TPREC) 1205 IFINCR = MOD(IFHR,ITPREC) 1206 ID(19) = IFHR 1207 ID(20) = 4 1208 IF (IFINCR.EQ.0) THEN 1209 ID(18) = IFHR-ITPREC 1210 ELSE 1211 ID(18) = IFHR-IFINCR 1212 ENDIF 1213 IF (ID(18).LT.0) ID(18) = 0 1214 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1215 CALL OUTPUT(IOUTYP,IGET(122),LVLS(1,IGET(122)), 1216 X GRID1,IMOUT,JMOUT) 1217 ENDIF 1218 C 1219 C ACCUMULATED BASEFLOW-GROUNDWATER RUNOFF. 1220 IF (IGET(123).GT.0) THEN 1221 CALL E2OUT(123,000,BGROFF,EGRID2, 1222 x GRID1,GRID2,IMOUT,JMOUT) 1223 ID(1:25) = 0 1224 IFHR = NTSD/TSPH + 0.5 1225 ITPREC = INT(TPREC) 1226 IFINCR = MOD(IFHR,ITPREC) 1227 ID(19) = IFHR 1228 ID(20) = 4 1229 IF (IFINCR.EQ.0) THEN 1230 ID(18) = IFHR-ITPREC 1231 ELSE 1232 ID(18) = IFHR-IFINCR 1233 ENDIF 1234 IF (ID(18).LT.0) ID(18) = 0 1235 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) Page 18 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1236 CALL OUTPUT(IOUTYP,IGET(123),LVLS(1,IGET(123)), 1237 X GRID1,IMOUT,JMOUT) 1238 ENDIF 1239 C 1240 C INSTANTANEOUS PRECIPITATION TYPE. 1241 IF (IGET(160).GT.0) THEN 1242 CALL CALWXT(T,Q,RES,PD,HTM,LMH,PREC,PT,AETA,ETA,IWX1) 1243 ID(1:25) = 0 1244 C 1245 C DECOMPOSE IWX1 ARRAY 1246 C 1247 DO J=JSTA,JEND 1248 DO I=1,IM 1249 IWX=IWX1(I,J) 1250 ISNO=MOD(IWX,2) 1251 IIP=MOD(IWX,4)/2 1252 IZR=MOD(IWX,8)/4 1253 IRAIN=IWX/8 1254 SNOW(I,J) = ISNO*1.0 1255 SLEET(I,J) = IIP*1.0 1256 FREEZR(I,J) = IZR*1.0 1257 RAIN(I,J) = IRAIN*1.0 1258 ENDDO 1259 ENDDO 1260 C 1261 C INTERPOLATE/OUTPUT REQUESTED SURFACE FIELDS. 1262 C 1263 C SNOW. 1264 ID(8) = 143 1265 CALL E2OUT(160,000,SNOW,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1266 CALL OUTPUT(IOUTYP,IGET(160),LVLS(1,IGET(160)), 1267 X GRID1,IMOUT,JMOUT) 1268 C ICE PELLETS. 1269 ID(8) = 142 1270 CALL E2OUT(160,000,SLEET,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1271 CALL OUTPUT(IOUTYP,IGET(160),LVLS(1,IGET(160)), 1272 X GRID1,IMOUT,JMOUT) 1273 C FREEZING RAIN. 1274 ID(8) = 141 1275 CALL E2OUT(160,000,FREEZR,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1276 CALL OUTPUT(IOUTYP,IGET(160),LVLS(1,IGET(160)), 1277 X GRID1,IMOUT,JMOUT) 1278 C RAIN. 1279 ID(8) = 140 1280 CALL E2OUT(160,000,RAIN,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1281 CALL OUTPUT(IOUTYP,IGET(160),LVLS(1,IGET(160)), 1282 X GRID1,IMOUT,JMOUT) 1283 ENDIF 1284 C 1285 C 1286 C 1287 C*** BLOCK 5. SURFACE EXCHANGE FIELDS. 1288 C 1289 C TIME AVERAGED SURFACE LATENT HEAT FLUX. 1290 IF (IGET(042).GT.0) THEN 1291 IF(ASRFC.GT.0.)THEN 1292 RRNUM=1./ASRFC Page 19 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1293 ELSE 1294 RRNUM=0. 1295 ENDIF 1296 DO J=JSTA,JEND 1297 DO I=1,IM 1298 EGRID1(I,J) = SFCLHX(I,J)*RRNUM 1299 ENDDO 1300 ENDDO 1301 CALL E2OUT(042,000,EGRID1,EGRID2, 1302 X GRID1,GRID2,IMOUT,JMOUT) 1303 ID(1:25) = 0 1304 IFHR = NTSD/TSPH + 0.5 1305 ITSRFC = INT(TSRFC) 1306 IFINCR = MOD(IFHR,ITSRFC) 1307 ID(19) = IFHR 1308 ID(20) = 3 1309 IF (IFINCR.EQ.0) THEN 1310 ID(18) = IFHR-ITSRFC 1311 ELSE 1312 ID(18) = IFHR-IFINCR 1313 ENDIF 1314 IF (ID(18).LT.0) ID(18) = 0 1315 CALL OUTPUT(IOUTYP,IGET(042),LVLS(1,IGET(042)), 1316 X GRID1,IMOUT,JMOUT) 1317 ENDIF 1318 C 1319 C TIME AVERAGED SURFACE SENSIBLE HEAT FLUX. 1320 IF (IGET(043).GT.0) THEN 1321 IF(ASRFC.GT.0.)THEN 1322 RRNUM=1./ASRFC 1323 ELSE 1324 RRNUM=0. 1325 ENDIF 1326 DO J=JSTA,JEND 1327 DO I=1,IM 1328 EGRID1(I,J) = SFCSHX(I,J)*RRNUM 1329 ENDDO 1330 ENDDO 1331 CALL E2OUT(043,000,EGRID1,EGRID2, 1332 X GRID1,GRID2,IMOUT,JMOUT) 1333 ID(1:25) = 0 1334 IFHR = NTSD/TSPH + 0.5 1335 ITSRFC = INT(TSRFC) 1336 IFINCR = MOD(IFHR,ITSRFC) 1337 ID(19) = IFHR 1338 ID(20) = 3 1339 IF (IFINCR.EQ.0) THEN 1340 ID(18) = IFHR-ITSRFC 1341 ELSE 1342 ID(18) = IFHR-IFINCR 1343 ENDIF 1344 IF (ID(18).LT.0) ID(18) = 0 1345 CALL OUTPUT(IOUTYP,IGET(043),LVLS(1,IGET(043)), 1346 X GRID1,IMOUT,JMOUT) 1347 ENDIF 1348 C 1349 C TIME AVERAGED SUB-SURFACE SENSIBLE HEAT FLUX. Page 20 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1350 IF (IGET(135).GT.0) THEN 1351 IF(ASRFC.GT.0.)THEN 1352 RRNUM=1./ASRFC 1353 ELSE 1354 RRNUM=0. 1355 ENDIF 1356 DO J=JSTA,JEND 1357 DO I=1,IM 1358 EGRID1(I,J) = SUBSHX(I,J)*RRNUM 1359 ENDDO 1360 ENDDO 1361 CALL E2OUT(135,000,EGRID1,EGRID2, 1362 X GRID1,GRID2,IMOUT,JMOUT) 1363 ID(1:25) = 0 1364 IFHR = NTSD/TSPH + 0.5 1365 ITSRFC = INT(TSRFC) 1366 IFINCR = MOD(IFHR,ITSRFC) 1367 ID(19) = IFHR 1368 ID(20) = 3 1369 IF (IFINCR.EQ.0) THEN 1370 ID(18) = IFHR-ITSRFC 1371 ELSE 1372 ID(18) = IFHR-IFINCR 1373 ENDIF 1374 IF (ID(18).LT.0) ID(18) = 0 1375 CALL OUTPUT(IOUTYP,IGET(135),LVLS(1,IGET(135)), 1376 X GRID1,IMOUT,JMOUT) 1377 ENDIF 1378 C 1379 C TIME AVERAGED SNOW PHASE CHANGE HEAT FLUX. 1380 IF (IGET(136).GT.0) THEN 1381 IF(ASRFC.GT.0.)THEN 1382 RRNUM=1./ASRFC 1383 ELSE 1384 RRNUM=0. 1385 ENDIF 1386 DO J=JSTA,JEND 1387 DO I=1,IM 1388 EGRID1(I,J) = SNOPCX(I,J)*RRNUM 1389 ENDDO 1390 ENDDO 1391 CALL E2OUT(136,000,EGRID1,EGRID2, 1392 X GRID1,GRID2,IMOUT,JMOUT) 1393 ID(1:25) = 0 1394 IFHR = NTSD/TSPH + 0.5 1395 ITSRFC = INT(TSRFC) 1396 IFINCR = MOD(IFHR,ITSRFC) 1397 ID(19) = IFHR 1398 ID(20) = 3 1399 IF (IFINCR.EQ.0) THEN 1400 ID(18) = IFHR-ITSRFC 1401 ELSE 1402 ID(18) = IFHR-IFINCR 1403 ENDIF 1404 IF (ID(18).LT.0) ID(18) = 0 1405 CALL OUTPUT(IOUTYP,IGET(136),LVLS(1,IGET(136)), 1406 X GRID1,IMOUT,JMOUT) Page 21 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1407 ENDIF 1408 C 1409 C TIME AVERAGED SURFACE MOMENTUM FLUX. 1410 IF (IGET(046).GT.0) THEN 1411 IF(ASRFC.GT.0.)THEN 1412 RRNUM=1./ASRFC 1413 ELSE 1414 RRNUM=0. 1415 ENDIF 1416 DO J=JSTA,JEND 1417 DO I=1,IM 1418 EGRID1(I,J) = SFCUVX(I,J)*RRNUM 1419 ENDDO 1420 ENDDO 1421 CALL E2OUT(046,000,EGRID1,EGRID2, 1422 X GRID1,GRID2,IMOUT,JMOUT) 1423 ID(1:25) = 0 1424 IFHR = NTSD/TSPH + 0.5 1425 ITSRFC = INT(TSRFC) 1426 IFINCR = MOD(IFHR,ITSRFC) 1427 ID(19) = IFHR 1428 ID(20) = 3 1429 IF (IFINCR.EQ.0) THEN 1430 ID(18) = IFHR-ITSRFC 1431 ELSE 1432 ID(18) = IFHR-IFINCR 1433 ENDIF 1434 IF (ID(18).LT.0) ID(18) = 0 1435 CALL OUTPUT(IOUTYP,IGET(046),LVLS(1,IGET(046)), 1436 X GRID1,IMOUT,JMOUT) 1437 ENDIF 1438 C 1439 C ACCUMULATED SURFACE EVAPORATION 1440 IF (IGET(047).GT.0) THEN 1441 DO J=JSTA,JEND 1442 DO I=1,IM 1443 EGRID1(I,J) = SFCEVP(I,J) 1444 ENDDO 1445 ENDDO 1446 CALL E2OUT(047,000,EGRID1,EGRID2, 1447 X GRID1,GRID2,IMOUT,JMOUT) 1448 ID(1:25) = 0 1449 IFHR = NTSD/TSPH + 0.5 1450 ITPREC = INT(TPREC) 1451 IFINCR = MOD(IFHR,ITPREC) 1452 ID(19) = IFHR 1453 ID(20) = 4 1454 IF (IFINCR.EQ.0) THEN 1455 ID(18) = IFHR-ITPREC 1456 ELSE 1457 ID(18) = IFHR-IFINCR 1458 ENDIF 1459 IF (ID(18).LT.0) ID(18) = 0 1460 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1461 CALL OUTPUT(IOUTYP,IGET(047),LVLS(1,IGET(047)), 1462 X GRID1,IMOUT,JMOUT) 1463 ENDIF Page 22 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1464 C 1465 C ACCUMULATED POTENTIAL EVAPORATION 1466 IF (IGET(137).GT.0) THEN 1467 DO J=JSTA,JEND 1468 DO I=1,IM 1469 EGRID1(I,J) = POTEVP(I,J) 1470 ENDDO 1471 ENDDO 1472 CALL E2OUT(137,000,EGRID1,EGRID2, 1473 X GRID1,GRID2,IMOUT,JMOUT) 1474 ID(1:25) = 0 1475 IFHR = NTSD/TSPH + 0.5 1476 ITPREC = INT(TPREC) 1477 IFINCR = MOD(IFHR,ITPREC) 1478 ID(19) = IFHR 1479 ID(20) = 4 1480 IF (IFINCR.EQ.0) THEN 1481 ID(18) = IFHR-ITPREC 1482 ELSE 1483 ID(18) = IFHR-IFINCR 1484 ENDIF 1485 IF (ID(18).LT.0) ID(18) = 0 1486 CALL SCLFLD(GRID1,1000.,IMOUT,JMOUT) 1487 CALL OUTPUT(IOUTYP,IGET(137),LVLS(1,IGET(137)), 1488 X GRID1,IMOUT,JMOUT) 1489 ENDIF 1490 C 1491 C ROUGHNESS LENGTH. 1492 IF (IGET(044).GT.0) THEN 1493 CALL E2OUT(044,000,Z0,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1494 ID(1:25) = 0 1495 CALL OUTPUT(IOUTYP,IGET(044),LVLS(1,IGET(044)), 1496 X GRID1,IMOUT,JMOUT) 1497 ENDIF 1498 C 1499 C FRICTION VELOCITY. 1500 IF (IGET(045).GT.0) THEN 1501 CALL E2OUT(045,000,USTAR,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1502 ID(1:25) = 0 1503 CALL OUTPUT(IOUTYP,IGET(045),LVLS(1,IGET(045)), 1504 X GRID1,IMOUT,JMOUT) 1505 ENDIF 1506 C 1507 C SURFACE DRAG COEFFICIENT. 1508 IF (IGET(132).GT.0) THEN 1509 CALL CALDRG(EGRID1) 1510 CALL E2OUT(132,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1511 ID(1:25) = 0 1512 CALL OUTPUT(IOUTYP,IGET(132),LVLS(1,IGET(132)), 1513 X GRID1,IMOUT,JMOUT) 1514 ENDIF 1515 C 1516 C SURFACE U AND/OR V COMPONENT WIND STRESS 1517 IF ( (IGET(133).GT.0) .OR. (IGET(134).GT.0) ) THEN 1518 CALL CALTAU(EGRID1,EGRID2) 1519 C 1520 C SURFACE U COMPONENT WIND STRESS. Page 23 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1521 IF (IGET(133).GT.0) THEN 1522 CALL E2OUT(133,000,EGRID1,EGRID2, 1523 X GRID1,GRID2,IMOUT,JMOUT) 1524 ID(1:25) = 0 1525 CALL OUTPUT(IOUTYP,IGET(133),LVLS(1,IGET(133)), 1526 X GRID1,IMOUT,JMOUT) 1527 ENDIF 1528 C 1529 C SURFACE V COMPONENT WIND STRESS 1530 IF (IGET(134).GT.0) THEN 1531 CALL E2OUT(134,000,EGRID2,EGRID1, 1532 X GRID1,GRID2,IMOUT,JMOUT) 1533 ID(1:25) = 0 1534 CALL OUTPUT(IOUTYP,IGET(134),LVLS(1,IGET(134)), 1535 X GRID1,IMOUT,JMOUT) 1536 ENDIF 1537 ENDIF 1538 C 1539 C INSTANTANEOUS SENSIBLE HEAT FLUX 1540 IF (IGET(154).GT.0) THEN 1541 CALL E2OUT(154,000,TWBS,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1542 ID(1:25) = 0 1543 CALL OUTPUT(IOUTYP,IGET(154),LVLS(1,IGET(154)), 1544 X GRID1,IMOUT,JMOUT) 1545 ENDIF 1546 C 1547 C INSTANTANEOUS LATENT HEAT FLUX 1548 IF (IGET(155).GT.0) THEN 1549 CALL E2OUT(155,000,QWBS,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1550 ID(1:25) = 0 1551 CALL OUTPUT(IOUTYP,IGET(155),LVLS(1,IGET(155)), 1552 X GRID1,IMOUT,JMOUT) 1553 ENDIF 1554 C 1555 C SURFACE EXCHANGE COEFF 1556 IF (IGET(169).GT.0) THEN 1557 CALL E2OUT(169,000,SFCEXC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1558 ID(1:25) = 0 1559 CALL OUTPUT(IOUTYP,IGET(169),LVLS(1,IGET(169)), 1560 X GRID1,IMOUT,JMOUT) 1561 ENDIF 1562 C 1563 C GREEN VEG FRACTION 1564 IF (IGET(170).GT.0) THEN 1565 CALL E2OUT(170,000,VEGFRC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1566 ID(1:25) = 0 1567 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 1568 CALL OUTPUT(IOUTYP,IGET(170),LVLS(1,IGET(170)), 1569 X GRID1,IMOUT,JMOUT) 1570 ENDIF 1571 C 1572 C INSTANTANEOUS GROUND HEAT FLUX 1573 IF (IGET(152).GT.0) THEN 1574 CALL E2OUT(152,000,GRNFLX,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1575 ID(1:25) = 0 1576 CALL OUTPUT(IOUTYP,IGET(152),LVLS(1,IGET(152)), 1577 X GRID1,IMOUT,JMOUT) Page 24 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1578 ENDIF 1579 ! 1580 ! VEGETATION TYPE 1581 IF (IGET(218).GT.0) THEN 1582 DO J=JSTA,JEND 1583 DO I=1,IM 1584 EGRID1(I,J) = FLOAT(IVGTYP(I,J)) 1585 ENDDO 1586 ENDDO 1587 CALL E2OUT(218,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1588 ID(1:25) = 0 1589 CALL OUTPUT(IOUTYP,IGET(218),LVLS(1,IGET(218)), 1590 X GRID1,IMOUT,JMOUT) 1591 ENDIF 1592 ! 1593 ! SOIL TYPE 1594 IF (IGET(219).GT.0) THEN 1595 DO J=JSTA,JEND 1596 DO I=1,IM 1597 EGRID1(I,J) = FLOAT(ISLTYP(I,J)) 1598 ENDDO 1599 ENDDO 1600 CALL E2OUT(219,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1601 ID(1:25) = 0 1602 CALL OUTPUT(IOUTYP,IGET(219),LVLS(1,IGET(219)), 1603 X GRID1,IMOUT,JMOUT) 1604 ENDIF 1605 ! 1606 ! SLOPE TYPE 1607 IF (IGET(223).GT.0) THEN 1608 DO J=JSTA,JEND 1609 DO I=1,IM 1610 EGRID1(I,J) = FLOAT(ISLOPE(I,J)) 1611 ENDDO 1612 ENDDO 1613 CALL E2OUT(223,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1614 ID(1:25) = 0 1615 ID(02)= 130 1616 CALL OUTPUT(IOUTYP,IGET(223),LVLS(1,IGET(223)), 1617 X GRID1,IMOUT,JMOUT) 1618 ENDIF 1619 ! 1620 ! CANOPY CONDUCTANCE 1621 IF (IGET(220).GT.0 .OR. IGET(234).GT.0 1622 & .OR. IGET(235).GT.0 .OR. IGET(236).GT.0 1623 & .OR. IGET(237).GT.0 .OR. IGET(238).GT.0 1624 & .OR. IGET(239).GT.0 .OR. IGET(240).GT.0 1625 & .OR. IGET(241).GT.0) THEN 1626 DO J=JSTA,JEND 1627 DO I=1,IM 1628 ! IF(abs(SM(I,J)-0.).lt.1.0E-5)THEN 1629 IF( (abs(SM(I,J)-0.) .lt. 1.0E-5) .AND. 1630 & (abs(SICE(I,J)-0.) .lt. 1.0E-5) ) THEN 1631 IF(CZMEAN(I,J).GT.1.E-6) THEN 1632 FACTRS=CZEN(I,J)/CZMEAN(I,J) 1633 ELSE 1634 FACTRS=0.0 Page 25 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1635 ENDIF 1636 SOLAR=HBM2(I,J)*RSWIN(I,J)*FACTRS 1637 SFCTMP=T(I,J,LMH(I,J)) 1638 SFCQ=Q(I,J,LMH(I,J)) 1639 SFCPRS=PT+PDSL(I,J)*AETA(LMH(I,J)) 1640 ! IF(IVGTYP(I,J).EQ.0)PRINT*,'IVGTYP ZERO AT ',I,J 1641 ! & ,SM(I,J) 1642 IVG=IVGTYP(I,J) 1643 ! IF(IVGTYP(I,J).EQ.0)IVG=7 1644 ! CALL CANRES(SOLAR,SFCTMP,SFCQ,SFCPRS 1645 ! & ,SMC(I,J,1:NSOIL),GC(I,J),RC,IVG,ISLTYP(I,J)) 1646 ! 1647 CALL CANRES(SOLAR,SFCTMP,SFCQ,SFCPRS 1648 & ,SH2O(I,J,1:NSOIL),GC(I,J),RC,IVG,ISLTYP(I,J) 1649 & ,RSMIN(I,J),NROOTS(I,J),SMCWLT(I,J),SMCREF(I,J) 1650 & ,RCS(I,J),RCQ(I,J),RCT(I,J),RCSOIL(I,J)) 1651 IF(abs(SMCWLT(I,J)-0.5).lt.1.e-5)print*, 1652 & 'LARGE SMCWLT',i,j,SM(I,J),ISLTYP(I,J),SMCWLT(I,J) 1653 ELSE 1654 GC(I,J)=0. 1655 RSMIN(I,J)=0. 1656 NROOTS(I,J)=0 1657 SMCWLT(I,J)=0. 1658 SMCREF(I,J)=0. 1659 RCS(I,J)=0. 1660 RCQ(I,J)=0. 1661 RCT(I,J)=0. 1662 RCSOIL(I,J)=0. 1663 END IF 1664 ENDDO 1665 ENDDO 1666 IF (IGET(220).GT.0 )THEN 1667 CALL E2OUT(220,000,GC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1668 ID(1:25) = 0 1669 ID(02)= 130 1670 CALL OUTPUT(IOUTYP,IGET(220),LVLS(1,IGET(220)), 1671 X GRID1,IMOUT,JMOUT) 1672 END IF 1673 ! 1674 IF (IGET(234).GT.0 )THEN 1675 CALL E2OUT(234,000,RSMIN,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1676 ID(1:25) = 0 1677 ID(02)= 130 1678 CALL OUTPUT(IOUTYP,IGET(234),LVLS(1,IGET(234)), 1679 X GRID1,IMOUT,JMOUT) 1680 END IF 1681 ! 1682 ! NO OF ROOT LAYERS 1683 IF (IGET(235).GT.0) THEN 1684 DO J=JSTA,JEND 1685 DO I=1,IM 1686 EGRID1(I,J) = FLOAT(NROOTS(I,J)) 1687 ENDDO 1688 ENDDO 1689 CALL E2OUT(235,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1690 ID(1:25) = 0 1691 ID(02)= 130 Page 26 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1692 CALL OUTPUT(IOUTYP,IGET(235),LVLS(1,IGET(235)), 1693 X GRID1,IMOUT,JMOUT) 1694 ENDIF 1695 ! 1696 IF (IGET(236).GT.0 )THEN 1697 CALL E2OUT(236,000,SMCWLT,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1698 ID(1:25) = 0 1699 ID(02)= 130 1700 CALL OUTPUT(IOUTYP,IGET(236),LVLS(1,IGET(236)), 1701 X GRID1,IMOUT,JMOUT) 1702 END IF 1703 ! 1704 IF (IGET(237).GT.0 )THEN 1705 CALL E2OUT(237,000,SMCREF,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1706 ID(1:25) = 0 1707 ID(02)= 130 1708 CALL OUTPUT(IOUTYP,IGET(237),LVLS(1,IGET(237)), 1709 X GRID1,IMOUT,JMOUT) 1710 END IF 1711 ! 1712 IF (IGET(238).GT.0 )THEN 1713 CALL E2OUT(238,000,RCS,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1714 ID(1:25) = 0 1715 ID(02)= 130 1716 CALL OUTPUT(IOUTYP,IGET(238),LVLS(1,IGET(238)), 1717 X GRID1,IMOUT,JMOUT) 1718 END IF 1719 ! 1720 IF (IGET(239).GT.0 )THEN 1721 CALL E2OUT(239,000,RCT,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1722 ID(1:25) = 0 1723 ID(02)= 130 1724 CALL OUTPUT(IOUTYP,IGET(239),LVLS(1,IGET(239)), 1725 X GRID1,IMOUT,JMOUT) 1726 END IF 1727 ! 1728 IF (IGET(240).GT.0 )THEN 1729 CALL E2OUT(240,000,RCQ,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1730 ID(1:25) = 0 1731 ID(02)= 130 1732 CALL OUTPUT(IOUTYP,IGET(240),LVLS(1,IGET(240)), 1733 X GRID1,IMOUT,JMOUT) 1734 END IF 1735 ! 1736 IF (IGET(241).GT.0 )THEN 1737 CALL E2OUT(241,000,RCSOIL,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1738 ID(1:25) = 0 1739 ID(02)= 130 1740 CALL OUTPUT(IOUTYP,IGET(241),LVLS(1,IGET(241)), 1741 X GRID1,IMOUT,JMOUT) 1742 END IF 1743 ! 1744 ENDIF 1745 C 1746 C DQADV 1747 IF (IGET(200).GT.0) THEN 1748 CALL E2OUT(200,000,DQADV,EGRID2,GRID1,GRID2,IMOUT,JMOUT) Page 27 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1749 ID(1:25) = 0 1750 CALL OUTPUT(IOUTYP,IGET(200),LVLS(1,IGET(200)), 1751 X GRID1,IMOUT,JMOUT) 1752 ENDIF 1753 C 1754 C DQFLX / WVCONV 1755 IF (IGET(201).GT.0) THEN 1756 CALL E2OUT(201,000,DQFLX,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1757 ID(1:25) = 0 1758 IFHR = NTSD/TSPH + 0.5 1759 ITPREC = INT(TPREC) 1760 IFINCR = MOD(IFHR,ITPREC) 1761 ID(19) = IFHR 1762 ID(20) = 4 1763 IF (IFINCR.EQ.0) THEN 1764 ID(18) = IFHR-ITPREC 1765 ELSE 1766 ID(18) = IFHR-IFINCR 1767 ENDIF 1768 IF (ID(18).LT.0) ID(18) = 0 1769 CALL OUTPUT(IOUTYP,IGET(201),LVLS(1,IGET(201)), 1770 X GRID1,IMOUT,JMOUT) 1771 1772 DO J=JSTA,JEND 1773 DO I=1,IM 1774 EGRID1(I,J)=DQFLX(I,J)-DQFLX7(I,J) 1775 ENDDO 1776 ENDDO 1777 1778 CALL E2OUT(201,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1779 ID(1:25) = 0 1780 IFHR = NTSD/TSPH + 0.5 1781 ITPREC = INT(TPREC) 1782 IFINCR = MOD(IFHR,ITPREC) 1783 ID(19) = IFHR 1784 ID(20) = 4 1785 IF (IFINCR.EQ.0) THEN 1786 ID(18) = IFHR-ITPREC 1787 ELSE 1788 ID(18) = IFHR-IFINCR 1789 ENDIF 1790 IF (ID(18).LT.0) ID(18) = 0 1791 ID(9) = 101 1792 ID(10) = 0 1793 ID(11) = 70 1794 CALL OUTPUT(IOUTYP,IGET(201),LVLS(1,IGET(201)), 1795 X GRID1,IMOUT,JMOUT) 1796 ENDIF 1797 C 1798 C DCFLX / WCCONV 1799 IF (IGET(202).GT.0) THEN 1800 CALL E2OUT(202,000,DCFLX,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1801 ID(1:25) = 0 1802 IFHR = NTSD/TSPH + 0.5 1803 ITPREC = INT(TPREC) 1804 IFINCR = MOD(IFHR,ITPREC) 1805 ID(19) = IFHR Page 28 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1806 ID(20) = 4 1807 IF (IFINCR.EQ.0) THEN 1808 ID(18) = IFHR-ITPREC 1809 ELSE 1810 ID(18) = IFHR-IFINCR 1811 ENDIF 1812 IF (ID(18).LT.0) ID(18) = 0 1813 CALL OUTPUT(IOUTYP,IGET(202),LVLS(1,IGET(202)), 1814 X GRID1,IMOUT,JMOUT) 1815 1816 DO J=JSTA,JEND 1817 DO I=1,IM 1818 EGRID1(I,J)=DCFLX(I,J)-DCFLX7(I,J) 1819 ENDDO 1820 ENDDO 1821 1822 CALL E2OUT(202,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1823 ID(1:25) = 0 1824 IFHR = NTSD/TSPH + 0.5 1825 ITPREC = INT(TPREC) 1826 IFINCR = MOD(IFHR,ITPREC) 1827 ID(19) = IFHR 1828 ID(20) = 4 1829 IF (IFINCR.EQ.0) THEN 1830 ID(18) = IFHR-ITPREC 1831 ELSE 1832 ID(18) = IFHR-IFINCR 1833 ENDIF 1834 IF (ID(18).LT.0) ID(18) = 0 1835 ID(9) = 101 1836 ID(10) = 0 1837 ID(11) = 70 1838 CALL OUTPUT(IOUTYP,IGET(202),LVLS(1,IGET(202)), 1839 X GRID1,IMOUT,JMOUT) 1840 1841 ENDIF 1842 C 1843 C FQU / WVUFLX 1844 IF (IGET(203).GT.0) THEN 1845 CALL E2OUT(203,000,FQU,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1846 ID(1:25) = 0 1847 IFHR = NTSD/TSPH + 0.5 1848 ITPREC = INT(TPREC) 1849 IFINCR = MOD(IFHR,ITPREC) 1850 ID(19) = IFHR 1851 ID(20) = 4 1852 IF (IFINCR.EQ.0) THEN 1853 ID(18) = IFHR-ITPREC 1854 ELSE 1855 ID(18) = IFHR-IFINCR 1856 ENDIF 1857 IF (ID(18).LT.0) ID(18) = 0 1858 CALL OUTPUT(IOUTYP,IGET(203),LVLS(1,IGET(203)), 1859 X GRID1,IMOUT,JMOUT) 1860 1861 DO J=JSTA,JEND 1862 DO I=1,IM Page 29 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1863 EGRID1(I,J)=FQU(I,J)-FQU7(I,J) 1864 ENDDO 1865 ENDDO 1866 1867 CALL E2OUT(203,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1868 ID(1:25) = 0 1869 IFHR = NTSD/TSPH + 0.5 1870 ITPREC = INT(TPREC) 1871 IFINCR = MOD(IFHR,ITPREC) 1872 ID(19) = IFHR 1873 ID(20) = 4 1874 IF (IFINCR.EQ.0) THEN 1875 ID(18) = IFHR-ITPREC 1876 ELSE 1877 ID(18) = IFHR-IFINCR 1878 ENDIF 1879 IF (ID(18).LT.0) ID(18) = 0 1880 ID(9) = 101 1881 ID(10) = 0 1882 ID(11) = 70 1883 CALL OUTPUT(IOUTYP,IGET(203),LVLS(1,IGET(203)), 1884 X GRID1,IMOUT,JMOUT) 1885 1886 ENDIF 1887 C 1888 C FQV / WVVFLX 1889 IF (IGET(204).GT.0) THEN 1890 CALL E2OUT(204,000,FQV,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1891 ID(1:25) = 0 1892 IFHR = NTSD/TSPH + 0.5 1893 ITPREC = INT(TPREC) 1894 IFINCR = MOD(IFHR,ITPREC) 1895 ID(19) = IFHR 1896 ID(20) = 4 1897 IF (IFINCR.EQ.0) THEN 1898 ID(18) = IFHR-ITPREC 1899 ELSE 1900 ID(18) = IFHR-IFINCR 1901 ENDIF 1902 IF (ID(18).LT.0) ID(18) = 0 1903 CALL OUTPUT(IOUTYP,IGET(204),LVLS(1,IGET(204)), 1904 X GRID1,IMOUT,JMOUT) 1905 1906 DO J=JSTA,JEND 1907 DO I=1,IM 1908 EGRID1(I,J)=FQV(I,J)-FQV7(I,J) 1909 ENDDO 1910 ENDDO 1911 1912 CALL E2OUT(204,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1913 ID(1:25) = 0 1914 IFHR = NTSD/TSPH + 0.5 1915 ITPREC = INT(TPREC) 1916 IFINCR = MOD(IFHR,ITPREC) 1917 ID(19) = IFHR 1918 ID(20) = 4 1919 IF (IFINCR.EQ.0) THEN Page 30 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1920 ID(18) = IFHR-ITPREC 1921 ELSE 1922 ID(18) = IFHR-IFINCR 1923 ENDIF 1924 IF (ID(18).LT.0) ID(18) = 0 1925 ID(9) = 101 1926 ID(10) = 0 1927 ID(11) = 70 1928 CALL OUTPUT(IOUTYP,IGET(204),LVLS(1,IGET(204)), 1929 X GRID1,IMOUT,JMOUT) 1930 1931 ENDIF 1932 C 1933 C FCU / WCUFLX 1934 IF (IGET(205).GT.0) THEN 1935 CALL E2OUT(205,000,FCU,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1936 ID(1:25) = 0 1937 IFHR = NTSD/TSPH + 0.5 1938 ITPREC = INT(TPREC) 1939 IFINCR = MOD(IFHR,ITPREC) 1940 ID(19) = IFHR 1941 ID(20) = 4 1942 IF (IFINCR.EQ.0) THEN 1943 ID(18) = IFHR-ITPREC 1944 ELSE 1945 ID(18) = IFHR-IFINCR 1946 ENDIF 1947 IF (ID(18).LT.0) ID(18) = 0 1948 CALL OUTPUT(IOUTYP,IGET(205),LVLS(1,IGET(205)), 1949 X GRID1,IMOUT,JMOUT) 1950 1951 DO J=JSTA,JEND 1952 DO I=1,IM 1953 EGRID1(I,J)=FCU(I,J)-FCU7(I,J) 1954 ENDDO 1955 ENDDO 1956 1957 CALL E2OUT(205,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1958 ID(1:25) = 0 1959 IFHR = NTSD/TSPH + 0.5 1960 ITPREC = INT(TPREC) 1961 IFINCR = MOD(IFHR,ITPREC) 1962 ID(19) = IFHR 1963 ID(20) = 4 1964 IF (IFINCR.EQ.0) THEN 1965 ID(18) = IFHR-ITPREC 1966 ELSE 1967 ID(18) = IFHR-IFINCR 1968 ENDIF 1969 IF (ID(18).LT.0) ID(18) = 0 1970 ID(9) = 101 1971 ID(10) = 0 1972 ID(11) = 70 1973 CALL OUTPUT(IOUTYP,IGET(205),LVLS(1,IGET(205)), 1974 X GRID1,IMOUT,JMOUT) 1975 1976 ENDIF Page 31 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 1977 C 1978 C FCV / WCVFLX 1979 IF (IGET(206).GT.0) THEN 1980 CALL E2OUT(206,000,FCV,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1981 ID(1:25) = 0 1982 IFHR = NTSD/TSPH + 0.5 1983 ITPREC = INT(TPREC) 1984 IFINCR = MOD(IFHR,ITPREC) 1985 ID(19) = IFHR 1986 ID(20) = 4 1987 IF (IFINCR.EQ.0) THEN 1988 ID(18) = IFHR-ITPREC 1989 ELSE 1990 ID(18) = IFHR-IFINCR 1991 ENDIF 1992 IF (ID(18).LT.0) ID(18) = 0 1993 CALL OUTPUT(IOUTYP,IGET(206),LVLS(1,IGET(206)), 1994 X GRID1,IMOUT,JMOUT) 1995 1996 DO J=JSTA,JEND 1997 DO I=1,IM 1998 EGRID1(I,J)=FCV(I,J)-FCV7(I,J) 1999 ENDDO 2000 ENDDO 2001 2002 CALL E2OUT(206,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 2003 ID(1:25) = 0 2004 IFHR = NTSD/TSPH + 0.5 2005 ITPREC = INT(TPREC) 2006 IFINCR = MOD(IFHR,ITPREC) 2007 ID(19) = IFHR 2008 ID(20) = 4 2009 IF (IFINCR.EQ.0) THEN 2010 ID(18) = IFHR-ITPREC 2011 ELSE 2012 ID(18) = IFHR-IFINCR 2013 ENDIF 2014 IF (ID(18).LT.0) ID(18) = 0 2015 ID(9) = 101 2016 ID(10) = 0 2017 ID(11) = 70 2018 CALL OUTPUT(IOUTYP,IGET(206),LVLS(1,IGET(206)), 2019 X GRID1,IMOUT,JMOUT) 2020 2021 ENDIF 2022 2023 2024 C 2025 C VAPINC / WVINC 2026 IF (IGET(207).GT.0) THEN 2027 CALL E2OUT(207,000,VAPINC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 2028 ID(1:25) = 0 2029 IFHR = NTSD/TSPH + 0.5 2030 ITPREC = INT(TPREC) 2031 IFINCR = MOD(IFHR,ITPREC) 2032 ID(19) = IFHR 2033 ID(20) = 4 Page 32 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 2034 IF (IFINCR.EQ.0) THEN 2035 ID(18) = IFHR-ITPREC 2036 ELSE 2037 ID(18) = IFHR-IFINCR 2038 ENDIF 2039 IF (ID(18).LT.0) ID(18) = 0 2040 CALL OUTPUT(IOUTYP,IGET(207),LVLS(1,IGET(207)), 2041 X GRID1,IMOUT,JMOUT) 2042 2043 DO J=JSTA,JEND 2044 DO I=1,IM 2045 EGRID1(I,J)=VAPINC(I,J)-VAPINC7(I,J) 2046 ENDDO 2047 ENDDO 2048 2049 CALL E2OUT(207,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 2050 ID(1:25) = 0 2051 IFHR = NTSD/TSPH + 0.5 2052 ITPREC = INT(TPREC) 2053 IFINCR = MOD(IFHR,ITPREC) 2054 ID(19) = IFHR 2055 ID(20) = 4 2056 IF (IFINCR.EQ.0) THEN 2057 ID(18) = IFHR-ITPREC 2058 ELSE 2059 ID(18) = IFHR-IFINCR 2060 ENDIF 2061 IF (ID(18).LT.0) ID(18) = 0 2062 ID(9) = 101 2063 ID(10) = 0 2064 ID(11) = 70 2065 CALL OUTPUT(IOUTYP,IGET(207),LVLS(1,IGET(207)), 2066 X GRID1,IMOUT,JMOUT) 2067 2068 ENDIF 2069 C 2070 C CLDINC / WCINC 2071 IF (IGET(208).GT.0) THEN 2072 CALL E2OUT(208,000,CLDINC,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 2073 ID(1:25) = 0 2074 IFHR = NTSD/TSPH + 0.5 2075 ITPREC = INT(TPREC) 2076 IFINCR = MOD(IFHR,ITPREC) 2077 ID(19) = IFHR 2078 ID(20) = 4 2079 IF (IFINCR.EQ.0) THEN 2080 ID(18) = IFHR-ITPREC 2081 ELSE 2082 ID(18) = IFHR-IFINCR 2083 ENDIF 2084 IF (ID(18).LT.0) ID(18) = 0 2085 CALL OUTPUT(IOUTYP,IGET(208),LVLS(1,IGET(208)), 2086 X GRID1,IMOUT,JMOUT) 2087 2088 DO J=JSTA,JEND 2089 DO I=1,IM 2090 EGRID1(I,J)=CLDINC(I,J)-CLDINC7(I,J) Page 33 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f 2091 ENDDO 2092 ENDDO 2093 2094 CALL E2OUT(208,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 2095 ID(1:25) = 0 2096 IFHR = NTSD/TSPH + 0.5 2097 ITPREC = INT(TPREC) 2098 IFINCR = MOD(IFHR,ITPREC) 2099 ID(19) = IFHR 2100 ID(20) = 4 2101 IF (IFINCR.EQ.0) THEN 2102 ID(18) = IFHR-ITPREC 2103 ELSE 2104 ID(18) = IFHR-IFINCR 2105 ENDIF 2106 IF (ID(18).LT.0) ID(18) = 0 2107 ID(9) = 101 2108 ID(10) = 0 2109 ID(11) = 70 2110 CALL OUTPUT(IOUTYP,IGET(208),LVLS(1,IGET(208)), 2111 X GRID1,IMOUT,JMOUT) 2112 2113 ENDIF 2114 2115 2116 2117 C END OF ROUTINE 2118 print*,'End of SURFCE2' 2119 C 2120 RETURN 2121 END Page 34 Source Listing SURFCE2 2014-12-17 20:48 Entry Points SURFCE2.f ENTRY POINTS Name surfce2_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 40 Label 253 182,183 A2 Param 41 R(4) 4 scalar 219 A3 Param 41 R(4) 4 scalar 219 A4 Param 41 R(4) 4 scalar 219 ABS Func 495 scalar 495,496,1363,1364,1385 ACCLIQ Scalar 7 R(4) 4 2 91719 COM ACMPRE Common 2 1467508 ACMRDS Common 2 2201268 ACMSFC Common 2 2568144 ACPREC Scalar 7 R(4) 4 2 91719 COM 249,812,832 ACSNOM Scalar 3 R(4) 4 2 91719 COM 913 ACSNOW Scalar 3 R(4) 4 2 91719 COM 892 AD05 Param 33 R(4) 4 scalar AETA Scalar 8 R(4) 4 1 45 COM 976,1373 AETA2 Scalar 7 R(4) 4 1 45 COM AI Local 212 R(4) 4 scalar 212,215,220 AKHS Scalar 3 R(4) 4 2 91719 COM AKMS Scalar 3 R(4) 4 2 91719 COM ALBASE Scalar 14 R(4) 4 2 91719 COM ALBEDO Scalar 15 R(4) 4 2 91719 COM ALPINT Scalar 4 R(4) 4 3 4219074 COM ALSL Scalar 7 R(4) 4 1 39 COM AMAX1 Func 207 scalar 207 AMIN1 Func 211 scalar 211,220 ANCPRC Local 105 R(4) 4 2 91719 249,872 APREC Scalar 3 R(4) 4 2 91719 COM ARDSW Scalar 3 R(4) 4 scalar COM ASRFC Scalar 3 R(4) 4 scalar COM 1025,1026,1055,1056,1085,1086,1115 ,1116,1145,1146 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 BETA Param 37 R(4) 4 scalar 37 BGROFF Scalar 4 R(4) 4 2 91719 COM 955 BI Local 213 R(4) 4 scalar 213,216,220 BOUND Subr 276 276,285,310,329,467,584,616 BTG Param 37 R(4) 4 scalar C2K Param 96 R(4) 4 scalar CALDRG Subr 1243 1243 CALRH2 Subr 613 613 CALTAU Subr 1252 1252 CALWXT Subr 976 976 CANRES Subr 1381 1381 Page 35 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References CAPA Param 28 R(4) 4 scalar 35,189,190,200 CFRACH Scalar 3 R(4) 4 2 91719 COM CFRACL Scalar 3 R(4) 4 2 91719 COM CFRACM Scalar 3 R(4) 4 2 91719 COM CFT0 Param 33 R(4) 4 scalar CI Scalar 4 R(4) 4 scalar COM CLDEFI Scalar 6 R(4) 4 2 91719 COM CLDINC Scalar 6 R(4) 4 2 91719 COM 1806,1824 CLDINC7 Scalar 6 R(4) 4 2 91719 COM 1824 CLDWTR Common 2 17610408 CLIMIT Local 145 R(4) 4 scalar 145,153 CM1 Param 23 R(4) 4 scalar CM2 Param 23 R(4) 4 scalar CM3 Param 23 R(4) 4 scalar CMC Scalar 6 R(4) 4 2 91719 COM 440,498 CMLD Scalar 6 R(4) 4 scalar COM CNVCLD Common 2 1467504 CP Param 34 R(4) 4 scalar 36,40 CPDR Param 36 R(4) 4 scalar CPHI0 Scalar 9 R(4) 4 scalar COM CS Scalar 4 R(4) 4 scalar COM CTLBLK Common 2 96 CUPPT Scalar 3 R(4) 4 2 91719 COM CUPREC Scalar 7 R(4) 4 2 91719 COM 249,852 CWM Scalar 3 R(4) 4 3 4127355 COM 153 CZEN Scalar 12 R(4) 4 2 91719 COM 1366 CZMEAN Scalar 16 R(4) 4 2 91719 COM 1365,1366 D00 Param 25 R(4) 4 scalar 241,276,285,467 D001 Param 25 R(4) 4 scalar 244,598 D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar 241 D05 Param 29 R(4) 4 scalar 33 D1 Param 25 R(4) 4 scalar D115 Param 29 R(4) 4 scalar D125 Param 27 R(4) 4 scalar D25 Param 27 R(4) 4 scalar D50 Param 25 R(4) 4 scalar D608 Param 29 R(4) 4 scalar DATSET Scalar 3 CHAR 6 scalar COM DBOT Local 355 R(4) 4 scalar 355,357,376,378,396,398 DCFLX Scalar 15 R(4) 4 2 91719 COM 1534,1552 DCFLX7 Scalar 17 R(4) 4 2 91719 COM 1552 DEC Scalar 8 R(4) 4 1 250 COM DETA Scalar 8 R(4) 4 1 45 COM DETA2 Scalar 7 R(4) 4 1 45 COM DEWPOINT Subr 318 318,601 DFRLG Scalar 7 R(4) 4 1 46 COM DI Scalar 4 R(4) 4 scalar 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 DQADV Scalar 19 R(4) 4 2 91719 COM 1482 DQFLX Scalar 14 R(4) 4 2 91719 COM 1490,1508 Page 36 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References DQFLX7 Scalar 16 R(4) 4 2 91719 COM 1508 DS Scalar 4 R(4) 4 scalar COM DT Scalar 4 R(4) 4 scalar COM DTD Scalar 3 R(4) 4 scalar COM DTOP Local 351 R(4) 4 scalar 351,353,355,356,372,374,376,377,39 2,394,396,397 DTQ2 Scalar 3 R(4) 4 scalar COM 794 DTR Param 34 R(4) 4 scalar 34 DWPSFC Local 104 R(4) 4 2 91719 318,319 E2OUT Subr 259 259,267,275,284,293,301,309,319,32 7,349,370,390,408,418,429,440,449, 466,474,512,519,526,533,541,549,57 2,583,602,614,627,650,665,677,689, 701,718,733,745,757,769,787,803,81 2,832,852,872,892,913,934,955,999, 1004,1009,1014,1035,1065,1095,1125 ,1155,1180,1206,1227,1235,1244,125 6,1265,1275,1283,1291,1299,1308,13 21,1334,1347,1401,1409,1423,1431,1 439,1447,1455,1463,1471,1482,1490, 1512,1534,1556,1579,1601,1624,1646 ,1669,1691,1714,1736,1761,1783,180 6,1828 ECAN Local 110 R(4) 4 2 91719 499,502,512 EDIR Local 110 R(4) 4 2 91719 499,503,519 EGRDI Scalar 9 R(4) 4 2 91719 COM EGRDJ Scalar 9 R(4) 4 2 91719 COM EGRID1 Local 106 R(4) 4 2 91719 346,349,367,370,387,390,583,601,60 2,613,614,627,784,787,798,799,803, 1032,1035,1062,1065,1092,1095,1122 ,1125,1152,1155,1177,1180,1203,120 6,1243,1244,1252,1256,1265,1318,13 21,1331,1334,1344,1347,1420,1423,1 508,1512,1552,1556,1597,1601,1642, 1646,1687,1691,1732,1736,1779,1783 ,1824,1828 EGRID2 Local 106 R(4) 4 2 91719 259,267,275,284,293,301,309,319,32 7,349,370,390,408,418,429,440,449, 466,474,512,519,526,533,541,549,57 2,602,614,665,677,689,701,733,745, 757,769,787,803,812,832,852,872,89 2,913,934,955,999,1004,1009,1014,1 035,1065,1095,1125,1155,1180,1206, 1227,1235,1244,1252,1256,1265,1275 ,1283,1291,1299,1308,1321,1334,134 7,1401,1409,1423,1431,1439,1447,14 55,1463,1471,1482,1490,1512,1534,1 556,1579,1601,1624,1646,1669,1691, 1714,1736,1761,1783,1806,1828 ELDR Param 36 R(4) 4 scalar ELWDWN Scalar 13 R(4) 4 2 91719 COM ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39,243,597 EPSQ2 Param 40 R(4) 4 scalar EPSR Scalar 11 R(4) 4 2 91719 COM Page 37 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References ERAD Param 34 R(4) 4 scalar ERLAM0 Scalar 9 R(4) 4 scalar COM ESNOW Local 110 R(4) 4 2 91719 499,505,533 ETA Scalar 8 R(4) 4 1 46 COM 976 ETCALC Subr 497 497 ETRANS Local 110 R(4) 4 2 91719 499,504,526 EVLAT Scalar 8 R(4) 4 2 91762 COM EVLON Scalar 8 R(4) 4 2 91762 COM EVP Local 104 R(4) 4 2 91719 243,244,318,597,598,601 EXP Func 187 scalar 187,188,219 EXTRA Common 3 71506268 FACTRS Local 1366 R(4) 4 scalar 1366,1368,1370 FCNEV Scalar 5 R(4) 4 2 91719 COM FCNEV7 Scalar 6 R(4) 4 2 91719 COM FCSEV Scalar 5 R(4) 4 2 91719 COM FCSEV7 Scalar 6 R(4) 4 2 91719 COM FCU Scalar 10 R(4) 4 2 91719 COM 1669,1687 FCU7 Scalar 12 R(4) 4 2 91719 COM 1687 FCV Scalar 10 R(4) 4 2 91719 COM 1714,1732 FCV7 Scalar 12 R(4) 4 2 91719 COM 1732 FIELD Scalar 2 CHAR 20 1 250 COM FIRST Scalar 101 L(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM 193 FLOAT Func 1318 scalar 1318,1331,1344,1420 FMD Param 39 R(4) 4 scalar 39 FMW Param 39 R(4) 4 scalar 39 FQNEV Scalar 2 R(4) 4 2 91719 COM FQNEV1 Scalar 20 R(4) 4 2 91719 COM FQNEV7 Scalar 3 R(4) 4 2 91719 COM FQSEV Scalar 2 R(4) 4 2 91719 COM FQSEV1 Scalar 20 R(4) 4 2 91719 COM FQSEV7 Scalar 3 R(4) 4 2 91719 COM FQU Scalar 9 R(4) 4 2 91719 COM 1579,1597 FQU7 Scalar 11 R(4) 4 2 91719 COM 1597 FQV Scalar 9 R(4) 4 2 91719 COM 1624,1642 FQV7 Scalar 11 R(4) 4 2 91719 COM 1642 FREEZR Local 109 R(4) 4 2 91719 990,1009 FSL Scalar 8 R(4) 4 2 91719 COM FVTLON Scalar 7 R(4) 4 2 91762 COM G Param 32 R(4) 4 scalar 32,187,188 GC Local 107 R(4) 4 2 91719 1382,1388,1401 GDLAT Scalar 3 R(4) 4 2 91719 COM GDLON Scalar 3 R(4) 4 2 91719 COM GDTLAT Scalar 4 R(4) 4 2 91719 COM GDTLON Scalar 4 R(4) 4 2 91719 COM GFFC Scalar 14 R(4) 4 2 91719 COM GI Param 32 R(4) 4 scalar 36,193 GLAT Scalar 12 R(4) 4 2 91719 COM GLON Scalar 12 R(4) 4 2 91719 COM GRID1 Local 108 R(4) 4 2 0 259,262,267,270,275,276,279,284,28 5,288,293,296,301,304,309,310,313, 319,322,327,328,329,332,349,358,37 0,379,390,400,408,413,418,422,424, 429,433,435,440,442,444,449,451,45 3,466,467,470,474,477,512,515,519, Page 38 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References 522,526,529,533,537,541,545,549,55 3,572,578,583,584,590,602,608,614, 615,616,622,627,633,650,652,665,66 7,677,679,689,691,701,703,718,720, 733,735,745,747,757,759,769,771,78 7,790,803,805,807,812,825,827,832, 845,847,852,865,867,872,885,887,89 3,906,908,914,927,929,935,948,950, 956,969,971,999,1001,1004,1006,100 9,1011,1014,1016,1036,1050,1066,10 80,1096,1110,1126,1140,1156,1170,1 181,1194,1196,1207,1220,1222,1227, 1230,1235,1238,1244,1247,1257,1260 ,1266,1269,1275,1278,1283,1286,129 1,1294,1299,1301,1303,1308,1311,13 21,1324,1334,1337,1347,1351,1401,1 405,1409,1413,1423,1427,1431,1435, 1439,1443,1447,1451,1455,1459,1463 ,1467,1471,1475,1482,1485,1490,150 4,1512,1529,1534,1548,1556,1573,15 79,1593,1601,1618,1624,1638,1646,1 663,1669,1683,1691,1708,1714,1728, 1736,1753,1761,1775,1783,1800,1806 ,1820,1828,1845 GRID2 Local 108 R(4) 4 2 0 259,267,275,284,293,301,309,319,32 7,349,370,390,408,418,429,440,449, 466,474,512,519,526,533,541,549,57 2,583,602,614,627,650,654,665,677, 689,701,718,722,733,745,757,769,78 7,803,812,832,852,872,893,914,935, 956,999,1004,1009,1014,1036,1066,1 096,1126,1156,1181,1207,1227,1235, 1244,1257,1266,1275,1283,1291,1299 ,1308,1321,1334,1347,1401,1409,142 3,1431,1439,1447,1455,1463,1471,14 82,1490,1512,1534,1556,1579,1601,1 624,1646,1669,1691,1714,1736,1761, 1783,1806,1828 GRNFLX Scalar 4 R(4) 4 2 91719 COM 1308 H1 Param 26 R(4) 4 scalar 33,240,329,467,616 H100 Param 27 R(4) 4 scalar 328,329,615,616 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 207,310,584 H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar H2 Param 26 R(4) 4 scalar H4 Param 26 R(4) 4 scalar H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar 276,285,310,584 HBM2 Scalar 3 R(4) 4 2 91719 COM 1370 HBOT Scalar 13 R(4) 4 2 91719 COM Page 39 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References HDAC Scalar 16 R(4) 4 2 91719 COM HDACV Scalar 16 R(4) 4 2 91719 COM HTLAT Scalar 5 R(4) 4 2 91719 COM HTLON Scalar 5 R(4) 4 2 91719 COM HTM Scalar 5 R(4) 4 3 4127355 COM 976 HTOP Scalar 13 R(4) 4 2 91719 COM I Local 150 I(4) 4 scalar 150,151,152,153,154,155,156,157,15 9,162,183,186,187,188,189,190,193, 196,199,200,206,207,208,218,227,23 8,240,241,242,243,244,249,345,346, 366,367,386,387,459,460,461,463,49 1,495,496,497,498,499,500,502,503, 504,505,506,507,596,597,598,783,78 4,797,798,799,982,983,988,989,990, 991,1031,1032,1061,1062,1091,1092, 1121,1122,1151,1152,1176,1177,1202 ,1203,1317,1318,1330,1331,1343,134 4,1361,1363,1364,1365,1366,1370,13 71,1372,1373,1376,1382,1383,1384,1 385,1386,1388,1389,1390,1391,1392, 1393,1394,1395,1396,1419,1420,1507 ,1508,1551,1552,1596,1597,1641,164 2,1686,1687,1731,1732,1778,1779,18 23,1824 ICNT Scalar 7 I(4) 4 1 1024 COM ID Scalar 11 I(4) 4 1 25 COM 260,268,277,286,294,302,311,320,33 0,350,356,357,371,377,378,391,397, 398,399,409,411,419,420,421,430,43 1,432,441,450,468,475,513,520,527, 534,535,542,543,550,551,573,575,57 6,585,587,588,603,605,606,617,619, 620,628,630,631,643,645,646,661,66 3,664,673,675,676,685,687,688,697, 699,700,711,713,714,729,731,732,74 1,743,744,753,755,756,765,767,768, 788,804,813,817,818,820,822,824,83 3,837,838,840,842,844,853,857,858, 860,862,864,873,877,878,880,882,88 4,894,898,899,901,903,905,915,919, 920,922,924,926,936,940,941,943,94 5,947,957,961,962,964,966,968,977, 998,1003,1008,1013,1037,1041,1042, 1044,1046,1048,1067,1071,1072,1074 ,1076,1078,1097,1101,1102,1104,110 6,1108,1127,1131,1132,1134,1136,11 38,1157,1161,1162,1164,1166,1168,1 182,1186,1187,1189,1191,1193,1208, 1212,1213,1215,1217,1219,1228,1236 ,1245,1258,1267,1276,1284,1292,130 0,1309,1322,1335,1348,1349,1402,14 03,1410,1411,1424,1425,1432,1433,1 440,1441,1448,1449,1456,1457,1464, 1465,1472,1473,1483,1491,1495,1496 ,1498,1500,1502,1513,1517,1518,152 0,1522,1524,1525,1526,1527,1535,15 Page 40 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References 39,1540,1542,1544,1546,1557,1561,1 562,1564,1566,1568,1569,1570,1571, 1580,1584,1585,1587,1589,1591,1602 ,1606,1607,1609,1611,1613,1614,161 5,1616,1625,1629,1630,1632,1634,16 36,1647,1651,1652,1654,1656,1658,1 659,1660,1661,1670,1674,1675,1677, 1679,1681,1692,1696,1697,1699,1701 ,1703,1704,1705,1706,1715,1719,172 0,1722,1724,1726,1737,1741,1742,17 44,1746,1748,1749,1750,1751,1762,1 766,1767,1769,1771,1773,1784,1788, 1789,1791,1793,1795,1796,1797,1798 ,1807,1811,1812,1814,1816,1818,182 9,1833,1834,1836,1838,1840,1841,18 42,1843 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 IEGRD Scalar 10 I(4) 4 2 91719 COM IEGRDK Scalar 6 I(4) 4 2 91719 COM IFHR Local 814 I(4) 4 scalar 814,816,817,820,822,834,836,837,84 0,842,854,856,857,860,862,874,876, 877,880,882,895,897,898,901,903,91 6,918,919,922,924,937,939,940,943, 945,958,960,961,964,966,1038,1040, 1041,1044,1046,1068,1070,1071,1074 ,1076,1098,1100,1101,1104,1106,112 8,1130,1131,1134,1136,1158,1160,11 61,1164,1166,1183,1185,1186,1189,1 191,1209,1211,1212,1215,1217,1492, 1494,1495,1498,1500,1514,1516,1517 ,1520,1522,1536,1538,1539,1542,154 4,1558,1560,1561,1564,1566,1581,15 83,1584,1587,1589,1603,1605,1606,1 609,1611,1626,1628,1629,1632,1634, 1648,1650,1651,1654,1656,1671,1673 ,1674,1677,1679,1693,1695,1696,169 9,1701,1716,1718,1719,1722,1724,17 38,1740,1741,1744,1746,1763,1765,1 766,1769,1771,1785,1787,1788,1791, 1793,1808,1810,1811,1814,1816,1830 ,1832,1833,1836,1838 IFILV Scalar 11 I(4) 4 1 250 COM IFINCR Local 816 I(4) 4 scalar 816,819,822,836,839,842,856,859,86 2,876,879,882,897,900,903,918,921, 924,939,942,945,960,963,966,1040,1 043,1046,1070,1073,1076,1100,1103, 1106,1130,1133,1136,1160,1163,1166 ,1185,1188,1191,1211,1214,1217,149 4,1497,1500,1516,1519,1522,1538,15 Page 41 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References 41,1544,1560,1563,1566,1583,1586,1 589,1605,1608,1611,1628,1631,1634, 1650,1653,1656,1673,1676,1679,1695 ,1698,1701,1718,1721,1724,1740,174 3,1746,1765,1768,1771,1787,1790,17 93,1810,1813,1816,1832,1835,1838 IGET Scalar 8 I(4) 4 1 250 COM 144,176,177,178,179,180,247,248,25 8,261,266,269,274,278,283,287,292, 295,300,303,308,312,317,321,326,33 1,342,343,358,363,364,379,383,384, 400,407,412,417,423,428,434,439,44 3,448,452,457,469,473,476,487,488, 489,511,514,518,521,525,528,532,53 6,540,544,548,552,564,565,566,571, 577,582,589,594,607,612,621,626,63 2,641,649,651,652,653,654,660,666, 667,672,678,679,684,690,691,696,70 2,703,709,717,719,720,721,722,728, 734,735,740,746,747,752,758,759,76 4,770,771,781,789,793,806,811,826, 831,846,851,866,871,886,891,907,91 2,928,933,949,954,970,975,1000,100 5,1010,1015,1024,1049,1054,1079,10 84,1109,1114,1139,1144,1169,1174,1 195,1200,1221,1226,1229,1234,1237, 1242,1246,1251,1255,1259,1264,1268 ,1274,1277,1282,1285,1290,1293,129 8,1302,1307,1310,1315,1323,1328,13 36,1341,1350,1355,1356,1357,1358,1 359,1400,1404,1408,1412,1417,1426, 1430,1434,1438,1442,1446,1450,1454 ,1458,1462,1466,1470,1474,1481,148 4,1489,1503,1528,1533,1547,1572,15 78,1592,1617,1623,1637,1662,1668,1 682,1707,1713,1727,1752,1760,1774, 1799,1805,1819,1844 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 IIP Local 985 I(4) 4 scalar 985,989 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,10,11,12,13,14,15, 16,17,19,20,102,103,104,105,106,10 7,108,109,110,111,112,113,114,150, 183,318,345,366,386,459,491,596,60 1,613,783,797,982,1031,1061,1091,1 121,1151,1176,1202,1317,1330,1343, 1361,1419,1507,1551,1596,1641,1686 ,1731,1778,1823 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,259,262,267,270,275,276,279,28 4,285,288,293,296,301,304,309,310, 313,319,322,327,328,329,332,349,35 Page 42 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References 8,370,379,390,400,408,413,418,422, 424,429,433,435,440,442,444,449,45 1,453,466,467,470,474,477,512,515, 519,522,526,529,533,537,541,545,54 9,553,572,578,583,584,590,602,608, 614,615,616,622,627,633,650,652,65 4,665,667,677,679,689,691,701,703, 718,720,722,733,735,745,747,757,75 9,769,771,787,790,803,805,807,812, 825,827,832,845,847,852,865,867,87 2,885,887,893,906,908,914,927,929, 935,948,950,956,969,971,999,1001,1 004,1006,1009,1011,1014,1016,1036, 1050,1066,1080,1096,1110,1126,1140 ,1156,1170,1181,1194,1196,1207,122 0,1222,1227,1230,1235,1238,1244,12 47,1257,1260,1266,1269,1275,1278,1 283,1286,1291,1294,1299,1301,1303, 1308,1311,1321,1324,1334,1337,1347 ,1351,1401,1405,1409,1413,1423,142 7,1431,1435,1439,1443,1447,1451,14 55,1459,1463,1467,1471,1475,1482,1 485,1490,1504,1512,1529,1534,1548, 1556,1573,1579,1593,1601,1618,1624 ,1638,1646,1663,1669,1683,1691,170 8,1714,1728,1736,1753,1761,1775,17 83,1800,1806,1820,1828,1845 IMT Param 4 I(4) 4 scalar 7,8 IMX Param 3 I(4) 4 scalar 3,4,9,10,11,12 INPES Param 6 I(4) 4 scalar 11,14 INT Func 815 scalar 815,835,855,875,896,917,938,959,10 39,1069,1099,1129,1159,1184,1210,1 493,1515,1537,1559,1582,1604,1627, 1649,1672,1694,1717,1739,1764,1786 ,1809,1831 IOUNIT Common 4 20 IOUT Scalar 3 I(4) 4 scalar COM IOUTYP Scalar 8 I(4) 4 scalar COM 261,269,278,287,295,303,312,321,33 1,358,379,400,412,423,434,443,452, 469,476,514,521,528,536,544,552,57 7,589,607,621,632,651,653,666,678, 690,702,719,721,734,746,758,770,78 9,806,826,846,866,886,907,928,949, 970,1000,1005,1010,1015,1049,1079, 1109,1139,1169,1195,1221,1229,1237 ,1246,1259,1268,1277,1285,1293,130 2,1310,1323,1336,1350,1404,1412,14 26,1434,1442,1450,1458,1466,1474,1 484,1503,1528,1547,1572,1592,1617, 1637,1662,1682,1707,1727,1752,1774 ,1799,1819,1844 IQ Scalar 9 I(4) 4 1 250 COM IRAIN Local 987 I(4) 4 scalar 987,991 IS Scalar 9 I(4) 4 1 250 COM ISHDE Scalar 4 I(4) 4 1 61 COM Page 43 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References ISLOPE Scalar 5 I(4) 4 2 91719 COM 1344 ISLTYP Scalar 5 I(4) 4 2 91719 COM 498,1331,1382,1386 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 ISNO Local 984 I(4) 4 scalar 984,988 ISVALUE Local 410 I(4) 4 scalar 410,411,574,575,576,586,587,588,60 4,605,606,618,619,620,629,630,631, 644,645,646,662,663,664,674,675,67 6,686,687,688,698,699,700,712,713, 714,730,731,732,742,743,744,754,75 5,756,766,767,768 ITAIL Param 11 I(4) 4 scalar ITB Param 1 I(4) 4 scalar 8,18 ITBQ Param 1 I(4) 4 scalar 9,18 ITPREC Local 815 I(4) 4 scalar 815,816,820,835,836,840,855,856,86 0,875,876,880,896,897,901,917,918, 922,938,939,943,959,960,964,1184,1 185,1189,1210,1211,1215,1493,1494, 1498,1515,1516,1520,1537,1538,1542 ,1559,1560,1564,1582,1583,1587,160 4,1605,1609,1627,1628,1632,1649,16 50,1654,1672,1673,1677,1694,1695,1 699,1717,1718,1722,1739,1740,1744, 1764,1765,1769,1786,1787,1791,1809 ,1810,1814,1831,1832,1836 ITSRFC Local 1039 I(4) 4 scalar 1039,1040,1044,1069,1070,1074,1099 ,1100,1104,1129,1130,1134,1159,116 0,1164 ITSTLOC Scalar 4 I(4) 4 scalar COM IUNIT Scalar 10 I(4) 4 scalar COM IUP Scalar 7 I(4) 4 scalar COM IVEG Local 461 I(4) 4 scalar 461,462 IVG Local 1376 I(4) 4 scalar 1376,1382 IVGTYP Scalar 5 I(4) 4 2 91719 COM 461,1318,1376 IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IW Local 108 R(4) 4 2 91719 146,152,155,157,159,162,227,613 IWGT Scalar 11 I(4) 4 2 91719 COM IWM1 Local 108 R(4) 4 scalar 152,159 IWX Local 983 I(4) 4 scalar 983,984,985,986,987 IWX1 Local 102 I(4) 4 2 91719 976,983 IXM Scalar 3 I(4) 4 scalar COM IYM Scalar 3 I(4) 4 scalar COM IZR Local 986 I(4) 4 scalar 986,990 J Local 149 I(4) 4 scalar 149,151,152,153,154,155,156,157,15 9,162,182,186,187,188,189,190,193, 196,199,200,206,207,208,218,227,23 8,240,241,242,243,244,249,344,346, 365,367,385,387,458,460,461,463,49 0,495,496,497,498,499,500,502,503, 504,505,506,507,595,597,598,782,78 4,796,798,799,981,983,988,989,990, 991,1030,1032,1060,1062,1090,1092, 1120,1122,1150,1152,1175,1177,1201 Page 44 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References ,1203,1316,1318,1329,1331,1342,134 4,1360,1363,1364,1365,1366,1370,13 71,1372,1373,1376,1382,1383,1384,1 385,1386,1388,1389,1390,1391,1392, 1393,1394,1395,1396,1418,1420,1506 ,1508,1550,1552,1595,1597,1640,164 2,1685,1687,1730,1732,1777,1779,18 22,1824 JAM Param 6 I(4) 4 scalar 3 JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEGRD Scalar 10 I(4) 4 2 91719 COM JEGRDK Scalar 6 I(4) 4 2 91719 COM JEND Scalar 6 I(4) 4 scalar COM 149,182,344,365,385,458,490,595,78 2,796,981,1030,1060,1090,1120,1150 ,1175,1201,1316,1329,1342,1360,141 8,1506,1550,1595,1640,1685,1730,17 77,1822 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 JM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,10,11,12,13,14,15, 16,17,19,20,102,103,104,105,106,10 7,108,109,110,111,112,113,114,318, 601,613 JM2 Param 6 I(4) 4 scalar JMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 108,259,262,267,270,275,276,279,28 4,285,288,293,296,301,304,309,310, 313,319,322,327,328,329,332,349,35 8,370,379,390,400,408,413,418,422, 424,429,433,435,440,442,444,449,45 1,453,466,467,470,474,477,512,515, 519,522,526,529,533,537,541,545,54 9,553,572,578,583,584,590,602,608, 614,615,616,622,627,633,650,652,65 4,665,667,677,679,689,691,701,703, 718,720,722,733,735,745,747,757,75 9,769,771,787,790,803,805,807,812, 825,827,832,845,847,852,865,867,87 2,885,887,893,906,908,914,927,929, 935,948,950,956,969,971,999,1001,1 004,1006,1009,1011,1014,1016,1036, 1050,1066,1080,1096,1110,1126,1140 ,1156,1170,1181,1194,1196,1207,122 0,1222,1227,1230,1235,1238,1244,12 47,1257,1260,1266,1269,1275,1278,1 283,1286,1291,1294,1299,1301,1303, 1308,1311,1321,1324,1334,1337,1347 ,1351,1401,1405,1409,1413,1423,142 7,1431,1435,1439,1443,1447,1451,14 55,1459,1463,1467,1471,1475,1482,1 485,1490,1504,1512,1529,1534,1548, 1556,1573,1579,1593,1601,1618,1624 Page 45 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References ,1638,1646,1663,1669,1683,1691,170 8,1714,1728,1736,1753,1761,1775,17 83,1800,1806,1820,1828,1845 JMT Param 4 I(4) 4 scalar 7,8 JMX Param 3 I(4) 4 scalar 3,4,9,10,11,12 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 149,182,344,365,385,458,490,595,78 2,796,981,1030,1060,1090,1120,1150 ,1175,1201,1316,1329,1342,1360,141 8,1506,1550,1595,1640,1685,1730,17 77,1822 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 8,18 JTBQ Param 1 I(4) 4 scalar 18 JTSTLOC Scalar 4 I(4) 4 scalar COM 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 KTM Scalar 3 I(4) 4 scalar COM KUNIT Scalar 10 I(4) 4 scalar COM L Local 148 I(4) 4 scalar 148,151,153,154,156,340,343,346,35 2,355,358,364,367,373,376,379,384, 387,393,396,400 LB Param 9 I(4) 4 scalar LC Scalar 3 I(4) 4 2 91719 COM LCNTRL Scalar 5 I(4) 4 scalar COM LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LIST Scalar 3 I(4) 4 scalar COM LLGRDS Common 2 8071788 LM Param 3 I(4) 4 scalar 2,3,5,7,8,9,19,148 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 151,186,976,1371,1372,1373 LMHK Local 186 I(4) 4 scalar 186,187,188 LMV Scalar 5 I(4) 4 2 91719 COM LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4,7,8 LS Local 352 I(4) 4 scalar 352,353,373,374,393,394 LSL Scalar 3 I(4) 4 scalar COM LSM Param 3 I(4) 4 scalar 7 LUNCO Scalar 5 I(4) 4 scalar COM Page 46 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 248,261,269,278,287,295,303,312,32 1,331,343,364,384,412,423,434,443, 452,469,476,514,521,528,536,544,55 2,577,589,607,621,632,652,654,667, 679,691,703,720,722,735,747,759,77 1,789,806,826,846,866,886,907,928, 949,970,1000,1005,1010,1015,1049,1 079,1109,1139,1169,1195,1221,1229, 1237,1246,1259,1268,1277,1285,1293 ,1302,1310,1323,1336,1350,1404,141 2,1426,1434,1442,1450,1458,1466,14 74,1484,1503,1528,1547,1572,1592,1 617,1637,1662,1682,1707,1727,1752, 1774,1799,1819,1844 MAPOT Common 2 1448 MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM MOD Func 575 scalar 575,576,587,588,605,606,619,620,63 0,631,645,646,663,664,675,676,687, 688,699,700,713,714,731,732,743,74 4,755,756,767,768,816,836,856,876, 897,918,939,960,984,985,986,1040,1 070,1100,1130,1160,1185,1211,1494, 1516,1538,1560,1583,1605,1628,1650 ,1673,1695,1718,1740,1765,1787,181 0,1832 MTSTPE Scalar 4 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 356 scalar 356,357,377,378,397,398 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 NRDSW Scalar 3 I(4) 4 scalar COM NROOT Param 1 I(4) 4 scalar NROOTS Local 102 I(4) 4 2 91719 1383,1390,1420 NSHDE Scalar 4 I(4) 4 scalar COM NSOIL Param 1 I(4) 4 scalar 6,7,340,1382 NSRFC Scalar 3 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 814,834,854,874,895,916,937,958,10 38,1068,1098,1128,1158,1183,1209,1 Page 47 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References 492,1514,1536,1558,1581,1603,1626, 1648,1671,1693,1716,1738,1763,1785 ,1808,1830 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 101 L(4) 4 scalar COM ONEPS Param 39 R(4) 4 scalar 243,597 OUTPUT Subr 261 261,269,278,287,295,303,312,321,33 1,358,379,400,412,423,434,443,452, 469,476,514,521,528,536,544,552,57 7,589,607,621,632,651,653,666,678, 690,702,719,721,734,746,758,770,78 9,806,826,846,866,886,907,928,949, 970,1000,1005,1010,1015,1049,1079, 1109,1139,1169,1195,1221,1229,1237 ,1246,1259,1268,1277,1285,1293,130 2,1310,1323,1336,1350,1404,1412,14 26,1434,1442,1450,1458,1466,1474,1 484,1503,1528,1547,1572,1592,1617, 1637,1662,1682,1707,1727,1752,1774 ,1799,1819,1844 OVERRC Param 33 R(4) 4 scalar 33 P10 Local 114 R(4) 4 2 91719 187,189,689 P1000 Param 35 R(4) 4 scalar 189,190,200 P1D Local 105 R(4) 4 2 91719 P30 Local 114 R(4) 4 2 91719 188,190,757 PARALLEL Common 9 8240 PCTSNO Scalar 4 R(4) 4 2 91719 COM 463,466 PD Scalar 3 R(4) 4 2 91719 COM 187,188,196,976 PDSL Scalar 6 R(4) 4 2 91719 COM 1373 PDVP1 Scalar 6 R(4) 4 2 91719 COM PFSL Scalar 7 R(4) 4 2 91719 COM PHOUR Scalar 2 R(4) 4 2 91719 COM PHYS Common 2 7690036 PINT Scalar 4 R(4) 4 3 4219074 COM PL Scalar 5 R(4) 4 scalar COM PLQ Scalar 6 R(4) 4 scalar COM POTEVP Scalar 5 R(4) 4 2 91719 COM 497,1203 PPTASM Common 1 19811316 PPTDAT Scalar 3 R(4) 4 3 275157 COM PQ0 Param 41 R(4) 4 scalar 218 PREC Scalar 7 R(4) 4 2 91719 COM 798,976 PSFC Local 103 R(4) 4 2 91719 196,200,218,243,259,267 PSHLTR Scalar 4 R(4) 4 2 91719 COM 597,613,627 PSLP Scalar 6 R(4) 4 2 91719 COM PT Scalar 5 R(4) 4 scalar COM 187,188,196,976,1373 PT2 Scalar 3 R(4) 4 scalar COM PTBL Scalar 18 R(4) 4 2 10184 COM PTRACE Param 93 R(4) 4 scalar PTSL Scalar 7 R(4) 4 2 91719 COM PVRBLS Common 2 23480064 Page 48 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References Q Scalar 5 R(4) 4 3 4127355 COM 976,1372 Q10 Scalar 3 R(4) 4 2 91719 COM 701 Q1D Local 105 R(4) 4 2 91719 Q2 Scalar 9 R(4) 4 3 4127355 COM Q30 Scalar 5 R(4) 4 2 91719 COM 769 QFLX Common 1 8438148 QI Local 220 R(4) 4 scalar 220,223,228 QINT Local 221 R(4) 4 scalar 221,225,230 QS Scalar 3 R(4) 4 2 91719 COM 206 QS0 Scalar 8 R(4) 4 1 134 COM QSAT Local 223 R(4) 4 scalar 223,225,228,230,235,238,242 QSFC Local 103 R(4) 4 2 91719 206,207,238,242,243,309 QSHLTR Scalar 4 R(4) 4 2 91719 COM 583,597,613 QSL Scalar 9 R(4) 4 2 91719 COM QW Local 218 R(4) 4 scalar 218,220,221,235 QWBS Scalar 5 R(4) 4 2 91719 COM 497,1283 QZ0 Scalar 4 R(4) 4 2 91719 COM R Scalar 5 R(4) 4 scalar COM RADIN Scalar 11 R(4) 4 2 91719 COM RADOT Scalar 11 R(4) 4 2 91719 COM RAIN Local 109 R(4) 4 2 91719 991,1014 RC Local 1382 R(4) 4 scalar 1382 RCAPA Param 35 R(4) 4 scalar RCQ Local 113 R(4) 4 2 91719 1384,1394,1463 RCS Local 113 R(4) 4 2 91719 1384,1393,1447 RCSOIL Local 113 R(4) 4 2 91719 1384,1396,1471 RCT Local 113 R(4) 4 2 91719 1384,1395,1455 RD Param 32 R(4) 4 scalar 32,36,40 RDOCP Param 40 R(4) 4 scalar RDP Scalar 5 R(4) 4 scalar COM RDPQ Scalar 6 R(4) 4 scalar COM RDQ Scalar 5 R(4) 4 scalar COM RDTH Scalar 5 R(4) 4 scalar COM RDTHE Scalar 5 R(4) 4 scalar COM RDTHEQ Scalar 6 R(4) 4 scalar COM RDTPHS Local 794 R(4) 4 scalar 794,798 RES Scalar 3 R(4) 4 2 91719 COM 976 RESTRT Scalar 101 L(4) 4 scalar COM RF Scalar 5 R(4) 4 2 91719 COM RG Param 36 R(4) 4 scalar RHSFC Local 103 R(4) 4 2 91719 238,240,241,242,327 RITE2 Scalar 5 L(4) 4 scalar COM RITEHD Scalar 5 L(4) 4 scalar COM RMXSNAL Scalar 11 R(4) 4 2 91719 COM ROG Param 32 R(4) 4 scalar ROI Scalar 4 R(4) 4 scalar COM ROS Scalar 4 R(4) 4 scalar COM RQSTFLD Common 7 64130 RRNUM Local 1026 R(4) 4 scalar 1026,1028,1032,1056,1058,1062,1086 ,1088,1092,1116,1118,1122,1146,114 8,1152 RSMIN Local 112 R(4) 4 2 91719 1383,1389,1409 RSWIN Scalar 3 R(4) 4 2 91719 COM 1370 RSWOUT Scalar 3 R(4) 4 2 91719 COM RSWTOA Scalar 3 R(4) 4 2 91719 COM Page 49 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References RTD Param 34 R(4) 4 scalar RTDPTH Scalar 7 R(4) 4 1 4 COM RUN Scalar 101 L(4) 4 scalar COM SBD Scalar 6 R(4) 4 scalar COM SCLFLD Subr 328 328,422,433,442,451,615,805,825,84 5,865,885,906,927,948,969,1194,122 0,1301 SEC2HR Param 97 R(4) 4 scalar SFCEVP Scalar 4 R(4) 4 2 91719 COM 1177 SFCEXC Scalar 3 R(4) 4 2 91719 COM 1291 SFCLHX Scalar 3 R(4) 4 2 91719 COM 1032 SFCPRS Local 1373 R(4) 4 scalar 1373,1381 SFCQ Local 1372 R(4) 4 scalar 1372,1381 SFCSHX Scalar 3 R(4) 4 2 91719 COM 1062 SFCTMP Local 1371 R(4) 4 scalar 1371,1381 SFCUVX Scalar 4 R(4) 4 2 91719 COM 1152 SH2O Scalar 7 R(4) 4 3 366876 COM 387,498,1382 SI Scalar 6 R(4) 4 2 91719 COM 474 SICE Scalar 3 R(4) 4 2 91719 COM 496,1364 SIGMA Scalar 101 L(4) 4 scalar COM SIGT4 Scalar 16 R(4) 4 2 91719 COM SLDPTH Scalar 7 R(4) 4 1 4 COM 353,355,374,376,394,396 SLEET Local 109 R(4) 4 2 91719 989,1004 SLP Scalar 5 R(4) 4 2 91719 COM SM Scalar 3 R(4) 4 2 91719 COM 495,1363,1386 SMC Scalar 6 R(4) 4 3 366876 COM 367 SMCDRY Local 111 R(4) 4 2 91719 499,506,541 SMCMAX Local 111 R(4) 4 2 91719 500,507,549 SMCREF Local 112 R(4) 4 2 91719 1383,1392,1439 SMCWLT Local 112 R(4) 4 2 91719 1383,1385,1386,1391,1431 SMSTAV Scalar 3 R(4) 4 2 91719 COM 418 SMSTOT Scalar 3 R(4) 4 2 91719 COM 429 SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar SNCOVR Local 462 R(4) 4 scalar 462,463 SNEQV Local 460 R(4) 4 scalar 460,462 SNFRAC Subr 462 462 SNO Scalar 6 R(4) 4 2 91719 COM 449,460,497 SNOPCX Scalar 4 R(4) 4 2 91719 COM 1122 SNOW Local 109 R(4) 4 2 91719 988,999 SOIL Common 2 8438180 SOILTB Scalar 3 R(4) 4 2 91719 COM 408 SOLAR Local 1370 R(4) 4 scalar 1370,1381 SPHI0 Scalar 9 R(4) 4 scalar COM SPL Scalar 7 R(4) 4 1 39 COM SQS Scalar 8 R(4) 4 1 134 COM SR Scalar 3 R(4) 4 2 91719 COM 784 SRFDSP Common 2 4402512 SSROFF Scalar 4 R(4) 4 2 91719 COM 934 SST Scalar 14 R(4) 4 2 91719 COM STC Scalar 6 R(4) 4 3 366876 COM 346 STDOUT Scalar 2 I(4) 4 scalar COM STDRD Scalar 101 L(4) 4 scalar COM STHE Scalar 8 R(4) 4 1 76 COM STHEQ Scalar 9 R(4) 4 1 152 COM Page 50 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References SUBSHX Scalar 4 R(4) 4 2 91719 COM 1092 SURFCE2 Subr 1 SVALUE Scalar 8 R(4) 4 scalar COM SWNET Scalar 12 R(4) 4 2 91719 COM T Scalar 5 R(4) 4 3 4127355 COM 154,156,187,188,976,1371 T10 Local 114 R(4) 4 2 91719 189,677 T1D Local 105 R(4) 4 2 91719 T30 Local 114 R(4) 4 2 91719 190,745 T500 Scalar 6 R(4) 4 2 91719 COM TDTD Scalar 3 R(4) 4 scalar COM TDTQ2 Scalar 3 R(4) 4 scalar COM TFRZ Param 39 R(4) 4 scalar TG Scalar 14 R(4) 4 2 91719 COM TH10 Scalar 3 R(4) 4 2 91719 COM 189,665 TH30 Scalar 5 R(4) 4 2 91719 COM 190,733 THE0 Scalar 8 R(4) 4 1 76 COM THE0Q Scalar 9 R(4) 4 1 152 COM THL Scalar 5 R(4) 4 scalar COM THS Scalar 3 R(4) 4 2 91719 COM 199 THSFC Local 104 R(4) 4 2 91719 199,200,301 THZ0 Scalar 4 R(4) 4 2 91719 COM TLAT Scalar 2 R(4) 4 3 4127355 COM TLM0D Scalar 6 R(4) 4 scalar COM TMT0 Local 210 R(4) 4 scalar 210,211,214,220,222,224 TMT15 Local 211 R(4) 4 scalar 211,221 TPH0D Scalar 6 R(4) 4 scalar COM TPREC Scalar 3 R(4) 4 scalar COM 815,835,855,875,896,917,938,959,11 84,1210,1493,1515,1537,1559,1582,1 604,1627,1649,1672,1694,1717,1739, 1764,1786,1809,1831 TRDSW Scalar 3 R(4) 4 scalar COM TSFC Local 103 R(4) 4 2 91719 200,208,293 TSFCK Local 208 R(4) 4 scalar 208,210,219 TSHDE Scalar 8 R(4) 4 1 61 COM TSHLTR Scalar 4 R(4) 4 2 91719 COM 572,613 TSL Scalar 8 R(4) 4 2 91719 COM TSPH Scalar 5 R(4) 4 scalar COM 814,834,854,874,895,916,937,958,10 38,1068,1098,1128,1158,1183,1209,1 492,1514,1536,1558,1581,1603,1626, 1648,1671,1693,1716,1738,1763,1785 ,1808,1830 TSRFC Scalar 3 R(4) 4 scalar COM 1039,1069,1099,1129,1159 TTBL Scalar 18 R(4) 4 2 10184 COM TTBLQ Scalar 18 R(4) 4 2 66880 COM TTV Scalar 8 R(4) 4 2 91719 COM TWBS Scalar 5 R(4) 4 2 91719 COM 1275 U Scalar 5 R(4) 4 3 4127355 COM U00 Scalar 3 R(4) 4 2 91719 COM U10 Scalar 3 R(4) 4 2 91719 COM 650 U30 Scalar 5 R(4) 4 2 91719 COM 718 UA Local 106 R(4) 4 2 91719 UL Scalar 3 R(4) 4 1 90 COM USTAR Scalar 3 R(4) 4 2 91719 COM 1235 UZ0 Scalar 4 R(4) 4 2 91719 COM V Scalar 5 R(4) 4 3 4127355 COM Page 51 Source Listing SURFCE2 2014-12-17 20:48 Symbol Table SURFCE2.f Name Object Declared Type Bytes Dimen Elements Attributes References V10 Scalar 3 R(4) 4 2 91719 COM 650 V30 Scalar 5 R(4) 4 2 91719 COM 718 VA Local 106 R(4) 4 2 91719 VAPINC Scalar 5 R(4) 4 2 91719 COM 1761,1779 VAPINC7 Scalar 5 R(4) 4 2 91719 COM 1779 VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VEGFRC Scalar 5 R(4) 4 2 91719 COM 497,1299 VRBLS Common 2 67138308 VTLAT Scalar 7 R(4) 4 2 91719 COM VTLON Scalar 7 R(4) 4 2 91719 COM 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 WIJ Scalar 11 R(4) 4 2 91719 COM WIJP Scalar 12 R(4) 4 2 91719 COM WIPJ Scalar 11 R(4) 4 2 91719 COM WIPJP Scalar 12 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 1227 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 ZSFC Local 104 R(4) 4 2 91719 193,275,284 Page 52 Source Listing SURFCE2 2014-12-17 20:48 Subprograms/Common Blocks SURFCE2.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References ACMPRE Common 2 1467508 ACMRDS Common 2 2201268 ACMSFC Common 2 2568144 CLDWTR Common 2 17610408 CNVCLD Common 2 1467504 CTLBLK Common 2 96 EXTRA Common 3 71506268 IOUNIT Common 4 20 LLGRDS Common 2 8071788 LOOPS Common 2 748952 MAPOT Common 2 1448 MASKS Common 2 34853220 PARALLEL Common 9 8240 PHYS Common 2 7690036 PPTASM Common 1 19811316 PVRBLS Common 2 23480064 QFLX Common 1 8438148 RQSTFLD Common 7 64130 SOIL Common 2 8438180 SRFDSP Common 2 4402512 SURFCE2 Subr 1 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 Page 53 Source Listing SURFCE2 2014-12-17 20:48 SURFCE2.f no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -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 : SURFCE2.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100