Page 1 Source Listing PRECPD 2014-12-17 20:47 /tmpnwprd/ifortFw1aCX.i 1 # 1 "PRECPD.F" 2 SUBROUTINE PRECPD 3 C ****************************************************************** 4 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 5 C . . . 6 C SUBPROGRAM: PRECPD LARGE SCALE PRECIPITATION 7 C PRGRMMR: ZHAO ORG: W/NP22 DATE: ??-??-?? 8 C 9 C ABSTRACT: 10 C PRECPD COMPUTES THE GRID SCALE PRECIPITATION. 11 C 12 C PROGRAM HISTORY LOG: 13 C 94-??-?? ZHAO - ORIGINATOR 14 C 95-03-25 BLACK - CONVERSION FROM 1-D TO 2-D IN HORIZONTAL 15 C 95-11-20 ABELES - PARALLEL OPTIMIZATION 16 C 96-03-29 BLACK - REMOVED SCRCH COMMON 17 C 96-07-18 ZHAO - NEW WMIN CALCULATION 18 C 96-09-25 BALDWIN - NEW SR CALCULATION 19 C 98-11-02 BLACK - MODIFICATION FOR DISTRIBUTED MEMORY 20 C 21 C USAGE: CALL PRECPD FROM MAIN PROGRAM EBU 22 C 23 C INPUT ARGUMENT LIST: 24 C NONE 25 C 26 C OUTPUT ARGUMENT LIST: 27 C NONE 28 C 29 C OUTPUT FILES: 30 C NONE 31 C 32 C SUBPROGRAMS CALLED: 33 C 34 C UNIQUE: NONE 35 C 36 C LIBRARY: NONE 37 C 38 C COMMON BLOCKS: CTLBLK 39 C LOOPS 40 C MASKS 41 C PHYS 42 C VRBLS 43 C CLDWTR 44 C PVRBLS 45 C ACMCLH 46 C 47 C ATTRIBUTES: 48 C LANGUAGE: FORTRAN 90 49 C MACHINE : IBM SP 50 C$$$ 51 C ****************************************************************** 52 P A R A M E T E R 53 & (A1=610.78,A2=17.2693882,A3=273.16,A4=35.86 54 &, PQ0=379.90516,TRESH=.95,R=287.04,C0=0.15 55 &, CP=1004.6,ELWV=2.50E6,ELIV=2.834E6,ROW=1.E3,G=9.8 56 &, EPSQ=2.E-12,DLDT=2274.0,ELIW=ELIV-ELWV) 57 C Page 2 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 58 P A R A M E T E R 59 & (ARCP=A2*(A3-A4)/CP,RCP=1./CP,PQ0C=PQ0*TRESH,RROG=1./(ROW*G) 60 &, RROW=1./ROW) 61 C---------------------------------------------------------------------- 62 INCLUDE "parmeta" 85 INCLUDE "parm.tbl" 87 INCLUDE "mpp.h" 88 # 1 "./sp.h" 1 89 # 4 90 91 # 65 "PRECPD.F" 2 92 C---------------------------------------------------------------------- 151 P A R A M E T E R 152 & (IMJM=IM*JM-JM/2,JAM=6+2*(JM-10) 153 &, LP1=LM+1,LTOP=1,LBOT=LM) 154 C 155 P A R A M E T E R 156 & (IMJM_LOC=IDIM2*JDIM2,LDA=(IDIM2-IDIM1+1)*(JDIM2-JDIM1+1)) 157 C----------------------------------------------------------------------- 158 L O G I C A L 159 & RUN,FIRST,RESTRT,SIGMA,NOZ 160 C---------------------------------------------------------------------- 161 INCLUDE "CTLBLK.comm" 162 C----------------------------------------------------------------------- 172 INCLUDE "LOOPS.comm" 173 C----------------------------------------------------------------------- 180 INCLUDE "MASKS.comm" 181 C----------------------------------------------------------------------- 192 INCLUDE "PHYS.comm" 193 C----------------------------------------------------------------------- 221 INCLUDE "VRBLS.comm" 222 C----------------------------------------------------------------------- 232 INCLUDE "CLDWTR.comm" 233 C----------------------------------------------------------------------- 241 INCLUDE "PVRBLS.comm" 242 C----------------------------------------------------------------------- 272 INCLUDE "ACMCLH.comm" 273 C----------------------------------------------------------------------- 280 INCLUDE "PPTASM.comm" 281 C----------------------------------------------------------------------- 305 D I M E N S I O N 306 & IPREC(IMJM_LOC),JPREC(IMJM_LOC) 307 &,T_T(LM,IDIM1:IDIM2,JDIM1:JDIM2) 308 &,Q_T(LM,IDIM1:IDIM2,JDIM1:JDIM2) 309 &,TRAIN_T(LM,IDIM1:IDIM2,JDIM1:JDIM2) 310 &,HTM_T(LM,IDIM1:IDIM2,JDIM1:JDIM2) 311 &,CWM_T(LM,IDIM1:IDIM2,JDIM1:JDIM2) 312 &,TLAT_T(LM,IDIM1:IDIM2,JDIM1:JDIM2) 313 C----------------------------------------------------------------------- 314 D I M E N S I O N 315 & PRECRL1(IMJM_LOC),PRECSL1(IMJM_LOC),IWL1(IMJM_LOC) 316 C----------------------------------------------------------------------- 317 R E A L 318 & KE,INIT,MI0 319 C----------------------------------------------------------------------- 320 C*********************************************************************** 321 C--------------PREPARATORY CALCULATIONS--------------------------------- Page 3 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 322 DTPH = NPHS * DT 323 RDTPH = 1. / DTPH 324 TWODT= DTPH 325 RTWODT=1./TWODT 326 KE=2.0E-5 327 US=1. 328 EPS=0.622E0 329 CCLIMIT=1.0E-3 330 CLIMIT=1.0E-20 331 CWS=0.025 332 CSM1=5.0000E-8 333 CRS1=5.00000E-6 334 CRS2=6.66600E-10 335 CR=5.0E-4 336 MI0=5.0E-4 337 AA2=1.25E-3 338 AVRAIN=AVRAIN+1. 339 ARATIM=ARATIM+1 340 C-------------------PADDING CLOUD MIXING RATIO IF TOO SMALL------------- 341 !$omp parallel do 342 DO 20 L=1,LM 343 DO J=MYJS,MYJE 344 DO I=MYIS,MYIE 345 CWM(I,J,L)=CWM(I,J,L)*HTM(I,J,L)*HBM2(I,J) 346 IF(CWM(I,J,L).LT.0.)CWM(I,J,L)=0. 347 C------------------PADDING SPECIFIC HUMIDITY IF TOO SMALL--------------- 348 IF(Q(I,J,L).LT.EPSQ)Q(I,J,L)=EPSQ*HTM(I,J,L) 349 ENDDO 350 ENDDO 351 20 CONTINUE 352 C 353 UTIM=1. 354 C----------------------------------------------------------------------- 355 !$omp parallel do 356 DO N=1,IMJM_LOC 357 IWL1(N)=0 358 PRECRL1(N)=0. 359 PRECSL1(N)=0. 360 ENDDO 361 C------------CHOOSE THE COLUMNS WHERE PREC CAN BE PRODUCED-------------- 362 NPRE=0 363 DO 35 J=MYJS2,MYJE2 364 DO 35 I=MYIS,MYIE 365 C 366 DO L=2,LM 367 TTEMP=0.025*(T(I,J,L)-273.16) 368 WFIX=0.9814*EXP(0.01873*L) 369 WMIN=0.1E-3*EXP(TTEMP)*WFIX 370 IF(CWM(I,J,L).GT.WMIN)GO TO 33 371 ENDDO 372 C 373 GO TO 35 374 33 NPRE=NPRE+1 375 IPREC(NPRE)=I 376 JPREC(NPRE)=J 377 35 CONTINUE 378 C------------------------------------------------------------------------ Page 4 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 379 C*** 380 C*** TRANSPOSE ARRAYS 381 C*** 382 !$omp parallel sections 383 !$omp section 384 CALL SGETMO(T,LDA,LDA,LM,T_T,LM) 385 !$omp section 386 CALL SGETMO(Q,LDA,LDA,LM,Q_T,LM) 387 !$omp section 388 CALL SGETMO(CWM,LDA,LDA,LM,CWM_T,LM) 389 !$omp section 390 CALL SGETMO(HTM,LDA,LDA,LM,HTM_T,LM) 391 !$omp section 392 CALL SGETMO(TRAIN,LDA,LDA,LM,TRAIN_T,LM) 393 !$omp section 394 CALL SGETMO(TLAT,LDA,LDA,LM,TLAT_T,LM) 395 !$omp end parallel sections 396 C----------------------------------------------------------------------- 397 C-----------------BEGINING OF PRECIPITATION CALCULATION----------------- 398 C----------------------------------------------------------------------- 399 C*** 400 C*** LOOP OVER ALL POSSIBLE PRECIPITATION POINTS 401 C*** 402 !$omp parallel do 403 !$omp& private(aai,aetal,ai,amaxcm,amaxps,amaxrq,bi,ccr,conde, 404 !$omp& const,cpdr,cs,cwmk,detal,erk,err,errt,ers,erst,expf,fi, 405 !$omp& fiw,hbm2k,hh,htmk,i,iwl,j,lml,mi0,pdsl,pid,ppr,pps, 406 !$omp& pracw,praut,precrk,precrl,precsk,precsl,precss,psaci, 407 !$omp& psaut,psm,psm1,psm2,qc,qct,qi,qint,qintt,qit,qk,qq, 408 !$omp& qtemp,qw,qwt,rconde,rprs,rq,rqkll,rqt,rqtt,tk,tmt0, 409 !$omp& tmt0k,tmt0t,tmt15,tmt15t,totppt,tt,ttemp,ttlat,u00ij, 410 !$omp& u00kl,u00klt,ull,wfix,wmink,ww) 411 C 412 DO 300 N=1,NPRE 413 C 414 I=IPREC(N) 415 J=JPREC(N) 416 HBM2K=HBM2(I,J) 417 PDSL=RES(I,J)*PD(I,J) 418 CONST=PDSL/G*TWODT 419 LML=LM-LMH(I,J) 420 U00IJ=U00(I,J) 421 C 422 DO 180 L=2,LM 423 C 424 DETAL=DETA(L) 425 ULL=UL(L) 426 AETAL=AETA(L) 427 WFIX=0.9814*EXP(0.01873*L) 428 C 429 PRECRL=0. 430 PRECSL=0. 431 PRAUT=0. 432 PSAUT=0. 433 PRACW=0. 434 PSACI=0. 435 ERR =0. Page 5 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 436 ERS =0. 437 PSM =0. 438 PSM1 =0. 439 PSM2 =0. 440 PPR =0. 441 PPS =0. 442 CPDR =0. 443 HH =0. 444 PID =0. 445 IWL =0. 446 CONDE=0. 447 RCONDE=0. 448 C----------------------------------------------------------------------- 449 TT=T_T(L,I,J) 450 QQ=Q_T(L,I,J) 451 WW=CWM_T(L,I,J) 452 HTMK=HTM_T(L,I,J) 453 TTLAT=TLAT_T(L,I,J) 454 C----------------------------------------------------------------------- 455 U00KL=U00IJ+UL(L+LML)*(0.95-U00IJ)*UTIM 456 TTEMP=0.025*(TT-273.16) 457 WMINK=0.1E-3*EXP(TTEMP)*WFIX 458 C----------------------------------------------------------------------- 459 C----------CHOOSE THE POINTS WHERE PRECIPITATION CAN BE PRODUCED-------- 460 C----------------------------------------------------------------------- 461 PRECRK=AMAX1(0.,PRECRL1(N)) 462 PRECSK=AMAX1(0.,PRECSL1(N)) 463 HH=HTMK*HBM2K 464 IF(WW.LT.WMINK.AND.(PRECRK+PRECSK).EQ.0.)THEN 465 PID=0. 466 ELSE 467 PID=HH 468 ENDIF 469 C----------------------------------------------------------------------- 470 C-------------------QW, QI AND QINT------------------------------------- 471 C----------------------------------------------------------------------- 472 IF(PID.EQ.1.)THEN 473 CONDE=CONST*DETAL 474 RCONDE=1./CONDE 475 TK=TT 476 QK=QQ 477 TMT0=(TK-273.16)*HH 478 TMT15=AMIN1(TMT0,-15.)*HH 479 AI=0.008855 480 BI=1. 481 C 482 IF(TMT0.LT.-20.)THEN 483 AI=0.007225 484 BI=0.9674 485 ENDIF 486 C 487 QW=HH*PQ0/(PDSL*AETAL+PT) 488 1 *EXP(HH*A2*(TK-A3)/(TK-A4)) 489 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 490 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 491 IF(TMT0.LE.-40.)QINT=QI 492 C-------------------ICE-WATER ID NUMBER IW------------------------------ Page 6 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 493 IF(TMT0.LT.-15.)THEN 494 FI=QK-U00KL*QI 495 IF(FI.GT.0..OR.WW.GT.CLIMIT) THEN 496 IWL=1 497 ELSE 498 IWL=0 499 ENDIF 500 ENDIF 501 C 502 IF(TMT0.LT.0.0.AND.TMT0.GE.-15.0)THEN 503 IWL=0 504 IF(IWL1(N).EQ.1.AND.WW.GT.CLIMIT)IWL=1 505 ENDIF 506 C 507 IF(TMT0.GE.0.)THEN 508 IWL=0 509 ENDIF 510 C----------------THE SATUATION SPECIFIC HUMIDITY------------------------ 511 FIW=FLOAT(IWL) 512 QC=(1.-FIW)*QINT+FIW*QI 513 C----------------THE RELATIVE HUMIDITY---------------------------------- 514 IF(QC.LE.0.)THEN 515 RQ=1.E-10 516 ELSE 517 RQ=QK/QC 518 ENDIF 519 C----------------CLOUD COVER RATIO CCR---------------------------------- 520 IF(RQ.LE.U00KL)THEN 521 CCR=0. 522 ELSE 523 RQKLL=AMIN1(US,RQ) 524 CCR=1.-SQRT((US-RQKLL)/(US-U00KL)) 525 ENDIF 526 C-----------CORRECT CCR IF IT IS TOO SMALL IN LARGE CWM REGIONS-------- 527 IF(CCR.GE.0.01.AND.CCR.LE.0.2.AND 528 1 .WW.GE.0.2E-3)THEN 529 CCR=AMIN1(1.,WW*1.0E3) 530 ENDIF 531 ENDIF 532 60 CONTINUE 533 C----------------------------------------------------------------------- 534 C------------------PRECIPITATION PRODUCTION RATES----------------------- 535 C------------------AUTO-CONVERT RATES----------------------------------- 536 C----------------------------------------------------------------------- 537 IF(PID.EQ.1.)THEN 538 IWK=IWL 539 CWMK=AMAX1(0.,WW-CLIMIT) 540 MI0=WMINK 541 C 542 IF(IWK.EQ.1)THEN 543 EXPF=EXP(0.025*TMT0) 544 AA1=1.E-3*EXPF 545 PSAUT=AA1*AMAX1(0.,CWMK-MI0) 546 CPDR=-PSAUT*TWODT 547 IF(-CPDR.GE.CWMK)THEN 548 CPDR=-CWMK 549 PSAUT=-CPDR*RTWODT Page 7 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 550 ENDIF 551 ELSE 552 AMAXCM=AMAX1(0.,CWMK-MI0) 553 PRAUT=C0*AMAXCM*AMAXCM 554 CPDR=-PRAUT*TWODT 555 IF(-CPDR.GE.CWMK)THEN 556 CPDR=-CWMK 557 PRAUT=-CPDR*RTWODT 558 ENDIF 559 ENDIF 560 PPR=PRAUT*CONDE 561 PPS=PSAUT*CONDE 562 ENDIF 563 C 564 IF(PID.EQ.1.)THEN 565 WW=CPDR*HH+WW 566 PRECRL=PRECRL1(N)+PPR*HH 567 PRECSL=PRECSL1(N)+PPS*HH 568 ENDIF 569 C----------------------------------------------------------------------- 570 C-----------------------ACCRETIONS-------------------------------------- 571 C----------------------------------------------------------------------- 572 IF(PID.EQ.1.)THEN 573 IWK=IWL 574 CWMK=WW 575 PRECRK=AMAX1(0.,PRECRL1(N)) 576 PRECSK=AMAX1(0.,PRECSL1(N)) 577 IF(IWK.EQ.1)THEN 578 EXPF=EXP(0.025*TMT0) 579 CS=AA2*EXPF 580 PSACI=CS*AMAX1(0.,CWMK)*PRECSK 581 CPDR=-PSACI*TWODT 582 IF(-CPDR.GE.CWMK)THEN 583 CPDR=-CWMK 584 PSACI=-CPDR*RTWODT 585 ENDIF 586 ELSE 587 PRACW=CR*AMAX1(0.,CWMK)*(PRECRK+PRECSK) 588 CPDR=-PRACW*TWODT 589 IF(-CPDR.GE.CWMK)THEN 590 CPDR=-CWMK 591 PRACW=-CPDR*RTWODT 592 ENDIF 593 ENDIF 594 PPR=PRACW*CONDE 595 PPS=PSACI*CONDE 596 ENDIF 597 C 598 IF(PID.EQ.1.)THEN 599 WW=CPDR*HH+WW 600 PRECRL=PRECRL+PPR*HH 601 PRECSL=PRECSL+PPS*HH 602 ENDIF 603 C----------------------------------------------------------------------- 604 C-----EVAPORATION/CONDENSATION OF PRECIPITATION------------------------- 605 C***** ERR & ERS POSITIVE--EVAPORATION 606 C***** ERR & ERS NEGTIVE---CONDENSATION Page 8 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 607 C----------------------------------------------------------------------- 608 IF(PID.EQ.1.0)THEN 609 QK=QQ 610 TMT0K=TMT0 611 IF(TMT0K.LT.-30.)TMT0K=-30. 612 PRECRK=AMAX1(0.,PRECRL) 613 PRECSK=AMAX1(0.,PRECSL) 614 C---------------------------------------------------------------------- 615 C INCREASE THE EVAPORATION/CONDENSATION FOR STRONG/LIGHT PREC 616 C---------------------------------------------------------------------- 617 U00KLT=U00KL 618 AMAXRQ=AMAX1(0.,U00KL-RQ) 619 ERR=KE*AMAXRQ*PRECRK**0.5 620 C 621 IF(TMT0.GE.0.)THEN 622 ERS=0. 623 ELSE 624 ERS=(CRS1+CRS2*TMT0K)*AMAXRQ*PRECSK/U00KLT 625 ENDIF 626 C 627 IF(ERR+ERS.LE.1.E-20) GO TO 125 628 C---------------CORRECT IF OVER-EVAPO./COND. OCCURS-------------------- 629 HHT=HH*TWODT 630 TTEMP=TT-RCP*(ELWV*ERR+ELIV*ERS)*HHT 631 QTEMP=QQ+HHT*(ERR+ERS) 632 TMT0T=(TTEMP-273.16)*HH 633 IF(TMT0T.LT.-30.)TMT0T=-30. 634 TMT15T=AMIN1(TMT0T,-15.)*HH 635 AI=0.008855 636 BI=1. 637 C 638 IF(TMT0T.LT.-20.)THEN 639 AI=0.007225 640 BI=0.9674 641 ENDIF 642 C 643 QWT=HH*PQ0/(PDSL*AETAL+PT) 644 1 *EXP(HH*A2*(TTEMP-A3)/(TTEMP-A4)) 645 QIT=QWT*(BI+AI*AMIN1(TMT0T,0.)) 646 QINTT=QWT*(1.-0.00032*TMT15T*(TMT15T+15.)) 647 IF(TMT0T.LE.-40.)QINTT=QIT 648 FIW=FLOAT(IWL) 649 QCT=(1.-FIW)*QINTT+FIW*QIT 650 C 651 IF(QCT.LE.1.E-10) THEN 652 RQT=1.E-10 653 RQTT=1.E-10 654 ELSE 655 RQT=QTEMP/QCT 656 RQTT=QQ/QCT 657 ENDIF 658 C 659 IF(RQT.LE.U00KL) GO TO 125 660 C 661 ERK=(U00KL-RQTT)*QCT*RTWODT 662 RPRS=ERK/(PRECRK+PRECSK) 663 ERRT=PRECRK*RPRS Page 9 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 664 ERST=PRECSK*RPRS 665 ERR=AMAX1(0.,0.5*(ERR+ERRT)) 666 ERS=AMAX1(0.,0.5*(ERS+ERST)) 667 C 668 125 CONTINUE 669 C 670 PPR=-ERR*CONDE 671 PPS=-ERS*CONDE 672 C 673 IF(-PPR.GE.PRECRK)THEN 674 PPR=-PRECRK 675 ERR=-PPR*RCONDE 676 ENDIF 677 C 678 IF(-PPS.GE.PRECSK)THEN 679 PPS=-PRECSK 680 ERS=-PPS*RCONDE 681 ENDIF 682 C 683 ENDIF 684 C 685 IF(PID.EQ.1.)THEN 686 PRECRL=PRECRL+PPR*HH 687 PRECSL=PRECSL+PPS*HH 688 ENDIF 689 C----------------------------------------------------------------------- 690 C--------------------MELTING OF THE SNOW-------------------------------- 691 C----------------------------------------------------------------------- 692 IF(PID.EQ.1.)THEN 693 CWMK=WW 694 AMAXPS=AMAX1(0.,PRECSL) 695 C 696 IF(TMT0.GT.0.)THEN 697 PSM1=CSM1*TMT0*TMT0*AMAXPS 698 PSM2=CWS*CR*CWMK*AMAXPS 699 PSM=PSM1+PSM2 700 ELSE 701 PSM1=0. 702 PSM2=0. 703 PSM=0. 704 ENDIF 705 C 706 PPR=PSM*CONDE 707 PPS=-PSM*CONDE 708 C 709 IF(-PPS.GE.AMAXPS)THEN 710 PPS=-AMAXPS 711 PPR=AMAXPS 712 PSM1=-PPS*RCONDE 713 PSM2=0. 714 PSM=PSM1 715 ENDIF 716 C 717 ENDIF 718 C 719 IF(PID.EQ.1.)THEN 720 PRECRL=PRECRL+PPR*HH Page 10 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 721 PRECSL=PRECSL+PPS*HH 722 ENDIF 723 C----------------------------------------------------------------------- 724 C---------------UPDATE T AND Q------------------------------------------ 725 C----------------------------------------------------------------------- 726 IF(PID.EQ.1.)THEN 727 HHT=HH*TWODT 728 TT=-RCP*(ELWV*ERR+ELIV*ERS+ELIW*PSM1) 729 1 *HHT+TT 730 QQ=(ERR+ERS)*HHT+QQ 731 TTLAT=-RCP*(ELWV*ERR+ELIV*ERS+ELIW*PSM1 732 & -ELWV*(PRAUT+PRACW) 733 & -ELIV*(PSAUT+PSACI))*HHT+TTLAT 734 ENDIF 735 C 736 IF(HH.EQ.1.)THEN 737 IWL1(N)=IWL 738 PRECRL1(N)=PRECRL 739 PRECSL1(N)=PRECSL 740 ENDIF 741 C 742 C ACCUMULATE LATENT HEATING DUE TO GRID-SCALE PRECIP/EVAP. 743 C SCALE BY THE RECIPROCAL OF THE PERIOD AT WHICH THIS ROUTINE 744 C IS CALLED. THIS PERIOD IS THE PHYSICS TIMESTEP. 745 C 746 TRAIN_T(L,I,J)=TRAIN_T(L,I,J)+(TT-T_T(L,I,J))*RDTPH 747 T_T(L,I,J)=TT 748 Q_T(L,I,J)=QQ 749 CWM_T(L,I,J)=WW 750 TLAT_T(L,I,J)=TTLAT 751 180 CONTINUE 752 C----------------------------------------------------------------------- 753 C-------------------THE PRECIPITATION ON SFC---------------------------- 754 C----------------------------------------------------------------------- 755 PRECRS=PRECRL1(N)*RROW 756 PRECSS=PRECSL1(N)*RROW 757 C 758 APREC(I,J)=PRECRS+PRECSS 759 PREC(I,J)=PREC(I,J)+PRECRS+PRECSS 760 ACPREC(I,J)=ACPREC(I,J)+APREC(I,J) 761 C----------------------------------------------------------------------- 762 C---------------THE SNOW AND RAIN RATIO OF SFC PREC--------------------- 763 C----SR IS THE RATIO OF SNOW TO THE TOTAL PRECIP------------------------ 764 C----IF TOTAL PRECIP IS ZERO, SR IS ZERO-------------------------------- 765 C----------------------------------------------------------------------- 766 TOTPPT=PRECRS+PRECSS 767 IF (TOTPPT.GT.1.E-8) THEN 768 SR(I,J)=PRECSS/TOTPPT 769 ELSE 770 SR(I,J)=0. 771 ENDIF 772 C----------------------------------------------------------------------- 773 300 CONTINUE 774 C----------------------------------------------------------------------- 775 C*** 776 C*** TRANSPOSE BACK 777 C*** Page 11 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F 778 !$omp parallel sections 779 !$omp section 780 CALL SGETMO(T_T,LM,LM,LDA,T,LDA) 781 !$omp section 782 CALL SGETMO(Q_T,LM,LM,LDA,Q,LDA) 783 !$omp section 784 CALL SGETMO(CWM_T,LM,LM,LDA,CWM,LDA) 785 !$omp section 786 CALL SGETMO(TRAIN_T,LM,LM,LDA,TRAIN,LDA) 787 !$omp section 788 CALL SGETMO(TLAT_T,LM,LM,LDA,TLAT,LDA) 789 !$omp end parallel sections 790 C----------------------------------------------------------------------- 791 RETURN 792 END ENTRY POINTS Name precpd_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 125 Label 457 416,448 180 Label 540 211 20 Label 140 131 300 Label 562 201 33 Label 163 159 35 Label 166 152,153,162 60 Label 321 A1 Param 52 R(4) 4 scalar A2 Param 52 R(4) 4 scalar 58,277,433 A3 Param 52 R(4) 4 scalar 58,277,433 A4 Param 52 R(4) 4 scalar 58,277,433 AA1 Local 333 R(4) 4 scalar 333,334 AA2 Local 126 R(4) 4 scalar 126,368 ACCLIQ Scalar 13 R(4) 4 2 17955 COM ACMCLH Common 2 6463824 ACPREC Scalar 12 R(4) 4 2 17955 COM 549 ACUTIM Scalar 3 R(4) 4 scalar COM AETA Scalar 9 R(4) 4 1 45 COM 215 AETAL Local 215 R(4) 4 scalar 215,276,432 AFSI Scalar 26 R(4) 4 2 17955 COM AI Local 268 R(4) 4 scalar 268,272,278,424,428,434 AKHS Scalar 7 R(4) 4 2 17955 COM AKMS Scalar 7 R(4) 4 2 17955 COM ALBASE Scalar 20 R(4) 4 2 17955 COM ALBEDO Scalar 21 R(4) 4 2 17955 COM AMAX1 Func 250 scalar 250,251,328,334,341,364,365,369,37 6,401,402,407,454,455,483 AMAXCM Local 341 R(4) 4 scalar 341,342 AMAXPS Local 483 R(4) 4 scalar 483,486,487,498,499,500 Page 12 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References AMAXRQ Local 407 R(4) 4 scalar 407,408,413 AMIN1 Func 267 scalar 267,278,312,318,423,434 APREC Scalar 3 R(4) 4 2 17955 COM 547,549 ARATIM Scalar 3 R(4) 4 scalar COM 128 ARCP Param 58 R(4) 4 scalar AVCNVC Scalar 3 R(4) 4 scalar COM AVRAIN Scalar 3 R(4) 4 scalar COM 127 BI Local 269 R(4) 4 scalar 269,273,278,425,429,434 C0 Param 53 R(4) 4 scalar 342 CCLIMIT Local 118 R(4) 4 scalar 118 CCR Local 310 R(4) 4 scalar 310,313,316,318 CI Scalar 6 R(4) 4 scalar COM CLDEFI Scalar 11 R(4) 4 2 17955 COM CLDINC Scalar 6 R(4) 4 2 17955 COM CLDINC7 Scalar 6 R(4) 4 2 17955 COM CLDWTR Common 2 3447720 CLIMIT Local 119 R(4) 4 scalar 119,284,293,328 CNVBOT Scalar 18 R(4) 4 2 17955 COM CNVTOP Scalar 18 R(4) 4 2 17955 COM CONDE Local 235 R(4) 4 scalar 235,262,263,349,350,383,384,459,46 0,495,496 CONST Local 207 R(4) 4 scalar 207,262 CP Param 54 R(4) 4 scalar 58 CPDR Local 231 R(4) 4 scalar 231,335,336,337,338,343,344,345,34 6,354,370,371,372,373,377,378,379, 380,388 CR Local 124 R(4) 4 scalar 124,376,487 CRS1 Local 122 R(4) 4 scalar 122,413 CRS2 Local 123 R(4) 4 scalar 123,413 CS Scalar 6 R(4) 4 scalar COM 368,369 CSM1 Local 121 R(4) 4 scalar 121,486 CTLBLK Common 4 112 CUPREC Scalar 13 R(4) 4 2 17955 COM CWM Scalar 3 R(4) 4 3 807975 COM 134,135,159,177,573 CWMK Local 328 R(4) 4 scalar 328,334,336,337,341,344,345,363,36 9,371,372,376,378,379,482,487 CWM_T Local 100 R(4) 4 3 807975 177,240,538,573 CWS Local 120 R(4) 4 scalar 120,487 CZEN Scalar 16 R(4) 4 2 17955 COM CZMEAN Scalar 23 R(4) 4 2 17955 COM DETA Scalar 9 R(4) 4 1 45 COM 213 DETAL Local 213 R(4) 4 scalar 213,262 DFRLG Scalar 9 R(4) 4 1 46 COM DI Scalar 6 R(4) 4 scalar COM DLDT Param 55 R(4) 4 scalar DS Scalar 6 R(4) 4 scalar COM DT Scalar 6 R(4) 4 scalar COM 111 DTD Scalar 5 R(4) 4 scalar COM DTPH Local 111 R(4) 4 scalar 111,112,113 DTQ2 Scalar 5 R(4) 4 scalar COM ELIV Param 54 R(4) 4 scalar 55,419,517,520,522 ELIW Param 55 R(4) 4 scalar 517,520 ELWV Param 54 R(4) 4 scalar 55,419,517,520,521 EPS Local 117 R(4) 4 scalar 117 EPSQ Param 55 R(4) 4 scalar 137 Page 13 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References EPSR Scalar 13 R(4) 4 2 17955 COM ERK Local 450 R(4) 4 scalar 450,451 ERR Local 224 R(4) 4 scalar 224,408,416,419,420,454,459,464,51 7,519,520 ERRT Local 452 R(4) 4 scalar 452,454 ERS Local 225 R(4) 4 scalar 225,411,413,416,419,420,455,460,46 9,517,519,520 ERST Local 453 R(4) 4 scalar 453,455 EXP Func 157 scalar 157,158,216,246,277,332,367,433 EXPF Local 332 R(4) 4 scalar 332,333,367,368 FI Local 283 R(4) 4 scalar 283,284 FIRST Scalar 74 L(4) 4 scalar COM FIS Scalar 7 R(4) 4 2 17955 COM FIW Local 300 R(4) 4 scalar 300,301,437,438 FLOAT Func 300 scalar 300,437 G Param 54 R(4) 4 scalar 58,207 G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GFFC Scalar 19 R(4) 4 2 17955 COM GLAT Scalar 15 R(4) 4 2 17955 COM GLB_TABLE Common 37 128 GLON Scalar 15 R(4) 4 2 17955 COM HBM2 Scalar 8 R(4) 4 2 17955 COM 134,205 HBM2K Local 205 R(4) 4 scalar 205,252 HBM3 Scalar 8 R(4) 4 2 17955 COM HBOT Scalar 17 R(4) 4 2 17955 COM HDAC Scalar 22 R(4) 4 2 17955 COM HDACV Scalar 22 R(4) 4 2 17955 COM HH Local 232 R(4) 4 scalar 232,252,256,266,267,276,277,354,35 5,356,388,389,390,418,421,423,432, 433,475,476,509,510,516,525 HHT Local 418 R(4) 4 scalar 418,419,420,516,518,519,522 HTM Scalar 6 R(4) 4 3 807975 COM 134,137,179 HTMG Scalar 52 R(4) 4 3 4183695 COM HTMK Local 241 R(4) 4 scalar 241,252 HTM_T Local 99 R(4) 4 3 807975 179,241 HTOP Scalar 17 R(4) 4 2 17955 COM I Local 133 I(4) 4 scalar 133,134,135,137,153,156,159,164,20 3,205,206,208,209,238,239,240,241, 242,535,536,537,538,539,547,548,54 9,557,559 IBROW Scalar 22 I(4) 4 scalar COM ICHUNKTAB Scalar 29 I(4) 4 1 8 COM IDAT Scalar 5 I(4) 4 1 3 COM IDIM1 Param 19 I(4) 4 scalar 2,3,4,5,6,7,8,9,10,11,12,13,14,15, 16,17,18,19,20,21,22,23,26,71,96,9 7,98,99,100,101 IDIM2 Param 19 I(4) 4 scalar 2,3,4,5,6,7,8,9,10,11,12,13,14,15, 16,17,18,19,20,21,22,23,26,71,96,9 7,98,99,100,101 IDTAD Scalar 6 I(4) 4 scalar COM IE_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 Page 14 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References IHHA Scalar 3 I(4) 4 1 760 COM IHLA Scalar 3 I(4) 4 1 760 COM IHRST Scalar 5 I(4) 4 scalar COM ILCOL Scalar 22 I(4) 4 scalar COM ILPAD1 Scalar 23 I(4) 4 scalar COM ILPAD2 Scalar 23 I(4) 4 scalar COM ILPAD3 Scalar 23 I(4) 4 scalar COM ILPAD4 Scalar 23 I(4) 4 scalar COM ILPAD5 Scalar 23 I(4) 4 scalar COM IM Param 3 I(4) 4 scalar 16,19,43,44,45,46,47,51,52,57,67 IMJM Param 67 I(4) 4 scalar IMJM_LOC Param 71 I(4) 4 scalar 95,104,145 INIT Local 107 R(4) 4 scalar INPES Param 11 I(4) 4 scalar 16,19,27,28,29,38,39 INUMQ Scalar 34 I(4) 4 1 100 COM IOUT Scalar 5 I(4) 4 scalar COM IPREC Local 95 I(4) 4 1 15327 164,203 IQUILT_GROUP Scalar 33 I(4) 4 scalar COM IRCOL Scalar 22 I(4) 4 scalar COM IRPAD1 Scalar 24 I(4) 4 scalar COM IRPAD2 Scalar 24 I(4) 4 scalar COM IRPAD3 Scalar 24 I(4) 4 scalar COM IRPAD4 Scalar 24 I(4) 4 scalar COM IRPAD5 Scalar 24 I(4) 4 scalar COM IS_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM ITAIL Param 16 I(4) 4 scalar ITB Param 1 I(4) 4 scalar 10,25 ITBQ Param 1 I(4) 4 scalar 11,25 ITEMP Scalar 47 I(4) 4 2 91719 COM ITEMP2 Scalar 47 I(4) 4 2 91719 COM ITROW Scalar 22 I(4) 4 scalar COM ITSTLOC Scalar 4 I(4) 4 scalar COM IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IWK Local 327 I(4) 4 scalar 327,331,362,366 IWL Local 234 I(4) 4 scalar 234,285,287,292,293,297,300,327,36 2,437,526 IWL1 Local 104 I(4) 4 1 15327 146,293,526 J Local 132 I(4) 4 scalar 132,134,135,137,152,156,159,165,20 4,205,206,208,209,238,239,240,241, 242,535,536,537,538,539,547,548,54 9,557,559 JAM Param 67 I(4) 4 scalar 3 JBPAD1 Scalar 25 I(4) 4 scalar COM JBPAD2 Scalar 25 I(4) 4 scalar COM JBPAD3 Scalar 25 I(4) 4 scalar COM JBPAD4 Scalar 25 I(4) 4 scalar COM JBPAD5 Scalar 25 I(4) 4 scalar COM JDIM1 Param 20 I(4) 4 scalar 2,3,4,5,6,7,8,9,10,11,12,13,14,15, 16,17,18,19,20,21,22,23,26,71,96,9 7,98,99,100,101 JDIM2 Param 20 I(4) 4 scalar 2,3,4,5,6,7,8,9,10,11,12,13,14,15, 16,17,18,19,20,21,22,23,26,71,96,9 7,98,99,100,101 Page 15 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References JE_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 JM Param 3 I(4) 4 scalar 17,20,43,44,45,46,47,51,52,57,67 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 JPREC Local 95 I(4) 4 1 15327 165,204 JRA Scalar 3 I(4) 4 1 760 COM JS_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM JTAIL Param 17 I(4) 4 scalar JTB Param 1 I(4) 4 scalar 10,25 JTBQ Param 1 I(4) 4 scalar 25 JTPAD1 Scalar 26 I(4) 4 scalar COM JTPAD2 Scalar 26 I(4) 4 scalar COM JTPAD3 Scalar 26 I(4) 4 scalar COM JTPAD4 Scalar 26 I(4) 4 scalar COM JTPAD5 Scalar 26 I(4) 4 scalar COM JTSTLOC Scalar 4 I(4) 4 scalar COM KE Local 107 R(4) 4 scalar 115,408 KTM Scalar 5 I(4) 4 scalar COM L Local 131 I(4) 4 scalar 131,134,135,137,155,156,157,159,21 1,213,214,215,216,238,239,240,241, 242,244,535,536,537,538,539 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LBOT Param 68 I(4) 4 scalar LC Scalar 5 I(4) 4 2 17955 COM LDA Param 71 I(4) 4 scalar 173,175,177,179,181,183,569,571,57 3,575,577 LIST Scalar 5 I(4) 4 scalar COM LM Param 3 I(4) 4 scalar 2,3,4,5,6,7,9,21,52,68,96,97,98,99 ,100,101,131,155,173,175,177,179,1 81,183,208,211,569,571,573,575,577 LMH Scalar 4 I(4) 4 2 17955 COM 208 LML Local 208 I(4) 4 scalar 208,244 LMV Scalar 4 I(4) 4 2 17955 COM LOOPS Common 2 158840 LP1 Param 68 I(4) 4 scalar 9 LSM Param 3 I(4) 4 scalar LTOP Param 68 I(4) 4 scalar MAPPINGS Common 56 5024 MASKS Common 2 6894720 MI0 Local 107 R(4) 4 scalar 125,329,334,341 MPI_COMM_COMP Scalar 30 I(4) 4 scalar COM MPI_COMM_INTER Scalar 31 I(4) 4 scalar COM MPI_COMM_INTER_ARRAY Scalar 32 I(4) 4 1 100 COM MPPCOM Common 1 1464 MTSTPE Scalar 4 I(4) 4 scalar COM MXSNAL Scalar 2 R(4) 4 2 17955 COM MYIE Scalar 4 I(4) 4 scalar COM 133,153 MYIE1 Scalar 4 I(4) 4 scalar COM MYIE1_P1 Scalar 8 I(4) 4 scalar COM MYIE1_P2 Scalar 8 I(4) 4 scalar COM MYIE1_P3 Scalar 8 I(4) 4 scalar COM Page 16 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References MYIE1_P4 Scalar 8 I(4) 4 scalar COM MYIE2 Scalar 4 I(4) 4 scalar COM MYIE2_P1 Scalar 9 I(4) 4 scalar COM MYIE_P1 Scalar 7 I(4) 4 scalar COM MYIE_P2 Scalar 7 I(4) 4 scalar COM MYIE_P3 Scalar 7 I(4) 4 scalar COM MYIE_P4 Scalar 7 I(4) 4 scalar COM MYIE_P5 Scalar 7 I(4) 4 scalar COM MYIS Scalar 4 I(4) 4 scalar COM 133,153 MYIS1 Scalar 4 I(4) 4 scalar COM MYIS1_P1 Scalar 6 I(4) 4 scalar COM MYIS1_P2 Scalar 6 I(4) 4 scalar COM MYIS1_P3 Scalar 6 I(4) 4 scalar COM MYIS1_P4 Scalar 6 I(4) 4 scalar COM MYIS2 Scalar 4 I(4) 4 scalar COM MYIS_P1 Scalar 5 I(4) 4 scalar COM MYIS_P2 Scalar 5 I(4) 4 scalar COM MYIS_P3 Scalar 5 I(4) 4 scalar COM MYIS_P4 Scalar 5 I(4) 4 scalar COM MYIS_P5 Scalar 5 I(4) 4 scalar COM MYJE Scalar 15 I(4) 4 scalar COM 132 MYJE1 Scalar 15 I(4) 4 scalar COM MYJE1_P1 Scalar 17 I(4) 4 scalar COM MYJE1_P2 Scalar 17 I(4) 4 scalar COM MYJE1_P3 Scalar 17 I(4) 4 scalar COM MYJE1_P4 Scalar 17 I(4) 4 scalar COM MYJE2 Scalar 15 I(4) 4 scalar COM 152 MYJE2_P1 Scalar 18 I(4) 4 scalar COM MYJE2_P2 Scalar 18 I(4) 4 scalar COM MYJE2_P3 Scalar 18 I(4) 4 scalar COM MYJE2_P4 Scalar 18 I(4) 4 scalar COM MYJE3 Scalar 15 I(4) 4 scalar COM MYJE3_P4 Scalar 18 I(4) 4 scalar COM MYJE4 Scalar 15 I(4) 4 scalar COM MYJE4_P1 Scalar 19 I(4) 4 scalar COM MYJE4_P4 Scalar 19 I(4) 4 scalar COM MYJE5 Scalar 15 I(4) 4 scalar COM MYJE5_P1 Scalar 19 I(4) 4 scalar COM MYJE5_P2 Scalar 19 I(4) 4 scalar COM MYJE_P1 Scalar 16 I(4) 4 scalar COM MYJE_P2 Scalar 16 I(4) 4 scalar COM MYJE_P3 Scalar 16 I(4) 4 scalar COM MYJE_P4 Scalar 16 I(4) 4 scalar COM MYJE_P5 Scalar 16 I(4) 4 scalar COM MYJS Scalar 10 I(4) 4 scalar COM 132 MYJS1 Scalar 10 I(4) 4 scalar COM MYJS1_P1 Scalar 12 I(4) 4 scalar COM MYJS1_P2 Scalar 12 I(4) 4 scalar COM MYJS1_P3 Scalar 12 I(4) 4 scalar COM MYJS1_P4 Scalar 12 I(4) 4 scalar COM MYJS2 Scalar 10 I(4) 4 scalar COM 152 MYJS2_P1 Scalar 13 I(4) 4 scalar COM MYJS2_P2 Scalar 13 I(4) 4 scalar COM MYJS2_P3 Scalar 13 I(4) 4 scalar COM MYJS2_P4 Scalar 13 I(4) 4 scalar COM Page 17 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References MYJS3 Scalar 10 I(4) 4 scalar COM MYJS3_P4 Scalar 13 I(4) 4 scalar COM MYJS4 Scalar 10 I(4) 4 scalar COM MYJS4_P1 Scalar 14 I(4) 4 scalar COM MYJS4_P4 Scalar 14 I(4) 4 scalar COM MYJS5 Scalar 10 I(4) 4 scalar COM MYJS5_P1 Scalar 14 I(4) 4 scalar COM MYJS5_P2 Scalar 14 I(4) 4 scalar COM MYJS_P1 Scalar 11 I(4) 4 scalar COM MYJS_P2 Scalar 11 I(4) 4 scalar COM MYJS_P3 Scalar 11 I(4) 4 scalar COM MYJS_P4 Scalar 11 I(4) 4 scalar COM MYJS_P5 Scalar 11 I(4) 4 scalar COM MYPE Scalar 2 I(4) 4 scalar COM MY_E Scalar 20 I(4) 4 scalar COM MY_IE_GLB Scalar 2 I(4) 4 scalar COM MY_IE_LOC Scalar 3 I(4) 4 scalar COM MY_IS_GLB Scalar 2 I(4) 4 scalar COM MY_IS_LOC Scalar 3 I(4) 4 scalar COM MY_JE_GLB Scalar 2 I(4) 4 scalar COM MY_JE_LOC Scalar 3 I(4) 4 scalar COM MY_JS_GLB Scalar 2 I(4) 4 scalar COM MY_JS_LOC Scalar 3 I(4) 4 scalar COM MY_N Scalar 20 I(4) 4 scalar COM MY_NE Scalar 21 I(4) 4 scalar COM MY_NEB Scalar 21 I(4) 4 1 8 COM MY_NW Scalar 21 I(4) 4 scalar COM MY_S Scalar 20 I(4) 4 scalar COM MY_SE Scalar 21 I(4) 4 scalar COM MY_SW Scalar 21 I(4) 4 scalar COM MY_W Scalar 20 I(4) 4 scalar COM N Local 145 I(4) 4 scalar 145,146,147,148,201,203,204,250,25 1,293,355,356,364,365,526,527,528, 544,545 NBC Scalar 5 I(4) 4 scalar COM NBOCO Scalar 6 I(4) 4 scalar COM NCNVC Scalar 7 I(4) 4 scalar COM NCP Scalar 6 I(4) 4 scalar COM NEST Scalar 7 I(4) 4 scalar COM NFCST Scalar 5 I(4) 4 scalar COM NHEAT Scalar 3 I(4) 4 scalar COM NOZ Local 74 L(4) 4 scalar NPES Scalar 2 I(4) 4 scalar COM NPHS Scalar 7 I(4) 4 scalar COM 111 NPRE Local 151 I(4) 4 scalar 151,163,164,165,201 NPREC Scalar 6 I(4) 4 scalar COM NRADL Scalar 7 I(4) 4 scalar COM NRADS Scalar 7 I(4) 4 scalar COM NSHDE Scalar 6 I(4) 4 scalar COM NSTART Scalar 6 I(4) 4 scalar COM NTDDMP Scalar 6 I(4) 4 scalar COM NTSD Scalar 6 I(4) 4 scalar COM NTSTM Scalar 6 I(4) 4 scalar COM PD Scalar 3 R(4) 4 2 17955 COM 206 PDSL Local 206 R(4) 4 scalar 206,207,276,432 Page 18 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References PHOUR Scalar 2 R(4) 4 2 17955 COM PHYS Common 4 1788916 PID Local 233 R(4) 4 scalar 233,254,256,261,326,353,361,387,39 7,474,481,508,515 PL Scalar 7 R(4) 4 scalar COM PLQ Scalar 8 R(4) 4 scalar COM PPR Local 229 R(4) 4 scalar 229,349,355,383,389,459,462,463,46 4,475,495,500,509 PPS Local 230 R(4) 4 scalar 230,350,356,384,390,460,467,468,46 9,476,496,498,499,501,510 PPTASM Common 1 3878292 PPTDAT Scalar 3 R(4) 4 3 53865 COM PQ0 Param 53 R(4) 4 scalar 58,276,432 PQ0C Param 58 R(4) 4 scalar PRACW Local 222 R(4) 4 scalar 222,376,377,380,383,521 PRAUT Local 220 R(4) 4 scalar 220,342,343,346,349,521 PREC Scalar 12 R(4) 4 2 17955 COM 548 PRECPD Subr 1 PRECRK Local 250 R(4) 4 scalar 250,253,364,376,401,408,451,452,46 2,463 PRECRL Local 218 R(4) 4 scalar 218,355,389,401,475,509,527 PRECRL1 Local 104 R(4) 4 1 15327 147,250,355,364,527,544 PRECRS Local 544 R(4) 4 scalar 544,547,548,555 PRECSK Local 251 R(4) 4 scalar 251,253,365,369,376,402,413,451,45 3,467,468 PRECSL Local 219 R(4) 4 scalar 219,356,390,402,476,483,510,528 PRECSL1 Local 104 R(4) 4 1 15327 148,251,356,365,528,545 PRECSS Local 545 R(4) 4 scalar 545,547,548,555,557 PSACI Local 223 R(4) 4 scalar 223,369,370,373,384,522 PSAUT Local 221 R(4) 4 scalar 221,334,335,338,350,522 PSHLTR Scalar 19 R(4) 4 2 17955 COM PSM Local 226 R(4) 4 scalar 226,488,492,495,496,503 PSM1 Local 227 R(4) 4 scalar 227,486,488,490,501,503,517,520 PSM2 Local 228 R(4) 4 scalar 228,487,488,491,502 PT Scalar 5 R(4) 4 scalar COM 276,432 PTBL Scalar 25 R(4) 4 2 10184 COM PVRBLS Common 2 5530140 Q Scalar 6 R(4) 4 3 807975 COM 137,175,571 Q10 Scalar 16 R(4) 4 2 17955 COM Q2 Scalar 21 R(4) 4 3 807975 COM Q30 Scalar 14 R(4) 4 2 17955 COM QC Local 301 R(4) 4 scalar 301,303,306 QCT Local 438 R(4) 4 scalar 438,440,444,445,450 QI Local 278 R(4) 4 scalar 278,280,283,301 QINT Local 279 R(4) 4 scalar 279,280,301 QINTT Local 435 R(4) 4 scalar 435,436,438 QIT Local 434 R(4) 4 scalar 434,436,438 QK Local 265 R(4) 4 scalar 265,283,306,398 QQ Local 239 R(4) 4 scalar 239,265,398,420,445,519,537 QS Scalar 6 R(4) 4 2 17955 COM QS0 Scalar 10 R(4) 4 1 134 COM QSHLTR Scalar 18 R(4) 4 2 17955 COM QTEMP Local 420 R(4) 4 scalar 420,444 QW Local 276 R(4) 4 scalar 276,278,279 QWBS Scalar 9 R(4) 4 2 17955 COM Page 19 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References QWT Local 432 R(4) 4 scalar 432,434,435 QZ0 Scalar 5 R(4) 4 2 17955 COM Q_T Local 97 R(4) 4 3 807975 175,239,537,571 R Param 53 R(4) 4 scalar RADIN Scalar 14 R(4) 4 2 17955 COM RADOT Scalar 14 R(4) 4 2 17955 COM RCONDE Local 236 R(4) 4 scalar 236,263,464,469,501 RCP Param 58 R(4) 4 scalar 419,517,520 RDP Scalar 7 R(4) 4 scalar COM RDPQ Scalar 8 R(4) 4 scalar COM RDQ Scalar 7 R(4) 4 scalar COM RDTH Scalar 7 R(4) 4 scalar COM RDTHE Scalar 7 R(4) 4 scalar COM RDTHEQ Scalar 8 R(4) 4 scalar COM RDTPH Local 112 R(4) 4 scalar 112,535 RES Scalar 7 R(4) 4 2 17955 COM 206 RESTRT Scalar 74 L(4) 4 scalar COM RF Scalar 8 R(4) 4 2 17955 COM ROI Scalar 6 R(4) 4 scalar COM ROS Scalar 6 R(4) 4 scalar COM ROW Param 54 R(4) 4 scalar 58,59 RPRS Local 451 R(4) 4 scalar 451,452,453 RQ Local 304 R(4) 4 scalar 304,306,309,312,407 RQKLL Local 312 R(4) 4 scalar 312,313 RQT Local 441 R(4) 4 scalar 441,444,448 RQTT Local 442 R(4) 4 scalar 442,445,450 RROG Param 58 R(4) 4 scalar RROW Param 59 R(4) 4 scalar 544,545 RTWODT Local 114 R(4) 4 scalar 114,338,346,373,380,450 RUN Scalar 74 L(4) 4 scalar COM SGETMO Subr 173 173,175,177,179,181,183,569,571,57 3,575,577 SI Scalar 10 R(4) 4 2 17955 COM SICE Scalar 4 R(4) 4 2 17955 COM SIGMA Scalar 74 L(4) 4 scalar COM SIGT4 Scalar 23 R(4) 4 2 17955 COM SINGLRST Local 2 L(4) 4 scalar SM Scalar 4 R(4) 4 2 17955 COM SNO Scalar 10 R(4) 4 2 17955 COM SQRT Func 313 scalar 313 SQS Scalar 10 R(4) 4 1 134 COM SR Scalar 5 R(4) 4 2 17955 COM 557,559 SST Scalar 20 R(4) 4 2 17955 COM STHE Scalar 10 R(4) 4 1 76 COM STHEQ Scalar 11 R(4) 4 1 152 COM SUBPOST Scalar 2 L(4) 4 scalar COM T Scalar 5 R(4) 4 3 807975 COM 156,173,569 TCUCN Scalar 5 R(4) 4 3 807975 COM TDTD Scalar 5 R(4) 4 scalar COM TDTQ2 Scalar 5 R(4) 4 scalar COM TEMP1 Scalar 43 R(4) 4 2 91719 COM TEMP10 Scalar 45 R(4) 4 2 91719 COM TEMP11 Scalar 45 R(4) 4 2 91719 COM TEMP12 Scalar 45 R(4) 4 2 91719 COM TEMP13 Scalar 46 R(4) 4 2 91719 COM Page 20 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References TEMP14 Scalar 46 R(4) 4 2 91719 COM TEMP15 Scalar 46 R(4) 4 2 91719 COM TEMP16 Scalar 46 R(4) 4 2 91719 COM TEMP2 Scalar 43 R(4) 4 2 91719 COM TEMP2X Scalar 51 R(4) 4 2 92971 COM TEMP3 Scalar 43 R(4) 4 2 91719 COM TEMP4 Scalar 43 R(4) 4 2 91719 COM TEMP5 Scalar 44 R(4) 4 2 91719 COM TEMP6 Scalar 44 R(4) 4 2 91719 COM TEMP7 Scalar 44 R(4) 4 2 91719 COM TEMP8 Scalar 44 R(4) 4 2 91719 COM TEMP9 Scalar 45 R(4) 4 2 91719 COM TEMPCOM Common 42 6603768 TG Scalar 19 R(4) 4 2 17955 COM TH10 Scalar 16 R(4) 4 2 17955 COM TH30 Scalar 14 R(4) 4 2 17955 COM THE0 Scalar 10 R(4) 4 1 76 COM THE0Q Scalar 11 R(4) 4 1 152 COM THEAT Scalar 3 R(4) 4 scalar COM THL Scalar 7 R(4) 4 scalar COM THS Scalar 6 R(4) 4 2 17955 COM THZ0 Scalar 5 R(4) 4 2 17955 COM TK Local 264 R(4) 4 scalar 264,266,277 TLAT Scalar 2 R(4) 4 3 807975 COM 183,577 TLAT_T Local 101 R(4) 4 3 807975 183,242,539,577 TMT0 Local 266 R(4) 4 scalar 266,267,271,278,280,282,291,296,33 2,367,399,410,485,486 TMT0K Local 399 R(4) 4 scalar 399,400,413 TMT0T Local 421 R(4) 4 scalar 421,422,423,427,434,436 TMT15 Local 267 R(4) 4 scalar 267,279 TMT15T Local 423 R(4) 4 scalar 423,435 TOPO Common 50 17478548 TOTPPT Local 555 R(4) 4 scalar 555,556,557 TRAIN Scalar 4 R(4) 4 3 807975 COM 181,575 TRAIN_T Local 98 R(4) 4 3 807975 181,535,575 TRESH Param 53 R(4) 4 scalar 58 TSHLTR Scalar 18 R(4) 4 2 17955 COM TT Local 238 R(4) 4 scalar 238,245,264,419,517,518,535,536 TTBL Scalar 25 R(4) 4 2 10184 COM TTBLQ Scalar 25 R(4) 4 2 66880 COM TTEMP Local 156 R(4) 4 scalar 156,158,245,246,419,421,433 TTLAT Local 242 R(4) 4 scalar 242,520,522,539 TTVG Scalar 51 R(4) 4 2 92971 COM TWBS Scalar 9 R(4) 4 2 17955 COM TWODT Local 113 R(4) 4 scalar 113,114,207,335,343,370,377,418,51 6 T_T Local 96 R(4) 4 3 807975 173,238,535,536,569 U Scalar 5 R(4) 4 3 807975 COM U00 Scalar 3 R(4) 4 2 17955 COM 209 U00IJ Local 209 R(4) 4 scalar 209,244 U00KL Local 244 R(4) 4 scalar 244,283,309,313,406,407,448,450 U00KLT Local 406 R(4) 4 scalar 406,413 U10 Scalar 17 R(4) 4 2 17955 COM U30 Scalar 15 R(4) 4 2 17955 COM UL Scalar 4 R(4) 4 1 90 COM 214,244 Page 21 Source Listing PRECPD 2014-12-17 20:47 Symbol Table PRECPD.F Name Object Declared Type Bytes Dimen Elements Attributes References ULL Local 214 R(4) 4 scalar 214 US Local 116 R(4) 4 scalar 116,312,313 USTAR Scalar 3 R(4) 4 2 17955 COM UTIM Local 142 R(4) 4 scalar 142,244 UZ0 Scalar 4 R(4) 4 2 17955 COM V Scalar 6 R(4) 4 3 807975 COM V10 Scalar 17 R(4) 4 2 17955 COM V30 Scalar 15 R(4) 4 2 17955 COM VAPINC Scalar 5 R(4) 4 2 17955 COM VAPINC7 Scalar 5 R(4) 4 2 17955 COM VBM2 Scalar 3 R(4) 4 2 17955 COM VBM3 Scalar 3 R(4) 4 2 17955 COM VRBLS Common 2 13143060 VTM Scalar 7 R(4) 4 3 807975 COM VZ0 Scalar 4 R(4) 4 2 17955 COM WFIX Local 157 R(4) 4 scalar 157,158,216,246 WMIN Local 158 R(4) 4 scalar 158,159 WMINK Local 246 R(4) 4 scalar 246,253,329 WW Local 240 R(4) 4 scalar 240,253,284,293,317,318,328,354,36 3,388,482,538 Z0 Scalar 3 R(4) 4 2 17955 COM Page 22 Source Listing PRECPD 2014-12-17 20:47 Subprograms/Common Blocks PRECPD.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References ACMCLH Common 2 6463824 CLDWTR Common 2 3447720 CTLBLK Common 4 112 GLB_TABLE Common 37 128 LOOPS Common 2 158840 MAPPINGS Common 56 5024 MASKS Common 2 6894720 MPPCOM Common 1 1464 PHYS Common 4 1788916 PPTASM Common 1 3878292 PRECPD Subr 1 PVRBLS Common 2 5530140 TEMPCOM Common 42 6603768 TOPO Common 50 17478548 VRBLS Common 2 13143060 COMPILER OPTIONS BEING USED -align nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores no -auto -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ Page 23 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F -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 precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits nono_precise_libm_functions -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 : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, Page 24 Source Listing PRECPD 2014-12-17 20:47 PRECPD.F /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : PRECPD.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100