Page 1 Source Listing TURBL 2025-03-12 18:22 /tmp/ifortRYJu7d.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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 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 2025-03-12 18:22 Entry Points TURBL.F ENTRY POINTS Name turbl_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 384 369 100 Label 553 464,465 125 Label 657 649,650 200 Label 749 725,726 40 Label 396 390 60 Label 438 416 AKH Local 311 R(4) 4 3 790020 564,682 AKHCOL Local 319 R(4) 4 3 790020 682,719 AKH_T Local 326 R(4) 4 3 790020 352,488,564,719,736 AKM Local 312 R(4) 4 3 790020 350,566,671,678,679,680 AKMCOL Local 318 R(4) 4 3 790020 678,721 AKMSV Local 320 R(4) 4 2 17955 692,745 AKM_T Local 327 R(4) 4 3 790020 353,488,566,721,746 AMAX1 Func 393 scalar 393,773 APE Local 310 R(4) 4 3 807975 381,448,662 APECOL_T Local 328 R(4) 4 3 807975 448,473,519,737 APESTS Local 375 R(4) 4 scalar 375,377,381 BTIM Local 642 R(4) 4 scalar 642,644,670,672 CAPA Param 158 R(4) 4 scalar 381 CKLQ Local 308 R(4) 4 2 17955 663,735 CT Local 309 R(4) 4 2 17955 514,735 CTLBLK Common 169 112 SAVE DIFCOF Subr 487 487 DPLM Local 524 R(4) 4 scalar 524,525 DPLM1 Local 530 R(4) 4 scalar 530,531 DYNAM Common 198 725404 SAVE DZLM Local 525 R(4) 4 scalar 525,528,534,535 DZLM1 Local 531 R(4) 4 scalar 531,534,535 EL Local 338 R(4) 4 1 44 473,481,487,495 EPSQ2 Param 159 R(4) 4 scalar 393,773 EPSZ Param 159 R(4) 4 scalar 364 EXCH Local 64 scalar 643,671 EXCH0 Subr 65 EXCH00 Subr 65 643 EXCH0001111 Subr 65 EXCH01 Subr 65 EXCH011 Subr 65 671 EXCH1 Subr 66 EXCH11 Subr 66 EXCH111 Subr 66 EXCH1111 Subr 66 EXCH11111 Subr 66 EXCH111111 Subr 67 EXCHM Module 63 63 Page 12 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References EXP Func 545 scalar 545 FAC1 Local 534 R(4) 4 scalar 534,537,538,544,547,548,549 FAC2 Local 535 R(4) 4 scalar 535,540,541,544,547,548,549 G Param 158 R(4) 4 scalar 158,545 GH Local 338 R(4) 4 1 44 473,481,487 GLB_TABLE Common 135 128 SAVE GM Local 338 R(4) 4 1 44 473,481,487 HPBL Local 471 R(4) 4 scalar 471,514 I Local 363 I(4) 4 scalar 363,364,365,373,374,377,378,379,38 1,392,393,400,401,402,403,404,419, 420,421,424,425,426,428,429,430,43 3,434,465,467,471,472,473,481,487, 488,495,500,501,502,503,508,509,51 1,512,513,514,515,516,517,518,519, 520,523,525,531,538,541,544,545,54 6,547,548,549,599,600,601,622,623, 624,650,651,652,653,654,655,656,67 7,678,679,680,681,682,689,690,691, 692,693,694,695,702,703,704,726,72 8,729,734,735,736,737,744,745,746, 772,773 IDIM1 Param 89 I(4) 4 scalar 162,180,187,188,190,191,192,204,20 5,206,207,208,209,226,227,228,229, 230,231,232,233,234,235,236,245,24 7,248,249,256,257,258,259,260,261, 262,263,264,265,266,267,268,269,27 0,271,272,274,279,289,300,308,309, 310,311,312,313,314,315,318,319,32 0,321,322,323,326,327,328,329,330, 331,332,333,334,335,662,663 IDIM2 Param 89 I(4) 4 scalar 162,180,187,188,190,191,192,204,20 5,206,207,208,209,226,227,228,229, 230,231,232,233,234,235,236,245,24 7,248,249,256,257,258,259,260,261, 262,263,264,265,266,267,268,269,27 0,271,272,274,279,289,300,308,309, 310,311,312,313,314,315,318,319,32 0,321,322,323,326,327,328,329,330, 331,332,333,334,335 IEXCH Subr 67 IGSTL Param 84 I(4) 4 scalar 89 IGSTR Param 84 I(4) 4 scalar 89 II Local 601 I(4) 4 scalar 601,602,603,604,605,606,607,608,60 9,610,611,624,625,626,627,628,629, 630,631,632,633,634 IIM Local 572 I(4) 4 scalar 572,581,582,583,584,585,586,587,58 8,589,590 IM Param 73 I(4) 4 scalar 86,89,141,142,143,144,145,149,150, 155,160,293,572,577,599,622 IMJM Param 160 I(4) 4 scalar INDX Common 286 6564 SAVE INDXG Common 292 8084 SAVE INIT_TIM Local 345 R(4) 4 scalar INPES Param 81 I(4) 4 scalar 86,89,125,126,127,136,137 ITAIL Param 86 I(4) 4 scalar Page 13 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References ITB Param 161 I(4) 4 scalar 223,238 ITBQ Param 161 I(4) 4 scalar 224,238 J Local 361 I(4) 4 scalar 361,364,365,371,374,377,378,379,38 1,391,393,399,401,402,403,404,418, 420,421,424,425,426,428,429,430,43 3,434,464,467,471,472,473,481,487, 488,495,500,501,502,503,508,509,51 1,512,513,514,515,516,517,518,519, 520,523,525,531,538,541,544,545,54 6,547,548,549,578,579,580,598,602, 603,604,605,606,607,608,609,610,61 1,620,621,649,651,652,653,654,655, 656,676,678,679,680,681,682,688,69 0,691,692,693,694,695,701,703,704, 725,728,729,734,735,736,737,744,74 5,746,771,773 JAM Param 160 I(4) 4 scalar 179,202 JDIM1 Param 90 I(4) 4 scalar 162,180,187,188,190,191,192,204,20 5,206,207,208,209,226,227,228,229, 230,231,232,233,234,235,236,245,24 7,248,249,256,257,258,259,260,261, 262,263,264,265,266,267,268,269,27 0,271,272,274,279,287,288,300,308, 309,310,311,312,313,314,315,318,31 9,320,321,322,323,326,327,328,329, 330,331,332,333,334,335,662,663 JDIM2 Param 90 I(4) 4 scalar 162,180,187,188,190,191,192,204,20 5,206,207,208,209,226,227,228,229, 230,231,232,233,234,235,236,245,24 7,248,249,256,257,258,259,260,261, 262,263,264,265,266,267,268,269,27 0,271,272,274,279,287,288,300,308, 309,310,311,312,313,314,315,318,31 9,320,321,322,323,326,327,328,329, 330,331,332,333,334,335 JGSTL Param 85 I(4) 4 scalar 90 JGSTR Param 85 I(4) 4 scalar 90 JJ Local 580 I(4) 4 scalar 580,581,582,583,584,585,586,587,58 8,589,590,621,625,626,627,628,629, 630,631,632,633,634 JJM Local 573 I(4) 4 scalar 573,625,626,627,628,629,630,631,63 2,633,634 JM Param 73 I(4) 4 scalar 87,90,141,142,143,144,145,149,150, 155,160,293,573,578,619,620 JNPES Param 81 I(4) 4 scalar 87,90,125,126,127,136,137 JTAIL Param 87 I(4) 4 scalar JTB Param 161 I(4) 4 scalar 223,238 JTBQ Param 161 I(4) 4 scalar 238 KTMQ2 Param 158 I(4) 4 scalar 495,734,744 L Local 369 I(4) 4 scalar 369,375,377,378,379,381,390,393,41 6,420,421,424,425,426,428,429,430, 433,434,675,678,679,680,681,682,70 0,703,704,770,773 LM Param 73 I(4) 4 scalar 150,160,190,191,200,222,247,248,27 4,310,322,323,328,331,332,333,334, Page 14 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References 335,369,390,416,444,446,448,450,45 4,456,562,715,717,757,759,761,763, 770 LM1 Param 160 I(4) 4 scalar 311,312,318,319,326,327,338,350,35 2,353,564,566,671,675,719,721 LMHK Local 467 I(4) 4 scalar 467,468,469,471,481,487,495,508,50 9,511,520,524,525,526,530,531,532, 544,545,547,548,549,728,734 LMHM Local 469 I(4) 4 scalar 469 LMHP Local 468 I(4) 4 scalar 468 LMVK Local 729 I(4) 4 scalar 729,744 LOOPS Common 178 158840 SAVE LP1 Param 160 I(4) 4 scalar 200,222,313,321,329,330,351,364,36 5,452,671,700,713 LPBL Local 471 I(4) 4 scalar 471 LSM Param 73 I(4) 4 scalar MAPPINGS Common 154 5024 SAVE MASKS Common 186 6894720 SAVE MIXLEN Subr 471 471 MPPCOM Common 99 1464 SAVE MPP_TIM Local 345 R(4) 4 scalar NHRZ Param 162 I(4) 4 scalar 444,446,448,450,452,454,456,562,56 4,566,713,715,717,719,721,757,759, 761,763 P30 Local 545 R(4) 4 scalar 545,546 PDSL Local 374 R(4) 4 scalar 374,375,378,523,524,526,530,532 PHYS Common 217 1788916 SAVE PRODQ2 Subr 481 481 PVRBLS Common 255 5530140 SAVE Q2COL_T Local 335 R(4) 4 3 807975 450,472,481,487,495,562 QCOL_T Local 334 R(4) 4 3 807975 446,472,518,736,757 RG Param 158 R(4) 4 scalar 365 ROG Param 158 R(4) 4 scalar 378,525,531 RWMSK Local 423 R(4) 4 scalar 423,427,431 SFCDIF Subr 511 511 SGETMO Subr 444 444,446,448,450,452,454,456,562,56 4,566,713,715,717,719,721,757,759, 761,763 SINGLRST Local 167 L(4) 4 scalar SURFCE Subr 662 662 T30 Local 544 R(4) 4 scalar 544,546 TCOL_T Local 333 R(4) 4 3 807975 444,472,518,736,759 TEMPCOM Common 140 6603768 SAVE TIMEF Func 344 R(8) 8 scalar 642,644,670,672 TIMING Common 346 16 SAVE TOPO Common 148 17478548 SAVE TURBL Subr 2 UCOL Local 322 R(4) 4 3 807975 424,433,454,508,548 UCOL_T Local 331 R(4) 4 3 807975 454,471,715,745,761 ULM Local 508 R(4) 4 scalar 508,518 UZ0H Local 314 R(4) 4 2 17955 354,401,512,643,651,652,653 VCOL Local 323 R(4) 4 3 807975 428,434,456,509,549 VCOL_T Local 332 R(4) 4 3 807975 456,471,717,745,763 VDIFH Subr 734 734 VDIFQ Subr 495 495 Page 15 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Object Declared Type Bytes Dimen Elements Attributes References VDIFV Subr 744 744 VLM Local 509 R(4) 4 scalar 509,518 VRBLS Common 244 13143060 SAVE VZ0H Local 315 R(4) 4 2 17955 355,403,512,643,654,655,656 WMSK Local 420 R(4) 4 scalar 420,422,423 WSTAR Local 513 R(4) 4 scalar 513 Z0EFFT Common 299 287280 SAVE ZCOL Local 321 R(4) 4 3 825930 703,713 ZCOL_T Local 329 R(4) 4 3 825930 452,473,488,495,519,737 ZCOL_T2 Local 330 R(4) 4 3 825930 713,746 ZEFF Local 338 R(4) 4 1 4 500,501,502,503,514 ZERO2 Subr 354 354,355 ZERO3 Subr 350 350,351 ZERO3_T Subr 352 352,353 ZINT Local 313 R(4) 4 3 825930 351,364,365,377,379,452,662,671,70 3,704 TYPE COMPONENTS/COMMON VARIABLES Name Type Bytes Offset Dimen Elements Attributes References ACCLIQ R(4) 4 1292760 2 17955 COM ACPREC R(4) 4 1220940 2 17955 COM AETA R(4) 4 396 1 45 COM 375,526,532 AETA2 R(4) 4 264 1 45 COM AFSI R(4) 4 5458320 2 17955 COM AKHS R(4) 4 646380 2 17955 COM 514,734 AKMS R(4) 4 574560 2 17955 COM 514,671,692,693,694 ALBASE R(4) 4 1009004 2 17955 COM ALBEDO R(4) 4 1080824 2 17955 COM CI R(4) 4 40 scalar COM CLDEFI R(4) 4 1077300 2 17955 COM CNVBOT R(4) 4 721724 2 17955 COM CNVTOP R(4) 4 649904 2 17955 COM CPGFU R(4) 4 150844 2 17955 COM CPGFV R(4) 4 4 scalar COM CS R(4) 4 28 scalar COM CUPREC R(4) 4 1364580 2 17955 COM CURV R(4) 4 222664 2 17955 COM CZEN R(4) 4 434444 2 17955 COM CZMEAN R(4) 4 1296284 2 17955 COM DDMPU R(4) 4 509944 2 17955 COM DDMPV R(4) 4 581764 2 17955 COM DETA R(4) 4 36 1 45 COM 378,524,530 DETA2 R(4) 4 84 1 45 COM DFL R(4) 4 940 1 46 COM DFRLG R(4) 4 444 1 46 COM 379 DI R(4) 4 44 scalar COM DS R(4) 4 32 scalar COM DT R(4) 4 48 scalar COM DTD R(4) 4 16 scalar COM DTQ2 R(4) 4 8 scalar COM 481,495,734,744 DX R(4) 4 7204 2 17955 COM DY R(4) 4 0 scalar COM Page 16 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Type Bytes Offset Dimen Elements Attributes References EF4T R(4) 4 32 scalar COM EM R(4) 4 1124 1 760 COM EMT R(4) 4 4164 1 760 COM EN R(4) 4 8 scalar COM ENT R(4) 4 12 scalar COM EPSR R(4) 4 75344 2 17955 COM ETA R(4) 4 756 1 46 COM EXCH_TIM R(4) 4 12 scalar COM 644,672 F R(4) 4 438124 2 17955 COM F4D R(4) 4 24 scalar COM F4Q R(4) 4 28 scalar COM F4Q2 R(4) 4 576 1 45 COM FAD R(4) 4 653584 2 17955 COM FCP R(4) 4 294484 2 17955 COM FDIV R(4) 4 366304 2 17955 COM FIRST L(4) 4 24 scalar COM FIS R(4) 4 12999420 2 17955 COM 365 G2LI I(4) 4 0 1 239 COM G2LJ I(4) 4 1912 1 389 COM GFFC R(4) 4 865364 2 17955 COM GLAT R(4) 4 290804 2 17955 COM GLON R(4) 4 362624 2 17955 COM HBM2 R(4) 4 6751080 2 17955 COM 393,402,404,651,652,653,654,655,65 6,678,679,680,682,692,693,694,773 HBM3 R(4) 4 6822900 2 17955 COM HBOT R(4) 4 578084 2 17955 COM HDAC R(4) 4 1152644 2 17955 COM HDACV R(4) 4 1224464 2 17955 COM HTM R(4) 4 287280 3 807975 COM 379,682 HTMG R(4) 4 743768 3 4183695 COM HTOP R(4) 4 506264 2 17955 COM IBROW I(4) 4 404 scalar COM ICHUNKTAB I(4) 4 620 1 8 COM IDAT I(4) 4 8 1 3 COM IDTAD I(4) 4 72 scalar COM IE_GLB_TABLE I(4) 4 32 1 8 COM IE_LOC_TABLE I(4) 4 556 1 8 COM IHE I(4) 4 0 1 399 COM 401,403,420,424,428 IHEG I(4) 4 0 1 387 COM IHHA I(4) 4 3040 1 760 COM IHLA I(4) 4 0 1 760 COM IHRST I(4) 4 20 scalar COM IHW I(4) 4 1596 1 399 COM 401,403,420,425,429 IHWG I(4) 4 1548 1 387 COM ILCOL I(4) 4 396 scalar COM ILPAD1 I(4) 4 412 scalar COM ILPAD2 I(4) 4 416 scalar COM ILPAD3 I(4) 4 420 scalar COM ILPAD4 I(4) 4 424 scalar COM ILPAD5 I(4) 4 428 scalar COM INUMQ I(4) 4 1064 1 100 COM IOUT I(4) 4 44 scalar COM IQUILT_GROUP I(4) 4 1060 scalar COM IRAD I(4) 4 6384 1 45 COM IRADG I(4) 4 6192 1 473 COM Page 17 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Type Bytes Offset Dimen Elements Attributes References IRCOL I(4) 4 400 scalar COM IRPAD1 I(4) 4 432 scalar COM IRPAD2 I(4) 4 436 scalar COM IRPAD3 I(4) 4 440 scalar COM IRPAD4 I(4) 4 444 scalar COM IRPAD5 I(4) 4 448 scalar COM IS_GLB_TABLE I(4) 4 0 1 8 COM IS_LOC_TABLE I(4) 4 492 1 8 COM ITEMP I(4) 4 5870016 2 91719 COM ITEMP2 I(4) 4 6236892 2 91719 COM ITROW I(4) 4 408 scalar COM IVE I(4) 4 3192 1 399 COM 651,654,678,692,703 IVEG I(4) 4 3096 1 387 COM IVHA I(4) 4 9120 1 760 COM IVLA I(4) 4 6080 1 760 COM IVW I(4) 4 4788 1 399 COM 652,655,679,693,703 IVWG I(4) 4 4644 1 387 COM JBPAD1 I(4) 4 452 scalar COM JBPAD2 I(4) 4 456 scalar COM JBPAD3 I(4) 4 460 scalar COM JBPAD4 I(4) 4 464 scalar COM JBPAD5 I(4) 4 468 scalar COM JE_GLB_TABLE I(4) 4 96 1 8 COM JE_LOC_TABLE I(4) 4 588 1 8 COM JRA I(4) 4 12160 1 760 COM JS_GLB_TABLE I(4) 4 64 1 8 COM JS_LOC_TABLE I(4) 4 524 1 8 COM JTPAD1 I(4) 4 472 scalar COM JTPAD2 I(4) 4 476 scalar COM JTPAD3 I(4) 4 480 scalar COM JTPAD4 I(4) 4 484 scalar COM JTPAD5 I(4) 4 488 scalar COM KTM I(4) 4 0 scalar COM L2GI I(4) 4 956 1 239 COM L2GJ I(4) 4 3468 1 389 COM LIST I(4) 4 40 scalar COM LMH I(4) 4 15200 2 17955 COM 467,728 LMV I(4) 4 87020 2 17955 COM 729 MPI_COMM_COMP I(4) 4 652 scalar COM MPI_COMM_INTER I(4) 4 656 scalar COM MPI_COMM_INTER_ARRAY I(4) 4 660 1 100 COM MXSNAL R(4) 4 3524 2 17955 COM MYIE I(4) 4 52 scalar COM 650,677,689,702,772 MYIE1 I(4) 4 56 scalar COM 726 MYIE1_P1 I(4) 4 120 scalar COM 465 MYIE1_P2 I(4) 4 124 scalar COM MYIE1_P3 I(4) 4 128 scalar COM MYIE1_P4 I(4) 4 132 scalar COM MYIE2 I(4) 4 60 scalar COM MYIE2_P1 I(4) 4 136 scalar COM MYIE_P1 I(4) 4 100 scalar COM 363,373,392,400,419 MYIE_P2 I(4) 4 104 scalar COM MYIE_P3 I(4) 4 108 scalar COM MYIE_P4 I(4) 4 112 scalar COM MYIE_P5 I(4) 4 116 scalar COM Page 18 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Type Bytes Offset Dimen Elements Attributes References MYIS I(4) 4 40 scalar COM 650,677,689,702,726,772 MYIS1 I(4) 4 44 scalar COM MYIS1_P1 I(4) 4 84 scalar COM MYIS1_P2 I(4) 4 88 scalar COM MYIS1_P3 I(4) 4 92 scalar COM MYIS1_P4 I(4) 4 96 scalar COM MYIS2 I(4) 4 48 scalar COM MYIS_P1 I(4) 4 64 scalar COM 363,373,392,400,419,465 MYIS_P2 I(4) 4 68 scalar COM MYIS_P3 I(4) 4 72 scalar COM MYIS_P4 I(4) 4 76 scalar COM MYIS_P5 I(4) 4 80 scalar COM MYJE I(4) 4 236 scalar COM 771 MYJE1 I(4) 4 240 scalar COM MYJE1_P1 I(4) 4 280 scalar COM 371 MYJE1_P2 I(4) 4 284 scalar COM MYJE1_P3 I(4) 4 288 scalar COM MYJE1_P4 I(4) 4 292 scalar COM MYJE2 I(4) 4 244 scalar COM 649,676,688,701,725 MYJE2_P1 I(4) 4 296 scalar COM 399,418,464 MYJE2_P2 I(4) 4 300 scalar COM MYJE2_P3 I(4) 4 304 scalar COM MYJE2_P4 I(4) 4 308 scalar COM MYJE3 I(4) 4 248 scalar COM MYJE3_P4 I(4) 4 312 scalar COM MYJE4 I(4) 4 252 scalar COM MYJE4_P1 I(4) 4 316 scalar COM MYJE4_P4 I(4) 4 320 scalar COM MYJE5 I(4) 4 256 scalar COM MYJE5_P1 I(4) 4 324 scalar COM MYJE5_P2 I(4) 4 328 scalar COM MYJE_P1 I(4) 4 260 scalar COM 361,391 MYJE_P2 I(4) 4 264 scalar COM MYJE_P3 I(4) 4 268 scalar COM MYJE_P4 I(4) 4 272 scalar COM MYJE_P5 I(4) 4 276 scalar COM MYJS I(4) 4 140 scalar COM 771 MYJS1 I(4) 4 144 scalar COM MYJS1_P1 I(4) 4 184 scalar COM 371 MYJS1_P2 I(4) 4 188 scalar COM MYJS1_P3 I(4) 4 192 scalar COM MYJS1_P4 I(4) 4 196 scalar COM MYJS2 I(4) 4 148 scalar COM 649,676,688,701,725 MYJS2_P1 I(4) 4 200 scalar COM 399,418,464 MYJS2_P2 I(4) 4 204 scalar COM MYJS2_P3 I(4) 4 208 scalar COM MYJS2_P4 I(4) 4 212 scalar COM MYJS3 I(4) 4 152 scalar COM MYJS3_P4 I(4) 4 216 scalar COM MYJS4 I(4) 4 156 scalar COM MYJS4_P1 I(4) 4 220 scalar COM MYJS4_P4 I(4) 4 224 scalar COM MYJS5 I(4) 4 160 scalar COM MYJS5_P1 I(4) 4 228 scalar COM MYJS5_P2 I(4) 4 232 scalar COM Page 19 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Type Bytes Offset Dimen Elements Attributes References MYJS_P1 I(4) 4 164 scalar COM 361,391 MYJS_P2 I(4) 4 168 scalar COM MYJS_P3 I(4) 4 172 scalar COM MYJS_P4 I(4) 4 176 scalar COM MYJS_P5 I(4) 4 180 scalar COM MYPE I(4) 4 0 scalar COM MY_E I(4) 4 336 scalar COM MY_IE_GLB I(4) 4 12 scalar COM 577,600,623 MY_IE_LOC I(4) 4 28 scalar COM MY_IS_GLB I(4) 4 8 scalar COM 572,600,601,623,624 MY_IS_LOC I(4) 4 24 scalar COM MY_JE_GLB I(4) 4 20 scalar COM 579,619 MY_JE_LOC I(4) 4 36 scalar COM MY_JS_GLB I(4) 4 16 scalar COM 573,579,580,597,621 MY_JS_LOC I(4) 4 32 scalar COM MY_N I(4) 4 332 scalar COM MY_NE I(4) 4 348 scalar COM MY_NEB I(4) 4 364 1 8 COM MY_NW I(4) 4 360 scalar COM MY_S I(4) 4 340 scalar COM MY_SE I(4) 4 352 scalar COM MY_SW I(4) 4 356 scalar COM MY_W I(4) 4 344 scalar COM NBC I(4) 4 36 scalar COM NBOCO I(4) 4 76 scalar COM NCNVC I(4) 4 92 scalar COM NCP I(4) 4 84 scalar COM NEST I(4) 4 108 scalar COM NFCST I(4) 4 32 scalar COM NHB_TIM R(4) 4 4 scalar COM NPES I(4) 4 4 scalar COM NPHS I(4) 4 88 scalar COM NPREC I(4) 4 68 scalar COM NRADL I(4) 4 100 scalar COM NRADS I(4) 4 96 scalar COM NSHDE I(4) 4 80 scalar COM NSTART I(4) 4 60 scalar COM NTDDMP I(4) 4 64 scalar COM NTSD I(4) 4 52 scalar COM NTSTM I(4) 4 56 scalar COM PD R(4) 4 0 2 17955 COM 374,519,523,545 PL R(4) 4 48 scalar COM PLQ R(4) 4 72 scalar COM PREC R(4) 4 1149120 2 17955 COM PSHLTR R(4) 4 2154600 2 17955 COM PT R(4) 4 20 scalar COM 375,519,526,532,545 PT2 R(4) 4 4 scalar COM PTBL R(4) 4 1439924 2 10184 COM Q R(4) 4 9767520 3 807975 COM 378,446,525,531,547,757 Q10 R(4) 4 1795500 2 17955 COM 516,582,603,626 Q2 R(4) 4 2226420 3 807975 COM 393,450,562,773 Q30 R(4) 4 1508220 2 17955 COM 517,547,588,609,632 QS R(4) 4 502740 2 17955 COM 511,691 QS0 R(4) 4 628 1 134 COM QSHLTR R(4) 4 2082780 2 17955 COM 516,586,607,630 Page 20 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Type Bytes Offset Dimen Elements Attributes References QWBS R(4) 4 861840 2 17955 COM QZ0 R(4) 4 359100 2 17955 COM 512,691,734 R R(4) 4 16 scalar COM RADIN R(4) 4 147164 2 17955 COM RADOT R(4) 4 218984 2 17955 COM RDETA R(4) 4 216 1 45 COM RDP R(4) 4 64 scalar COM RDPQ R(4) 4 76 scalar COM RDQ R(4) 4 56 scalar COM RDTH R(4) 4 60 scalar COM RDTHE R(4) 4 68 scalar COM RDTHEQ R(4) 4 80 scalar COM RES R(4) 4 13071240 2 17955 COM 374,523 RESTRT L(4) 4 28 scalar COM RES_TIM R(4) 4 8 scalar COM RF R(4) 4 718200 2 17955 COM ROI R(4) 4 36 scalar COM ROS R(4) 4 24 scalar COM RUN L(4) 4 4 scalar COM SI R(4) 4 1005480 2 17955 COM SICE R(4) 4 215460 2 17955 COM SIGMA L(4) 4 0 scalar COM 365 SIGT4 R(4) 4 1368104 2 17955 COM SM R(4) 4 143640 2 17955 COM 511,690,691 SNO R(4) 4 933660 2 17955 COM SQS R(4) 4 1164 1 134 COM SST R(4) 4 937184 2 17955 COM STHE R(4) 4 2004 1 76 COM STHEQ R(4) 4 2916 1 152 COM SUBPOST L(4) 4 104 scalar COM SURFCE_TIM R(4) 4 0 scalar COM T R(4) 4 71820 3 807975 COM 377,444,520,525,531,544,545,759 TDTD R(4) 4 20 scalar COM TDTQ2 R(4) 4 12 scalar COM TEMP1 R(4) 4 0 2 91719 COM TEMP10 R(4) 4 3301884 2 91719 COM TEMP11 R(4) 4 3668760 2 91719 COM TEMP12 R(4) 4 4035636 2 91719 COM TEMP13 R(4) 4 4402512 2 91719 COM TEMP14 R(4) 4 4769388 2 91719 COM TEMP15 R(4) 4 5136264 2 91719 COM TEMP16 R(4) 4 5503140 2 91719 COM TEMP2 R(4) 4 366876 2 91719 COM TEMP2X R(4) 4 0 2 92971 COM TEMP3 R(4) 4 733752 2 91719 COM TEMP4 R(4) 4 1100628 2 91719 COM TEMP5 R(4) 4 1467504 2 91719 COM TEMP6 R(4) 4 1834380 2 91719 COM TEMP7 R(4) 4 2201256 2 91719 COM TEMP8 R(4) 4 2568132 2 91719 COM TEMP9 R(4) 4 2935008 2 91719 COM TG R(4) 4 793544 2 17955 COM TH10 R(4) 4 1723680 2 17955 COM 515,581,602,625 TH30 R(4) 4 1436400 2 17955 COM 517,546,587,608,631 THE0 R(4) 4 1700 1 76 COM Page 21 Source Listing TURBL 2025-03-12 18:22 Symbol Table TURBL.F Name Type Bytes Offset Dimen Elements Attributes References THE0Q R(4) 4 2308 1 152 COM THL R(4) 4 52 scalar COM THS R(4) 4 430920 2 17955 COM 511,690 THZ0 R(4) 4 287280 2 17955 COM 512,690,734 TSHLTR R(4) 4 2010960 2 17955 COM 515,585,606,629 TTBL R(4) 4 1480660 2 10184 COM TTBLQ R(4) 4 1521396 2 66880 COM TTVG R(4) 4 371884 2 92971 COM TWBS R(4) 4 790020 2 17955 COM U R(4) 4 3303720 3 807975 COM 424,425,426,715,761 U10 R(4) 4 1867320 2 17955 COM 515,583,604,627 U30 R(4) 4 1580040 2 17955 COM 517,548,589,610,633 USTAR R(4) 4 71820 2 17955 COM 481,513 UZ0 R(4) 4 143640 2 17955 COM 401,402,651,744 V R(4) 4 6535620 3 807975 COM 428,429,430,717,763 V10 R(4) 4 1939140 2 17955 COM 515,584,605,628 V30 R(4) 4 1651860 2 17955 COM 517,549,590,611,634 VBM2 R(4) 4 0 2 17955 COM 681,695 VBM3 R(4) 4 71820 2 17955 COM VTM R(4) 4 3519180 3 807975 COM 420,421,424,425,426,428,429,430,68 1 VZ0 R(4) 4 215460 2 17955 COM 403,404,654,744 WPDAR R(4) 4 79024 2 17955 COM Z0 R(4) 4 0 2 17955 COM 514 ZEFFIJ R(4) 4 0 3 71820 COM 500,501,502,503 Page 22 Source Listing TURBL 2025-03-12 18:22 Subprograms/Common Blocks TURBL.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CTLBLK Common 169 112 SAVE DYNAM Common 198 725404 SAVE GLB_TABLE Common 135 128 SAVE INDX Common 286 6564 SAVE INDXG Common 292 8084 SAVE LOOPS Common 178 158840 SAVE MAPPINGS Common 154 5024 SAVE MASKS Common 186 6894720 SAVE MPPCOM Common 99 1464 SAVE PHYS Common 217 1788916 SAVE PVRBLS Common 255 5530140 SAVE TEMPCOM Common 140 6603768 SAVE TIMING Common 346 16 SAVE TOPO Common 148 17478548 SAVE TURBL Subr 2 VRBLS Common 244 13143060 SAVE Z0EFFT Common 299 287280 SAVE COMPILER OPTIONS BEING USED -align noall -align nonone -align nocommons -align nodcommons -align noqcommons -align nozcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -align norec32byte -align norec64byte -align noarray8byte -align noarray16byte -align noarray32byte -align noarray64byte -align noarray128byte -align noarray256byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume nobuffered_stdout -assume byterecl -assume nocontiguous_assumed_shape -assume nocontiguous_pointer -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_complex_align -assume old_unit_star -assume old_inquire_recl -assume old_ldout_format -assume old_ldout_zero -assume noold_logical_assign -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume noprotect_allocates -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume std_minus0_rounding -assume nostd_mod_proc_name -assume std_value -assume realloc_lhs -assume underscore -assume no2underscores -assume norecursion no -auto -auto_scalar no -bintext -ccdefault default Page 23 Source Listing TURBL 2025-03-12 18:22 TURBL.F -check noarg_temp_created -check noassume -check nobounds -check nocontiguous -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check noshape -check nostack -check nouninitialized -check noudio_iostat -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1910 -D __INTEL_COMPILER_UPDATE=3 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D __amd64 -D __amd64__ -D __INTEL_COMPILER_BUILD_DATE=20200925 -D __INTEL_OFFLOAD -D __MMX__ -D __SSE__ -D __SSE_MATH__ -D __SSE2__ -D __SSE2_MATH__ -D __SSE3__ -D __SSSE3__ -D __SSE4_1__ -D __SSE4_2__ -D __POPCNT__ -D __PCLMUL__ -D __AES__ -D __AVX__ -D __F16C__ -D __AVX_I__ -D __RDRND__ -D __FMA__ -D __FP_FAST_FMA -D __FP_FAST_FMAF -D __BMI__ -D __LZCNT__ -D __AVX2__ -D __haswell -D __haswell__ -D __tune_haswell__ -D __core_avx2 -D __core_avx2__ -D __tune_core_avx2__ -D __CRAY_X86_ROME -D __CRAYXT_COMPUTE_LINUX_TARGET -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 Page 24 Source Listing TURBL 2025-03-12 18:22 TURBL.F -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 -init noarrays -init nohuge -init noinfinity -init nominus_huge -init nominus_infinity -init nominus_tiny -init nonan -init nosnan -init notiny -init nozero no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude no -o -offload-build=host -openmp-simd -O2 no -pad_source -real_size 32 no -recursive -reentrancy threaded -vec=simd -show nofullpath -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w nodeclarations -w noexternals -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w noshape -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage no -wrap-margins -includepath : /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/, .f90,./.f90,/opt/cray/pe/mpich/8.1.12/ofi/intel/19.0/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/.f90,/usr/lib64/gcc/x86_64-suse-linux/7/include/.f90, /usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/.f90,/usr/include/.f90,/usr/include/.f90,/usr/include/.f90 -list filename : TURBL.lst no -o COMPILER: Intel(R) Fortran 19.1-1655