Page 1 Source Listing TURBL 2014-12-17 20:47 /tmpnwprd/ifort3KJebo.i 1 # 1 "TURBL.F" 2 SUBROUTINE TURBL 3 C ****************************************************************** 4 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 5 C . . . 6 C SUBPROGRAM: TURBL VERTICAL TURBULENT EXCHANGE 7 C PRGRMMR: JANJIC ORG: W/NP2 DATE: 95-03-20 8 C 9 C ABSTRACT: 10 C TURBL UPDATES THE TURBULENT KINETIC ENERGY WITH THE PROD- 11 C UCTION/DISSIPATION TERM AND THE VERTICAL DIFFUSION TERM 12 C DIFFUSION TERM (USING AN IMPLICIT FORMULATION). EXCHANGE 13 C COEFFICIENTS FOR THE SURFACE AND FOR ALL LAYER INTERFACES 14 C ARE THEN COMPUTED AND THE EXCHANGE IS EXECUTED. 15 C 16 C PROGRAM HISTORY LOG: 17 C 95-03-15 JANJIC - ORIGINATOR 18 C 95-03-28 BLACK - CONVERSION FROM 1-D TO 2-D IN HORIZONTAL 19 C 96-03-29 BLACK - ADDED EXTERNAL EDGE; REMOVED SCRCH COMMON 20 C 96-07-19 MESINGER - ADDED Z0 EFFECTIVE 21 C 98-??-?? TUCCILLO - MODIFIED FOR CLASS VIII PARALLELISM 22 C 98-10-27 BLACK - PARALLEL CHANGES INTO MOST RECENT CODE 23 C 24 C USAGE: CALL TURBL FROM MAIN PROGRAM EBU 25 C INPUT ARGUMENT LIST: 26 C NONE 27 C 28 C OUTPUT ARGUMENT LIST: 29 C NONE 30 C 31 C OUTPUT FILES: 32 C NONE 33 C 34 C SUBPROGRAMS CALLED: 35 C UNIQUE: MIXLEN 36 C PRODQ2 37 C DIFCOF 38 C SFCDIF 39 C VDIFH 40 C VDIFQ 41 C VDIFV 42 C 43 C LIBRARY: NONE 44 C 45 C COMMON BLOCKS: CTLBLK 46 C LOOPS 47 C MASKS 48 C DYNAM 49 C PHYS2 50 C VRBLS 51 C PVRBLS 52 C INDX 53 C Z0EFFT 54 C 55 C ATTRIBUTES: 56 C LANGUAGE: FORTRAN 90 57 C MACHINE : IBM SP Page 2 Source Listing TURBL 2014-12-17 20:47 TURBL.F 58 C$$$ 59 C*********************************************************************** 60 C----------------------------------------------------------------------- 61 C 62 INCLUDE "EXCHM.h" 70 INCLUDE "parmeta" 93 INCLUDE "mpp.h" 94 # 1 "./sp.h" 1 95 # 4 96 97 # 65 "TURBL.F" 2 98 C----------------------------------------------------------------------- 157 P A R A M E T E R 158 & (KTMQ2=1,CAPA=0.28589641,G=9.8,RG=1./G,ROG=287.04/G 159 &, EPSZ=1.E-4,EPSQ2=0.2 160 &, IMJM=IM*JM-JM/2,LM1=LM-1,LP1=LM+1,JAM=6+2*(JM-10) 161 &, ITB=76,JTB=134,ITBQ=152,JTBQ=440 162 &, NHRZ=(IDIM2-IDIM1+1)*(JDIM2-JDIM1+1)) 163 C----------------------------------------------------------------------- 164 INCLUDE "CTLBLK.comm" 165 C----------------------------------------------------------------------- 175 INCLUDE "LOOPS.comm" 176 C----------------------------------------------------------------------- 183 INCLUDE "MASKS.comm" 184 C----------------------------------------------------------------------- 195 INCLUDE "DYNAM.comm" 196 C----------------------------------------------------------------------- 212 INCLUDE "PHYS2.comm" 213 C----------------------------------------------------------------------- 241 INCLUDE "VRBLS.comm" 242 C----------------------------------------------------------------------- 252 INCLUDE "PVRBLS.comm" 253 C----------------------------------------------------------------------- 283 INCLUDE "INDX.comm" 284 C----------------------------------------------------------------------- 296 INCLUDE "Z0EFFT.comm" 297 C----------------------------------------------------------------------- 304 L O G I C A L 305 & RUN,FIRST,RESTRT,SIGMA 306 C----------------------------------------------------------------------- 307 R E A L 308 & CKLQ(IDIM1:IDIM2,JDIM1:JDIM2) 309 &,CT (IDIM1:IDIM2,JDIM1:JDIM2) 310 &,APE (IDIM1:IDIM2,JDIM1:JDIM2,LM) 311 &,AKH (IDIM1:IDIM2,JDIM1:JDIM2,LM1) 312 &,AKM (IDIM1:IDIM2,JDIM1:JDIM2,LM1) 313 &,ZINT(IDIM1:IDIM2,JDIM1:JDIM2,LP1) 314 &,UZ0H(IDIM1:IDIM2,JDIM1:JDIM2) 315 &,VZ0H(IDIM1:IDIM2,JDIM1:JDIM2) 316 C 317 R E A L 318 & AKMCOL(IDIM1:IDIM2,JDIM1:JDIM2,LM1) 319 &,AKHCOL(IDIM1:IDIM2,JDIM1:JDIM2,LM1) 320 &,AKMSV (IDIM1:IDIM2,JDIM1:JDIM2) 321 &,ZCOL (IDIM1:IDIM2,JDIM1:JDIM2,LP1) 322 &,UCOL (IDIM1:IDIM2,JDIM1:JDIM2,LM) 323 &,VCOL (IDIM1:IDIM2,JDIM1:JDIM2,LM) Page 3 Source Listing TURBL 2014-12-17 20:47 TURBL.F 324 C 325 R E A L 326 & AKH_T (LM1,IDIM1:IDIM2,JDIM1:JDIM2) 327 &,AKM_T (LM1,IDIM1:IDIM2,JDIM1:JDIM2) 328 &,APECOL_T(LM,IDIM1:IDIM2,JDIM1:JDIM2) 329 &,ZCOL_T (LP1,IDIM1:IDIM2,JDIM1:JDIM2) 330 &,ZCOL_T2 (LP1,IDIM1:IDIM2,JDIM1:JDIM2) 331 &,UCOL_T (LM,IDIM1:IDIM2,JDIM1:JDIM2) 332 &,VCOL_T (LM,IDIM1:IDIM2,JDIM1:JDIM2) 333 &,TCOL_T (LM,IDIM1:IDIM2,JDIM1:JDIM2) 334 &,QCOL_T (LM,IDIM1:IDIM2,JDIM1:JDIM2) 335 &,Q2COL_T (LM,IDIM1:IDIM2,JDIM1:JDIM2) 336 C 337 R E A L 338 & GM(LM1),GH(LM1),EL(LM1),ZEFF(4) 339 C 340 C----------------------------------------------------------------------- 341 C*** 342 C*** THE FOLLOWING ARE USED FOR TIMIMG PURPOSES ONLY 343 C*** 344 real*8 timef 345 real nhb_tim,mpp_tim,init_tim 346 common/timing/surfce_tim,nhb_tim,res_tim,exch_tim 347 save /timing/ 348 C*********************************************************************** 349 C----------------------------------------------------------------------- 350 CALL ZERO3(AKM,LM1) 351 CALL ZERO3(ZINT,LP1) 352 CALL ZERO3_T(AKH_T,LM1) 353 CALL ZERO3_T(AKM_T,LM1) 354 CALL ZERO2(UZ0H) 355 CALL ZERO2(VZ0H) 356 C----------------------------------------------------------------------- 357 C*** 358 C*** COMPUTE THE HEIGHTS OF THE LAYER INTERFACES AND THE EXNER FUNCTION 359 C*** 360 !$omp parallel do 361 DO J=MYJS_P1,MYJE_P1 ! This line is correct 362 c DO J=MYJS2_P1,MYJE2_P1 ! This line matches operations 363 DO I=MYIS_P1,MYIE_P1 364 ZINT(I,J,LP1)=EPSZ 365 IF(SIGMA)ZINT(I,J,LP1)=RG*FIS(I,J) 366 ENDDO 367 ENDDO 368 C 369 DO 10 L=LM,1,-1 370 !$omp parallel do private (apests,pdsl) 371 DO J=MYJS1_P1,MYJE1_P1 ! This line is correct 372 c DO J=MYJS2_P1,MYJE2_P1 ! This line matches operations 373 DO I=MYIS_P1,MYIE_P1 374 PDSL=PD(I,J)*RES(I,J) 375 APESTS=PDSL*AETA(L)+PT 376 C 377 ZINT(I,J,L)=ZINT(I,J,L+1)+T(I,J,L)/APESTS 378 1 *PDSL*DETA(L)*ROG*(Q(I,J,L)*0.608+1.) 379 ZINT(I,J,L)=(ZINT(I,J,L)-DFRLG(L))*HTM(I,J,L)+DFRLG(L) 380 C Page 4 Source Listing TURBL 2014-12-17 20:47 TURBL.F 381 APE(I,J,L)=(1.E5/APESTS)**CAPA 382 ENDDO 383 ENDDO 384 10 CONTINUE 385 C----------------------------------------------------------------------- 386 C*** 387 C*** REMOVE NEGATIVE Q2 388 C*** 389 !$omp parallel do 390 DO 40 L=1,LM 391 DO J=MYJS_P1,MYJE_P1 392 DO I=MYIS_P1,MYIE_P1 393 Q2(I,J,L)=AMAX1(Q2(I,J,L)*HBM2(I,J),EPSQ2) 394 ENDDO 395 ENDDO 396 40 CONTINUE 397 C----------------------------------------------------------------------- 398 !$omp parallel do 399 DO J=MYJS2_P1,MYJE2_P1 400 DO I=MYIS_P1,MYIE_P1 401 UZ0H(I,J)=(UZ0(I+IHE(J),J)+UZ0(I+IHW(J),J) 402 1 +UZ0(I,J+1)+UZ0(I,J-1))*HBM2(I,J)*0.25 403 VZ0H(I,J)=(VZ0(I+IHE(J),J)+VZ0(I+IHW(J),J) 404 1 +VZ0(I,J+1)+VZ0(I,J-1))*HBM2(I,J)*0.25 405 ENDDO 406 ENDDO 407 C----------------------------------------------------------------------- 408 C----------------------------------------------------------------------- 409 C*** PREPARE THE EXCHANGE COEFFICIENTS 410 C----------------------------------------------------------------------- 411 C----------------------------------------------------------------------- 412 C*** 413 C*** COMPUTE VELOCITY COMPONENTS AT H POINTS 414 C*** 415 !$omp parallel do private(rwmsk,wmsk) 416 DO 60 L=1,LM 417 C 418 DO J=MYJS2_P1,MYJE2_P1 419 DO I=MYIS_P1,MYIE_P1 420 WMSK=VTM(I+IHE(J),J,L)+VTM(I+IHW(J),J,L) 421 1 +VTM(I,J+1,L)+VTM(I,J-1,L) 422 IF(WMSK.GT.0.)THEN 423 RWMSK=1./WMSK 424 UCOL(I,J,L)=(U(I+IHE(J),J,L)*VTM(I+IHE(J),J,L) 425 1 +U(I+IHW(J),J,L)*VTM(I+IHW(J),J,L) 426 2 +U(I,J+1,L)*VTM(I,J+1,L)+U(I,J-1,L)*VTM(I,J-1,L)) 427 3 *RWMSK 428 VCOL(I,J,L)=(V(I+IHE(J),J,L)*VTM(I+IHE(J),J,L) 429 1 +V(I+IHW(J),J,L)*VTM(I+IHW(J),J,L) 430 2 +V(I,J+1,L)*VTM(I,J+1,L)+V(I,J-1,L)*VTM(I,J-1,L)) 431 3 *RWMSK 432 ELSE 433 UCOL(I,J,L)=0. 434 VCOL(I,J,L)=0. 435 ENDIF 436 ENDDO 437 ENDDO Page 5 Source Listing TURBL 2014-12-17 20:47 TURBL.F 438 60 CONTINUE 439 C*** 440 C*** FILL TRANSPOSED ARRAYS 441 C*** 442 !$omp parallel sections 443 !$omp section 444 CALL SGETMO(T,NHRZ,NHRZ,LM,TCOL_T,LM) 445 !$omp section 446 CALL SGETMO(Q,NHRZ,NHRZ,LM,QCOL_T,LM) 447 !$omp section 448 CALL SGETMO(APE,NHRZ,NHRZ,LM,APECOL_T,LM) 449 !$omp section 450 CALL SGETMO(Q2,NHRZ,NHRZ,LM,Q2COL_T,LM) 451 !$omp section 452 CALL SGETMO(ZINT,NHRZ,NHRZ,LP1,ZCOL_T,LP1) 453 !$omp section 454 CALL SGETMO(UCOL,NHRZ,NHRZ,LM,UCOL_T,LM) 455 !$omp section 456 CALL SGETMO(VCOL,NHRZ,NHRZ,LM,VCOL_T,LM) 457 !$omp end parallel sections 458 C---------------------------------------------------------------------- 459 C*** 460 C*** FIND THE MIXING LENGTH 461 C*** 462 !$omp parallel do private(el,gh,gm,hpbl,lmhk,lmhm,lmhp,lpbl) 463 !$omp& private(ulm,vlm,wstar,zeff) 464 DO 100 J=MYJS2_P1,MYJE2_P1 465 DO 100 I=MYIS_P1,MYIE1_P1 466 467 LMHK=LMH(I,J) 468 LMHP=LMHK+1 469 LMHM=LMHK-1 470 C 471 CALL MIXLEN(LMHK,LPBL,HPBL,UCOL_T(1,I,J),VCOL_T(1,I,J) 472 1, TCOL_T(1,I,J),QCOL_T(1,I,J),Q2COL_T(1,I,J) 473 2, APECOL_T(1,I,J),ZCOL_T(1,I,J),GM,GH,EL) 474 C 475 C----------------------------------------------------------------------- 476 C*** 477 C*** SOLVE FOR THE PRODUCTION/DISSIPATION OF 478 C*** THE TURBULENT KINETIC ENERGY 479 C*** 480 C 481 CALL PRODQ2(LMHK,DTQ2,USTAR(I,J),GM,GH,EL,Q2COL_T(1,I,J)) 482 C 483 C----------------------------------------------------------------------- 484 C*** 485 C*** FIND THE EXCHANGE COEFFICIENTS IN THE FREE ATMOSPHERE 486 C*** 487 CALL DIFCOF(LMHK,GM,GH,EL,Q2COL_T(1,I,J) 488 1, ZCOL_T(1,I,J),AKM_T(1,I,J),AKH_T(1,I,J)) 489 C----------------------------------------------------------------------- 490 C*** 491 C*** CARRY OUT THE VERTICAL DIFFUSION OF 492 C*** TURBULENT KINETIC ENERGY 493 C*** 494 C Page 6 Source Listing TURBL 2014-12-17 20:47 TURBL.F 495 CALL VDIFQ(LMHK,KTMQ2,DTQ2,Q2COL_T(1,I,J),EL,ZCOL_T(1,I,J)) 496 C----------------------------------------------------------------------- 497 C*** 498 C*** FIND THE Z0 EFFECTIVE 499 C*** 500 ZEFF(1)=ZEFFIJ(I,J,1) 501 ZEFF(2)=ZEFFIJ(I,J,2) 502 ZEFF(3)=ZEFFIJ(I,J,3) 503 ZEFF(4)=ZEFFIJ(I,J,4) 504 C----------------------------------------------------------------------- 505 C*** 506 C*** FIND THE SURFACE EXCHANGE COEFFICIENTS 507 C*** 508 ULM=UCOL(I,J,LMHK) 509 VLM=VCOL(I,J,LMHK) 510 C 511 CALL SFCDIF(LMHK,SM(I,J),THS(I,J),QS(I,J) 512 1, UZ0H(I,J),VZ0H(I,J),THZ0(I,J),QZ0(I,J) 513 2, USTAR(I,J),WSTAR 514 3, Z0(I,J),ZEFF,AKMS(I,J),AKHS(I,J),HPBL,CT(I,J) 515 4, U10(I,J),V10(I,J),TSHLTR(I,J),TH10(I,J) 516 5, QSHLTR(I,J),Q10(I,J) 517 5, TH30(I,J),Q30(I,J),U30(I,J),V30(I,J) 518 6, ULM,VLM,TCOL_T(1,I,J),QCOL_T(1,I,J) 519 7, APECOL_T(1,I,J),ZCOL_T(1,I,J),PD(I,J),PT 520 8, T(I,J,LMHK)) 521 C 522 523 PDSL=PD(I,J)*RES(I,J) 524 DPLM=PDSL*DETA(LMHK)*0.5 525 DZLM=ROG*DPLM*T(I,J,LMHK)*(1.+0.608*Q(I,J,LMHK)) 526 & /(PDSL*AETA(LMHK)+PT) 527 528 IF (DZLM.LT.30.0) THEN 529 530 DPLM1=PDSL*DETA(LMHK-1)*0.5 531 DZLM1=ROG*DPLM1*T(I,J,LMHK-1)*(1.+0.608*Q(I,J,LMHK-1)) 532 & /(PDSL*AETA(LMHK-1)+PT) 533 534 FAC1=(30.-DZLM)/(DZLM+DZLM1) 535 FAC2=((2.*DZLM+DZLM1)-30.)/(DZLM+DZLM1) 536 537 if (FAC1.GT.1.0 .or. FAC1.LT.0.0) then 538 print *, " *****************************", i,j,FAC1 539 end if 540 if (FAC2.GT.1.0 .or. FAC2.LT.0.0) then 541 print *, " *****************************", i,j,FAC2 542 end if 543 544 T30=FAC2*T(I,J,LMHK)+FAC1*T(I,J,LMHK-1) 545 P30=(PD(I,J)+PT)*EXP(-30.0*G/(287.04*T(I,J,LMHK))) 546 TH30(I,J)=T30*(1.E5/P30)**0.28589641 547 Q30(I,J)=FAC2*Q(I,J,LMHK)+FAC1*Q(I,J,LMHK-1) 548 U30(I,J)=FAC2*UCOL(I,J,LMHK)+FAC1*UCOL(I,J,LMHK-1) 549 V30(I,J)=FAC2*VCOL(I,J,LMHK)+FAC1*VCOL(I,J,LMHK-1) 550 551 END IF Page 7 Source Listing TURBL 2014-12-17 20:47 TURBL.F 552 553 100 CONTINUE 554 555 556 C------------------------------------------------------------------------ 557 C*** 558 C*** FILL STANDARD ARRAYS FROM TRANSPOSED ARRAYS 559 C*** 560 !$omp parallel sections 561 !$omp section 562 CALL SGETMO(Q2COL_T,LM,LM,NHRZ,Q2,NHRZ) 563 !$omp section 564 CALL SGETMO(AKH_T,LM1,LM1,NHRZ,AKH,NHRZ) 565 !$omp section 566 CALL SGETMO(AKM_T,LM1,LM1,NHRZ,AKM,NHRZ) 567 !$omp end parallel sections 568 C----------------------------------------------------------------------- 569 C*** 570 C*** UNCOMPUTED LOCATIONS MUST BE FILLED IN FOR THE POST-PROCESSOR 571 C*** 572 IIM=IM-MY_IS_GLB+1 573 JJM=JM-MY_JS_GLB+1 574 C 575 C*** EASTERN GLOBAL BOUNDARY 576 C 577 IF(MY_IE_GLB.EQ.IM)THEN 578 DO J=1,JM 579 IF(J.GE.MY_JS_GLB.AND.J.LE.MY_JE_GLB)THEN 580 JJ=J-MY_JS_GLB+1 581 TH10(IIM,JJ)=TH10(IIM-1,JJ) 582 Q10(IIM,JJ)=Q10(IIM-1,JJ) 583 U10(IIM,JJ)=U10(IIM-1,JJ) 584 V10(IIM,JJ)=V10(IIM-1,JJ) 585 TSHLTR(IIM,JJ)=TSHLTR(IIM-1,JJ) 586 QSHLTR(IIM,JJ)=QSHLTR(IIM-1,JJ) 587 TH30(IIM,JJ)=TH30(IIM-1,JJ) 588 Q30(IIM,JJ)=Q30(IIM-1,JJ) 589 U30(IIM,JJ)=U30(IIM-1,JJ) 590 V30(IIM,JJ)=V30(IIM-1,JJ) 591 ENDIF 592 ENDDO 593 ENDIF 594 C 595 C*** SOUTHERN GLOBAL BOUNDARY 596 C 597 IF(MY_JS_GLB.EQ.1)THEN 598 DO J=1,2 599 DO I=1,IM 600 IF(I.GE.MY_IS_GLB.AND.I.LE.MY_IE_GLB)THEN 601 II=I-MY_IS_GLB+1 602 TH10(II,J)=TH10(II,3) 603 Q10(II,J)=Q10(II,3) 604 U10(II,J)=U10(II,3) 605 V10(II,J)=V10(II,3) 606 TSHLTR(II,J)=TSHLTR(II,3) 607 QSHLTR(II,J)=QSHLTR(II,3) 608 TH30(II,J)=TH30(II,3) Page 8 Source Listing TURBL 2014-12-17 20:47 TURBL.F 609 Q30(II,J)=Q30(II,3) 610 U30(II,J)=U30(II,3) 611 V30(II,J)=V30(II,3) 612 ENDIF 613 ENDDO 614 ENDDO 615 ENDIF 616 C 617 C*** NORTHERN GLOBAL BOUNDARY 618 C 619 IF(MY_JE_GLB.EQ.JM)THEN 620 DO J=JM-1,JM 621 JJ=J-MY_JS_GLB+1 622 DO I=1,IM 623 IF(I.GE.MY_IS_GLB.AND.I.LE.MY_IE_GLB)THEN 624 II=I-MY_IS_GLB+1 625 TH10(II,JJ)=TH10(II,JJM-2) 626 Q10(II,JJ)=Q10(II,JJM-2) 627 U10(II,JJ)=U10(II,JJM-2) 628 V10(II,JJ)=V10(II,JJM-2) 629 TSHLTR(II,JJ)=TSHLTR(II,JJM-2) 630 QSHLTR(II,JJ)=QSHLTR(II,JJM-2) 631 TH30(II,JJ)=TH30(II,JJM-2) 632 Q30(II,JJ)=Q30(II,JJM-2) 633 U30(II,JJ)=U30(II,JJM-2) 634 V30(II,JJ)=V30(II,JJM-2) 635 ENDIF 636 ENDDO 637 ENDDO 638 ENDIF 639 C----------------------------------------------------------------------- 640 C----------------------------------------------------------------------- 641 C----------------------------------------------------------------------- 642 btim=timef() 643 CALL EXCH(UZ0H,1,VZ0H,1,1,1) 644 exch_tim=exch_tim+timef()-btim 645 C*** 646 C*** AVERAGE UZ0 AND VZ0 BACK TO V POINTS 647 C*** 648 !$omp parallel do 649 DO 125 J=MYJS2,MYJE2 650 DO 125 I=MYIS,MYIE 651 UZ0(I,J)=(UZ0H(I+IVE(J),J)*HBM2(I+IVE(J),J) 652 1 +UZ0H(I+IVW(J),J)*HBM2(I+IVW(J),J) 653 2 +UZ0H(I,J+1)*HBM2(I,J+1)+UZ0H(I,J-1)*HBM2(I,J-1))*0.25 654 VZ0(I,J)=(VZ0H(I+IVE(J),J)*HBM2(I+IVE(J),J) 655 1 +VZ0H(I+IVW(J),J)*HBM2(I+IVW(J),J) 656 2 +VZ0H(I,J+1)*HBM2(I,J+1)+VZ0H(I,J-1)*HBM2(I,J-1))*0.25 657 125 CONTINUE 658 C----------------------------------------------------------------------- 659 C*** 660 C*** EXECUTE THE GROUND PROCESSES 661 C*** 662 CALL SURFCE(APE(IDIM1,JDIM1,1),ZINT(IDIM1,JDIM1,1) 663 1, CKLQ(IDIM1,JDIM1)) 664 C 665 C----------------------------------------------------------------------- Page 9 Source Listing TURBL 2014-12-17 20:47 TURBL.F 666 C----------------------------------------------------------------------- 667 C*** EXECUTE THE VERTICAL EXCHANGE 668 C----------------------------------------------------------------------- 669 C----------------------------------------------------------------------- 670 btim=timef() 671 CALL EXCH(AKMS,1,AKM,LM1,ZINT,LP1,1,1) 672 exch_tim=exch_tim+timef()-btim 673 C 674 !$omp parallel do 675 DO L=1,LM1 676 DO J=MYJS2,MYJE2 677 DO I=MYIS,MYIE 678 AKMCOL(I,J,L)=(AKM(I+IVE(J),J,L)*HBM2(I+IVE(J),J) 679 1 +AKM(I+IVW(J),J,L)*HBM2(I+IVW(J),J) 680 2 +AKM(I,J+1,L)*HBM2(I,J+1)+AKM(I,J-1,L)*HBM2(I,J-1)) 681 3 *VTM(I,J,L)*VBM2(I,J)*0.25 682 AKHCOL(I,J,L)=AKH(I,J,L)*HTM(I,J,L)*HBM2(I,J) 683 ENDDO 684 ENDDO 685 ENDDO 686 C 687 !$omp parallel do 688 DO J=MYJS2,MYJE2 689 DO I=MYIS,MYIE 690 THZ0(I,J)=(1.-SM(I,J))*THS(I,J)+SM(I,J)*THZ0(I,J) 691 QZ0 (I,J)=(1.-SM(I,J))*QS (I,J)+SM(I,J)*QZ0 (I,J) 692 AKMSV(I,J)=(AKMS(I+IVE(J),J)*HBM2(I+IVE(J),J) 693 1 +AKMS(I+IVW(J),J)*HBM2(I+IVW(J),J) 694 2 +AKMS(I,J+1)*HBM2(I,J+1)+AKMS(I,J-1)*HBM2(I,J-1)) 695 3 *VBM2(I,J)*0.25 696 ENDDO 697 ENDDO 698 C 699 !$omp parallel do 700 DO L=1,LP1 701 DO J=MYJS2,MYJE2 702 DO I=MYIS,MYIE 703 ZCOL(I,J,L)=0.25*(ZINT(I+IVE(J),J,L)+ZINT(I+IVW(J),J,L) 704 1 +ZINT(I,J+1,L)+ZINT(I,J-1,L)) 705 ENDDO 706 ENDDO 707 ENDDO 708 C*** 709 C*** TRANSPOSE ARRAYS 710 C*** 711 !$omp parallel sections 712 !$omp section 713 CALL SGETMO(ZCOL,NHRZ,NHRZ,LP1,ZCOL_T2,LP1) 714 !$omp section 715 CALL SGETMO(U,NHRZ,NHRZ,LM,UCOL_T,LM) 716 !$omp section 717 CALL SGETMO(V,NHRZ,NHRZ,LM,VCOL_T,LM) 718 !$omp section 719 CALL SGETMO(AKHCOL,NHRZ,NHRZ,LM1,AKH_T,LM1) 720 !$omp section 721 CALL SGETMO(AKMCOL,NHRZ,NHRZ,LM1,AKM_T,LM1) 722 !$omp end parallel sections Page 10 Source Listing TURBL 2014-12-17 20:47 TURBL.F 723 C----------------------------------------------------------------------- 724 !$omp parallel do private(lmhk,lmvk) 725 DO 200 J=MYJS2,MYJE2 726 DO 200 I=MYIS,MYIE1 727 C 728 LMHK=LMH(I,J) 729 LMVK=LMV(I,J) 730 C*** 731 C*** CARRY OUT THE VERTICAL DIFFUSION OF 732 C*** TEMPERATURE AND WATER VAPOR 733 C*** 734 CALL VDIFH(LMHK,KTMQ2,DTQ2,THZ0(I,J),QZ0(I,J),AKHS(I,J) 735 1, CT(I,J),CKLQ(I,J) 736 2, TCOL_T(1,I,J),QCOL_T(1,I,J),AKH_T(1,I,J) 737 3, APECOL_T(1,I,J),ZCOL_T(1,I,J)) 738 C 739 C----------------------------------------------------------------------- 740 C*** 741 C*** CARRY OUT THE VERTICAL DIFFUSION OF 742 C*** VELOCITY COMPONENTS 743 C*** 744 CALL VDIFV(LMVK,KTMQ2,DTQ2,UZ0(I,J),VZ0(I,J) 745 1, AKMSV(I,J),UCOL_T(1,I,J),VCOL_T(1,I,J) 746 2, AKM_T(1,I,J),ZCOL_T2(1,I,J)) 747 C 748 C----------------------------------------------------------------------- 749 200 CONTINUE 750 C----------------------------------------------------------------------- 751 C----------------------------------------------------------------------- 752 C*** 753 C*** FILL STANDARD ARRAYS FROM TRANSPOSED ARRAYS 754 C*** 755 !$omp parallel sections 756 !$omp section 757 CALL SGETMO(QCOL_t,LM,LM,NHRZ,Q,NHRZ) 758 !$omp section 759 CALL SGETMO(TCOL_t,LM,LM,NHRZ,T,NHRZ) 760 !$omp section 761 CALL SGETMO(UCOL_t,LM,LM,NHRZ,U,NHRZ) 762 !$omp section 763 CALL SGETMO(VCOL_t,LM,LM,NHRZ,V,NHRZ) 764 !$omp end parallel sections 765 C 766 C*** 767 C*** REMOVE NEGATIVE Q2 768 C*** 769 !$omp parallel do 770 DO L=1,LM 771 DO J=MYJS,MYJE 772 DO I=MYIS,MYIE 773 Q2(I,J,L)=AMAX1(Q2(I,J,L)*HBM2(I,J),EPSQ2) 774 ENDDO 775 ENDDO 776 ENDDO 777 C---------------------------------------------------------------------- 778 RETURN 779 END Page 11 Source Listing TURBL 2014-12-17 20:47 Entry Points TURBL.F ENTRY POINTS Name turbl_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 171 156 100 Label 340 251,252 125 Label 444 436,437 200 Label 536 512,513 40 Label 183 177 60 Label 225 203 ACCLIQ Scalar 13 R(4) 4 2 17955 COM ACPREC Scalar 12 R(4) 4 2 17955 COM AETA Scalar 4 R(4) 4 1 45 COM 162,313,319 AETA2 Scalar 9 R(4) 4 1 45 COM AFSI Scalar 26 R(4) 4 2 17955 COM AKH Local 98 R(4) 4 3 790020 351,469 AKHCOL Local 106 R(4) 4 3 790020 469,506 AKHS Scalar 7 R(4) 4 2 17955 COM 301,521 AKH_T Local 113 R(4) 4 3 790020 139,275,351,506,523 AKM Local 99 R(4) 4 3 790020 137,353,458,465,466,467 AKMCOL Local 105 R(4) 4 3 790020 465,508 AKMS Scalar 7 R(4) 4 2 17955 COM 301,458,479,480,481 AKMSV Local 107 R(4) 4 2 17955 479,532 AKM_T Local 114 R(4) 4 3 790020 140,275,353,508,533 ALBASE Scalar 20 R(4) 4 2 17955 COM ALBEDO Scalar 21 R(4) 4 2 17955 COM AMAX1 Func 180 scalar 180,560 APE Local 97 R(4) 4 3 807975 168,235,449 APECOL_T Local 115 R(4) 4 3 807975 235,260,306,524 APESTS Local 162 R(4) 4 scalar 162,164,168 BTIM Local 429 R(4) 4 scalar 429,431,457,459 CAPA Param 67 R(4) 4 scalar 168 CI Scalar 6 R(4) 4 scalar COM CKLQ Local 95 R(4) 4 2 17955 450,522 CLDEFI Scalar 11 R(4) 4 2 17955 COM CNVBOT Scalar 18 R(4) 4 2 17955 COM CNVTOP Scalar 18 R(4) 4 2 17955 COM CPGFU Scalar 9 R(4) 4 2 17955 COM CPGFV Scalar 3 R(4) 4 scalar COM CS Scalar 6 R(4) 4 scalar COM CT Local 96 R(4) 4 2 17955 301,522 CTLBLK Common 4 112 CUPREC Scalar 13 R(4) 4 2 17955 COM CURV Scalar 9 R(4) 4 2 17955 COM CZEN Scalar 16 R(4) 4 2 17955 COM CZMEAN Scalar 23 R(4) 4 2 17955 COM DDMPU Scalar 12 R(4) 4 2 17955 COM DDMPV Scalar 12 R(4) 4 2 17955 COM Page 12 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References DETA Scalar 4 R(4) 4 1 45 COM 165,311,317 DETA2 Scalar 9 R(4) 4 1 45 COM DFL Scalar 4 R(4) 4 1 46 COM DFRLG Scalar 9 R(4) 4 1 46 COM 166 DI Scalar 6 R(4) 4 scalar COM DIFCOF Subr 274 274 DPLM Local 311 R(4) 4 scalar 311,312 DPLM1 Local 317 R(4) 4 scalar 317,318 DS Scalar 6 R(4) 4 scalar COM DT Scalar 6 R(4) 4 scalar COM DTD Scalar 5 R(4) 4 scalar COM DTQ2 Scalar 5 R(4) 4 scalar COM 268,282,521,531 DX Scalar 8 R(4) 4 2 17955 COM DY Scalar 3 R(4) 4 scalar COM DYNAM Common 2 725404 DZLM Local 312 R(4) 4 scalar 312,315,321,322 DZLM1 Local 318 R(4) 4 scalar 318,321,322 EF4T Scalar 3 R(4) 4 scalar COM EL Local 125 R(4) 4 1 44 260,268,274,282 EM Scalar 6 R(4) 4 1 760 COM EMT Scalar 6 R(4) 4 1 760 COM EN Scalar 3 R(4) 4 scalar COM ENT Scalar 3 R(4) 4 scalar COM EPSQ2 Param 68 R(4) 4 scalar 180,560 EPSR Scalar 13 R(4) 4 2 17955 COM EPSZ Param 68 R(4) 4 scalar 151 ETA Scalar 4 R(4) 4 1 46 COM EXCH Local 2 scalar 430,458 EXCH0 Subr 3 EXCH00 Subr 3 430 EXCH0001111 Subr 3 EXCH01 Subr 3 EXCH011 Subr 3 458 EXCH1 Subr 4 EXCH11 Subr 4 EXCH111 Subr 4 EXCH1111 Subr 4 EXCH11111 Subr 4 EXCH111111 Subr 5 EXCHM Module 1 1 EXCH_TIM Scalar 133 R(4) 4 scalar COM 431,459 EXP Func 332 scalar 332 F Scalar 11 R(4) 4 2 17955 COM F4D Scalar 3 R(4) 4 scalar COM F4Q Scalar 3 R(4) 4 scalar COM F4Q2 Scalar 4 R(4) 4 1 45 COM FAC1 Local 321 R(4) 4 scalar 321,324,325,331,334,335,336 FAC2 Local 322 R(4) 4 scalar 322,327,328,331,334,335,336 FAD Scalar 13 R(4) 4 2 17955 COM FCP Scalar 10 R(4) 4 2 17955 COM FDIV Scalar 10 R(4) 4 2 17955 COM FIRST Scalar 5 L(4) 4 scalar COM FIS Scalar 7 R(4) 4 2 17955 COM 152 G Param 67 R(4) 4 scalar 67,332 G2LI Scalar 55 I(4) 4 1 239 COM Page 13 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References G2LJ Scalar 55 I(4) 4 1 389 COM GFFC Scalar 19 R(4) 4 2 17955 COM GH Local 125 R(4) 4 1 44 260,268,274 GLAT Scalar 15 R(4) 4 2 17955 COM GLB_TABLE Common 37 128 GLON Scalar 15 R(4) 4 2 17955 COM GM Local 125 R(4) 4 1 44 260,268,274 HBM2 Scalar 8 R(4) 4 2 17955 COM 180,189,191,438,439,440,441,442,44 3,465,466,467,469,479,480,481,560 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 HPBL Local 258 R(4) 4 scalar 258,301 HTM Scalar 6 R(4) 4 3 807975 COM 166,469 HTMG Scalar 52 R(4) 4 3 4183695 COM HTOP Scalar 17 R(4) 4 2 17955 COM I Local 150 I(4) 4 scalar 150,151,152,160,161,164,165,166,16 8,179,180,187,188,189,190,191,206, 207,208,211,212,213,215,216,217,22 0,221,252,254,258,259,260,268,274, 275,282,287,288,289,290,295,296,29 8,299,300,301,302,303,304,305,306, 307,310,312,318,325,328,331,332,33 3,334,335,336,386,387,388,409,410, 411,437,438,439,440,441,442,443,46 4,465,466,467,468,469,476,477,478, 479,480,481,482,489,490,491,513,51 5,516,521,522,523,524,531,532,533, 559,560 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 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,20,21,22,23,26,71,95,96, 97,98,99,100,101,102,105,106,107,1 08,109,110,113,114,115,116,117,118 ,119,120,121,122,449,450 IDIM2 Param 19 I(4) 4 scalar 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,20,21,22,23,26,71,95,96, 97,98,99,100,101,102,105,106,107,1 08,109,110,113,114,115,116,117,118 ,119,120,121,122 IDTAD Scalar 6 I(4) 4 scalar COM IEXCH Subr 5 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 188,190,207,211,215 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 IHRST Scalar 5 I(4) 4 scalar COM IHW Scalar 3 I(4) 4 1 399 COM 188,190,207,212,216 Page 14 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References IHWG Scalar 9 I(4) 4 1 387 COM II Local 388 I(4) 4 scalar 388,389,390,391,392,393,394,395,39 6,397,398,411,412,413,414,415,416, 417,418,419,420,421 IIM Local 359 I(4) 4 scalar 359,368,369,370,371,372,373,374,37 5,376,377 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 9,16,19,43,44,45,46,47,51,52,57,69 ,359,364,386,409 IMJM Param 69 I(4) 4 scalar INDX Common 2 6564 INDXG Common 8 8084 INIT_TIM Local 132 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 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 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 70 I(4) 4 scalar 10,25 ITBQ Param 70 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 IVE Scalar 4 I(4) 4 1 399 COM 438,441,465,479,490 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 439,442,466,480,490 IVWG Scalar 9 I(4) 4 1 387 COM J Local 148 I(4) 4 scalar 148,151,152,158,161,164,165,166,16 8,178,180,186,188,189,190,191,205, 207,208,211,212,213,215,216,217,22 0,221,251,254,258,259,260,268,274, 275,282,287,288,289,290,295,296,29 8,299,300,301,302,303,304,305,306, 307,310,312,318,325,328,331,332,33 3,334,335,336,365,366,367,385,389, 390,391,392,393,394,395,396,397,39 8,407,408,436,438,439,440,441,442, 443,463,465,466,467,468,469,475,47 Page 15 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References 7,478,479,480,481,482,488,490,491, 512,515,516,521,522,523,524,531,53 2,533,558,560 JAM Param 69 I(4) 4 scalar 3,6 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 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,20,21,22,23,26,71,95,96, 97,98,99,100,101,102,105,106,107,1 08,109,110,113,114,115,116,117,118 ,119,120,121,122,449,450 JDIM2 Param 20 I(4) 4 scalar 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,20,21,22,23,26,71,95,96, 97,98,99,100,101,102,105,106,107,1 08,109,110,113,114,115,116,117,118 ,119,120,121,122 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 JJ Local 367 I(4) 4 scalar 367,368,369,370,371,372,373,374,37 5,376,377,408,412,413,414,415,416, 417,418,419,420,421 JJM Local 360 I(4) 4 scalar 360,412,413,414,415,416,417,418,41 9,420,421 JM Param 3 I(4) 4 scalar 9,17,20,43,44,45,46,47,51,52,57,69 ,360,365,406,407 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 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 70 I(4) 4 scalar 10,25 JTBQ Param 70 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 KTM Scalar 5 I(4) 4 scalar COM KTMQ2 Param 67 I(4) 4 scalar 282,521,531 L Local 156 I(4) 4 scalar 156,162,164,165,166,168,177,180,20 3,207,208,211,212,213,215,216,217, 220,221,462,465,466,467,468,469,48 7,490,491,557,560 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LIST Scalar 5 I(4) 4 scalar COM LM Param 3 I(4) 4 scalar 4,5,6,7,9,21,52,69,97,109,110,115, 118,119,120,121,122,156,177,203,23 1,233,235,237,241,243,349,502,504, 544,546,548,550,557 Page 16 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References LM1 Param 69 I(4) 4 scalar 98,99,105,106,113,114,125,137,139, 140,351,353,458,462,506,508 LMH Scalar 4 I(4) 4 2 17955 COM 254,515 LMHK Local 254 I(4) 4 scalar 254,255,256,258,268,274,282,295,29 6,298,307,311,312,313,317,318,319, 331,332,334,335,336,515,521 LMHM Local 256 I(4) 4 scalar 256 LMHP Local 255 I(4) 4 scalar 255 LMV Scalar 4 I(4) 4 2 17955 COM 516 LMVK Local 516 I(4) 4 scalar 516,531 LOOPS Common 2 158840 LP1 Param 69 I(4) 4 scalar 4,9,100,108,116,117,138,151,152,23 9,458,487,500 LPBL Local 258 I(4) 4 scalar 258 LSM Param 3 I(4) 4 scalar MAPPINGS Common 56 5024 MASKS Common 2 6894720 MIXLEN Subr 258 258 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 MPP_TIM Local 132 R(4) 4 scalar MXSNAL Scalar 2 R(4) 4 2 17955 COM MYIE Scalar 4 I(4) 4 scalar COM 437,464,476,489,559 MYIE1 Scalar 4 I(4) 4 scalar COM 513 MYIE1_P1 Scalar 8 I(4) 4 scalar COM 252 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 150,160,179,187,206 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 437,464,476,489,513,559 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 150,160,179,187,206,252 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 558 MYJE1 Scalar 15 I(4) 4 scalar COM MYJE1_P1 Scalar 17 I(4) 4 scalar COM 158 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 Page 17 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References MYJE2 Scalar 15 I(4) 4 scalar COM 436,463,475,488,512 MYJE2_P1 Scalar 18 I(4) 4 scalar COM 186,205,251 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 148,178 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 558 MYJS1 Scalar 10 I(4) 4 scalar COM MYJS1_P1 Scalar 12 I(4) 4 scalar COM 158 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 436,463,475,488,512 MYJS2_P1 Scalar 13 I(4) 4 scalar COM 186,205,251 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 148,178 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 364,387,410 MY_IE_LOC Scalar 3 I(4) 4 scalar COM MY_IS_GLB Scalar 2 I(4) 4 scalar COM 359,387,388,410,411 MY_IS_LOC Scalar 3 I(4) 4 scalar COM MY_JE_GLB Scalar 2 I(4) 4 scalar COM 366,406 MY_JE_LOC Scalar 3 I(4) 4 scalar COM MY_JS_GLB Scalar 2 I(4) 4 scalar COM 360,366,367,384,408 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 Page 18 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 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 NHB_TIM Scalar 132 R(4) 4 scalar COM NHRZ Param 71 I(4) 4 scalar 231,233,235,237,239,241,243,349,35 1,353,500,502,504,506,508,544,546, 548,550 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 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 P30 Local 332 R(4) 4 scalar 332,333 PD Scalar 3 R(4) 4 2 17955 COM 161,306,310,332 PDSL Local 161 R(4) 4 scalar 161,162,165,310,311,313,317,319 PHYS Common 4 1788916 PL Scalar 7 R(4) 4 scalar COM PLQ Scalar 8 R(4) 4 scalar COM PREC Scalar 12 R(4) 4 2 17955 COM PRODQ2 Subr 268 268 PSHLTR Scalar 19 R(4) 4 2 17955 COM PT Scalar 3 R(4) 4 scalar COM 162,306,313,319,332 PT2 Scalar 5 R(4) 4 scalar COM PTBL Scalar 25 R(4) 4 2 10184 COM PVRBLS Common 2 5530140 Q Scalar 6 R(4) 4 3 807975 COM 165,233,312,318,334,544 Q10 Scalar 16 R(4) 4 2 17955 COM 303,369,390,413 Q2 Scalar 21 R(4) 4 3 807975 COM 180,237,349,560 Q2COL_T Local 122 R(4) 4 3 807975 237,259,268,274,282,349 Q30 Scalar 14 R(4) 4 2 17955 COM 304,334,375,396,419 QCOL_T Local 121 R(4) 4 3 807975 233,259,305,523,544 QS Scalar 6 R(4) 4 2 17955 COM 298,478 QS0 Scalar 10 R(4) 4 1 134 COM QSHLTR Scalar 18 R(4) 4 2 17955 COM 303,373,394,417 QWBS Scalar 9 R(4) 4 2 17955 COM QZ0 Scalar 5 R(4) 4 2 17955 COM 299,478,521 R Scalar 3 R(4) 4 scalar COM RADIN Scalar 14 R(4) 4 2 17955 COM RADOT Scalar 14 R(4) 4 2 17955 COM RDETA Scalar 4 R(4) 4 1 45 COM RDP Scalar 7 R(4) 4 scalar COM RDPQ Scalar 8 R(4) 4 scalar COM Page 19 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 RES Scalar 7 R(4) 4 2 17955 COM 161,310 RESTRT Scalar 5 L(4) 4 scalar COM RES_TIM Scalar 133 R(4) 4 scalar COM RF Scalar 8 R(4) 4 2 17955 COM RG Param 67 R(4) 4 scalar 152 ROG Param 67 R(4) 4 scalar 165,312,318 ROI Scalar 6 R(4) 4 scalar COM ROS Scalar 6 R(4) 4 scalar COM RUN Scalar 5 L(4) 4 scalar COM RWMSK Local 210 R(4) 4 scalar 210,214,218 SFCDIF Subr 298 298 SGETMO Subr 231 231,233,235,237,239,241,243,349,35 1,353,500,502,504,506,508,544,546, 548,550 SI Scalar 10 R(4) 4 2 17955 COM SICE Scalar 4 R(4) 4 2 17955 COM SIGMA Scalar 5 L(4) 4 scalar COM 152 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 298,477,478 SNO Scalar 10 R(4) 4 2 17955 COM SQS Scalar 10 R(4) 4 1 134 COM 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 SURFCE Subr 449 449 SURFCE_TIM Scalar 133 R(4) 4 scalar COM T Scalar 5 R(4) 4 3 807975 COM 164,231,307,312,318,331,332,546 T30 Local 331 R(4) 4 scalar 331,333 TCOL_T Local 120 R(4) 4 3 807975 231,259,305,523,546 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 Page 20 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References TG Scalar 19 R(4) 4 2 17955 COM TH10 Scalar 16 R(4) 4 2 17955 COM 302,368,389,412 TH30 Scalar 14 R(4) 4 2 17955 COM 304,333,374,395,418 THE0 Scalar 10 R(4) 4 1 76 COM THE0Q Scalar 11 R(4) 4 1 152 COM THL Scalar 7 R(4) 4 scalar COM THS Scalar 6 R(4) 4 2 17955 COM 298,477 THZ0 Scalar 5 R(4) 4 2 17955 COM 299,477,521 TIMEF Func 131 R(8) 8 scalar 429,431,457,459 TIMING Common 133 16 TOPO Common 50 17478548 TSHLTR Scalar 18 R(4) 4 2 17955 COM 302,372,393,416 TTBL Scalar 25 R(4) 4 2 10184 COM TTBLQ Scalar 25 R(4) 4 2 66880 COM TTVG Scalar 51 R(4) 4 2 92971 COM TURBL Subr 1 TWBS Scalar 9 R(4) 4 2 17955 COM U Scalar 5 R(4) 4 3 807975 COM 211,212,213,502,548 U10 Scalar 17 R(4) 4 2 17955 COM 302,370,391,414 U30 Scalar 15 R(4) 4 2 17955 COM 304,335,376,397,420 UCOL Local 109 R(4) 4 3 807975 211,220,241,295,335 UCOL_T Local 118 R(4) 4 3 807975 241,258,502,532,548 ULM Local 295 R(4) 4 scalar 295,305 USTAR Scalar 3 R(4) 4 2 17955 COM 268,300 UZ0 Scalar 4 R(4) 4 2 17955 COM 188,189,438,531 UZ0H Local 101 R(4) 4 2 17955 141,188,299,430,438,439,440 V Scalar 6 R(4) 4 3 807975 COM 215,216,217,504,550 V10 Scalar 17 R(4) 4 2 17955 COM 302,371,392,415 V30 Scalar 15 R(4) 4 2 17955 COM 304,336,377,398,421 VBM2 Scalar 3 R(4) 4 2 17955 COM 468,482 VBM3 Scalar 3 R(4) 4 2 17955 COM VCOL Local 110 R(4) 4 3 807975 215,221,243,296,336 VCOL_T Local 119 R(4) 4 3 807975 243,258,504,532,550 VDIFH Subr 521 521 VDIFQ Subr 282 282 VDIFV Subr 531 531 VLM Local 296 R(4) 4 scalar 296,305 VRBLS Common 2 13143060 VTM Scalar 7 R(4) 4 3 807975 COM 207,208,211,212,213,215,216,217,46 8 VZ0 Scalar 4 R(4) 4 2 17955 COM 190,191,441,531 VZ0H Local 102 R(4) 4 2 17955 142,190,299,430,441,442,443 WMSK Local 207 R(4) 4 scalar 207,209,210 WPDAR Scalar 8 R(4) 4 2 17955 COM WSTAR Local 300 R(4) 4 scalar 300 Z0 Scalar 3 R(4) 4 2 17955 COM 301 Z0EFFT Common 2 287280 ZCOL Local 108 R(4) 4 3 825930 490,500 ZCOL_T Local 116 R(4) 4 3 825930 239,260,275,282,306,524 ZCOL_T2 Local 117 R(4) 4 3 825930 500,533 ZEFF Local 125 R(4) 4 1 4 287,288,289,290,301 ZEFFIJ Scalar 3 R(4) 4 3 71820 COM 287,288,289,290 ZERO2 Subr 141 141,142 ZERO3 Subr 137 137,138 ZERO3_T Subr 139 139,140 Page 21 Source Listing TURBL 2014-12-17 20:47 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References ZINT Local 100 R(4) 4 3 825930 138,151,152,164,166,239,449,458,49 0,491 Page 22 Source Listing TURBL 2014-12-17 20:47 Subprograms/Common Blocks TURBL.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CTLBLK Common 4 112 DYNAM Common 2 725404 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 PVRBLS Common 2 5530140 TEMPCOM Common 42 6603768 TIMING Common 133 16 TOPO Common 50 17478548 TURBL Subr 1 VRBLS Common 2 13143060 Z0EFFT Common 2 287280 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 23 Source Listing TURBL 2014-12-17 20:47 TURBL.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 24 Source Listing TURBL 2014-12-17 20:47 TURBL.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 : TURBL.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100