Page 1 Source Listing CUCNVC 2014-12-17 20:46 /tmpnwprd/ifortCnZaaL.i 1 # 1 "CUCNVC.F" 2 SUBROUTINE CUCNVC 3 C ****************************************************************** 4 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 5 C . . . 6 C SUBPROGRAM: CUCNVC CONVECTIVE PRECIPITATION PARAMETERIZATION 7 C PRGRMMR: JANJIC ORG: W/NP2 DATE: 93-11-02 8 C 9 C ABSTRACT: 10 C CUCNVC CALCULATES THE SUB-GRID SCALE CONVECTION INCLUDING 11 C DEEP AND SHALLOW CONVECTIVE CLOUDS FOLLOWING THE SCHEME DESCRIBED 12 C BY JANJIC (1994) BUT WITH SIGNIFICANT MODIFICATIONS. IN ADDITION, 13 C THE LATENT HEAT RELEASE AND MOISTURE CHANGE DUE TO PRECIPITATING 14 C AND NON-PRECIPITATING CLOUDS ARE COMPUTED. 15 C 16 C 17 C PROGRAM HISTORY LOG: 18 C 87-09-?? JANJIC - ORIGINATOR 19 C 90-11-21 JANJIC - TWO SETS OF DSP PROFILES (FAST AND SLOW) 20 C REPLACE THE ORIGINAL ONE SET 21 C 95-03-25 BLACK - CONVERSION FROM 1-D TO 2-D IN HORIZONTAL 22 C 96-03-28 BLACK - ADDED EXTERNAL EDGE 23 C 98-11-02 BLACK - MODIFIED FOR DISTRIBUTED MEMORY 24 C 25 C USAGE: CALL CUCNVC FROM MAIN PROGRAM EBU 26 C 27 C INPUT ARGUMENT LIST: 28 C NONE 29 C 30 C OUTPUT ARGUMENT LIST: 31 C NONE 32 C 33 C OUTPUT FILES: 34 C NONE 35 C 36 C SUBPROGRAMS CALLED: 37 C 38 C UNIQUE: 39 C TTBLEX 40 C 41 C LIBRARY: 42 C NONE 43 C 44 C COMMON BLOCKS: CTLBLK 45 C LOOPS 46 C MASKS 47 C PHYS 48 C VRBLS 49 C CNVCLD 50 C PVRBLS 51 C ACMCLH 52 C INDX 53 C 54 C ATTRIBUTES: 55 C LANGUAGE: FORTRAN 90 56 C MACHINE : IBM SP 57 C$$$ Page 2 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 58 C ****************************************************************** 59 C * * 60 C * REFERENCES: * 61 C * * 62 C * JANJIC, Z.I., 1994: THE STEP-MOUNTAIN ETA COORDINATE MODEL: * 63 C * FURTHER DEVELOPMENTS OF THE CONVECTION, VISCOUS SUBLAYER * 64 C * AND TURBULENCE CLOSURE SCHEMES. MONTHLY WEATHER REVIEW, * 65 C * VOL. 122, 927-945. * 66 C * * 67 C ****************************************************************** 68 C *** WARNING: THIS SUBROUTINE WILL NOT WORK IF LM.LT.12; 69 C MUST BE CALLED IN THE SAME STEP WITH PROFQ2 BECAUSE PROFQ 70 C DEFINES APE; 71 C---------------------------------------------------------------------- 72 INCLUDE "cuparm" 73 C---------------------------------------------------------------------- 164 INCLUDE "parmeta" 187 INCLUDE "parm.tbl" 189 INCLUDE "mpp.h" 190 C---------------------------------------------------------------------- 249 P A R A M E T E R 250 CVVVVVVVVVV INSTABILITY FOR TOO LARGE LSH VVVVVVVVVVVVVVVVVVVVVVVVVVVVV 251 & (KSMUD=0,NROW= 0) 252 C & (KSMUD=0,NROW= 5) 253 CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 254 P A R A M E T E R 255 & (IMJM=IM*JM-JM/2,JAM=6+2*(JM-10) 256 &, IMJM_LOC=IDIM2*JDIM2 257 CVVVVVVVVVV INSTABILITY FOR TOO LARGE LSH VVVVVVVVVVVVVVVVVVVVVVVVVVVVV 258 C &, LP1=LM+1,LM1=LM-1,LNO=1,LSH=LM/3-1,LSHU=LM/2-1,LQM=LM/5,KBM=3 259 C &, LP1=LM+1,LM1=LM-1,LNO=1,LSH=LM/3 ,LSHU=LM/2-1,LQM=LM/5,KBM=3 260 C &, LP1=LM+1,LM1=LM-1,LNO=3,LSH=LM/3 ,LSHU=LM/2-1,LQM=LM/5,KBM=3 261 C &, LP1=LM+1,LM1=LM-1,LNO=2,LSH=LM/3-1,LSHU=LM/2-1,LQM=LM/5,KBM=3 262 C &, LP1=LM+1,LM1=LM-1,LNO=2,LSH=LM/3-2,LSHU=LM/2-1,LQM=LM/5,KBM=3 263 &, LP1=LM+1,LM1=LM-1) 264 CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 265 C---------------------------------------------------------------------- 266 L O G I C A L 267 & RUN,FIRST,RESTRT,SIGMA 268 C---------------------------------------------------------------------- 269 INCLUDE "CTLBLK.comm" 270 C----------------------------------------------------------------------- 280 INCLUDE "LOOPS.comm" 281 C----------------------------------------------------------------------- 288 INCLUDE "MASKS.comm" 289 C----------------------------------------------------------------------- 300 INCLUDE "PHYS.comm" 301 C----------------------------------------------------------------------- 329 INCLUDE "VRBLS.comm" 330 C----------------------------------------------------------------------- 340 INCLUDE "CNVCLD.comm" 341 C----------------------------------------------------------------------- 348 INCLUDE "PVRBLS.comm" 349 C----------------------------------------------------------------------- 379 INCLUDE "ACMCLH.comm" 380 C----------------------------------------------------------------------- 387 INCLUDE "INDX.comm" Page 3 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 388 CYL--------------------------------------------------------------------- 400 INCLUDE "PPTASM.comm" 401 CYL--------------------------------------------------------------------- 425 D I M E N S I O N 426 & TREFK (LM),QREFK (LM),PK (LM),APEK (LM),TK (LM) 427 &,THSK (LM),PSK (LM),APESK (LM),QK (LM),THERK (LM) 428 &,THVREF(LM),THEVRF(LM),THVMOD(LM),DIFT (LM),DIFQ (LM) 429 &,QSATK (LM),FPK (LM) 430 &,NTOPD (LM),NBOTD (LM),NTOPS (LM),NBOTS (LM) 431 &,NDPTHD(LM),NDPTHS(LM) 432 C 433 D I M E N S I O N 434 & LTOP (IDIM1:IDIM2,JDIM1:JDIM2),LBOT (IDIM1:IDIM2,JDIM1:JDIM2) 435 &,PTOP (IDIM1:IDIM2,JDIM1:JDIM2),PBOT (IDIM1:IDIM2,JDIM1:JDIM2) 436 &,IPTB (IDIM1:IDIM2,JDIM1:JDIM2),ITHTB (IDIM1:IDIM2,JDIM1:JDIM2) 437 &,PDSL (IDIM1:IDIM2,JDIM1:JDIM2),APEBT (IDIM1:IDIM2,JDIM1:JDIM2) 438 &,TBT (IDIM1:IDIM2,JDIM1:JDIM2),Q2BT (IDIM1:IDIM2,JDIM1:JDIM2) 439 &,QQ (IDIM1:IDIM2,JDIM1:JDIM2),PP (IDIM1:IDIM2,JDIM1:JDIM2) 440 &,PSP (IDIM1:IDIM2,JDIM1:JDIM2),THBT (IDIM1:IDIM2,JDIM1:JDIM2) 441 &,THESP (IDIM1:IDIM2,JDIM1:JDIM2),P (IDIM1:IDIM2,JDIM1:JDIM2) 442 &,BTH (IDIM1:IDIM2,JDIM1:JDIM2),STH (IDIM1:IDIM2,JDIM1:JDIM2) 443 &,T00 (IDIM1:IDIM2,JDIM1:JDIM2),T10 (IDIM1:IDIM2,JDIM1:JDIM2) 444 &,T01 (IDIM1:IDIM2,JDIM1:JDIM2),T11 (IDIM1:IDIM2,JDIM1:JDIM2) 445 &,WF1 (IDIM1:IDIM2,JDIM1:JDIM2),WF2 (IDIM1:IDIM2,JDIM1:JDIM2) 446 &,WF3 (IDIM1:IDIM2,JDIM1:JDIM2),WF4 (IDIM1:IDIM2,JDIM1:JDIM2) 447 &,PRECOL(IDIM1:IDIM2,JDIM1:JDIM2) 448 C 449 &,IBUOY (IMJM_LOC),JBUOY (IMJM_LOC) 450 &,IDEEP (IMJM_LOC),JDEEP (IMJM_LOC) 451 &,ISHAL (IMJM_LOC),JSHAL (IMJM_LOC) 452 &,ILRES (IMJM_LOC),JLRES (IMJM_LOC) 453 &,IHRES (IMJM_LOC),JHRES (IMJM_LOC) 454 C 455 &,APE (IDIM1:IDIM2,JDIM1:JDIM2,LM) 456 &,TREF (IDIM1:IDIM2,JDIM1:JDIM2,LM) 457 &,TMOD (IDIM1:IDIM2,JDIM1:JDIM2,LM) 458 &,QMOD (IDIM1:IDIM2,JDIM1:JDIM2,LM) 459 C 460 D I M E N S I O N 461 & DSPB (IDIM1:IDIM2,JDIM1:JDIM2),DSP0 (IDIM1:IDIM2,JDIM1:JDIM2) 462 &,DSPT (IDIM1:IDIM2,JDIM1:JDIM2) 463 &,TL (IDIM1:IDIM2,JDIM1:JDIM2),QL (IDIM1:IDIM2,JDIM1:JDIM2) 464 &,TNE (IDIM1:IDIM2,JDIM1:JDIM2),TSE (IDIM1:IDIM2,JDIM1:JDIM2) 465 &,QNE (IDIM1:IDIM2,JDIM1:JDIM2),QSE (IDIM1:IDIM2,JDIM1:JDIM2) 466 C----------------------------------------------------------------------- 467 CALL ZERO2(DSP0) 468 CALL ZERO2(DSPB) 469 CALL ZERO2(DSPT) 470 CALL ZERO2(PSP) 471 C---------------------------------------------------------------------- 472 AVCNVC=AVCNVC+1. 473 ACUTIM=ACUTIM+1. 474 DTCNVC=NCNVC*DT 475 RDTCNVC=1./DTCNVC 476 TAUK=DTCNVC/TREL 477 c possible future change for shallow convection 478 c TAUKSC=DTCNVC/(5.*TREL) Page 4 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 479 C----------------------------------------------------------------------- 480 C************************** DIAGNOSTICS ******************************** 481 C----------------------------------------------------------------------- 482 DO L=1,LM 483 NTOPD (L)=0 484 NBOTD (L)=0 485 NTOPS (L)=0 486 NBOTS (L)=0 487 NDPTHS(L)=0 488 NDPTHD(L)=0 489 ENDDO 490 C*********************************************************************** 491 C---------------------------PREPARATIONS-------------------------------- 492 C----------------------------------------------------------------------- 493 !$omp parallel do 494 DO 120 J=MYJS,MYJE 495 DO 120 I=MYIS,MYIE 496 LBOT (I,J)=LMH(I,J) 497 THESP(I,J)=0. 498 PDSL (I,J)=RES(I,J)*PD(I,J) 499 PRECOL(I,J)=0. 500 TREF(I,J,1)=T(I,J,1) 501 120 CONTINUE 502 C----------------------------------------------------------------------- 503 C--------------PADDING SPECIFIC HUMIDITY IF TOO SMALL------------------- 504 C RESTORE APE TO SCRATCH ARRAY 505 C----------------------------------------------------------------------- 506 !$omp parallel do private(apests) 507 DO 130 L=1,LM 508 DO J=MYJS,MYJE 509 DO I=MYIS,MYIE 510 APESTS=PDSL(I,J)*AETA(L)+PT 511 APE(I,J,L)=(1.E5/APESTS)**CAPA 512 IF(Q(I,J,L).LT.EPSQ)Q(I,J,L)=HTM(I,J,L)*EPSQ 513 ENDDO 514 ENDDO 515 130 CONTINUE 516 C----------------------------------------------------------------------- 517 C--------------SEARCH FOR MAXIMUM BUOYANCY LEVEL------------------------ 518 C----------------------------------------------------------------------- 519 DO 170 KB=1,LM 520 C----------------------------------------------------------------------- 521 C--------------TRIAL MAXIMUM BUOYANCY LEVEL VARIABLES------------------- 522 C----------------------------------------------------------------------- 523 !$omp parallel do 524 !$omp& private(apesp,bq,bqs00k,bqs10k,iq,it,ittb,ittbk,iqtb, 525 !$omp& lmhk,p00k,p01k,p10k,p11k,pkl,pp1,psfck,qbt,qq1, 526 !$omp& sq,sqs00k,sqs10k,tpsp,tq,tth,tthbt,tthes) 527 DO 155 J=MYJS2,MYJE2 528 DO 150 I=MYIS1,MYIE1 529 C 530 PKL=AETA(KB)*PDSL(I,J)+PT 531 LMHK=LMH(I,J) 532 PSFCK=AETA(LMHK)*PDSL(I,J)+PT 533 c now searching over a scaled depth in finding the parcel 534 c with the max theta-e instead of the old 130 mb 535 c IF(PKL.LT.PSFCK-PBM)GO TO 150 Page 5 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 536 c IF(KB.GT.LMHK)GO TO 150 537 IF(KB.LE.LMHK .AND. PKL.GE.0.80*PSFCK) THEN 538 QBT=Q(I,J,KB) 539 TTHBT=T(I,J,KB)*APE(I,J,KB) 540 TTH=(TTHBT-THL)*RDTH 541 QQ1=TTH-AINT(TTH) 542 ITTB=INT(TTH)+1 543 C--------------KEEPING INDICES WITHIN THE TABLE------------------------- 544 IF(ITTB.LT.1)THEN 545 ITTB=1 546 QQ1=0. 547 ENDIF 548 IF(ITTB.GE.JTB)THEN 549 ITTB=JTB-1 550 QQ1=0. 551 ENDIF 552 CONTINUE 553 C--------------BASE AND SCALING FACTOR FOR SPEC. HUMIDITY--------------- 554 ITTBK=ITTB 555 BQS00K=QS0(ITTBK) 556 SQS00K=SQS(ITTBK) 557 BQS10K=QS0(ITTBK+1) 558 SQS10K=SQS(ITTBK+1) 559 C--------------SCALING SPEC. HUMIDITY & TABLE INDEX--------------------- 560 BQ = (BQS10K-BQS00K)*QQ1+BQS00K 561 SQ = (SQS10K-SQS00K)*QQ1+SQS00K 562 TQ = (QBT-BQ)/SQ*RDQ 563 PP1 =TQ - AINT(TQ) 564 IQTB=INT(TQ)+1 565 C--------------KEEPING INDICES WITHIN THE TABLE------------------------- 566 IF(IQTB.LT.1)THEN 567 IQTB=1 568 PP1=0. 569 ENDIF 570 IF(IQTB.GE.ITB)THEN 571 IQTB=ITB-1 572 PP1=0. 573 ENDIF 574 C--------------SATURATION PRESSURE AT FOUR SURROUNDING TABLE PTS.------- 575 IQ=IQTB 576 IT=ITTB 577 P00K=PTBL(IQ ,IT ) 578 P10K=PTBL(IQ+1,IT ) 579 P01K=PTBL(IQ ,IT+1) 580 P11K=PTBL(IQ+1,IT+1) 581 C--------------SATURATION POINT VARIABLES AT THE BOTTOM----------------- 582 TPSP=P00K+(P10K-P00K)*PP1+(P01K-P00K)*QQ1 583 1 +(P00K-P10K-P01K+P11K)*PP1*QQ1 584 APESP=(1.E5/TPSP)**CAPA 585 TTHES=TTHBT*EXP(ELOCP*QBT*APESP/TTHBT) 586 C--------------CHECK FOR MAXIMUM BUOYANCY------------------------------- 587 IF(TTHES.GT.THESP(I,J))THEN 588 PSP (I,J)=TPSP 589 THBT (I,J)=TTHBT 590 THESP(I,J)=TTHES 591 ENDIF 592 ENDIF Page 6 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 593 150 CONTINUE 594 155 CONTINUE 595 C----------------------------------------------------------------------- 596 170 CONTINUE 597 C----------------------------------------------------------------------- 598 C---------CHOOSE CLOUD BASE AS MODEL LEVEL JUST BELOW PSP-------------- 599 C----------------------------------------------------------------------- 600 DO 240 L=1,LM1 601 AETAL=AETA(L) 602 C 603 !$omp parallel do 604 DO J=MYJS2,MYJE2 605 DO I=MYIS,MYIE 606 P(I,J)=PDSL(I,J)*AETAL+PT 607 IF(P(I,J).LT.PSP(I,J).AND.P(I,J).GE.PQM)LBOT(I,J)=L+1 608 ENDDO 609 ENDDO 610 C 611 240 CONTINUE 612 C*** 613 C*** WARNING: LBOT MUST NOT BE GT LMH(I,J)-1 IN SHALLOW CONVECTION 614 C*** MAKE SURE CLOUD BASE IS AT LEAST PONE ABOVE THE SURFACE 615 C*** 616 !$omp parallel do private(lmhij,psfck) 617 DO 250 J=MYJS2,MYJE2 618 DO 250 I=MYIS,MYIE 619 LMHIJ=LMH(I,J) 620 PBOT(I,J)=AETA(LBOT(I,J))*PDSL(I,J)+PT 621 PSFCK=AETA(LMHIJ)*PDSL(I,J)+PT 622 IF(PBOT(I,J).GE.PSFCK-PONE.OR.LBOT(I,J).GE.LMHIJ)THEN 623 C 624 DO L=1,LMHIJ-1 625 P(I,J)=AETA(L)*PDSL(I,J)+PT 626 IF(P(I,J).LT.PSFCK-PONE)LBOT(I,J)=L 627 ENDDO 628 C 629 PBOT(I,J)=AETA(LBOT(I,J))*PDSL(I,J)+PT 630 ENDIF 631 250 CONTINUE 632 C--------------CLOUD TOP COMPUTATION------------------------------------ 633 !$omp parallel do 634 DO J=MYJS,MYJE 635 DO I=MYIS,MYIE 636 LTOP(I,J)=LBOT(I,J) 637 PTOP(I,J)=PBOT(I,J) 638 ENDDO 639 ENDDO 640 C----------------------------------------------------------------------- 641 !$omp parallel do 642 !$omp& private(ihres,ilres,iptb,ithtb,jhres,jlres, 643 !$omp& knumh,knuml,pp,presk,qq) 644 C 645 DO 290 L=LM,1,-1 646 C 647 C--------------SCALING PRESSURE & TT TABLE INDEX------------------------ 648 KNUML=0 649 KNUMH=0 Page 7 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 650 C 651 DO 270 J=MYJS2,MYJE2 652 DO 270 I=MYIS1,MYIE1 653 PRESK=PDSL(I,J)*AETA(L)+PT 654 IF(PRESK.LT.PLQ)THEN 655 KNUML=KNUML+1 656 ILRES(KNUML)=I 657 JLRES(KNUML)=J 658 ELSE 659 KNUMH=KNUMH+1 660 IHRES(KNUMH)=I 661 JHRES(KNUMH)=J 662 ENDIF 663 270 CONTINUE 664 C*** 665 C*** COMPUTE PARCEL TEMPERATURE ALONG MOIST ADIABAT FOR PRESSUREPL 677 C** 678 IF(KNUMH.GT.0)THEN 679 CALL TTBLEX(TREF(IDIM1,JDIM1,L),TTBLQ,ITBQ,JTBQ,KNUMH 680 1, IHRES,JHRES,PDSL,AETA(L),HTM(IDIM1,JDIM1,L) 681 2, PT,PLQ,QQ(IDIM1,JDIM1),PP(IDIM1,JDIM1) 682 3, RDPQ,THE0Q,STHEQ,RDTHEQ 683 4, THESP(IDIM1,JDIM1),IPTB(IDIM1,JDIM1) 684 5, ITHTB(IDIM1,JDIM1)) 685 ENDIF 686 290 CONTINUE 687 C--------------BUOYANCY CHECK------------------------------------------- 688 DO 295 L=LM,1,-1 689 !$omp parallel do 690 DO J=MYJS2,MYJE2 691 DO I=MYIS1,MYIE1 692 IF(TREF(I,J,L).GT.T(I,J,L)-DTTOP)LTOP(I,J)=L 693 ENDDO 694 ENDDO 695 C----------------------------------------------------------------------- 696 295 CONTINUE 697 C-----------------CLOUD TOP PRESSURE------------------------------------ 698 !$omp parallel do 699 DO J=MYJS2,MYJE2 700 DO I=MYIS1,MYIE1 701 PTOP(I,J)=AETA(LTOP(I,J))*PDSL(I,J)+PT 702 ENDDO 703 ENDDO 704 C----------------------------------------------------------------------- 705 C--------------DEFINE AND SMOOTH DSPS AND CLDEFI------------------------ 706 C ************ UNIFIED OR SEPARATE LAND/SEA CONV *********************** Page 8 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 707 C----------------------------------------------------------------------- 708 C 709 IF(UNIS)THEN 710 !$omp parallel do private(efi) 711 DO J=MYJS,MYJE 712 DO I=MYIS,MYIE 713 EFI=CLDEFI(I,J) 714 DSPB(I,J)=(EFI-EFIMN)*SLOPBS+DSPBSS 715 DSP0(I,J)=(EFI-EFIMN)*SLOP0S+DSP0SS 716 DSPT(I,J)=(EFI-EFIMN)*SLOPTS+DSPTSS 717 ENDDO 718 ENDDO 719 ELSEIF(.NOT.UNIL)THEN 720 !$omp parallel do private(efi) 721 DO J=MYJS,MYJE 722 DO I=MYIS,MYIE 723 EFI=CLDEFI(I,J) 724 DSPB(I,J)=((EFI-EFIMN)*SLOPBS+DSPBSS)*SM(I,J) 725 1 +((EFI-EFIMN)*SLOPBL+DSPBSL)*(1.-SM(I,J)) 726 DSP0(I,J)=((EFI-EFIMN)*SLOP0S+DSP0SS)*SM(I,J) 727 1 +((EFI-EFIMN)*SLOP0L+DSP0SL)*(1.-SM(I,J)) 728 DSPT(I,J)=((EFI-EFIMN)*SLOPTS+DSPTSS)*SM(I,J) 729 1 +((EFI-EFIMN)*SLOPTL+DSPTSL)*(1.-SM(I,J)) 730 ENDDO 731 ENDDO 732 ELSE 733 !$omp parallel do private(efi) 734 DO J=MYJS,MYJE 735 DO I=MYIS,MYIE 736 EFI=CLDEFI(I,J) 737 DSPB(I,J)=((EFI-EFIMN)*SLOPBL+DSPBSL) 738 DSP0(I,J)=((EFI-EFIMN)*SLOP0L+DSP0SL) 739 DSPT(I,J)=((EFI-EFIMN)*SLOPTL+DSPTSL) 740 ENDDO 741 ENDDO 742 ENDIF 743 C 744 C--------------EXTENDING SEA STRUCTURES INLAND ALONG COASTLINE---------- 745 IF(NROW.GT.0.AND..NOT.UNIS.AND..NOT.UNIL)THEN 746 !$omp parallel do 747 DO J=MYJS,MYJE 748 DO I=MYIS,MYIE 749 WF1(I,J)=0. 750 WF2(I,J)=0. 751 WF3(I,J)=0. 752 WF4(I,J)=0. 753 ENDDO 754 ENDDO 755 C 756 KROW=NROW 757 C 758 DO 350 KOUNT=1,KROW 759 !$omp parallel do 760 DO 345 J=MYJS2,MYJE2 761 DO 345 I=MYIS1,MYIE1 762 WF1(I,J)=(DSPB(I+IHW(J),J-1)+DSPB(I+IHE(J),J-1) 763 1 +DSPB(I+IHW(J),J+1)+DSPB(I+IHE(J),J+1)+4.*DSPB(I,J)) Page 9 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 764 2 *HBM2(I,J)*0.125 765 WF2(I,J)=(DSP0(I+IHW(J),J-1)+DSP0(I+IHE(J),J-1) 766 1 +DSP0(I+IHW(J),J+1)+DSP0(I+IHE(J),J+1)+4.*DSP0(I,J)) 767 2 *HBM2(I,J)*0.125 768 WF3(I,J)=(DSPT(I+IHW(J),J-1)+DSPT(I+IHE(J),J-1) 769 1 +DSPT(I+IHW(J),J+1)+DSPT(I+IHE(J),J+1)+4.*DSPT(I,J)) 770 2 *HBM2(I,J)*0.125 771 WF4(I,J)=(CLDEFI(I+IHW(J),J-1)+CLDEFI(I+IHE(J),J-1) 772 1 +CLDEFI(I+IHW(J),J+1)+CLDEFI(I+IHE(J),J+1)+4.*CLDEFI(I,J)) 773 2 *HBM2(I,J)*0.125 774 345 CONTINUE 775 C 776 !$omp parallel do private(rsmk,smk) 777 DO J=MYJS,MYJE 778 DO I=MYIS,MYIE 779 SMK =SM(I,J) 780 RSMK=1.-SMK 781 DSPB (I,J)=DSPB (I,J)*SMK+WF1(I,J)*RSMK 782 DSP0 (I,J)=DSP0 (I,J)*SMK+WF2(I,J)*RSMK 783 DSPT (I,J)=DSPT (I,J)*SMK+WF3(I,J)*RSMK 784 CLDEFI(I,J)=CLDEFI(I,J)*SMK+WF4(I,J)*RSMK 785 ENDDO 786 ENDDO 787 C 788 350 CONTINUE 789 C----------------------------------------------------------------------- 790 ENDIF 791 C--------------INITIALIZE CHANGES OF T AND Q DUE TO CONVECTION---------- 792 !$omp parallel do 793 DO 360 L=1,LM 794 DO J=MYJS,MYJE 795 DO I=MYIS,MYIE 796 TMOD(I,J,L)=0. 797 QMOD(I,J,L)=0. 798 ENDDO 799 ENDDO 800 360 CONTINUE 801 C--------------CLEAN UP AND GATHER DEEP CONVECTION POINTS--------------- 802 !$omp parallel do 803 DO 380 J=MYJS2,MYJE2 804 DO 380 I=MYIS1,MYIE1 805 IF(LTOP(I,J).GE.LBOT(I,J))THEN 806 LBOT(I,J)=0 807 LTOP(I,J)=LBOT(I,J) 808 PTOP(I,J)=PBOT(I,J) 809 ENDIF 810 IF(HBM2(I,J).LT.0.90)THEN 811 LBOT(I,J)=0 812 LTOP(I,J)=LBOT(I,J) 813 PTOP(I,J)=PBOT(I,J) 814 ENDIF 815 IF(PTOP(I,J).GT.PBOT(I,J)-PNO.OR.LTOP(I,J).GT.LBOT(I,J)-2) 816 1 CLDEFI(I,J)=AVGEFI*SM(I,J)+STEFI*(1.-SM(I,J)) 817 380 CONTINUE 818 C 819 KHDEEP=0 820 PSHNEW=20000. Page 10 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 821 DO J=MYJS2,MYJE2 822 DO I=MYIS1,MYIE1 823 PSFCIJ=PD(I,J)+PT 824 c depth of cloud required to make the point a deep convection point 825 c is now a scaled value of the PSFC instead of 290 mb everywhere 826 DEPMIN=PSHNEW*PSFCIJ*1.E-5 827 DEPTH=PBOT(I,J)-PTOP(I,J) 828 IF(DEPTH .GE. DEPMIN) THEN 829 c IF(PTOP(I,J).LT.PBOT(I,J)-PSH)THEN 830 KHDEEP=KHDEEP+1 831 IDEEP(KHDEEP)=I 832 JDEEP(KHDEEP)=J 833 ENDIF 834 ENDDO 835 ENDDO 836 C************* HORIZONTAL LOOP FOR DEEP CONVECTION ********************* 837 !$omp parallel do 838 !$omp& private(apek,apekl,apekxx,apekxy,apesk,avrgt,avtgtl,dentpy, 839 !$omp& depmin,depth,depwl,dhdt,difq,difql,dift,diftl,drheat, 840 !$omp& drhdp,dsp,dsp0k,dspbk,dsptk,dthem,efi,fefi,hcorr, 841 !$omp& i,j,l0,l0m1,lb,lbm1,lbtk,lcor,lqm,lshu,ltp1,ltp2, 842 !$omp& ltpk,ltsh,pbtk,pk,pk0,pkb,pkl,pkt,preck,psfcij,psk, 843 !$omp& pthrs,ptpk,qk,qkl,qrefk,qsatk,rdp0t,rhh,rhl,rhmax, 844 !$omp& sumde,sumdp,therk,therkx,therky,thsk,thskl,tk,tkl, 845 !$omp& trefk,trefkx,tskl) 846 C*********************************************************************** 847 DO 600 N=1,KHDEEP 848 C*********************************************************************** 849 I=IDEEP(N) 850 J=JDEEP(N) 851 PSFCIJ=PD(I,J)+PT 852 LTPK=LTOP(I,J) 853 LBTK=LBOT(I,J) 854 CDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCD 855 CDCDCDCDCDCDC DEEP CONVECTION DCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCD 856 CDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCD 857 LB =LBTK 858 EFI =CLDEFI(I,J) 859 DSPBK=DSPB (I,J) 860 DSP0K=DSP0 (I,J) 861 DSPTK=DSPT (I,J) 862 C--------------INITIALIZE VARIABLES IN THE CONVECTIVE COLUMN------------ 863 C*** 864 C*** ONE SHOULD NOTE THAT THE VALUES ASSIGNED TO THE ARRAY TREFK 865 C*** IN THE 410 LOOP ARE REALLY ONLY RELEVANT IN ANCHORING THE 866 C*** REFERENCE TEMPERATURE PROFILE AT LEVEL LB. WHEN BUILDING THE 867 C*** REFERENCE PROFILE FROM CLOUD BASE, THEN ASSIGNING THE 868 C*** AMBIENT TEMPERATURE TO TREFK IS ACCEPTABLE. HOWEVER, WHEN 869 C*** BUILDING THE REFERENCE PROFILE FROM SOME OTHER LEVEL (SUCH AS 870 C*** ONE LEVEL ABOVE THE GROUND), THEN TREFK SHOULD BE FILLED WITH 871 C*** THE TEMPERATURES IN TREF(I,J,L) WHICH ARE THE TEMPERATURES OF 872 C*** THE MOIST ADIABAT THROUGH CLOUD BASE. BY THE TIME THE LINE 873 C*** NUMBERED 450 HAS BEEN REACHED, TREFK ACTUALLY DOES HOLD THE 874 C*** REFERENCE TEMPERATURE PROFILE. 875 C*** 876 DO 410 L=1,LM 877 DIFT (L)=0. Page 11 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 878 DIFQ (L)=0. 879 TKL =T(I,J,L) 880 TK (L)=TKL 881 TREFK (L)=TKL 882 QKL =Q(I,J,L) 883 QK (L)=QKL 884 QREFK (L)=QKL 885 PKL =AETA(L)*PDSL(I,J)+PT 886 PK (L)=PKL 887 PSK (L)=PKL 888 APEKL =APE(I,J,L) 889 APEK (L)=APEKL 890 THERK (L)=TREF(I,J,L)*APEKL 891 C 410 THVMOD(L)=(QKL*0.608+1.)*TKL*APEKL 892 410 CONTINUE 893 C--------------DEEP CONVECTION REFERENCE TEMPERATURE PROFILE------------ 894 LTP1=LTPK+1 895 LBM1=LB-1 896 PKB=PK(LB) 897 PKT=PK(LTPK) 898 C--------------TEMPERATURE REFERENCE PROFILE BELOW FREEZING LEVEL------- 899 L0=LB 900 PK0=PK(LB) 901 TREFKX=TREFK(LB) 902 THERKX=THERK(LB) 903 APEKXX=APEK(LB) 904 THERKY=THERK(LBM1) 905 APEKXY=APEK(LBM1) 906 C 907 DO 420 L=LBM1,LTPK,-1 908 IF(T(I,J,L+1).LT.TFRZ)GO TO 430 909 STABDL=STABD 910 TREFKX=((THERKY-THERKX)*STABDL 911 1 +TREFKX*APEKXX)/APEKXY 912 TREFK(L)=TREFKX 913 APEKXX=APEKXY 914 THERKX=THERKY 915 APEKXY=APEK(L-1) 916 THERKY=THERK(L-1) 917 L0=L 918 PK0=PK(L0) 919 420 CONTINUE 920 C--------------FREEZING LEVEL AT OR ABOVE THE CLOUD TOP----------------- 921 L0M1=L0-1 922 GO TO 450 923 C--------------TEMPERATURE REFERENCE PROFILE ABOVE FREEZING LEVEL------- 924 430 L0M1=L0-1 925 RDP0T=1./(PK0-PKT) 926 DTHEM=THERK(L0)-TREFK(L0)*APEK(L0) 927 C 928 DO L=LTPK,L0M1 929 TREFK(L)=(THERK(L)-(PK(L)-PKT)*DTHEM*RDP0T)/APEK(L) 930 ENDDO 931 C----------------------------------------------------------------------- 932 C--------------DEEP CONVECTION REFERENCE HUMIDITY PROFILE--------------- 933 C----------------------------------------------------------------------- 934 c reference profile had been too dry in the case where the cloud Page 12 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 935 c base was close to the freezing level - this is now changed 936 450 DEPTH=PFRZ*PSFCIJ*1.E-5 937 DEPWL=PKB-PK0 938 DO 460 L=LTPK,LB 939 C----------------------------------------------------------------------- 940 C--------------SATURATION PRESSURE DIFFERENCE--------------------------- 941 C----------------------------------------------------------------------- 942 c IF(PKB-PK0.GE.PFRZ)THEN 943 IF(DEPWL .GE. DEPTH) THEN 944 IF(L.LT.L0)THEN 945 DSP=((PK0-PK(L))*DSPTK+(PK(L)-PKT)*DSP0K)/(PK0-PKT) 946 ELSE 947 DSP=((PKB-PK(L))*DSP0K+(PK(L)-PK0)*DSPBK)/(PKB-PK0) 948 ENDIF 949 ELSE 950 c DSP=DSPC 951 DSP=DSP0K 952 IF(L.LT.L0) 953 1 DSP=( (PK0-PK(L))*DSPTK+(PK(L)-PKT)*DSP0K)/(PK0-PKT) 954 ENDIF 955 956 C--------------HUMIDITY PROFILE----------------------------------------- 957 IF(PK(L).GT.PQM)THEN 958 PSK(L)=PK(L)+DSP 959 APESK(L)=(1.E5/PSK(L))**CAPA 960 THSK(L)=TREFK(L)*APEK(L) 961 QREFK(L)=PQ0/PSK(L)*EXP(A2*(THSK(L)-A3*APESK(L)) 962 1 /(THSK(L)-A4*APESK(L))) 963 ELSE 964 QREFK(L)=Q(I,J,L) 965 ENDIF 966 460 CONTINUE 967 C--------------ENTHALPY CONSERVATION INTEGRAL-------------------------- 968 DO 520 ITER=1,2 969 C----------------------------------------------------------------------- 970 SUMDE=0. 971 SUMDP=0. 972 C 973 DO L=LTPK,LB 974 SUMDE=((TK(L)-TREFK(L))*CP+(QK(L)-QREFK(L))*ELWV)*DETA(L) 975 1 +SUMDE 976 SUMDP=SUMDP+DETA(L) 977 ENDDO 978 C 979 HCORR=SUMDE/(SUMDP-DETA(LTPK)) 980 LCOR=LTPK+1 981 C-----------------------FIND LQM---------------------------------------- 982 DO L=1,LB 983 IF(PK(L).LE.PQM)LQM=L 984 ENDDO 985 C--------------ABOVE LQM CORRECT TEMPERATURE ONLY----------------------- 986 IF(LCOR.LE.LQM)THEN 987 DO L=LCOR,LQM 988 TREFK(L)=TREFK(L)+HCORR*RCP 989 ENDDO 990 LCOR=LQM+1 991 ENDIF Page 13 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 992 C--------------BELOW LQM CORRECT BOTH TEMPERATURE AND MOISTURE---------- 993 C 994 DO 510 L=LCOR,LB 995 TSKL=TREFK(L)*APEK(L)/APESK(L) 996 DHDT=QREFK(L)*A23M4L/(TSKL-A4)**2+CP 997 TREFK(L)=HCORR/DHDT+TREFK(L) 998 THSKL=TREFK(L)*APEK(L) 999 QREFK(L)=PQ0/PSK(L)*EXP(A2*(THSKL-A3*APESK(L)) 1000 1 /(THSKL-A4*APESK(L))) 1001 510 CONTINUE 1002 C----------------------------------------------------------------------- 1003 520 CONTINUE 1004 C--------------HEATING, MOISTENING, PRECIPITATION----------------------- 1005 DENTPY=0. 1006 AVRGT =0. 1007 PRECK =0. 1008 C 1009 DO 530 L=LTPK,LB 1010 TKL =TK(L) 1011 DIFTL =(TREFK(L)-TKL )*TAUK 1012 DIFQL =(QREFK(L)-QK(L))*TAUK 1013 AVRGTL =(TKL+TKL+DIFTL) 1014 DENTPY =(DIFTL*CP+DIFQL*ELWV)*DETA(L)/AVRGTL+DENTPY 1015 AVRGT =AVRGTL*DETA(L)+AVRGT 1016 PRECK =DETA(L)*DIFTL+PRECK 1017 DIFT(L)=DIFTL 1018 DIFQ(L)=DIFQL 1019 530 CONTINUE 1020 C 1021 DENTPY=DENTPY+DENTPY 1022 AVRGT =AVRGT/(SUMDP+SUMDP) 1023 C-------------SWAP IF ENTROPY AND/OR PRECIP .LT. 0 ...------------------ 1024 IF(DENTPY.LT.EPSNTP.OR.PRECK.LT.0.)THEN 1025 IF(OCT90)THEN 1026 CLDEFI(I,J)=EFIMN 1027 ELSE 1028 CLDEFI(I,J)=EFIMN*SM(I,J)+STEFI*(1.-SM(I,J)) 1029 ENDIF 1030 C 1031 C--------------SEARCH FOR SHALLOW CLOUD TOP----------------------------- 1032 LBTK=LBOT(I,J) 1033 LTSH=LBTK 1034 LBM1=LBTK-1 1035 PBTK=PK(LBTK) 1036 c use new threshold for cloud depth 1037 c PTPK=PBTK-PSH 1038 PSFCIJ=PD(I,J)+PT 1039 DEPMIN=PSHNEW*PSFCIJ*1.E-5 1040 PTPK=PBTK-DEPMIN 1041 C-------------CLOUD TOP IS THE LEVEL JUST BELOW PBTK-PSH---------------- 1042 DO L=1,LM 1043 IF(PK(L).LE.PTPK)LTPK=L+1 1044 ENDDO 1045 PTPK=PK(LTPK) 1046 C--------------HIGHEST LEVEL ALLOWED IS LEVEL JUST BELOW PSHU----------- 1047 IF(PTPK.LE.PSHU)THEN 1048 DO L=1,LM Page 14 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1049 IF(PK(L).LE.PSHU)LSHU=L+1 1050 ENDDO 1051 LTPK=LSHU 1052 PTPK=PK(LTPK) 1053 ENDIF 1054 C 1055 LTP1=LTPK+1 1056 LTP2=LTPK+2 1057 C----------------------------------------------------------------------- 1058 DO L=LTPK,LBTK 1059 QSATK(L)=PQ0/PK(L)*EXP(A2*(TK(L)-A3)/(TK(L)-A4)) 1060 ENDDO 1061 C----------------------------------------------------------------------- 1062 RHH=QK(LTPK)/QSATK(LTPK) 1063 RHMAX=0. 1064 C 1065 DO 570 L=LTP1,LBM1 1066 RHL=QK(L)/QSATK(L) 1067 DRHDP=(RHH-RHL)/(PK(L-1)-PK(L)) 1068 IF(DRHDP.GT.RHMAX)THEN 1069 LTSH=L-1 1070 RHMAX=DRHDP 1071 ENDIF 1072 RHH=RHL 1073 570 CONTINUE 1074 C 1075 LTOP(I,J)=LTSH 1076 C---------------CLOUD MUST BE AT LEAST TWO LAYERS THICK----------------- 1077 IF(LBOT(I,J)-LTOP(I,J).LT.2)LTOP(I,J)=LBOT(I,J)-2 1078 C 1079 PTOP(I,J)=PK(LTOP(I,J)) 1080 GO TO 600 1081 ENDIF 1082 C----------------------------------------------------------------------- 1083 C--------------... DEEP CONVECTION OTHERWISE---------------------------- 1084 C----------------------------------------------------------------------- 1085 DRHEAT=(PRECK*SM(I,J)+AMAX1(EPSP,PRECK)*(1.-SM(I,J))) 1086 1 *CP/AVRGT 1087 EFI=EFIFC*DENTPY/DRHEAT 1088 C 1089 C************** UNIFIED OR SEPARATE LAND/SEA CONV. ************** 1090 C 1091 IF(OCT90)THEN 1092 IF(UNIS)THEN 1093 EFI=CLDEFI(I,J)*FCB+EFI*FCC 1094 ELSEIF(.NOT.UNIL)THEN 1095 EFI=(CLDEFI(I,J)*FCB+EFI*FCC)*SM(I,J)+1.-SM(I,J) 1096 ELSE 1097 EFI=1. 1098 ENDIF 1099 ELSE 1100 EFI=CLDEFI(I,J)*FCB+EFI*FCC 1101 ENDIF 1102 C 1103 IF(EFI.GT.1. ) EFI=1. 1104 IF(EFI.LT.EFIMN) EFI=EFIMN 1105 IF(PRECK.EQ.0.) EFI=1. Page 15 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1106 CLDEFI(I,J)=EFI 1107 C 1108 FEFI=EFMNT+SLOPE*(EFI-EFIMN) 1109 C FEFI=AMAX1(EFI,EFMNT) 1110 C 1111 PRECK=PRECK*FEFI 1112 C 1113 C--------------UPDATE PRECIPITATION, TEMPERATURE & MOISTURE------------- 1114 C 1115 PRECOL(I,J)=PDSL(I,J)*PRECK*CPRLG 1116 PREC (I,J)=PDSL(I,J)*PRECK*CPRLG+PREC (I,J) 1117 CUPREC(I,J)=PDSL(I,J)*PRECK*CPRLG+CUPREC(I,J) 1118 ACPREC(I,J)=PDSL(I,J)*PRECK*CPRLG+ACPREC(I,J) 1119 CUPPT(I,J) =PDSL(I,J)*PRECK*CPRLG+CUPPT(I,J) 1120 C 1121 DO L=LTPK,LB 1122 TMOD(I,J,L)=DIFT(L)*FEFI 1123 QMOD(I,J,L)=DIFQ(L)*FEFI 1124 CYL 1125 TLAT(I,J,L)=DIFT(L)*FEFI + TLAT(I,J,L) 1126 CYL 1127 ENDDO 1128 C 1129 CDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCD 1130 CDCDCDCDCDCDC END OF DEEP CONVECTION DCDCDCDCDCDCDCD 1131 CDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCDCD 1132 C 1133 C----------------------------------------------------------------------- 1134 600 CONTINUE 1135 C----------------------------------------------------------------------- 1136 NDEEP=0 1137 C 1138 DO 620 J=MYJS2,MYJE2 1139 DO 620 I=MYIS,MYIE 1140 LTPK=LTOP(I,J) 1141 LBTK=LBOT(I,J) 1142 LB =LMH(I,J)-1 1143 PSFCIJ=PD(I,J)+PT 1144 DEPMIN=PSHNEW*PSFCIJ*1.E-5 1145 IF(PTOP(I,J).LT.PBOT(I,J)-DEPMIN)THEN 1146 NDEEP=NDEEP+1 1147 NDEPTH=LB-LTPK 1148 NTOPD (LTPK )=NTOPD (LTPK )+1 1149 NBOTD (LB )=NBOTD (LB )+1 1150 NDPTHD(NDEPTH)=NDPTHD(NDEPTH)+1 1151 ENDIF 1152 620 CONTINUE 1153 NNEG=KHDEEP-NDEEP 1154 C 1155 C--------------GATHER SHALLOW CONVECTION POINTS------------------------- 1156 C 1157 KHSHAL=0 1158 NDSTN =0 1159 NDSTP =0 1160 C 1161 DO 630 J=MYJS2,MYJE2 1162 DO 630 I=MYIS,MYIE Page 16 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1163 IF(PTOP(I,J).GT.PBOT(I,J)-PNO.OR. 1164 1 LTOP(I,J).GT.LBOT(I,J)-2)GO TO 630 1165 PSFCIJ=PD(I,J)+PT 1166 DEPMIN=PSHNEW*PSFCIJ*1.E-5 1167 IF(PTOP(I,J)+1..GE.PBOT(I,J)-DEPMIN)THEN 1168 KHSHAL=KHSHAL+1 1169 ISHAL(KHSHAL)=I 1170 JSHAL(KHSHAL)=J 1171 ENDIF 1172 630 CONTINUE 1173 C 1174 C************* HORIZONTAL LOOP FOR SHALLOW CONVECTION ****************** 1175 !$omp parallel do 1176 !$omp& private(apek,apekl,apekxx,apekxy,bqk,bqs00k,bqs10k,den,dentpy, 1177 !$omp& dpkl,dpmix,dqref,dst,dstq,dtdeta,fpk,fptk,i,iq,it,j, 1178 !$omp& lbm1,lbtk,ltp1,ltpk,otsum,part1,part2,part3,pk,pkl, 1179 !$omp& pkxxxx,pkxxxy,potsum,ppk,psum,ptpk,pz0,qk,qkl,qnew, 1180 !$omp& qotsum,qqk,qrefk,qrfkl,qrftp,qsatk,qsum,rdpsum,rtbar, 1181 !$omp& smix,sqk,sqs00k,sqs10k,sumdp,sumdt,tcorr,thvmkl, 1182 !$omp& thvref,tk,tkl,tqk,trefk,trefkx,trfkl,tthk) 1183 C 1184 DO 800 N=1,KHSHAL 1185 C*********************************************************************** 1186 I=ISHAL(N) 1187 J=JSHAL(N) 1188 C----------------------------------------------------------------------- 1189 CSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCS 1190 CSCSCSCSCSCSC SHALLOW CONVECTION CSCSCSCSCSCSCSCSCSCSCS 1191 CSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCS 1192 C----------------------------------------------------------------------- 1193 PZ0=PD(I,J)+PT 1194 LLMH=LMH(I,J) 1195 C 1196 DO 650 L=1,LLMH 1197 TKL =T (I,J,L) 1198 TK (L) =TKL 1199 TREFK(L) =TKL 1200 QKL =Q (I,J,L) 1201 QK (L) =QKL 1202 QREFK(L) =QKL 1203 PKL =AETA(L)*PDSL(I,J)+PT 1204 PK (L) =PKL 1205 QSATK(L) =PQ0/PK(L)*EXP(A2*(TK(L)-A3)/(TK(L)-A4)) 1206 APEKL =APE(I,J,L) 1207 CVVVVVVVVVVVV CHOOSE THE PRESSURE FUNCTION VVVVVVVVVVVVVVVVVVVVVVVVVVVVV 1208 C FPK (L) =ALOG(PKL) 1209 C FPK (L) =PKL 1210 C FPK (L) =-1./PKL 1211 CAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAA 1212 APEK (L) =APEKL 1213 THVMKL =TKL*APEKL*(QKL*0.608+1.) 1214 C THVMOD(L)=THVMKL 1215 THVREF(L)=THVMKL 1216 650 CONTINUE 1217 C--------------SHALLOW CLOUD TOP----------------------------- 1218 LBTK=LBOT(I,J) 1219 LBM1=LBTK-1 Page 17 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1220 PTPK=PTOP(I,J) 1221 LTP1=LTOP(I,J) 1222 LTPK=LTOP(I,J)-1 1223 C----------------------------------------------------------------------- 1224 IF(PTOP(I,J).GT.PBOT(I,J)-PNO.OR.LTOP(I,J).GT.LBOT(I,J)-2)THEN 1225 LBOT(I,J)=0 1226 LTOP(I,J)=LBOT(I,J) 1227 PTOP(I,J)=PBOT(I,J) 1228 GO TO 800 1229 ENDIF 1230 C--------------SCALING POTENTIAL TEMPERATURE & TABLE INDEX AT TOP------- 1231 C 1232 C THTPK=T(I,J,LTP1)*APE(I,J,LTP1) 1233 THTPK=T(I,J,LTPK)*APE(I,J,LTPK) 1234 C 1235 TTHK =(THTPK-THL)*RDTH 1236 QQK =TTHK-AINT(TTHK) 1237 IT =INT(TTHK)+1 1238 IF(IT.LT.1)THEN 1239 IT=1 1240 QQK=0. 1241 ENDIF 1242 IF(IT.GE.JTB)THEN 1243 IT=JTB-1 1244 QQK=0. 1245 ENDIF 1246 C--------------BASE AND SCALING FACTOR FOR SPEC. HUMIDITY AT TOP-------- 1247 BQS00K=QS0(IT) 1248 SQS00K=SQS(IT) 1249 BQS10K=QS0(IT+1) 1250 SQS10K=SQS(IT+1) 1251 C--------------SCALING SPEC. HUMIDITY & TABLE INDEX AT TOP-------------- 1252 BQK=(BQS10K-BQS00K)*QQK+BQS00K 1253 SQK=(SQS10K-SQS00K)*QQK+SQS00K 1254 C 1255 C TQK=(Q(I,J,LTP1)-BQK)/SQK*RDQ 1256 TQK=(Q(I,J,LTPK)-BQK)/SQK*RDQ 1257 C 1258 PPK=TQK-AINT(TQK) 1259 IQ =INT(TQK)+1 1260 IF(IQ.LT.1)THEN 1261 IQ=1 1262 PPK=0. 1263 ENDIF 1264 IF(IQ.GE.ITB)THEN 1265 IQ=ITB-1 1266 PPK=0. 1267 ENDIF 1268 C--------------CLOUD TOP SATURATION POINT PRESSURE---------------------- 1269 PART1=(PTBL(IQ+1,IT)-PTBL(IQ,IT))*PPK 1270 PART2=(PTBL(IQ,IT+1)-PTBL(IQ,IT))*QQK 1271 PART3=(PTBL(IQ ,IT )-PTBL(IQ+1,IT ) 1272 1 -PTBL(IQ ,IT+1)+PTBL(IQ+1,IT+1))*PPK*QQK 1273 PTPK=PTBL(IQ,IT)+PART1+PART2+PART3 1274 C----------------------------------------------------------------------- 1275 DPMIX=PTPK-PSP(I,J) 1276 IF(ABS(DPMIX).LT.3000.)DPMIX=-3000. Page 18 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1277 C--------------TEMPERATURE PROFILE SLOPE-------------------------------- 1278 SMIX=(THTPK-THBT(I,J))/DPMIX*STABS 1279 C 1280 TREFKX=TREFK(LBTK+1) 1281 PKXXXX=PK(LBTK+1) 1282 PKXXXY=PK(LBTK) 1283 APEKXX=APEK(LBTK+1) 1284 APEKXY=APEK(LBTK) 1285 C 1286 DO 670 L=LBTK,LTP1,-1 1287 TREFKX=((PKXXXY-PKXXXX)*SMIX 1288 1 +TREFKX*APEKXX)/APEKXY 1289 TREFK(L)=TREFKX 1290 APEKXX=APEKXY 1291 PKXXXX=PKXXXY 1292 APEKXY=APEK(L-1) 1293 PKXXXY=PK(L-1) 1294 670 CONTINUE 1295 C--------------TEMPERATURE REFERENCE PROFILE CORRECTION----------------- 1296 SUMDT=0. 1297 SUMDP=0. 1298 C 1299 DO L=LTP1,LBTK 1300 SUMDT=(TK(L)-TREFK(L))*DETA(L)+SUMDT 1301 SUMDP=SUMDP+DETA(L) 1302 ENDDO 1303 C 1304 RDPSUM=1./SUMDP 1305 FPK(LBTK)=TREFK(LBTK) 1306 C 1307 TCORR=SUMDT*RDPSUM 1308 C 1309 DO L=LTP1,LBTK 1310 C TCORR=SUMDT/(SUMDP-DETA(LBTK)) 1311 C DO L=LTP1,LBM1 1312 C 1313 TRFKL =TREFK(L)+TCORR 1314 TREFK(L)=TRFKL 1315 FPK (L)=TRFKL 1316 ENDDO 1317 C--------------HUMIDITY PROFILE EQUATIONS------------------------------- 1318 PSUM =0. 1319 QSUM =0. 1320 POTSUM=0. 1321 QOTSUM=0. 1322 OTSUM =0. 1323 DST =0. 1324 FPTK =FPK(LTP1) 1325 C 1326 DO 700 L=LTP1,LBTK 1327 DPKL =FPK(L)-FPTK 1328 PSUM =DPKL *DETA(L)+PSUM 1329 QSUM =QK(L)*DETA(L)+QSUM 1330 RTBAR =2./(TREFK(L)+TK(L)) 1331 OTSUM =DETA(L)*RTBAR+OTSUM 1332 POTSUM=DPKL *RTBAR*DETA(L)+POTSUM 1333 QOTSUM=QK(L) *RTBAR*DETA(L)+QOTSUM Page 19 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1334 DST =(TREFK(L)-TK(L))*RTBAR*DETA(L)+DST 1335 700 CONTINUE 1336 C 1337 PSUM =PSUM*RDPSUM 1338 QSUM =QSUM*RDPSUM 1339 ROTSUM=1./OTSUM 1340 POTSUM=POTSUM*ROTSUM 1341 QOTSUM=QOTSUM*ROTSUM 1342 DST =DST *ROTSUM*CP/ELWV 1343 C--------------ENSURE POSITIVE ENTROPY CHANGE--------------------------- 1344 IF(DST.GT.0.)THEN 1345 C 1346 C DSTQ=DST*EPSUP 1347 LBOT(I,J)=0 1348 LTOP(I,J)=LBOT(I,J) 1349 PTOP(I,J)=PBOT(I,J) 1350 cccc NDSTP=NDSTP+1 1351 GO TO 800 1352 C 1353 ELSE 1354 DSTQ=DST*EPSDN 1355 ENDIF 1356 C--------------CHECK FOR ISOTHERMAL ATMOSPHERE-------------------------- 1357 DEN=POTSUM-PSUM 1358 C 1359 IF(-DEN/PSUM.LT.5.E-5)THEN 1360 LBOT(I,J)=0 1361 LTOP(I,J)=LBOT(I,J) 1362 PTOP(I,J)=PBOT(I,J) 1363 GO TO 800 1364 C 1365 C--------------SLOPE OF THE REFERENCE HUMIDITY PROFILE------------------ 1366 ELSE 1367 DQREF=(QOTSUM-DSTQ-QSUM)/DEN 1368 ENDIF 1369 C------------ HUMIDITY DOES NOT INCREASE WITH HEIGHT-------------------- 1370 IF(DQREF.LT.0.)THEN 1371 LBOT(I,J)=0 1372 LTOP(I,J)=LBOT(I,J) 1373 PTOP(I,J)=PBOT(I,J) 1374 GO TO 800 1375 ENDIF 1376 C--------------HUMIDITY AT THE CLOUD TOP-------------------------------- 1377 QRFTP=QSUM-DQREF*PSUM 1378 C--------------HUMIDITY PROFILE----------------------------------------- 1379 C 1380 DO 720 L=LTP1,LBTK 1381 QRFKL=(FPK(L)-FPTK)*DQREF+QRFTP 1382 C 1383 C--------------SUPERSATURATION OR NEGATIVE Q NOT ALLOWED---------------- 1384 C 1385 QNEW =(QRFKL-QK(L))*TAUK+QK(L) 1386 IF(QNEW.GT.QSATK(L)*STRESH.OR.QNEW.LT.0.)THEN 1387 LBOT(I,J)=0 1388 LTOP(I,J)=LBOT(I,J) 1389 PTOP(I,J)=PBOT(I,J) 1390 GO TO 800 Page 20 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1391 ENDIF 1392 C 1393 C----------------------------------------------------------------------- 1394 THVREF(L)=TREFK(L)*APEK(L)*(QRFKL*0.608+1.) 1395 QREFK(L)=QRFKL 1396 720 CONTINUE 1397 C--------------ELIMINATE IMPOSSIBLE SLOPES (BETTS, DTHETA/DQ)----------- 1398 DO 730 L=LTP1,LBTK 1399 DTDETA=(THVREF(L-1)-THVREF(L))/(AETA(L)-AETA(L-1)) 1400 IF(DTDETA.LT.EPSTH)THEN 1401 LBOT(I,J)=0 1402 LTOP(I,J)=LBOT(I,J) 1403 PTOP(I,J)=PBOT(I,J) 1404 GO TO 800 1405 ENDIF 1406 730 CONTINUE 1407 C*************************** DIAGNOSTICS ****************************** 1408 cccc IF(DST.GT.0.)THEN 1409 cccc NDSTP=NDSTP+1 1410 cccc ELSE 1411 cccc NDSTN=NDSTN+1 1412 cccc ENDIF 1413 C********************************************************************** 1414 DENTPY=0. 1415 C 1416 DO L=LTP1,LBTK 1417 DENTPY=((TREFK(L)-TK(L))*CP+(QREFK(L)-QK(L))*ELWV) 1418 1 /(TK(L)+TREFK(L))*DETA(L)+DENTPY 1419 ENDDO 1420 C 1421 C----------------------------------------------------------------------- 1422 C--------------RELAXATION TOWARDS REFERENCE PROFILES-------------------- 1423 C----------------------------------------------------------------------- 1424 C 1425 C 1426 DO 750 L=LTP1,LBTK 1427 TMOD(I,J,L)=(TREFK(L)-TK(L))*TAUK 1428 QMOD(I,J,L)=(QREFK(L)-QK(L))*TAUK 1429 750 CONTINUE 1430 C----------------------------------------------------------------------- 1431 CSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCS 1432 CSCSCSCSCSCSC END OF SHALLOW CONVECTION SCSCSCSCSCSCSCSCS 1433 CSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCSCS 1434 C----------------------------------------------------------------------- 1435 800 CONTINUE 1436 C----------------------------------------------------------------------- 1437 C 1438 C************************** DIAGNOSTICS ******************************** 1439 NSHAL=0 1440 C 1441 DO 820 J=MYJS2,MYJE2 1442 DO 820 I=MYIS,MYIE 1443 LTPK=LTOP(I,J) 1444 LBTK=LBOT(I,J) 1445 PTPK=PTOP(I,J) 1446 PBTK=PBOT(I,J) 1447 IF(PTPK.GT.PBTK-PNO.OR.LTPK.GT.LBTK-2)GO TO 820 Page 21 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1448 PSFCIJ=PD(I,J)+PT 1449 DEPMIN=PSHNEW*PSFCIJ*1.E-5 1450 C IF(PTPK.GE.PBTK-PSH)THEN 1451 IF(PTPK.GE.PBTK-DEPMIN)THEN 1452 NSHAL=NSHAL+1 1453 NTOPS(LTPK)=NTOPS(LTPK)+1 1454 NBOTS(LBTK)=NBOTS(LBTK)+1 1455 NDEPTH=LBTK-LTPK 1456 NDPTHS(NDEPTH)=NDPTHS(NDEPTH)+1 1457 ENDIF 1458 820 CONTINUE 1459 NEGDS=KHSHAL-NSHAL 1460 C*********************************************************************** 1461 C 1462 C--------------SMOOTHING TEMPERATURE & HUMIDITY CORRECTIONS------------- 1463 IF(KSMUD.EQ.0)THEN 1464 !$omp parallel do 1465 C 1466 DO 900 L=1,LM 1467 C*** 1468 C*** UPDATE THE FUNDAMENTAL PROGNOSTIC ARRAYS 1469 C*** 1470 DO 830 J=MYJS,MYJE 1471 DO 830 I=MYIS,MYIE 1472 T(I,J,L)=T(I,J,L)+TMOD(I,J,L) 1473 Q(I,J,L)=Q(I,J,L)+QMOD(I,J,L) 1474 C*** 1475 C*** ACCUMULATE LATENT HEATING DUE TO CONVECTION. 1476 C*** SCALE BY THE RECIPROCAL OF THE PERIOD AT WHICH THIS ROUTINE 1477 C*** IS CALLED. THIS PERIOD IS THE CONVECTION TIMESTEP. 1478 C*** 1479 TCUCN(I,J,L)=TCUCN(I,J,L)+TMOD(I,J,L)*RDTCNVC 1480 830 CONTINUE 1481 900 CONTINUE 1482 ELSE 1483 !$omp parallel do private(ql,qne,qse,tl,tne,tse) 1484 C 1485 DO 910 L=1,LM 1486 C 1487 CALL ZERO2(QL) 1488 CALL ZERO2(QNE) 1489 CALL ZERO2(QSE) 1490 CALL ZERO2(TL) 1491 CALL ZERO2(TNE) 1492 CALL ZERO2(TSE) 1493 C----------------------------------------------------------------------- 1494 DO J=MYJS,MYJE 1495 DO I=MYIS,MYIE 1496 TL(I,J)=TMOD(I,J,L) 1497 QL(I,J)=QMOD(I,J,L) 1498 ENDDO 1499 ENDDO 1500 C----------------------------------------------------------------------- 1501 C----------------------------------------------------------------------- 1502 NSMUD=KSMUD 1503 C 1504 DO 870 KS=1,NSMUD Page 22 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1505 C 1506 DO J=MYJS,MYJE1 1507 DO I=MYIS,MYIE 1508 TNE(I,J)=(TL(I+IHE(J),J+1)-TL(I,J)) 1509 1 *HTM(I,J,L)*HTM(I+IHE(J),J+1,L) 1510 QNE(I,J)=(QL(I+IHE(J),J+1)-QL(I,J)) 1511 1 *HTM(I,J,L)*HTM(I+IHE(J),J+1,L) 1512 ENDDO 1513 ENDDO 1514 C 1515 DO J=MYJS1,MYJE 1516 DO I=MYIS,MYIE 1517 TSE(I,J)=(TL(I+IHE(J),J-1)-TL(I,J)) 1518 1 *HTM(I+IHE(J),J-1,L)*HTM(I,J,L) 1519 QSE(I,J)=(QL(I+IHE(J),J-1)-QL(I,J)) 1520 1 *HTM(I+IHE(J),J-1,L)*HTM(I,J,L) 1521 ENDDO 1522 ENDDO 1523 C 1524 DO J=MYJS2,MYJE2 1525 DO I=MYIS,MYIE 1526 TL(I,J)=(TNE(I,J)-TNE(I+IHW(J),J-1)+TSE(I,J)-TSE(I+IHW(J),J+1)) 1527 1 *HBM2(I,J)*0.125+TL(I,J) 1528 QL(I,J)=(QNE(I,J)-QNE(I+IHW(J),J-1)+QSE(I,J)-QSE(I+IHW(J),J+1)) 1529 1 *HBM2(I,J)*0.125+QL(I,J) 1530 ENDDO 1531 ENDDO 1532 C 1533 870 CONTINUE 1534 C----------------------------------------------------------------------- 1535 C*** 1536 C*** UPDATE THE FUNDAMENTAL PROGNOSTIC ARRAYS 1537 C*** 1538 DO J=MYJS,MYJE 1539 DO I=MYIS,MYIE 1540 T(I,J,L)=T(I,J,L)+TL(I,J) 1541 Q(I,J,L)=Q(I,J,L)+QL(I,J) 1542 ENDDO 1543 ENDDO 1544 C*** 1545 C*** ACCUMULATE LATENT HEATING DUE TO CONVECTION. 1546 C*** SCALE BY THE RECIPROCAL OF THE PERIOD AT WHICH THIS ROUTINE 1547 C*** IS CALLED. THIS PERIOD IS THE CONVECTION TIMESTEP. 1548 C*** 1549 DO J=MYJS,MYJE 1550 DO I=MYIS,MYIE 1551 TCUCN(I,J,L)=TCUCN(I,J,L)+TL(I,J)*RDTCNVC 1552 ENDDO 1553 ENDDO 1554 C 1555 910 CONTINUE 1556 C 1557 ENDIF 1558 C--------------SAVE CLOUD TOP AND BOTTOM FOR RADIATION------------------ 1559 !$omp parallel do 1560 DO J=MYJS,MYJE 1561 DO I=MYIS,MYIE Page 23 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F 1562 IF (LTOP(I,J).GT.0 .AND. LTOP(I,J).LT.LBOT(I,J)) THEN 1563 CUtop=FLOAT(LTOP(I,J)) 1564 CUbot=FLOAT(LBOT(I,J)) 1565 HTOP(I,J)=MIN(CUtop,HTOP(I,J)) 1566 HBOT(I,J)=MAX(CUbot,HBOT(I,J)) 1567 CNVTOP(I,J)=MIN(CUtop,CNVTOP(I,J)) 1568 CNVBOT(I,J)=MAX(CUbot,CNVBOT(I,J)) 1569 ENDIF 1570 ENDDO 1571 ENDDO 1572 C----------------------------------------------------------------------- 1573 C************************* DIAGNOSTICS ********************************* 1574 C 1575 C WRITE(LIST,950)NTSD,NSHAL,NDEEP,NNEG,NEGDS,NDSTN,NDSTP 1576 C DO 940 L=1,LM 1577 C WRITE(LIST,952)L 1578 C WRITE(LIST,954)NBOTS(L),NTOPS(L),NDPTHS(L) 1579 C 1, NBOTD(L),NTOPD(L),NDPTHD(L) 1580 C 940 CONTINUE 1581 C 950 FORMAT(' NTSD=',I3,I8,' SHALLOW, ',I4,' DEEP, ', 1582 C 1 I4,' NEG., ',I4,' NEG. SHALL.,',I4,' DST.LT.0, ',I4,' DST.GT.0') 1583 C 952 FORMAT(' LAYER (FROM TOP),',I2) 1584 C 954 FORMAT(' NBOTS=',I4,' NTOPS=',I4,' NDPTHS=',I4, 1585 C 1 ' NBOTD=',I4,' NTOPD=',I4,' NDPTHD=',I4) 1586 C*********************************************************************** 1587 RETURN 1588 END ENTRY POINTS Name cucnvc_ Page 24 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 120 Label 193 186,187 130 Label 207 199 150 Label 285 220 155 Label 286 219 170 Label 288 211 240 Label 303 292 250 Label 323 309,310 270 Label 355 343,344 290 Label 378 337 295 Label 388 380 345 Label 466 452,453 350 Label 480 450 360 Label 492 485 380 Label 509 495,496 410 Label 584 568 420 Label 611 599 430 Label 616 600 450 Label 628 614 460 Label 658 630 510 Label 693 686 520 Label 695 660 530 Label 711 701 570 Label 765 757 600 Label 826 539,772 620 Label 844 830,831 630 Label 864 853,854,856 650 Label 908 888 670 Label 986 978 700 Label 1027 1018 720 Label 1088 1072 730 Label 1098 1090 750 Label 1121 1118 800 Label 1127 876,920,1043,1055,1066,1082,1096 820 Label 1150 1133,1134,1139 830 Label 1172 1162,1163 870 Label 1225 1196 900 Label 1173 1158 910 Label 1247 1177 A2 Param 7 R(4) 4 scalar 88,653,691,751,897 A23M4L Param 88 R(4) 4 scalar 688 A3 Param 7 R(4) 4 scalar 88,653,691,751,897 A4 Param 7 R(4) 4 scalar 88,654,688,692,751,897 ABS Func 968 scalar 968 ACCLIQ Scalar 13 R(4) 4 2 17955 COM ACMCLH Common 2 6463824 ACPREC Scalar 12 R(4) 4 2 17955 COM 810 ACUTIM Scalar 3 R(4) 4 scalar COM 165 AETA Scalar 9 R(4) 4 1 45 COM 202,222,224,293,312,313,317,321,34 5,361,372,393,577,895,1091 AETAL Local 293 R(4) 4 scalar 293,298 AFSI Scalar 26 R(4) 4 2 17955 COM Page 25 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References AINT Func 233 scalar 233,255,928,950 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 777 scalar 777 APE Local 147 R(4) 4 3 807975 203,231,580,898,925 APEBT Local 129 R(4) 4 2 17955 APEK Local 118 R(4) 4 1 45 581,595,597,607,618,621,652,687,69 0,904,975,976,984,1086 APEKL Local 580 R(4) 4 scalar 580,581,582,898,904,905 APEKXX Local 595 R(4) 4 scalar 595,603,605,975,980,982 APEKXY Local 597 R(4) 4 scalar 597,603,605,607,976,980,982,984 APESK Local 119 R(4) 4 1 45 651,653,654,687,691,692 APESP Local 276 R(4) 4 scalar 276,277 APESTS Local 202 R(4) 4 scalar 202,203 APREC Scalar 3 R(4) 4 2 17955 COM ARATIM Scalar 3 R(4) 4 scalar COM AVCNVC Scalar 3 R(4) 4 scalar COM 164 AVGEFI Param 75 R(4) 4 scalar 508 AVRAIN Scalar 3 R(4) 4 scalar COM AVRGT Local 698 R(4) 4 scalar 698,707,714,778 AVRGTL Local 705 R(4) 4 scalar 705,706,707 BQ Local 252 R(4) 4 scalar 252,254 BQK Local 944 R(4) 4 scalar 944,948 BQS00K Local 247 R(4) 4 scalar 247,252,939,944 BQS10K Local 249 R(4) 4 scalar 249,252,941,944 BTH Local 134 R(4) 4 2 17955 CAPA Param 5 R(4) 4 scalar 203,276,651 CFRACH Scalar 4 R(4) 4 2 17955 COM CFRACL Scalar 3 R(4) 4 2 17955 COM CFRACM Scalar 4 R(4) 4 2 17955 COM CI Scalar 6 R(4) 4 scalar COM CLDEFI Scalar 11 R(4) 4 2 17955 COM 405,415,428,463,464,476,508,550,71 8,720,785,787,792,798 CLDINC Scalar 6 R(4) 4 2 17955 COM CLDINC7 Scalar 6 R(4) 4 2 17955 COM CNVBOT Scalar 18 R(4) 4 2 17955 COM 1260 CNVCLD Common 2 287280 CNVTOP Scalar 18 R(4) 4 2 17955 COM 1259 CP Param 5 R(4) 4 scalar 89,666,688,706,778,1034,1109 CPRLG Param 89 R(4) 4 scalar 807,808,809,810,811 CS Scalar 6 R(4) 4 scalar COM CTLBLK Common 4 112 CUBOT Local 1256 R(4) 4 scalar 1256,1258,1260 CUCNVC Subr 1 CUPPT Scalar 3 R(4) 4 2 17955 COM 811 CUPREC Scalar 13 R(4) 4 2 17955 COM 809 CUTOP Local 1255 R(4) 4 scalar 1255,1257,1259 CZEN Scalar 16 R(4) 4 2 17955 COM CZMEAN Scalar 23 R(4) 4 2 17955 COM DEN Local 1049 R(4) 4 scalar 1049,1051,1059 DENTPY Local 697 R(4) 4 scalar 697,706,713,716,779,1106,1109,1110 DEPMIN Local 518 R(4) 4 scalar 518,520,731,732,836,837,858,859,11 41,1143 Page 26 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References DEPTH Local 519 R(4) 4 scalar 519,520,628,635 DEPWL Local 629 R(4) 4 scalar 629,635 DETA Scalar 9 R(4) 4 1 45 COM 666,668,671,706,707,708,992,993,10 20,1021,1023,1024,1025,1026,1110 DFRLG Scalar 9 R(4) 4 1 46 COM DHDT Local 688 R(4) 4 scalar 688,689 DI Scalar 6 R(4) 4 scalar COM DIFQ Local 120 R(4) 4 1 45 570,710,815 DIFQL Local 704 R(4) 4 scalar 704,706,710 DIFT Local 120 R(4) 4 1 45 569,709,814,817 DIFTL Local 703 R(4) 4 scalar 703,705,706,708,709 DPKL Local 1019 R(4) 4 scalar 1019,1020,1024 DPMIX Local 967 R(4) 4 scalar 967,968,970 DQREF Local 1059 R(4) 4 scalar 1059,1062,1069,1073 DRHDP Local 759 R(4) 4 scalar 759,760,762 DRHEAT Local 777 R(4) 4 scalar 777,779 DS Scalar 6 R(4) 4 scalar COM DSP Local 637 R(4) 4 scalar 637,639,643,645,650 DSP0 Local 153 R(4) 4 2 17955 159,407,418,430,457,458,474,552 DSP0FL Param 56 R(4) 4 scalar 61,81 DSP0FS Param 60 R(4) 4 scalar 62,84 DSP0K Local 552 R(4) 4 scalar 552,637,639,643,645 DSP0SL Param 61 R(4) 4 scalar 81,419,430 DSP0SS Param 62 R(4) 4 scalar 84,407,418 DSPB Local 153 R(4) 4 2 17955 160,406,416,429,454,455,473,551 DSPBFL Param 56 R(4) 4 scalar 61,80 DSPBFS Param 60 R(4) 4 scalar 62,83 DSPBK Local 551 R(4) 4 scalar 551,639 DSPBSL Param 61 R(4) 4 scalar 80,417,429 DSPBSS Param 62 R(4) 4 scalar 83,406,416 DSPC Param 75 R(4) 4 scalar DSPT Local 154 R(4) 4 2 17955 161,408,420,431,460,461,475,553 DSPTFL Param 56 R(4) 4 scalar 61,82 DSPTFS Param 60 R(4) 4 scalar 62,85 DSPTK Local 553 R(4) 4 scalar 553,637,645 DSPTSL Param 61 R(4) 4 scalar 82,421,431 DSPTSS Param 62 R(4) 4 scalar 85,408,420 DST Local 1015 R(4) 4 scalar 1015,1026,1034,1036,1046 DSTQ Local 1046 R(4) 4 scalar 1046,1059 DT Scalar 6 R(4) 4 scalar COM 166 DTCNVC Local 166 R(4) 4 scalar 166,167,168 DTD Scalar 5 R(4) 4 scalar COM DTDETA Local 1091 R(4) 4 scalar 1091,1092 DTHEM Local 618 R(4) 4 scalar 618,621 DTQ2 Scalar 5 R(4) 4 scalar COM DTTOP Param 9 R(4) 4 scalar 384 EFI Local 405 R(4) 4 scalar 405,406,407,408,415,416,417,418,41 9,420,421,428,429,430,431,550,779, 785,787,789,792,795,796,797,798,80 0 EFIFC Param 74 R(4) 4 scalar 779 EFIMN Param 46 R(4) 4 scalar 75,80,81,82,83,84,85,86,406,407,40 8,416,417,418,419,420,421,429,430, 431,718,720,796,800 EFMNT Param 46 R(4) 4 scalar 86,800 Page 27 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References ELIVW Param 6 R(4) 4 scalar 89 ELOCP Param 89 R(4) 4 scalar 277 ELWV Param 6 R(4) 4 scalar 88,89,666,706,1034,1109 EPSDN Param 11 R(4) 4 scalar 1046 EPSNTP Param 74 R(4) 4 scalar 716 EPSP Param 75 R(4) 4 scalar 777 EPSQ Param 6 R(4) 4 scalar 204 EPSR Scalar 13 R(4) 4 2 17955 COM EPSTH Param 11 R(4) 4 scalar 1092 EPSUP Param 11 R(4) 4 scalar EXP Func 277 scalar 277,653,691,751,897 FCB Param 46 R(4) 4 scalar 785,787,792 FCC Param 46 R(4) 4 scalar 46,785,787,792 FEFI Local 800 R(4) 4 scalar 800,803,814,815,817 FIRST Scalar 95 L(4) 4 scalar COM FIS Scalar 7 R(4) 4 2 17955 COM FLOAT Func 1255 scalar 1255,1256 FPK Local 121 R(4) 4 1 45 997,1007,1016,1019,1073 FPTK Local 1016 R(4) 4 scalar 1016,1019,1073 FSL Param 56 R(4) 4 scalar 61 FSS Param 46 R(4) 4 scalar 62 G Param 5 R(4) 4 scalar 89 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 456,459,462,465,502,1219,1221 HBM3 Scalar 8 R(4) 4 2 17955 COM HBOT Scalar 17 R(4) 4 2 17955 COM 1258 HCORR Local 671 R(4) 4 scalar 671,680,689 HDAC Scalar 22 R(4) 4 2 17955 COM HDACV Scalar 22 R(4) 4 2 17955 COM HTM Scalar 6 R(4) 4 3 807975 COM 204,361,372,1201,1203,1210,1212 HTMG Scalar 52 R(4) 4 3 4183695 COM HTOP Scalar 17 R(4) 4 2 17955 COM 1257 I Local 187 I(4) 4 scalar 187,188,189,190,191,192,201,202,20 3,204,220,222,223,224,230,231,279, 280,281,282,297,298,299,310,311,31 2,313,314,317,318,321,327,328,329, 344,345,348,352,383,384,392,393,40 4,405,406,407,408,414,415,416,417, 418,419,420,421,427,428,429,430,43 1,440,441,442,443,444,453,454,455, 456,457,458,459,460,461,462,463,46 4,465,470,471,473,474,475,476,487, 488,489,496,497,498,499,500,502,50 3,504,505,507,508,514,515,519,523, 541,543,544,545,550,551,552,553,57 1,574,577,580,582,600,656,718,720, 724,730,767,769,771,777,785,787,79 2,798,807,808,809,810,811,814,815, 817,831,832,833,834,835,837,854,85 5,856,857,859,861,878,885,886,889, Page 28 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References 892,895,898,910,912,913,914,916,91 7,918,919,925,948,967,970,1039,104 0,1041,1052,1053,1054,1063,1064,10 65,1079,1080,1081,1093,1094,1095,1 119,1120,1134,1135,1136,1137,1138, 1140,1163,1164,1165,1171,1187,1188 ,1189,1199,1200,1201,1202,1203,120 8,1209,1210,1211,1212,1217,1218,12 19,1220,1221,1231,1232,1233,1242,1 243,1253,1254,1255,1256,1257,1258, 1259,1260 IBROW Scalar 22 I(4) 4 scalar COM IBUOY Local 141 I(4) 4 1 15327 ICHUNKTAB Scalar 29 I(4) 4 1 8 COM IDAT Scalar 5 I(4) 4 1 3 COM IDEEP Local 142 I(4) 4 1 15327 523,541 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,126,127 ,128,129,130,131,132,133,134,135,1 36,137,138,139,147,148,149,150,153 ,154,155,156,157,360,361,362,364,3 65,371,372,373,375,376 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,84,126, 127,128,129,130,131,132,133,134,13 5,136,137,138,139,147,148,149,150, 153,154,155,156,157 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 IHE Scalar 3 I(4) 4 1 399 COM 454,455,457,458,460,461,463,464,12 00,1201,1202,1203,1209,1210,1211,1 212 IHEG Scalar 9 I(4) 4 1 387 COM IHHA Scalar 3 I(4) 4 1 760 COM IHLA Scalar 3 I(4) 4 1 760 COM IHRES Local 145 I(4) 4 1 15327 352,372 IHRST Scalar 5 I(4) 4 scalar COM IHW Scalar 3 I(4) 4 1 399 COM 454,455,457,458,460,461,463,464,12 18,1220 IHWG Scalar 9 I(4) 4 1 387 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 ILRES Local 144 I(4) 4 1 15327 348,361 IM Param 3 I(4) 4 scalar 9,16,19,43,44,45,46,47,51,52,57,83 IMJM Param 83 I(4) 4 scalar IMJM_LOC Param 84 I(4) 4 scalar 141,142,143,144,145 INDX Common 2 6564 INDXG Common 8 8084 Page 29 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References INPES Param 11 I(4) 4 scalar 16,19,27,28,29,38,39 INT Func 234 scalar 234,256,929,951 INUMQ Scalar 34 I(4) 4 1 100 COM IOUT Scalar 5 I(4) 4 scalar COM IPTB Local 128 I(4) 4 2 17955 364,375 IQ Local 267 I(4) 4 scalar 267,269,270,271,272,951,952,953,95 6,957,961,962,963,964,965 IQTB Local 256 I(4) 4 scalar 256,258,259,262,263,267 IQUILT_GROUP Scalar 33 I(4) 4 scalar COM IRAD Scalar 5 I(4) 4 1 45 COM IRADG Scalar 9 I(4) 4 1 473 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 ISHAL Local 143 I(4) 4 1 15327 861,878 IS_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM IT Local 268 I(4) 4 scalar 268,269,270,271,272,929,930,931,93 4,935,939,940,941,942,961,962,963, 964,965 ITAIL Param 16 I(4) 4 scalar ITB Param 1 I(4) 4 scalar 10,25,262,263,360,956,957 ITBQ Param 1 I(4) 4 scalar 11,25,371 ITEMP Scalar 47 I(4) 4 2 91719 COM ITEMP2 Scalar 47 I(4) 4 2 91719 COM ITER Local 660 I(4) 4 scalar 660 ITHTB Local 128 I(4) 4 2 17955 365,376 ITROW Scalar 22 I(4) 4 scalar COM ITSTLOC Scalar 4 I(4) 4 scalar COM ITTB Local 234 I(4) 4 scalar 234,236,237,240,241,246,268 ITTBK Local 246 I(4) 4 scalar 246,247,248,249,250 IVE Scalar 4 I(4) 4 1 399 COM IVEG Scalar 9 I(4) 4 1 387 COM IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IVW Scalar 4 I(4) 4 1 399 COM IVWG Scalar 9 I(4) 4 1 387 COM J Local 186 I(4) 4 scalar 186,188,189,190,191,192,200,202,20 3,204,219,222,223,224,230,231,279, 280,281,282,296,298,299,309,311,31 2,313,314,317,318,321,326,328,329, 343,345,349,353,382,384,391,393,40 3,405,406,407,408,413,415,416,417, 418,419,420,421,426,428,429,430,43 1,439,441,442,443,444,452,454,455, 456,457,458,459,460,461,462,463,46 4,465,469,471,473,474,475,476,486, 488,489,495,497,498,499,500,502,50 3,504,505,507,508,513,515,519,524, 542,543,544,545,550,551,552,553,57 1,574,577,580,582,600,656,718,720, 724,730,767,769,771,777,785,787,79 Page 30 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References 2,798,807,808,809,810,811,814,815, 817,830,832,833,834,835,837,853,85 5,856,857,859,862,879,885,886,889, 892,895,898,910,912,913,914,916,91 7,918,919,925,948,967,970,1039,104 0,1041,1052,1053,1054,1063,1064,10 65,1079,1080,1081,1093,1094,1095,1 119,1120,1133,1135,1136,1137,1138, 1140,1162,1164,1165,1171,1186,1188 ,1189,1198,1200,1201,1202,1203,120 7,1209,1210,1211,1212,1216,1218,12 19,1220,1221,1230,1232,1233,1241,1 243,1252,1254,1255,1256,1257,1258, 1259,1260 JAM Param 83 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 JBUOY Local 141 I(4) 4 1 15327 JDEEP Local 142 I(4) 4 1 15327 524,542 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,126,127 ,128,129,130,131,132,133,134,135,1 36,137,138,139,147,148,149,150,153 ,154,155,156,157,360,361,362,364,3 65,371,372,373,375,376 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,84,126, 127,128,129,130,131,132,133,134,13 5,136,137,138,139,147,148,149,150, 153,154,155,156,157 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 JHRES Local 145 I(4) 4 1 15327 353,372 JLRES Local 144 I(4) 4 1 15327 349,361 JM Param 3 I(4) 4 scalar 9,17,20,43,44,45,46,47,51,52,57,83 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 JRA Scalar 3 I(4) 4 1 760 COM JSHAL Local 143 I(4) 4 1 15327 862,879 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,240,241,360,934,935 JTBQ Param 1 I(4) 4 scalar 25,371 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 KB Local 211 I(4) 4 scalar 211,222,229,230,231 Page 31 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References KHDEEP Local 511 I(4) 4 scalar 511,522,523,524,539,845 KHSHAL Local 849 I(4) 4 scalar 849,860,861,862,876,1151 KNUMH Local 341 I(4) 4 scalar 341,351,352,353,370,371 KNUML Local 340 I(4) 4 scalar 340,347,348,349,359,360 KOUNT Local 450 I(4) 4 scalar 450 KROW Local 448 I(4) 4 scalar 448,450 KS Local 1196 I(4) 4 scalar 1196 KSMUD Param 79 I(4) 4 scalar 1155,1194 KTM Scalar 5 I(4) 4 scalar COM L Local 174 I(4) 4 scalar 174,175,176,177,178,179,180,199,20 2,203,204,292,293,299,316,317,318, 337,345,360,361,371,372,380,384,48 5,488,489,568,569,570,571,572,573, 574,575,576,577,578,579,580,581,58 2,599,600,604,607,608,609,620,621, 630,636,637,639,644,645,649,650,65 1,652,653,654,656,665,666,668,674, 675,679,680,686,687,688,689,690,69 1,692,701,702,703,704,706,707,708, 709,710,734,735,740,741,750,751,75 7,758,759,761,813,814,815,817,888, 889,890,891,892,893,894,895,896,89 7,898,904,907,978,981,984,985,991, 992,993,1001,1005,1006,1007,1018,1 019,1020,1021,1022,1023,1024,1025, 1026,1072,1073,1077,1078,1086,1087 ,1090,1091,1108,1109,1110,1118,111 9,1120,1158,1164,1165,1171,1177,11 88,1189,1201,1203,1210,1212,1232,1 233,1243 L0 Local 591 I(4) 4 scalar 591,609,610,613,616,618,636,644 L0M1 Local 613 I(4) 4 scalar 613,616,620 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LB Local 549 I(4) 4 scalar 549,587,588,591,592,593,594,595,63 0,665,674,686,701,813,834,839,841 LBM1 Local 587 I(4) 4 scalar 587,596,597,599,726,757,911 LBOT Local 126 I(4) 4 2 17955 188,299,312,314,318,321,328,497,49 8,499,503,504,507,545,724,769,833, 856,910,916,917,918,1039,1040,1052 ,1053,1063,1064,1079,1080,1093,109 4,1136,1254,1256 LBTK Local 545 I(4) 4 scalar 545,549,724,725,726,727,750,833,91 0,911,972,973,974,975,976,978,991, 997,1001,1018,1072,1090,1108,1118, 1136,1139,1146,1147 LCOR Local 672 I(4) 4 scalar 672,678,679,682,686 LIST Scalar 5 I(4) 4 scalar COM LLMH Local 886 I(4) 4 scalar 886,888 LM Param 3 I(4) 4 scalar 2,4,5,6,7,9,21,52,91,118,119,120,1 21,122,123,147,148,149,150,174,199 ,211,337,380,485,568,734,740,1158, 1177 LM1 Param 91 I(4) 4 scalar 292 LMH Scalar 4 I(4) 4 2 17955 COM 188,223,311,834,886 Page 32 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References LMHIJ Local 311 I(4) 4 scalar 311,313,314,316 LMHK Local 223 I(4) 4 scalar 223,224,229 LMV Scalar 4 I(4) 4 2 17955 COM LOOPS Common 2 158840 LP1 Param 91 I(4) 4 scalar 9 LQM Local 675 I(4) 4 scalar 675,678,679,682 LSHU Local 741 I(4) 4 scalar 741,743 LSM Param 3 I(4) 4 scalar LTOP Local 126 I(4) 4 2 17955 328,384,393,497,499,504,507,544,76 7,769,771,832,856,913,914,916,918, 1040,1053,1064,1080,1094,1135,1254 ,1255 LTP1 Local 586 I(4) 4 scalar 586,747,757,913,978,991,1001,1016, 1018,1072,1090,1108,1118 LTP2 Local 748 I(4) 4 scalar 748 LTPK Local 544 I(4) 4 scalar 544,586,589,599,620,630,665,671,67 2,701,735,737,743,744,747,748,750, 754,813,832,839,840,914,925,948,11 35,1139,1145,1147 LTSH Local 725 I(4) 4 scalar 725,761,767 MAPPINGS Common 56 5024 MASKS Common 2 6894720 MAX Func 1258 scalar 1258,1260 MIN Func 1257 scalar 1257,1259 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 187,201,297,310,327,404,414,427,44 0,470,487,831,854,1134,1163,1187,1 199,1208,1217,1231,1242,1253 MYIE1 Scalar 4 I(4) 4 scalar COM 220,344,383,392,453,496,514 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 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 187,201,297,310,327,404,414,427,44 0,470,487,831,854,1134,1163,1187,1 199,1208,1217,1231,1242,1253 MYIS1 Scalar 4 I(4) 4 scalar COM 220,344,383,392,453,496,514 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 Page 33 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 186,200,326,403,413,426,439,469,48 6,1162,1186,1207,1230,1241,1252 MYJE1 Scalar 15 I(4) 4 scalar COM 1198 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 219,296,309,343,382,391,452,495,51 3,830,853,1133,1216 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 186,200,326,403,413,426,439,469,48 6,1162,1186,1198,1230,1241,1252 MYJS1 Scalar 10 I(4) 4 scalar COM 1207 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 219,296,309,343,382,391,452,495,51 3,830,853,1133,1216 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 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 Page 34 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 539 I(4) 4 scalar 539,541,542,876,878,879 NBC Scalar 5 I(4) 4 scalar COM NBOCO Scalar 6 I(4) 4 scalar COM NBOTD Local 122 I(4) 4 1 45 176,841 NBOTS Local 122 I(4) 4 1 45 178,1146 NCNVC Scalar 7 I(4) 4 scalar COM 166 NCP Scalar 6 I(4) 4 scalar COM NDEEP Local 828 I(4) 4 scalar 828,838,845 NDEPTH Local 839 I(4) 4 scalar 839,842,1147,1148 NDPTHD Local 123 I(4) 4 1 45 180,842 NDPTHS Local 123 I(4) 4 1 45 179,1148 NDSTN Local 850 I(4) 4 scalar 850 NDSTP Local 851 I(4) 4 scalar 851 NEGDS Local 1151 I(4) 4 scalar 1151 NEST Scalar 7 I(4) 4 scalar COM NFCST Scalar 5 I(4) 4 scalar COM NHEAT Scalar 3 I(4) 4 scalar COM NNEG Local 845 I(4) 4 scalar 845 NPES Scalar 2 I(4) 4 scalar COM NPHS Scalar 7 I(4) 4 scalar COM NPREC Scalar 6 I(4) 4 scalar COM NRADL Scalar 7 I(4) 4 scalar COM NRADS Scalar 7 I(4) 4 scalar COM NROW Param 79 I(4) 4 scalar 437,448 NSHAL Local 1131 I(4) 4 scalar 1131,1144,1151 NSHDE Scalar 6 I(4) 4 scalar COM NSMUD Local 1194 I(4) 4 scalar 1194,1196 NSTART Scalar 6 I(4) 4 scalar COM NTDDMP Scalar 6 I(4) 4 scalar COM NTOPD Local 122 I(4) 4 1 45 175,840 NTOPS Local 122 I(4) 4 1 45 177,1145 NTSD Scalar 6 I(4) 4 scalar COM NTSTM Scalar 6 I(4) 4 scalar COM OCT90 Param 2 L(4) 4 scalar 717,783 OTSUM Local 1014 R(4) 4 scalar 1014,1023,1031 P Local 133 R(4) 4 2 17955 298,299,317,318 Page 35 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References P00K Local 269 R(4) 4 scalar 269,274,275 P01K Local 271 R(4) 4 scalar 271,274,275 P10K Local 270 R(4) 4 scalar 270,274,275 P11K Local 272 R(4) 4 scalar 272,275 PART1 Local 961 R(4) 4 scalar 961,965 PART2 Local 962 R(4) 4 scalar 962,965 PART3 Local 963 R(4) 4 scalar 963,965 PBM Param 12 R(4) 4 scalar PBOT Local 127 R(4) 4 2 17955 312,314,321,329,500,505,507,519,83 7,855,859,916,919,1041,1054,1065,1 081,1095,1138 PBTK Local 727 R(4) 4 scalar 727,732,1138,1139,1143 PD Scalar 3 R(4) 4 2 17955 COM 190,515,543,730,835,857,885,1140 PDSL Local 129 R(4) 4 2 17955 190,202,222,224,298,312,313,317,32 1,345,361,372,393,577,807,808,809, 810,811,895 PFRZ Param 13 R(4) 4 scalar 628 PHOUR Scalar 2 R(4) 4 2 17955 COM PHYS Common 4 1788916 PK Local 118 R(4) 4 1 45 578,588,589,592,610,621,637,639,64 5,649,650,675,727,735,737,741,744, 751,759,771,896,897,973,974,985 PK0 Local 592 R(4) 4 scalar 592,610,617,629,637,639,645 PKB Local 588 R(4) 4 scalar 588,629,639 PKL Local 222 R(4) 4 scalar 222,229,577,578,579,895,896 PKT Local 589 R(4) 4 scalar 589,617,621,637,645 PKXXXX Local 973 R(4) 4 scalar 973,979,983 PKXXXY Local 974 R(4) 4 scalar 974,979,983,985 PL Scalar 7 R(4) 4 scalar COM 362 PLQ Scalar 8 R(4) 4 scalar COM 346,373 PNO Param 12 R(4) 4 scalar 507,855,916,1139 PONE Param 12 R(4) 4 scalar 314,318 POTSUM Local 1012 R(4) 4 scalar 1012,1024,1032,1049 PP Local 131 R(4) 4 2 17955 362,373 PP1 Local 255 R(4) 4 scalar 255,260,264,274,275 PPK Local 950 R(4) 4 scalar 950,954,958,961,964 PPTASM Common 1 3878292 PPTDAT Scalar 3 R(4) 4 3 53865 COM PQ0 Param 8 R(4) 4 scalar 653,691,751,897 PQM Param 12 R(4) 4 scalar 299,649,675 PREC Scalar 12 R(4) 4 2 17955 COM 808 PRECK Local 699 R(4) 4 scalar 699,708,716,777,797,803,807,808,80 9,810,811 PRECOL Local 139 R(4) 4 2 17955 191,807 PRESK Local 345 R(4) 4 scalar 345,346 PSFCIJ Local 515 R(4) 4 scalar 515,518,543,628,730,731,835,836,85 7,858,1140,1141 PSFCK Local 224 R(4) 4 scalar 224,229,313,314,318 PSH Param 12 R(4) 4 scalar PSHLTR Scalar 19 R(4) 4 2 17955 COM PSHNEW Local 512 R(4) 4 scalar 512,518,731,836,858,1141 PSHU Param 13 R(4) 4 scalar 739,741 PSK Local 119 R(4) 4 1 45 579,650,651,653,691 PSP Local 132 R(4) 4 2 17955 162,280,299,967 PSUM Local 1010 R(4) 4 scalar 1010,1020,1029,1049,1051,1069 Page 36 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References PT Scalar 5 R(4) 4 scalar COM 202,222,224,298,312,313,317,321,34 5,362,373,393,515,543,577,730,835, 857,885,895,1140 PTBL Scalar 25 R(4) 4 2 10184 COM 269,270,271,272,961,962,963,964,96 5 PTOP Local 127 R(4) 4 2 17955 329,393,500,505,507,519,771,837,85 5,859,912,916,919,1041,1054,1065,1 081,1095,1137 PTPK Local 732 R(4) 4 scalar 732,735,737,739,744,912,965,967,11 37,1139,1143 PVRBLS Common 2 5530140 PZ0 Local 885 R(4) 4 scalar 885 Q Scalar 6 R(4) 4 3 807975 COM 204,230,574,656,892,948,1165,1233 Q10 Scalar 16 R(4) 4 2 17955 COM Q2 Scalar 21 R(4) 4 3 807975 COM Q2BT Local 130 R(4) 4 2 17955 Q30 Scalar 14 R(4) 4 2 17955 COM QBT Local 230 R(4) 4 scalar 230,254,277 QK Local 119 R(4) 4 1 45 575,666,704,754,758,893,1021,1025, 1077,1109,1120 QKL Local 574 R(4) 4 scalar 574,575,576,892,893,894,905 QL Local 155 R(4) 4 2 17955 1179,1189,1202,1211,1220,1221,1233 QMOD Local 150 R(4) 4 3 807975 489,815,1120,1165,1189 QNE Local 157 R(4) 4 2 17955 1180,1202,1220 QNEW Local 1077 R(4) 4 scalar 1077,1078 QOTSUM Local 1013 R(4) 4 scalar 1013,1025,1033,1059 QQ Local 131 R(4) 4 2 17955 362,373 QQ1 Local 233 R(4) 4 scalar 233,238,242,252,253,274,275 QQK Local 928 R(4) 4 scalar 928,932,936,944,945,962,964 QREFK Local 118 R(4) 4 1 45 576,653,656,666,688,691,704,894,10 87,1109,1120 QRFKL Local 1073 R(4) 4 scalar 1073,1077,1086,1087 QRFTP Local 1069 R(4) 4 scalar 1069,1073 QS Scalar 6 R(4) 4 2 17955 COM QS0 Scalar 10 R(4) 4 1 134 COM 247,249,939,941 QSATK Local 121 R(4) 4 1 45 751,754,758,897,1078 QSE Local 157 R(4) 4 2 17955 1181,1211,1220 QSHLTR Scalar 18 R(4) 4 2 17955 COM QSUM Local 1011 R(4) 4 scalar 1011,1021,1030,1059,1069 QWBS Scalar 9 R(4) 4 2 17955 COM QZ0 Scalar 5 R(4) 4 2 17955 COM RADIN Scalar 14 R(4) 4 2 17955 COM RADOT Scalar 14 R(4) 4 2 17955 COM RCP Param 89 R(4) 4 scalar 680 RDP Scalar 7 R(4) 4 scalar COM 363 RDP0T Local 617 R(4) 4 scalar 617,621 RDPQ Scalar 8 R(4) 4 scalar COM 374 RDPSUM Local 996 R(4) 4 scalar 996,999,1029,1030 RDQ Scalar 7 R(4) 4 scalar COM 254,948 RDTCNVC Local 167 R(4) 4 scalar 167,1171,1243 RDTH Scalar 7 R(4) 4 scalar COM 232,927 RDTHE Scalar 7 R(4) 4 scalar COM 363 RDTHEQ Scalar 8 R(4) 4 scalar COM 374 RES Scalar 7 R(4) 4 2 17955 COM 190 RESTRT Scalar 95 L(4) 4 scalar COM Page 37 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References RF Scalar 8 R(4) 4 2 17955 COM RHF Param 11 R(4) 4 scalar RHH Local 754 R(4) 4 scalar 754,759,764 RHL Local 758 R(4) 4 scalar 758,759,764 RHMAX Local 755 R(4) 4 scalar 755,760,762 ROI Scalar 6 R(4) 4 scalar COM ROS Scalar 6 R(4) 4 scalar COM ROTSUM Local 1031 R(4) 4 scalar 1031,1032,1033,1034 ROW Param 6 R(4) 4 scalar 89 RSMK Local 472 R(4) 4 scalar 472,473,474,475,476 RTBAR Local 1022 R(4) 4 scalar 1022,1023,1024,1025,1026 RUN Scalar 95 L(4) 4 scalar COM SI Scalar 10 R(4) 4 2 17955 COM SICE Scalar 4 R(4) 4 2 17955 COM SIGMA Scalar 95 L(4) 4 scalar COM SIGT4 Scalar 23 R(4) 4 2 17955 COM SINGLRST Local 2 L(4) 4 scalar SLOP0L Param 81 R(4) 4 scalar 419,430 SLOP0S Param 84 R(4) 4 scalar 407,418 SLOPBL Param 80 R(4) 4 scalar 417,429 SLOPBS Param 83 R(4) 4 scalar 406,416 SLOPE Param 86 R(4) 4 scalar 800 SLOPTL Param 82 R(4) 4 scalar 421,431 SLOPTS Param 85 R(4) 4 scalar 408,420 SM Scalar 4 R(4) 4 2 17955 COM 416,417,418,419,420,421,471,508,72 0,777,787 SMIX Local 970 R(4) 4 scalar 970,979 SMK Local 471 R(4) 4 scalar 471,472,473,474,475,476 SNO Scalar 10 R(4) 4 2 17955 COM SQ Local 253 R(4) 4 scalar 253,254 SQK Local 945 R(4) 4 scalar 945,948 SQS Scalar 10 R(4) 4 1 134 COM 248,250,940,942 SQS00K Local 248 R(4) 4 scalar 248,253,940,945 SQS10K Local 250 R(4) 4 scalar 250,253,942,945 SST Scalar 20 R(4) 4 2 17955 COM STABD Param 9 R(4) 4 scalar 601 STABDL Local 601 R(4) 4 scalar 601,602 STABFC Param 9 R(4) 4 scalar STABS Param 9 R(4) 4 scalar 970 STEFI Param 76 R(4) 4 scalar 508,720 STH Local 134 R(4) 4 2 17955 STHE Scalar 10 R(4) 4 1 76 COM 363 STHEQ Scalar 11 R(4) 4 1 152 COM 374 STRESH Param 8 R(4) 4 scalar 1078 SUBPOST Scalar 2 L(4) 4 scalar COM SUMDE Local 662 R(4) 4 scalar 662,666,667,671 SUMDP Local 663 R(4) 4 scalar 663,668,671,714,989,993,996 SUMDT Local 988 R(4) 4 scalar 988,992,999 T Scalar 5 R(4) 4 3 807975 COM 192,231,384,571,600,889,925,1164,1 232 T00 Local 135 R(4) 4 2 17955 T01 Local 136 R(4) 4 2 17955 T1 Param 8 R(4) 4 scalar T10 Local 135 R(4) 4 2 17955 T11 Local 136 R(4) 4 2 17955 Page 38 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References TAUK Local 168 R(4) 4 scalar 168,703,704,1077,1119,1120 TBT Local 130 R(4) 4 2 17955 TCORR Local 999 R(4) 4 scalar 999,1005 TCUCN Scalar 5 R(4) 4 3 807975 COM 1171,1243 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 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 TFRZ Param 8 R(4) 4 scalar 600 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 THBT Local 132 R(4) 4 2 17955 281,970 THE0 Scalar 10 R(4) 4 1 76 COM 363 THE0Q Scalar 11 R(4) 4 1 152 COM 374 THEAT Scalar 3 R(4) 4 scalar COM THERK Local 119 R(4) 4 1 45 582,594,596,608,618,621 THERKX Local 594 R(4) 4 scalar 594,602,606 THERKY Local 596 R(4) 4 scalar 596,602,606,608 THESP Local 133 R(4) 4 2 17955 189,279,282,364,375 THEVRF Local 120 R(4) 4 1 45 THL Scalar 7 R(4) 4 scalar COM 232,927 THS Scalar 6 R(4) 4 2 17955 COM THSK Local 119 R(4) 4 1 45 652,653,654 THSKL Local 690 R(4) 4 scalar 690,691,692 THTPK Local 925 R(4) 4 scalar 925,927,970 THVMKL Local 905 R(4) 4 scalar 905,907 THVMOD Local 120 R(4) 4 1 45 THVREF Local 120 R(4) 4 1 45 907,1086,1091 THZ0 Scalar 5 R(4) 4 2 17955 COM TK Local 118 R(4) 4 1 45 572,666,702,751,890,897,992,1022,1 026,1109,1110,1119 TKL Local 571 R(4) 4 scalar 571,572,573,702,703,705,889,890,89 1,905 TL Local 155 R(4) 4 2 17955 1182,1188,1200,1209,1218,1219,1232 ,1243 TLAT Scalar 2 R(4) 4 3 807975 COM 817 TMOD Local 149 R(4) 4 3 807975 488,814,1119,1164,1171,1188 TNE Local 156 R(4) 4 2 17955 1183,1200,1218 Page 39 Source Listing CUCNVC 2014-12-17 20:46 Symbol Table CUCNVC.F Name Object Declared Type Bytes Dimen Elements Attributes References TOPO Common 50 17478548 TPSP Local 274 R(4) 4 scalar 274,276,280 TQ Local 254 R(4) 4 scalar 254,255,256 TQK Local 948 R(4) 4 scalar 948,950,951 TRAIN Scalar 4 R(4) 4 3 807975 COM TREF Local 148 R(4) 4 3 807975 192,360,371,384,582 TREFK Local 118 R(4) 4 1 45 573,593,604,618,621,652,666,680,68 7,689,690,703,891,972,981,992,997, 1005,1006,1022,1026,1086,1109,1110 ,1119 TREFKX Local 593 R(4) 4 scalar 593,602,603,604,972,979,980,981 TREL Param 74 R(4) 4 scalar 168 TRFKL Local 1005 R(4) 4 scalar 1005,1006,1007 TSE Local 156 R(4) 4 2 17955 1184,1209,1218 TSHLTR Scalar 18 R(4) 4 2 17955 COM TSKL Local 687 R(4) 4 scalar 687,688 TTBL Scalar 25 R(4) 4 2 10184 COM 360 TTBLEX Subr 360 360,371 TTBLQ Scalar 25 R(4) 4 2 66880 COM 371 TTH Local 232 R(4) 4 scalar 232,233,234 TTHBT Local 231 R(4) 4 scalar 231,232,277,281 TTHES Local 277 R(4) 4 scalar 277,279,282 TTHK Local 927 R(4) 4 scalar 927,928,929 TTVG Scalar 51 R(4) 4 2 92971 COM TWBS Scalar 9 R(4) 4 2 17955 COM U Scalar 5 R(4) 4 3 807975 COM U10 Scalar 17 R(4) 4 2 17955 COM U30 Scalar 15 R(4) 4 2 17955 COM UNIL Param 2 L(4) 4 scalar 411,437,786 UNIS Param 2 L(4) 4 scalar 401,437,784 USTAR Scalar 3 R(4) 4 2 17955 COM 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 WF1 Local 137 R(4) 4 2 17955 441,454,473 WF2 Local 137 R(4) 4 2 17955 442,457,474 WF3 Local 138 R(4) 4 2 17955 443,460,475 WF4 Local 138 R(4) 4 2 17955 444,463,476 Z0 Scalar 3 R(4) 4 2 17955 COM ZERO2 Subr 159 159,160,161,162,1179,1180,1181,118 2,1183,1184 Page 40 Source Listing CUCNVC 2014-12-17 20:46 Subprograms/Common Blocks CUCNVC.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References ACMCLH Common 2 6463824 CNVCLD Common 2 287280 CTLBLK Common 4 112 CUCNVC Subr 1 GLB_TABLE Common 37 128 INDX Common 2 6564 INDXG Common 8 8084 LOOPS Common 2 158840 MAPPINGS Common 56 5024 MASKS Common 2 6894720 MPPCOM Common 1 1464 PHYS Common 4 1788916 PPTASM Common 1 3878292 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 Page 41 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F -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 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/, Page 42 Source Listing CUCNVC 2014-12-17 20:46 CUCNVC.F /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/, /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : CUCNVC.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100