Page 1 Source Listing VTADV 2014-12-17 20:47 /tmpnwprd/ifortvfhTSS.i 1 # 1 "VTADV.F" 2 c&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 3 SUBROUTINE VTADV 4 C ****************************************************************** 5 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 6 C . . . 7 C SUBPROGRAM: VTADV VERTICAL ADVECTION 8 C PRGRMMR: JANJIC ORG: W/NP22 DATE: 93-11-17 9 C 10 C ABSTRACT: 11 C VTADV CALCULATES THE CONTRIBUTION OF THE VERTICAL ADVECTION 12 C TO THE TENDENCIES OF TEMPERATURE, SPECIFIC HUMIDITY, WIND 13 C COMPONENTS, AND TURBULENT KINETIC ENERGY AND THEN UPDATES THOSE 14 C VARIABLES. FOR ALL VARIABLES EXCEPT SPECIFIC HUMIDITY A 15 C SIMPLE CENTERED DIFFERENCE SCHEME IN SPACE IS USED IN 16 C CONJUNCTION WITH THE PURE EULER-BACKWARD TIME SCHEME. 17 C A PIECEWISE LINEAR SCHEME IS USED TO CALCULATE THE VERTICAL 18 C ADVECTION OF SPECIFIC HUMIDITY SO THAT NO FALSE MAXIMA OR 19 C MINIMA ARE PRODUCED. 20 C 21 C PROGRAM HISTORY LOG: 22 C 87-06-?? JANJIC - ORIGINATOR 23 C 90-??-?? MESINGER - INSERTED PIECEWISE LINEAR SCHEME FOR 24 C SPECIFIC HUMIDITY 25 C 95-03-25 BLACK - CONVERSION FROM 1-D TO 2-D IN HORIZONTAL 26 C 95-11-20 ABELES - PARALLEL OPTIMIZATION 27 C 96-03-29 BLACK - ADDED EXTERNAL EDGE; REMOVED SCRCH COMMON 28 C 98-10-30 BLACK - MODIFIED FOR DISTRIBUTED MEMORY 29 C 00-02-04 BLACK - ADDED CLOUD WATER/ICE 30 C 01-12-11 BLACK - SMOOTHING FOR CFL VIOLATION 31 C 32 C USAGE: CALL VTADV FROM MAIN PROGRAM EBU 33 C INPUT ARGUMENT LIST: 34 C NONE 35 C 36 C OUTPUT ARGUMENT LIST: 37 C NONE 38 C 39 C OUTPUT FILES: 40 C NONE 41 C 42 C SUBPROGRAMS CALLED: 43 C 44 C UNIQUE: NONE 45 C 46 C LIBRARY: NONE 47 C 48 C COMMON BLOCKS: CTLBLK 49 C MASKS 50 C DYNAM 51 C VRBLS 52 C CONTIN 53 C PVRBLS 54 C CLDWTR 55 C INDX 56 C 57 C ATTRIBUTES: Page 2 Source Listing VTADV 2014-12-17 20:47 VTADV.F 58 C LANGUAGE: FORTRAN 90 59 C MACHINE : IBM SP 60 C$$$ 61 C*********************************************************************** 62 P A R A M E T E R 63 & (EDQMX=2.E-5,EDQMN=-2.E-5,EPSQ=1.E-12,EPSQ2=0.2,KSMUD=0) 64 C 65 P A R A M E T E R 66 & (CFL_MAX=0.97) 67 C----------------------------------------------------------------------- 68 INCLUDE "parmeta" 91 INCLUDE "mpp.h" 92 # 1 "./sp.h" 1 93 # 4 94 95 # 70 "VTADV.F" 2 96 C----------------------------------------------------------------------- 155 P A R A M E T E R 156 & (IMJM=IM*JM-JM/2,JAM=6+2*(JM-10) 157 &, LM1=LM-1,LM2=LM-2,LP1=LM+1) 158 C----------------------------------------------------------------------- 159 L O G I C A L 160 & RUN,FIRST,RESTRT,SIGMA,NOSLA 161 C---------------------------------------------------------------------- 162 INCLUDE "CTLBLK.comm" 163 C----------------------------------------------------------------------- 173 INCLUDE "MASKS.comm" 174 C----------------------------------------------------------------------- 185 INCLUDE "DYNAM.comm" 186 C----------------------------------------------------------------------- 202 INCLUDE "VRBLS.comm" 203 C----------------------------------------------------------------------- 213 INCLUDE "CONTIN.comm" 214 C----------------------------------------------------------------------- 224 INCLUDE "PVRBLS.comm" 225 C----------------------------------------------------------------------- 255 INCLUDE "CLDWTR.comm" 256 C----------------------------------------------------------------------- 264 INCLUDE "INDX.comm" 265 C----------------------------------------------------------------------- 277 D I M E N S I O N 278 & WFA ( LM1),WFB ( LM1) 279 C 280 D I M E N S I O N 281 & ETADTL(IDIM1:IDIM2,JDIM1:JDIM2) 282 &,TTA (IDIM1:IDIM2,JDIM1:JDIM2),TQ2A (IDIM1:IDIM2,JDIM1:JDIM2) 283 &,TUA (IDIM1:IDIM2,JDIM1:JDIM2),TVA (IDIM1:IDIM2,JDIM1:JDIM2) 284 &,TTB (IDIM1:IDIM2,JDIM1:JDIM2),TQ2B (IDIM1:IDIM2,JDIM1:JDIM2) 285 &,TUB (IDIM1:IDIM2,JDIM1:JDIM2),TVB (IDIM1:IDIM2,JDIM1:JDIM2) 286 &,ASTI (IDIM1:IDIM2,JDIM1:JDIM2),ASBI (IDIM1:IDIM2,JDIM1:JDIM2) 287 &,DQTI (IDIM1:IDIM2,JDIM1:JDIM2),DQBI (IDIM1:IDIM2,JDIM1:JDIM2) 288 &,VM (IDIM1:IDIM2,JDIM1:JDIM2) 289 &,RPDX (IDIM1:IDIM2,JDIM1:JDIM2),RPDY (IDIM1:IDIM2,JDIM1:JDIM2) 290 &,QDEDT (IDIM1:IDIM2,JDIM1:JDIM2),QDEUT (IDIM1:IDIM2,JDIM1:JDIM2) 291 &,QDEDB (IDIM1:IDIM2,JDIM1:JDIM2),QDEUB (IDIM1:IDIM2,JDIM1:JDIM2) 292 &,EDTD (IDIM1:IDIM2,JDIM1:JDIM2),EDBD (IDIM1:IDIM2,JDIM1:JDIM2) 293 &,EDBF (IDIM1:IDIM2,JDIM1:JDIM2) Page 3 Source Listing VTADV 2014-12-17 20:47 VTADV.F 294 &,DQDE (IDIM1:IDIM2,JDIM1:JDIM2),DQDEB (IDIM1:IDIM2,JDIM1:JDIM2) 295 &,SEDB (IDIM1:IDIM2,JDIM1:JDIM2) 296 C 297 D I M E N S I O N 298 & FNE (IDIM1:IDIM2,JDIM1:JDIM2),FSE (IDIM1:IDIM2,JDIM1:JDIM2) 299 C 300 D I M E N S I O N 301 & TSTL (IDIM1:IDIM2,JDIM1:JDIM2,LM) 302 &,USTL (IDIM1:IDIM2,JDIM1:JDIM2,LM) 303 &,VSTL (IDIM1:IDIM2,JDIM1:JDIM2,LM) 304 &,SAM (IDIM1:IDIM2,JDIM1:JDIM2,LM) 305 &,QBI (IDIM1:IDIM2,JDIM1:JDIM2,LM) 306 &,Q2ST (IDIM1:IDIM2,JDIM1:JDIM2,LM) 307 &,ARRAY1(IDIM1:IDIM2,JDIM1:JDIM2,LM1) 308 &,ARRAY2(IDIM1:IDIM2,JDIM1:JDIM2,LM1) 309 C----------------------------------------------------------------------- 310 R E A L 311 &,ALLOCATABLE,DIMENSION(:,:,:) :: WATER 312 C----------------------------------------------------------------------- 313 R E A L 314 & VAD_TEND1(IDIM1:IDIM2,JDIM1:JDIM2,LM) 315 &,VAD_TEND2(IDIM1:IDIM2,JDIM1:JDIM2,LM) 316 &,VAD_TNDX1(LM),VAD_TNDX2(LM) 317 ! 318 I N T E G E R 319 & LBOT_CFL_T(IDIM1:IDIM2,JDIM1:JDIM2) 320 &,LTOP_CFL_T(IDIM1:IDIM2,JDIM1:JDIM2) 321 &,LBOT_CFL_U(IDIM1:IDIM2,JDIM1:JDIM2) 322 &,LTOP_CFL_U(IDIM1:IDIM2,JDIM1:JDIM2) 323 &,LBOT_CFL_V(IDIM1:IDIM2,JDIM1:JDIM2) 324 &,LTOP_CFL_V(IDIM1:IDIM2,JDIM1:JDIM2) 325 C----------------------------------------------------------------------- 326 C----------------------------------------------------------------------- 327 DTAD=IDTAD*DT 328 C----------------------------------------------------------------------- 329 C--------------DEFINE ADDED UPSTREAM ADVECTION CONSTANTS---------------- 330 C----------------------------------------------------------------------- 331 DO 25 L=1,LM1 332 WFA(L)=DETA(L )/(DETA(L)+DETA(L+1)) 333 WFB(L)=DETA(L+1)/(DETA(L)+DETA(L+1)) 334 25 CONTINUE 335 C--------------NO MOISTURE SLOPE ADJUSTMENT IF NOT WANTED--------------- 336 NOSLA=.FALSE. 337 C IF FALSE, NUMBER OF MOISTURE SLOPE ADJUSTMENT PASSES 338 NMSAP=3 339 C--------------SMOOTHING VERTICAL VELOCITY AT H POINTS------------------ 340 IF(KSMUD.GT.0)THEN 341 !$omp parallel do 342 DO 90 L=1,LM1 343 DO 50 J=MYJS_P4,MYJE_P4 344 DO 50 I=MYIS_P4,MYIE_P4 345 ETADT(I,J,L)=ETADT(I,J,L)*HBM2(I,J) 346 50 CONTINUE 347 C----------------------------------------------------------------------- 348 NSMUD=KSMUD 349 C*** 350 C*** THE FNE, FSE, ETADTL, AND ETADT ARRAYS Page 4 Source Listing VTADV 2014-12-17 20:47 VTADV.F 351 C*** ARE ON OR ASSOCIATED WITH H POINTS 352 C*** 353 DO 90 KS=1,NSMUD 354 DO 80 J=MYJS_P3,MYJE1_P3 355 DO 80 I=MYIS_P3,MYIE_P3 356 FNE(I,J)=(ETADT(I+IHE(J),J+1,L)-ETADT(I,J,L)) 357 1 *HTM(I,J,L+1)*HTM(I+IHE(J),J+1,L+1) 358 80 CONTINUE 359 DO 82 J=MYJS1_P3,MYJE_P3 360 DO 82 I=MYIS_P3,MYIE_P3 361 FSE(I,J)=(ETADT(I+IHE(J),J-1,L)-ETADT(I,J,L)) 362 1 *HTM(I+IHE(J),J-1,L+1)*HTM(I,J,L+1) 363 82 CONTINUE 364 DO 84 J=MYJS2_P1,MYJE2_P1 365 DO 84 I=MYIS_P1,MYIE_P1 366 ETADTL(I,J)=(FNE(I,J)-FNE(I+IHW(J),J-1) 367 1 +FSE(I,J)-FSE(I+IHW(J),J+1))*HBM2(I,J) 368 84 CONTINUE 369 DO 86 J=MYJS2_P1,MYJE2_P1 370 DO 86 I=MYIS_P1,MYIE_P1 371 ETADT(I,J,L)=ETADTL(I,J)*0.125+ETADT(I,J,L) 372 86 CONTINUE 373 90 CONTINUE 374 C----------------------------------------------------------------------- 375 ENDIF 376 !----------------------------------------------------------------------- 377 ! 378 !*** IF THE CFL CRITERION IS VIOLATED THEN LOCATE VERTICAL LIMITS 379 !*** BETWEEN WHICH TO SMOOTH THE TENDENCIES 380 ! 381 !----------------------------------------------------------------------- 382 !$omp parallel do 383 DO J=MYJS,MYJE 384 DO I=MYIS,MYIE 385 LTOP_CFL_T(I,J)=0 386 LBOT_CFL_T(I,J)=0 387 LTOP_CFL_U(I,J)=0 388 LBOT_CFL_U(I,J)=0 389 LTOP_CFL_V(I,J)=0 390 LBOT_CFL_V(I,J)=0 391 ENDDO 392 ENDDO 393 ! 394 DO L=1,LM1 395 ! 396 !$omp parallel do private(cfl) 397 DO J=MYJS2,MYJE2 398 DO I=MYIS,MYIE 399 ! 400 !*** MASS POINTS 401 ! 402 CFL=ETADT(I,J,L)*DTAD*HBM2(I,J)/(0.5*(DETA(L)+DETA(L+1))) 403 IF(ABS(CFL).GT.CFL_MAX)THEN 404 IF(LTOP_CFL_T(I,J).EQ.0)LTOP_CFL_T(I,J)=MAX(L,2) 405 IF(LBOT_CFL_T(I,J).LT.L)LBOT_CFL_T(I,J)=MAX(L,2) 406 ENDIF 407 ! Page 5 Source Listing VTADV 2014-12-17 20:47 VTADV.F 408 !*** U COMPONENT 409 ! 410 CFL=(ETADT(I+IVW(J),J,L)+ETADT(I+IVE(J),J,L))*DTAD*VBM2(I,J) 411 1 /(DETA(L)+DETA(L+1)) 412 ! 413 IF(ABS(CFL).GT.CFL_MAX)THEN 414 IF(LTOP_CFL_U(I,J).EQ.0)LTOP_CFL_U(I,J)=MAX(L,2) 415 IF(LBOT_CFL_U(I,J).LT.L)LBOT_CFL_U(I,J)=MAX(L,2) 416 ENDIF 417 ! 418 !*** V COMPONENT 419 ! 420 CFL=(ETADT(I,J-1,L)+ETADT(I,J+1,L))*DTAD*VBM2(I,J) 421 1 /(DETA(L)+DETA(L+1)) 422 ! 423 IF(ABS(CFL).GT.CFL_MAX)THEN 424 IF(LTOP_CFL_V(I,J).EQ.0)LTOP_CFL_V(I,J)=MAX(L,2) 425 IF(LBOT_CFL_V(I,J).LT.L)LBOT_CFL_V(I,J)=MAX(L,2) 426 ENDIF 427 ! 428 ENDDO 429 ENDDO 430 ! 431 ENDDO 432 ! 433 C----------------------------------------------------------------------- 434 C--------------VERTICAL (MATSUNO) ADVECTION OF T------------------------ 435 C----------------------------------------------------------------------- 436 !$omp parallel do 437 DO 100 J=MYJS,MYJE 438 DO 100 I=MYIS,MYIE 439 TTB(I,J)=0. 440 100 CONTINUE 441 ! 442 DO L=1,LM1 443 !$omp parallel do private(ttak) 444 DO J=MYJS2,MYJE2 445 DO I=MYIS,MYIE 446 TTAK =(T(I,J,L+1)-T(I,J,L))*ETADT(I,J,L)*F4D 447 TSTL(I,J,L)=(TTAK +TTB(I,J))*RDETA(L)+T(I,J,L) 448 VAD_TEND1(I,J,L)=(TTAK +TTB(I,J))*RDETA(L) 449 TTB(I,J)=TTAK 450 ENDDO 451 ENDDO 452 ENDDO 453 ! 454 DO J=MYJS2,MYJE2 455 DO I=MYIS,MYIE 456 VAD_TEND1(I,J,LM)=TTB(I,J)*RDETA(LM) 457 ENDDO 458 ENDDO 459 ! 460 DO 110 L=1,LM1 461 !$omp parallel do private(ttak) 462 DO 110 J=MYJS2,MYJE2 463 DO 110 I=MYIS,MYIE 464 TSTL(I,J,L)=VAD_TEND1(I,J,L)+T(I,J,L) Page 6 Source Listing VTADV 2014-12-17 20:47 VTADV.F 465 110 CONTINUE 466 C 467 !$omp parallel do 468 DO 120 J=MYJS2,MYJE2 469 DO 120 I=MYIS,MYIE 470 TSTL(I,J,LM)=T(I,J,LM)+VAD_TEND1(I,J,LM) 471 120 CONTINUE 472 C----------------------------------------------------------------------- 473 C--------------SECOND (BACKWARD) MATSUNO STEP--------------------------- 474 C----------------------------------------------------------------------- 475 !$omp parallel do 476 DO 125 J=MYJS,MYJE 477 DO 125 I=MYIS,MYIE 478 TTB(I,J)=0. 479 125 CONTINUE 480 C 481 DO L=1,LM1 482 !$omp parallel do private(ttak) 483 DO J=MYJS2,MYJE2 484 DO I=MYIS,MYIE 485 TTAK =(TSTL(I,J,L+1)-TSTL(I,J,L))*ETADT(I,J,L)*F4D 486 VAD_TEND1(I,J,L)=(TTAK +TTB(I,J))*RDETA(L) 487 TTB(I,J)=TTAK 488 ENDDO 489 ENDDO 490 ENDDO 491 ! 492 DO J=MYJS2,MYJE2 493 DO I=MYIS,MYIE 494 VAD_TEND1(I,J,LM)=TTB(I,J)*RDETA(LM) 495 ENDDO 496 ENDDO 497 !----------------------------------------------------------------------- 498 ! 499 !*** IF THE CFL CRITERION IS VIOLATED THEN VERTICALLY SMOOTH 500 !*** THE TENDENCY 501 ! 502 !----------------------------------------------------------------------- 503 !$omp parallel do 504 !$omp& private(lstart,lstop,vad_tndx1) 505 DO J=MYJS2,MYJE2 506 DO I=MYIS,MYIE 507 ! 508 IF(LTOP_CFL_T(I,J).GT.0)THEN 509 LSTART=LTOP_CFL_T(I,J) 510 LSTOP =MIN(LBOT_CFL_T(I,J),LM-1) 511 ! 512 DO L=LSTART,LSTOP 513 VAD_TNDX1(L)=(VAD_TEND1(I,J,L-1)+VAD_TEND1(I,J,L+1) 514 1 +2.*VAD_TEND1(I,J,L))*0.25 515 ENDDO 516 DO L=LSTART,LSTOP 517 VAD_TEND1(I,J,L)=VAD_TNDX1(L) 518 ENDDO 519 ENDIF 520 ! 521 ENDDO Page 7 Source Listing VTADV 2014-12-17 20:47 VTADV.F 522 ENDDO 523 ! 524 DO 140 L=1,LM1 525 !$omp parallel do 526 DO 140 J=MYJS2,MYJE2 527 DO 140 I=MYIS,MYIE 528 T(I,J,L)=VAD_TEND1(I,J,L)+T(I,J,L) 529 ! 530 140 CONTINUE 531 C 532 !$omp parallel do 533 DO 150 J=MYJS2,MYJE2 534 DO 150 I=MYIS,MYIE 535 T(I,J,LM)=T(I,J,LM)+VAD_TEND1(I,J,LM) 536 150 CONTINUE 537 C----------------------------------------------------------------------- 538 C----------------------------------------------------------------------- 539 C--------------PIECEWISE LINEAR UPSTREAM VERTICAL ADVECTION ------------ 540 C-------------------------- OF Q AND CLOUD ----------------------------- 541 C----------------------------------------------------------------------- 542 ALLOCATE(WATER(IDIM1:IDIM2,JDIM1:JDIM2,LM),STAT=I) 543 C----------------------------------------------------------------------- 544 C INTIALIZE Q AT THE BOTTOM INTERFACE AND THE SLOPE ADJUSTMENT 545 C MASK (SAM=1 FOR SA PERMITTED, 0 FOR NOT PERMITTED) 546 C----------------------------------------------------------------------- 547 C 548 C*** LOOP OVER WATER VARIABLES 549 C 550 DO 400 NWTR=1,2 551 C----------------------------------------------------------------------- 552 IF(NWTR.EQ.1)THEN 553 !$omp parallel do 554 DO L=1,LM 555 DO J=JDIM1,JDIM2 556 DO I=IDIM1,IDIM2 557 WATER(I,J,L)=Q(I,J,L) 558 ENDDO 559 ENDDO 560 ENDDO 561 ELSE 562 !$omp parallel do 563 DO L=1,LM 564 DO J=JDIM1,JDIM2 565 DO I=IDIM1,IDIM2 566 WATER(I,J,L)=CWM(I,J,L) 567 ENDDO 568 ENDDO 569 ENDDO 570 ENDIF 571 C----------------------------------------------------------------------- 572 C 573 !$omp parallel do 574 DO 175 L=1,LM 575 DO 175 J=MYJS2,MYJE2 576 DO 175 I=MYIS,MYIE 577 QBI(I,J,L)=WATER(I,J,L) 578 SAM(I,J,L)=1. Page 8 Source Listing VTADV 2014-12-17 20:47 VTADV.F 579 175 CONTINUE 580 IF(NOSLA) GO TO 290 581 C--------------THE SLOPE ADJUSTMENT CODE-------------------------------- 582 C NO SLOPE PERMITTED AT THE TOP AND AT THE BOTTOM LAYER 583 C----------------------------------------------------------------------- 584 !$omp parallel do 585 DO 190 J=MYJS2,MYJE2 586 DO 190 I=MYIS,MYIE 587 SAM(I,J, 1)=0. 588 SAM(I,J,LM)=0. 589 190 CONTINUE 590 C 591 !$omp parallel do 592 DO 200 L=1,LM1 593 DO 200 J=MYJS2,MYJE2 594 DO 200 I=MYIS,MYIE 595 SAM(I,J,L)=SAM(I,J,L)*HTM(I,J,L+1) 596 200 CONTINUE 597 C----------------------------------------------------------------------- 598 C NOW, SEARCH FOR THE MAXIMA AND MINIMA OF Q (AT THE FIRST 599 C PASS) AND FOR LAYERS WHICH HAD OVERADJUSTED (AT SUBSEQUENT 600 C PASSES) DUE TO ROUND-OFF ERRORS 601 C----------------------------------------------------------------------- 602 !$omp parallel do private(dqbi,dqti,extrem) 603 DO 220 L=2,LM1 604 DO 220 J=MYJS2,MYJE2 605 DO 220 I=MYIS,MYIE 606 DQTI(I,J)=WATER(I,J,L)-WATER(I,J,L-1) 607 DQBI(I,J)=WATER(I,J,L+1)-WATER(I,J,L) 608 EXTREM=DQTI(I,J)*DQBI(I,J) 609 IF(EXTREM.LE.0.)SAM(I,J,L)=0. 610 220 CONTINUE 611 C 612 !$omp parallel do 613 DO 230 L=2,LM1 614 DO 230 J=MYJS2,MYJE2 615 DO 230 I=MYIS,MYIE 616 ARRAY1(I,J,L)=WFA(L-1)*(1.-SAM(I,J,L-1))+WFB(L-1) 617 ARRAY2(I,J,L)=WFA(L)+WFB(L)*(1.-SAM(I,J,L+1)) 618 230 CONTINUE 619 DO 260 MSA=1,NMSAP 620 C----------------------------------------------------------------------- 621 C CALCULATE DQ AT INTERFACES AND ADJUST THE SLOPES WHERE 622 C AND TO THE EXTENT PERMITTED OBSERVING THE MONOTONICITY 623 C CONDITION (E.G. VAN LEER, J. COMP. PHYS. 1977, 276-299) 624 C----------------------------------------------------------------------- 625 !$omp parallel do 626 DO 240 J=MYJS2,MYJE2 627 DO 240 I=MYIS,MYIE 628 DQBI(I,J)=2.*WATER(I,J,2)-QBI(I,J,2) -QBI(I,J,1) 629 240 CONTINUE 630 C 631 DO 250 L=2,LM1 632 !$omp parallel do private(asbik,astik,dqtik) 633 DO 250 J=MYJS2,MYJE2 634 DO 250 I=MYIS,MYIE 635 DQTIK =DQBI(I,J) Page 9 Source Listing VTADV 2014-12-17 20:47 VTADV.F 636 ASTIK =ARRAY1(I,J,L)*DQTIK 637 DQBI(I,J)=2.*WATER(I,J,L+1)-QBI(I,J,L+1)-QBI(I,J,L) 638 ASBIK =ARRAY2(I,J,L)*DQBI(I,J) 639 QBI(I,J,L)=QBI(I,J,L) 640 1 +(ASTIK-SIGN(1.,ASTIK) 641 2 *DIM(ABS(ASTIK),ABS(ASBIK)))*SAM(I,J,L) 642 250 CONTINUE 643 260 CONTINUE 644 C----------------------------------------------------------------------- 645 C SLOPE ADJUSTMENT OF THE LAYERS ABOVE THAT NEXT TO THE SURFACE 646 C IS DONE; NOW ADJUST THE LOWERMOST LAYER 647 C----------------------------------------------------------------------- 648 DO 270 L=9,LM1 649 !$omp parallel do 650 DO 270 J=MYJS2,MYJE2 651 DO 270 I=MYIS,MYIE 652 IF(HTM(I,J,L+1).EQ.0.)QBI(I,J,L)=2.*WATER(I,J,L)-QBI(I,J,L-1) 653 270 CONTINUE 654 C 655 !$omp parallel do 656 DO 280 J=MYJS2,MYJE2 657 DO 280 I=MYIS,MYIE 658 QBI(I,J,LM)=2.*WATER(I,J,LM)-QBI(I,J,LM1) 659 280 CONTINUE 660 C----------------------------------------------------------------------- 661 C--------------END OF THE SLOPE ADJUSTMENT CODE------------------------- 662 C----------------------------------------------------------------------- 663 290 CONTINUE 664 !$omp parallel do 665 DO 300 J=MYJS2,MYJE2 666 DO 300 I=MYIS,MYIE 667 QDEDB(I,J)=0. 668 QDEUB(I,J)=0. 669 DQDEB(I,J)=2.*(QBI(I,J,1)-WATER(I,J,1))*RDETA(1) 670 EDBD (I,J)=0. 671 300 CONTINUE 672 C 673 DO 320 L=1,LM1 674 !$omp parallel do private(dqdek,edbfk,edtdk,qdedtk,qdeutk,sedbk) 675 DO 320 J=MYJS2,MYJE2 676 DO 320 I=MYIS,MYIE 677 QDEDTK =QDEDB(I,J) 678 QDEUTK =QDEUB(I,J) 679 SEDBK =SIGN(1.,ETADT(I,J,L)) 680 DQDEK =DQDEB(I,J) 681 DQDEB(I,J)=2.*(QBI(I,J,L+1)-WATER(I,J,L+1))*RDETA(L+1) 682 EDBFK =ETADT(I,J,L)*F4D 683 QDEDB(I,J)=(1.+SEDBK)*(QBI(I,J,L)+DQDEK*EDBFK)*(-EDBFK) 684 QDEUB(I,J)=(1.-SEDBK)*(2.*WATER(I,J,L+1)-QBI(I,J,L+1) 685 1 +DQDEB(I,J)*EDBFK)*EDBFK 686 EDTDK =EDBD(I,J) 687 EDBD (I,J)=ETADT(I,J,L)*(-F4Q) 688 WATER(I,J,L)=WATER(I,J,L)+(QDEDTK-QDEUTK-QDEDB(I,J)+QDEUB(I,J) 689 1 +WATER(I,J,L)*(EDBD(I,J)-EDTDK))*RDETA(L) 690 320 CONTINUE 691 C 692 !$omp parallel do Page 10 Source Listing VTADV 2014-12-17 20:47 VTADV.F 693 DO 330 J=MYJS2,MYJE2 694 DO 330 I=MYIS,MYIE 695 WATER(I,J,LM)=WATER(I,J,LM)+(QDEDB(I,J)-QDEUB(I,J) 696 1 +WATER(I,J,LM)*(-EDBD(I,J)))*RDETA(LM) 697 330 CONTINUE 698 C-------NEGATIVE MOISTURE MAY OCCUR DUE TO VIOLATION OF THE CFL--------- 699 DO 350 L=1,LM1 700 !$omp parallel do 701 DO 350 J=MYJS2,MYJE2 702 DO 350 I=MYIS,MYIE 703 IF(WATER(I,J,L).LT.EPSQ)THEN 704 DQBI(I,J)=WATER(I,J,L) 705 WATER(I,J,L)=EPSQ 706 WATER(I,J,L+1)=WATER(I,J,L+1)+DETA(L)*RDETA(L+1)*DQBI(I,J) 707 ENDIF 708 350 CONTINUE 709 C 710 !$omp parallel do 711 DO 360 J=MYJS2,MYJE2 712 DO 360 I=MYIS,MYIE 713 IF(WATER(I,J,LM).LT.EPSQ)WATER(I,J,LM)=EPSQ 714 360 CONTINUE 715 C----------------------------------------------------------------------- 716 IF(NWTR.EQ.1)THEN 717 !$omp parallel do 718 DO L=1,LM 719 DO J=JDIM1,JDIM2 720 DO I=IDIM1,IDIM2 721 Q(I,J,L)=WATER(I,J,L) 722 ENDDO 723 ENDDO 724 ENDDO 725 ELSE 726 !$omp parallel do 727 DO L=1,LM 728 DO J=JDIM1,JDIM2 729 DO I=IDIM1,IDIM2 730 CWM(I,J,L)=WATER(I,J,L) 731 ENDDO 732 ENDDO 733 ENDDO 734 ENDIF 735 C----------------------------------------------------------------------- 736 400 CONTINUE 737 C 738 DEALLOCATE(WATER,STAT=IER) 739 C----------------------------------------------------------------------- 740 C--------------VERTICAL (MATSUNO) ADVECTION OF Q2----------------------- 741 C----------------------------------------------------------------------- 742 !$omp parallel do 743 DO 420 J=MYJS2,MYJE2 744 DO 420 I=MYIS,MYIE 745 TQ2B(I,J)=Q2(I,J,1)*ETADT(I,J,1)*F4Q2(1) 746 420 CONTINUE 747 C 748 DO 425 L=1,LM2 749 !$omp parallel do private(tq2ak) Page 11 Source Listing VTADV 2014-12-17 20:47 VTADV.F 750 DO 425 J=MYJS2,MYJE2 751 DO 425 I=MYIS,MYIE 752 TQ2AK=(Q2(I,J,L+1)-Q2(I,J,L))*(ETADT(I,J,L)+ETADT(I,J,L+1)) 753 1 *F4Q2(L+1) 754 Q2ST(I,J,L)=TQ2AK+TQ2B(I,J)+Q2(I,J,L) 755 TQ2B(I,J)=TQ2AK 756 425 CONTINUE 757 C 758 !$omp parallel do private(tq2ak) 759 DO 440 J=MYJS2,MYJE2 760 DO 440 I=MYIS,MYIE 761 TQ2AK=(Q2(I,J,LM)-Q2(I,J,LM1))*ETADT(I,J,LM1)*F4Q2(LM) 762 Q2ST(I,J,LM1)=TQ2AK+TQ2B(I,J)+Q2(I,J,LM1) 763 Q2ST(I,J,LM )=Q2(I,J,LM) 764 440 CONTINUE 765 C----------------------------------------------------------------------- 766 C--------------SECOND (BACKWARD) MATSUNO STEP--------------------------- 767 C----------------------------------------------------------------------- 768 !$omp parallel do 769 DO 450 J=MYJS2,MYJE2 770 DO 450 I=MYIS,MYIE 771 TQ2B(I,J)=Q2ST(I,J,1)*ETADT(I,J,1)*F4Q2(1) 772 450 CONTINUE 773 C 774 DO L=1,LM2 775 !$omp parallel do private(tq2ak) 776 DO J=MYJS2,MYJE2 777 DO I=MYIS,MYIE 778 TQ2AK =(Q2ST(I,J,L+1)-Q2ST(I,J,L)) 779 1 *(ETADT(I,J,L)+ETADT(I,J,L+1))*F4Q2(L+1) 780 VAD_TEND1(I,J,L)=TQ2AK+TQ2B(I,J) 781 TQ2B(I,J)=TQ2AK 782 ENDDO 783 ENDDO 784 ENDDO 785 ! 786 DO J=MYJS2,MYJE2 787 DO I=MYIS,MYIE 788 TQ2AK =(Q2ST(I,J,LM)-Q2ST(I,J,LM1))*ETADT(I,J,LM1)*F4Q2(LM) 789 VAD_TEND1(I,J,LM1)=TQ2AK+TQ2B(I,J) 790 ENDDO 791 ENDDO 792 !----------------------------------------------------------------------- 793 ! 794 !*** IF THE CFL CRITERION IS VIOLATED THEN VERTICALLY SMOOTH 795 !*** THE TENDENCY 796 ! 797 !----------------------------------------------------------------------- 798 ! 799 !$omp parallel do 800 !$omp& private(cfl,lbot_cfl,lstart,lstop,ltop_cfl,vad_tend1,vad_tndx1) 801 DO J=MYJS2,MYJE2 802 DO I=MYIS,MYIE 803 ! 804 IF(LTOP_CFL_T(I,J).GT.0)THEN 805 LSTART=LTOP_CFL_T(I,J) 806 LSTOP =MIN(LBOT_CFL_T(I,J),LM-2) Page 12 Source Listing VTADV 2014-12-17 20:47 VTADV.F 807 ! 808 DO L=LSTART,LSTOP 809 VAD_TNDX1(L)=(VAD_TEND1(I,J,L-1)+VAD_TEND1(I,J,L+1) 810 1 +2.*VAD_TEND1(I,J,L))*0.25 811 ENDDO 812 DO L=LSTART,LSTOP 813 VAD_TEND1(I,J,L)=VAD_TNDX1(L) 814 ENDDO 815 ENDIF 816 ! 817 ENDDO 818 ENDDO 819 C 820 DO 470 L=1,LM2 821 !$omp parallel do 822 DO 470 J=MYJS2,MYJE2 823 DO 470 I=MYIS,MYIE 824 Q2(I,J,L)=VAD_TEND1(I,J,L)+Q2(I,J,L) 825 Q2(I,J,L)=AMAX1(Q2(I,J,L),EPSQ2) 826 470 CONTINUE 827 C 828 !$omp parallel do 829 DO 480 J=MYJS2,MYJE2 830 DO 480 I=MYIS,MYIE 831 Q2(I,J,LM1)=VAD_TEND1(I,J,LM1)+Q2(I,J,LM1) 832 Q2(I,J,LM1)=AMAX1(Q2(I,J,LM1),EPSQ2) 833 480 CONTINUE 834 C----------------------------------------------------------------------- 835 C--------------DEFINITION OF VARIABLES NEEDED AT V POINTS--------------- 836 C----------------------------------------------------------------------- 837 !$omp parallel do 838 DO 500 L=1,LM1 839 DO 500 J=MYJS_P1,MYJE_P1 840 DO 500 I=MYIS_P1,MYIE_P1 841 ETADT(I,J,L)=ETADT(I,J,L)*PDSL(I,J)*HBM2(I,J) 842 500 CONTINUE 843 C 844 !$omp parallel do 845 DO 510 J=MYJS2,MYJE2 846 DO 510 I=MYIS,MYIE 847 RPDX(I,J)=1./(PDSL(I+IVW(J),J)+PDSL(I+IVE(J),J)) 848 RPDY(I,J)=1./(PDSL(I,J-1)+PDSL(I,J+1)) 849 510 CONTINUE 850 C----------------------------------------------------------------------- 851 C--------------VERTICAL (MATSUNO) ADVECTION OF U & V-------------------- 852 C----------------------------------------------------------------------- 853 !$omp parallel do 854 DO 520 J=MYJS,MYJE 855 DO 520 I=MYIS,MYIE 856 TUB(I,J)=0. 857 TVB(I,J)=0. 858 520 CONTINUE 859 C 860 DO L=1,LM1 861 !$omp parallel do private(tuak,tvak,vmk) 862 DO J=MYJS2,MYJE2 863 DO I=MYIS,MYIE Page 13 Source Listing VTADV 2014-12-17 20:47 VTADV.F 864 VMK =VTM(I,J,L+1)*VBM2(I,J) 865 TUAK =(ETADT(I+IVW(J),J,L)+ETADT(I+IVE(J),J,L)) 866 1 *(U(I,J,L+1)-U(I,J,L))*RPDX(I,J)*(VMK*F4D) 867 VAD_TEND1(I,J,L)=(TUAK+TUB(I,J))*RDETA(L) 868 TUB(I,J)=TUAK 869 TVAK =(ETADT(I,J-1,L)+ETADT(I,J+1,L))*(V(I,J,L+1)-V(I,J,L)) 870 1 *RPDY(I,J)*(VMK*F4D) 871 VAD_TEND2(I,J,L)=(TVAK+TVB(I,J))*RDETA(L) 872 TVB(I,J)=TVAK 873 ENDDO 874 ENDDO 875 ENDDO 876 ! 877 DO J=MYJS2,MYJE2 878 DO I=MYIS,MYIE 879 VAD_TEND1(I,J,LM)=TUB(I,J)*RDETA(LM) 880 VAD_TEND2(I,J,LM)=TVB(I,J)*RDETA(LM) 881 ENDDO 882 ENDDO 883 ! 884 DO 540 L=1,LM1 885 !$omp parallel do private(tuak,tvak,vmk) 886 DO 540 J=MYJS2,MYJE2 887 DO 540 I=MYIS,MYIE 888 USTL(I,J,L)=VAD_TEND1(I,J,L)+U(I,J,L) 889 VSTL(I,J,L)=VAD_TEND2(I,J,L)+V(I,J,L) 890 TVB(I,J)=TVAK 891 540 CONTINUE 892 C 893 !$omp parallel do 894 DO 550 J=MYJS2,MYJE2 895 DO 550 I=MYIS,MYIE 896 USTL(I,J,LM)=U(I,J,LM)+VAD_TEND1(I,J,LM) 897 VSTL(I,J,LM)=V(I,J,LM)+VAD_TEND2(I,J,LM) 898 550 CONTINUE 899 C----------------------------------------------------------------------- 900 C--------------SECOND (BACKWARD) MATSUNO STEP--------------------------- 901 C----------------------------------------------------------------------- 902 !$omp parallel do 903 DO 560 J=MYJS,MYJE 904 DO 560 I=MYIS,MYIE 905 TUB(I,J)=0. 906 TVB(I,J)=0. 907 560 CONTINUE 908 C 909 DO L=1,LM1 910 !$omp parallel do private(tuak,tvak,vmk) 911 DO J=MYJS2,MYJE2 912 DO I=MYIS,MYIE 913 VMK =VTM(I,J,L+1)*VBM2(I,J) 914 TUAK =(ETADT(I+IVW(J),J,L)+ETADT(I+IVE(J),J,L)) 915 1 *(USTL(I,J,L+1)-USTL(I,J,L))*RPDX(I,J)*(VMK*F4D) 916 VAD_TEND1(I,J,L)=(TUAK+TUB(I,J))*RDETA(L) 917 TUB(I,J)=TUAK 918 TVAK =(ETADT(I,J-1,L)+ETADT(I,J+1,L)) 919 1 *(VSTL(I,J,L+1)-VSTL(I,J,L))*RPDY(I,J)*(VMK*F4D) 920 VAD_TEND2(I,J,L)=(TVAK+TVB(I,J))*RDETA(L) Page 14 Source Listing VTADV 2014-12-17 20:47 VTADV.F 921 TVB(I,J)=TVAK 922 ENDDO 923 ENDDO 924 ENDDO 925 ! 926 DO J=MYJS2,MYJE2 927 DO I=MYIS,MYIE 928 VAD_TEND1(I,J,LM)=TUB(I,J)*RDETA(LM) 929 VAD_TEND2(I,J,LM)=TVB(I,J)*RDETA(LM) 930 ENDDO 931 ENDDO 932 ! 933 !----------------------------------------------------------------------- 934 ! 935 !*** IF THE CFL CRITERION IS VIOLATED THEN VERTICALLY SMOOTH 936 !*** THE TENDENCIES 937 ! 938 !----------------------------------------------------------------------- 939 ! 940 !$omp parallel do 941 !$omp& private(lstart,lstop,vad_tndx1,vad_tndx2) 942 DO J=MYJS2,MYJE2 943 DO I=MYIS,MYIE 944 ! 945 !*** U COMPONENT 946 ! 947 IF(LTOP_CFL_U(I,J).GT.0)THEN 948 LSTART=LTOP_CFL_U(I,J) 949 LSTOP =MIN(LBOT_CFL_U(I,J),LM-1) 950 ! 951 DO L=LSTART,LSTOP 952 VAD_TNDX1(L)=(VAD_TEND1(I,J,L-1)+VAD_TEND1(I,J,L+1) 953 1 +2.*VAD_TEND1(I,J,L))*0.25 954 ENDDO 955 DO L=LSTART,LSTOP 956 VAD_TEND1(I,J,L)=VAD_TNDX1(L) 957 ENDDO 958 ENDIF 959 ! 960 !*** V COMPONENT 961 ! 962 IF(LTOP_CFL_V(I,J).GT.0)THEN 963 LSTART=LTOP_CFL_V(I,J) 964 LSTOP =MIN(LBOT_CFL_V(I,J),LM-1) 965 ! 966 DO L=LSTART,LSTOP 967 VAD_TNDX2(L)=(VAD_TEND2(I,J,L-1)+VAD_TEND2(I,J,L+1) 968 1 +2.*VAD_TEND2(I,J,L))*0.25 969 ENDDO 970 DO L=LSTART,LSTOP 971 VAD_TEND2(I,J,L)=VAD_TNDX2(L) 972 ENDDO 973 ENDIF 974 ! 975 ENDDO 976 ENDDO 977 C Page 15 Source Listing VTADV 2014-12-17 20:47 VTADV.F 978 DO 580 L=1,LM1 979 !$omp parallel do 980 DO 580 J=MYJS2,MYJE2 981 DO 580 I=MYIS,MYIE 982 U(I,J,L)=VAD_TEND1(I,J,L)+U(I,J,L) 983 V(I,J,L)=VAD_TEND2(I,J,L)+V(I,J,L) 984 580 CONTINUE 985 C 986 !$omp parallel do 987 DO 590 J=MYJS2,MYJE2 988 DO 590 I=MYIS,MYIE 989 U(I,J,LM)=U(I,J,LM)+VAD_TEND1(I,J,LM) 990 V(I,J,LM)=V(I,J,LM)+VAD_TEND2(I,J,LM) 991 590 CONTINUE 992 C----------------------------------------------------------------------- 993 RETURN 994 END ENTRY POINTS Name vtadv_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 100 Label 257 254,255 110 Label 282 277,279,280 120 Label 288 285,286 125 Label 296 293,294 140 Label 347 341,343,344 150 Label 353 350,351 175 Label 396 391,392,393 190 Label 406 402,403 200 Label 413 409,410,411 220 Label 427 420,421,422 230 Label 435 430,431,432 240 Label 446 443,444 25 Label 151 148 250 Label 459 448,450,451 260 Label 460 436 270 Label 470 465,467,468 280 Label 476 473,474 290 Label 480 397 300 Label 488 482,483 320 Label 507 490,492,493 330 Label 514 510,511 350 Label 525 516,518,519 360 Label 531 528,529 400 Label 553 367 420 Label 563 560,561 425 Label 573 565,567,568 440 Label 581 576,577 Page 16 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References 450 Label 589 586,587 470 Label 643 637,639,640 480 Label 650 646,647 50 Label 163 160,161 500 Label 659 655,656,657 510 Label 666 662,663 520 Label 675 671,672 540 Label 708 701,703,704 550 Label 715 711,712 560 Label 724 720,721 580 Label 801 795,797,798 590 Label 808 804,805 80 Label 175 171,172 82 Label 180 176,177 84 Label 185 181,182 86 Label 189 186,187 90 Label 190 159,170 ABS Func 220 scalar 220,230,240,458 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 AFSI Scalar 26 R(4) 4 2 17955 COM AKHS Scalar 7 R(4) 4 2 17955 COM AKMS Scalar 7 R(4) 4 2 17955 COM AMAX1 Func 642 scalar 642,649 ARRAY1 Local 124 R(4) 4 3 790020 433,453 ARRAY2 Local 125 R(4) 4 3 790020 434,455 ASBI Local 103 R(4) 4 2 17955 ASBIK Local 455 R(4) 4 scalar 455,458 ASTI Local 103 R(4) 4 2 17955 ASTIK Local 453 R(4) 4 scalar 453,457,458 CFL Local 219 R(4) 4 scalar 219,220,227,230,237,240 CFL_MAX Param 65 R(4) 4 scalar 220,230,240 CLDEFI Scalar 11 R(4) 4 2 17955 COM CLDWTR Common 2 3447720 CONTIN Common 2 12999420 CPGFU Scalar 9 R(4) 4 2 17955 COM CPGFV Scalar 3 R(4) 4 scalar COM CTLBLK Common 4 112 CUPREC Scalar 13 R(4) 4 2 17955 COM CURV Scalar 9 R(4) 4 2 17955 COM CWM Scalar 3 R(4) 4 3 807975 COM 383,547 DDMPU Scalar 12 R(4) 4 2 17955 COM DDMPV Scalar 12 R(4) 4 2 17955 COM DETA Scalar 4 R(4) 4 1 45 COM 149,150,219,228,238,523 DFL Scalar 4 R(4) 4 1 46 COM DIM Func 458 scalar 458 DIV Scalar 6 R(4) 4 3 807975 COM DQBI Local 104 R(4) 4 2 17955 424,425,445,452,454,455,521,523 DQDE Local 111 R(4) 4 2 17955 DQDEB Local 111 R(4) 4 2 17955 486,497,498,502 DQDEK Local 497 R(4) 4 scalar 497,500 DQTI Local 104 R(4) 4 2 17955 423,425 DQTIK Local 452 R(4) 4 scalar 452,453 DT Scalar 6 R(4) 4 scalar COM 144 Page 17 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References DTAD Local 144 R(4) 4 scalar 144,219,227,237 DX Scalar 8 R(4) 4 2 17955 COM DY Scalar 3 R(4) 4 scalar COM DYNAM Common 2 725404 EDBD Local 109 R(4) 4 2 17955 487,503,504,506,513 EDBF Local 110 R(4) 4 2 17955 EDBFK Local 499 R(4) 4 scalar 499,500,502 EDQMN Param 62 R(4) 4 scalar EDQMX Param 62 R(4) 4 scalar EDTD Local 109 R(4) 4 2 17955 EDTDK Local 503 R(4) 4 scalar 503,506 EF4T Scalar 3 R(4) 4 scalar COM 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 EPSQ Param 62 R(4) 4 scalar 520,522,530 EPSQ2 Param 62 R(4) 4 scalar 642,649 ETA Scalar 4 R(4) 4 1 46 COM ETADT Scalar 7 R(4) 4 3 790020 COM 162,173,178,188,219,227,237,263,30 2,496,499,504,562,569,578,588,596, 605,658,682,686,731,735 ETADTL Local 98 R(4) 4 2 17955 183,188 EXTREM Local 425 R(4) 4 scalar 425,426 F Scalar 11 R(4) 4 2 17955 COM F4D Scalar 3 R(4) 4 scalar COM 263,302,499,683,687,732,736 F4Q Scalar 3 R(4) 4 scalar COM 504 F4Q2 Scalar 4 R(4) 4 1 45 COM 562,570,578,588,596,605 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 76 L(4) 4 scalar COM FIS Scalar 7 R(4) 4 2 17955 COM FNE Local 115 R(4) 4 2 17955 173,183 FSE Local 115 R(4) 4 2 17955 178,184 G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GLB_TABLE Common 37 128 HBM2 Scalar 8 R(4) 4 2 17955 COM 162,184,219,658 HBM3 Scalar 8 R(4) 4 2 17955 COM HTM Scalar 6 R(4) 4 3 807975 COM 174,179,412,469 HTMG Scalar 52 R(4) 4 3 4183695 COM I Local 161 I(4) 4 scalar 161,162,172,173,174,177,178,179,18 2,183,184,187,188,201,202,203,204, 205,206,207,215,219,221,222,227,23 1,232,237,241,242,255,256,262,263, 264,265,266,272,273,280,281,286,28 7,294,295,301,302,303,304,310,311, 323,325,326,327,330,331,334,344,34 5,351,352,359,373,374,382,383,393, 394,395,403,404,405,411,412,422,42 3,424,425,426,432,433,434,444,445, 451,452,453,454,455,456,458,468,46 9,474,475,483,484,485,486,487,493, 494,495,496,497,498,499,500,501,50 Page 18 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References 2,503,504,505,506,511,512,513,519, 520,521,522,523,529,530,537,538,54 6,547,561,562,568,569,571,572,577, 578,579,580,587,588,594,595,596,59 7,598,604,605,606,619,621,622,623, 626,627,630,640,641,642,647,648,64 9,657,658,663,664,665,672,673,674, 680,681,682,683,684,685,686,687,68 8,689,695,696,697,704,705,706,707, 712,713,714,721,722,723,729,730,73 1,732,733,734,735,736,737,738,744, 745,746,760,764,765,766,769,770,77 3,779,780,781,784,785,788,798,799, 800,805,806,807 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,21,26,98,99,100,101,102, 103,104,105,106,107,108,109,110,11 1,112,115,118,119,120,121,122,123, 124,125,131,132,136,137,138,139,14 0,141,359,373,382,537,546 IDIM2 Param 19 I(4) 4 scalar 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,21,26,98,99,100,101,102, 103,104,105,106,107,108,109,110,11 1,112,115,118,119,120,121,122,123, 124,125,131,132,136,137,138,139,14 0,141,359,373,382,537,546 IDTAD Scalar 6 I(4) 4 scalar COM 144 IER Local 555 I(4) 4 scalar 555 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 173,174,178,179 IHEG Scalar 9 I(4) 4 1 387 COM IHRST Scalar 5 I(4) 4 scalar COM IHW Scalar 3 I(4) 4 1 399 COM 183,184 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 IM Param 3 I(4) 4 scalar 9,16,19,43,44,45,46,47,51,52,57,72 IMJM Param 72 I(4) 4 scalar INDX Common 2 6564 INDXG Common 8 8084 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 Page 19 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References 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 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 227,664,682,731 IVEG Scalar 9 I(4) 4 1 387 COM IVW Scalar 4 I(4) 4 1 399 COM 227,664,682,731 IVWG Scalar 9 I(4) 4 1 387 COM J Local 160 I(4) 4 scalar 160,162,171,173,174,176,178,179,18 1,183,184,186,188,200,202,203,204, 205,206,207,214,219,221,222,227,23 1,232,237,241,242,254,256,261,263, 264,265,266,271,273,279,281,285,28 7,293,295,300,302,303,304,309,311, 322,325,326,327,330,331,334,343,34 5,350,352,372,374,381,383,392,394, 395,402,404,405,410,412,421,423,42 4,425,426,431,433,434,443,445,450, 452,453,454,455,456,458,467,469,47 3,475,482,484,485,486,487,492,494, 495,496,497,498,499,500,501,502,50 3,504,505,506,510,512,513,518,520, 521,522,523,528,530,536,538,545,54 7,560,562,567,569,571,572,576,578, 579,580,586,588,593,595,596,597,59 8,603,605,606,618,621,622,623,626, 627,630,639,641,642,646,648,649,65 6,658,662,664,665,671,673,674,679, 681,682,683,684,685,686,687,688,68 9,694,696,697,703,705,706,707,711, 713,714,720,722,723,728,730,731,73 2,733,734,735,736,737,738,743,745, 746,759,764,765,766,769,770,773,77 9,780,781,784,785,788,797,799,800, 804,806,807 JAM Param 72 I(4) 4 scalar 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,21,26,98,99,100,101,102, 103,104,105,106,107,108,109,110,11 1,112,115,118,119,120,121,122,123, 124,125,131,132,136,137,138,139,14 Page 20 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References 0,141,359,372,381,536,545 JDIM2 Param 20 I(4) 4 scalar 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,21,26,98,99,100,101,102, 103,104,105,106,107,108,109,110,11 1,112,115,118,119,120,121,122,123, 124,125,131,132,136,137,138,139,14 0,141,359,372,381,536,545 JE_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 JM Param 3 I(4) 4 scalar 9,17,20,43,44,45,46,47,51,52,57,72 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 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 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 KS Local 170 I(4) 4 scalar 170 KSMUD Param 62 I(4) 4 scalar 157,165 L Local 148 I(4) 4 scalar 148,149,150,159,162,173,174,178,17 9,188,211,219,221,222,227,228,231, 232,237,238,241,242,259,263,264,26 5,277,281,298,302,303,329,330,331, 333,334,341,345,371,374,380,383,39 1,394,395,409,412,420,423,424,426, 430,433,434,448,453,454,455,456,45 8,465,469,490,496,498,499,500,501, 504,505,506,516,520,521,522,523,53 5,538,544,547,565,569,570,571,591, 595,596,597,625,626,627,629,630,63 7,641,642,655,658,677,681,682,683, 684,686,688,701,705,706,726,730,73 1,732,733,735,736,737,768,769,770, 772,773,783,784,785,787,788,795,79 9,800 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LBOT_CFL_T Local 136 I(4) 4 2 17955 203,222,327,623 LBOT_CFL_U Local 138 I(4) 4 2 17955 205,232,766 LBOT_CFL_V Local 140 I(4) 4 2 17955 207,242,781 LC Scalar 5 I(4) 4 2 17955 COM LIST Scalar 5 I(4) 4 scalar COM LM Param 3 I(4) 4 scalar 3,4,5,6,7,21,52,73,118,119,120,121 ,122,123,131,132,133,273,287,311,3 27,352,359,371,380,391,405,475,512 ,513,530,535,544,578,580,605,623,6 96,697,713,714,745,746,766,781,806 ,807 LM1 Param 73 I(4) 4 scalar 95,124,125,148,159,211,259,277,298 ,341,409,420,430,448,465,475,490,5 16,578,579,605,606,648,649,655,677 Page 21 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References ,701,726,795 LM2 Param 73 I(4) 4 scalar 565,591,637 LP1 Param 73 I(4) 4 scalar 4 LSM Param 3 I(4) 4 scalar LSTART Local 326 I(4) 4 scalar 326,329,333,622,625,629,765,768,77 2,780,783,787 LSTOP Local 327 I(4) 4 scalar 327,329,333,623,625,629,766,768,77 2,781,783,787 LTOP_CFL_T Local 137 I(4) 4 2 17955 202,221,325,326,621,622 LTOP_CFL_U Local 139 I(4) 4 2 17955 204,231,764,765 LTOP_CFL_V Local 141 I(4) 4 2 17955 206,241,779,780 MAPPINGS Common 56 5024 MASKS Common 2 6894720 MAX Func 221 scalar 221,222,231,232,241,242 MIN Func 327 scalar 327,623,766,781 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 MSA Local 436 I(4) 4 scalar 436 MYIE Scalar 4 I(4) 4 scalar COM 201,215,255,262,272,280,286,294,30 1,310,323,344,351,393,403,411,422, 432,444,451,468,474,483,493,511,51 9,529,561,568,577,587,594,604,619, 640,647,663,672,680,695,704,712,72 1,729,744,760,798,805 MYIE1 Scalar 4 I(4) 4 scalar COM MYIE1_P1 Scalar 8 I(4) 4 scalar COM MYIE1_P2 Scalar 8 I(4) 4 scalar COM MYIE1_P3 Scalar 8 I(4) 4 scalar COM 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 182,187,657 MYIE_P2 Scalar 7 I(4) 4 scalar COM MYIE_P3 Scalar 7 I(4) 4 scalar COM 172,177 MYIE_P4 Scalar 7 I(4) 4 scalar COM 161 MYIE_P5 Scalar 7 I(4) 4 scalar COM MYIS Scalar 4 I(4) 4 scalar COM 201,215,255,262,272,280,286,294,30 1,310,323,344,351,393,403,411,422, 432,444,451,468,474,483,493,511,51 9,529,561,568,577,587,594,604,619, 640,647,663,672,680,695,704,712,72 1,729,744,760,798,805 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 182,187,657 MYIS_P2 Scalar 5 I(4) 4 scalar COM MYIS_P3 Scalar 5 I(4) 4 scalar COM 172,177 MYIS_P4 Scalar 5 I(4) 4 scalar COM 161 MYIS_P5 Scalar 5 I(4) 4 scalar COM Page 22 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References MYJE Scalar 15 I(4) 4 scalar COM 200,254,293,671,720 MYJE1 Scalar 15 I(4) 4 scalar COM MYJE1_P1 Scalar 17 I(4) 4 scalar COM MYJE1_P2 Scalar 17 I(4) 4 scalar COM MYJE1_P3 Scalar 17 I(4) 4 scalar COM 171 MYJE1_P4 Scalar 17 I(4) 4 scalar COM MYJE2 Scalar 15 I(4) 4 scalar COM 214,261,271,279,285,300,309,322,34 3,350,392,402,410,421,431,443,450, 467,473,482,492,510,518,528,560,56 7,576,586,593,603,618,639,646,662, 679,694,703,711,728,743,759,797,80 4 MYJE2_P1 Scalar 18 I(4) 4 scalar COM 181,186 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 656 MYJE_P2 Scalar 16 I(4) 4 scalar COM MYJE_P3 Scalar 16 I(4) 4 scalar COM 176 MYJE_P4 Scalar 16 I(4) 4 scalar COM 160 MYJE_P5 Scalar 16 I(4) 4 scalar COM MYJS Scalar 10 I(4) 4 scalar COM 200,254,293,671,720 MYJS1 Scalar 10 I(4) 4 scalar COM MYJS1_P1 Scalar 12 I(4) 4 scalar COM MYJS1_P2 Scalar 12 I(4) 4 scalar COM MYJS1_P3 Scalar 12 I(4) 4 scalar COM 176 MYJS1_P4 Scalar 12 I(4) 4 scalar COM MYJS2 Scalar 10 I(4) 4 scalar COM 214,261,271,279,285,300,309,322,34 3,350,392,402,410,421,431,443,450, 467,473,482,492,510,518,528,560,56 7,576,586,593,603,618,639,646,662, 679,694,703,711,728,743,759,797,80 4 MYJS2_P1 Scalar 13 I(4) 4 scalar COM 181,186 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 656 MYJS_P2 Scalar 11 I(4) 4 scalar COM Page 23 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References MYJS_P3 Scalar 11 I(4) 4 scalar COM 171 MYJS_P4 Scalar 11 I(4) 4 scalar COM 160 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 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 NMSAP Local 155 I(4) 4 scalar 155,436 NOSLA Local 76 L(4) 4 scalar 153,397 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 NSMUD Local 165 I(4) 4 scalar 165,170 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 NWTR Local 367 I(4) 4 scalar 367,369,533 OMGALF Scalar 5 R(4) 4 3 807975 COM PD Scalar 3 R(4) 4 2 17955 COM PDSL Scalar 3 R(4) 4 2 17955 COM 658,664,665 PREC Scalar 12 R(4) 4 2 17955 COM PSDT Scalar 3 R(4) 4 2 17955 COM PSHLTR Scalar 19 R(4) 4 2 17955 COM PT Scalar 3 R(4) 4 scalar COM PVRBLS Common 2 5530140 Q Scalar 6 R(4) 4 3 807975 COM 374,538 Q10 Scalar 16 R(4) 4 2 17955 COM Q2 Scalar 21 R(4) 4 3 807975 COM 562,569,571,578,579,580,641,642,64 8,649 Q2ST Local 123 R(4) 4 3 807975 571,579,580,588,595,605 Q30 Scalar 14 R(4) 4 2 17955 COM Page 24 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References QBI Local 122 R(4) 4 3 807975 394,445,454,456,469,475,486,498,50 0,501 QDEDB Local 108 R(4) 4 2 17955 484,494,500,505,512 QDEDT Local 107 R(4) 4 2 17955 QDEDTK Local 494 R(4) 4 scalar 494,505 QDEUB Local 108 R(4) 4 2 17955 485,495,501,505,512 QDEUT Local 107 R(4) 4 2 17955 QDEUTK Local 495 R(4) 4 scalar 495,505 QS Scalar 6 R(4) 4 2 17955 COM QSHLTR Scalar 18 R(4) 4 2 17955 COM QWBS Scalar 9 R(4) 4 2 17955 COM QZ0 Scalar 5 R(4) 4 2 17955 COM R Scalar 3 R(4) 4 scalar COM RDETA Scalar 4 R(4) 4 1 45 COM 264,265,273,303,311,486,498,506,51 3,523,684,688,696,697,733,737,745, 746 RES Scalar 7 R(4) 4 2 17955 COM RESTRT Scalar 76 L(4) 4 scalar COM RF Scalar 8 R(4) 4 2 17955 COM RPDX Local 106 R(4) 4 2 17955 664,683,732 RPDY Local 106 R(4) 4 2 17955 665,687,736 RTOP Scalar 4 R(4) 4 3 807975 COM RUN Scalar 76 L(4) 4 scalar COM SAM Local 121 R(4) 4 3 807975 395,404,405,412,426,433,434,458 SEDB Local 112 R(4) 4 2 17955 SEDBK Local 496 R(4) 4 scalar 496,500,501 SI Scalar 10 R(4) 4 2 17955 COM SICE Scalar 4 R(4) 4 2 17955 COM SIGMA Scalar 76 L(4) 4 scalar COM SIGN Func 457 scalar 457,496 SINGLRST Local 2 L(4) 4 scalar SM Scalar 4 R(4) 4 2 17955 COM SNO Scalar 10 R(4) 4 2 17955 COM SR Scalar 5 R(4) 4 2 17955 COM SUBPOST Scalar 2 L(4) 4 scalar COM T Scalar 5 R(4) 4 3 807975 COM 263,264,281,287,345,352 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 TH10 Scalar 16 R(4) 4 2 17955 COM Page 25 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References TH30 Scalar 14 R(4) 4 2 17955 COM THS Scalar 6 R(4) 4 2 17955 COM THZ0 Scalar 5 R(4) 4 2 17955 COM TOPO Common 50 17478548 TQ2A Local 99 R(4) 4 2 17955 TQ2AK Local 569 R(4) 4 scalar 569,571,572,578,579,595,597,598,60 5,606 TQ2B Local 101 R(4) 4 2 17955 562,571,572,579,588,597,598,606 TSHLTR Scalar 18 R(4) 4 2 17955 COM TSTL Local 118 R(4) 4 3 807975 264,281,287,302 TTA Local 99 R(4) 4 2 17955 TTAK Local 263 R(4) 4 scalar 263,264,265,266,302,303,304 TTB Local 101 R(4) 4 2 17955 256,264,265,266,273,295,303,304,31 1 TTVG Scalar 51 R(4) 4 2 92971 COM TUA Local 100 R(4) 4 2 17955 TUAK Local 682 R(4) 4 scalar 682,684,685,731,733,734 TUB Local 102 R(4) 4 2 17955 673,684,685,696,722,733,734,745 TVA Local 100 R(4) 4 2 17955 TVAK Local 686 R(4) 4 scalar 686,688,689,707,735,737,738 TVB Local 102 R(4) 4 2 17955 674,688,689,697,707,723,737,738,74 6 TWBS Scalar 9 R(4) 4 2 17955 COM U Scalar 5 R(4) 4 3 807975 COM 683,705,713,799,806 U00 Scalar 3 R(4) 4 2 17955 COM U10 Scalar 17 R(4) 4 2 17955 COM U30 Scalar 15 R(4) 4 2 17955 COM UL Scalar 4 R(4) 4 1 90 COM USTAR Scalar 3 R(4) 4 2 17955 COM USTL Local 119 R(4) 4 3 807975 705,713,732 UZ0 Scalar 4 R(4) 4 2 17955 COM V Scalar 6 R(4) 4 3 807975 COM 686,706,714,800,807 V10 Scalar 17 R(4) 4 2 17955 COM V30 Scalar 15 R(4) 4 2 17955 COM VAD_TEND1 Local 131 R(4) 4 3 807975 265,273,281,287,303,311,330,331,33 4,345,352,597,606,626,627,630,641, 648,684,696,705,713,733,745,769,77 0,773,799,806 VAD_TEND2 Local 132 R(4) 4 3 807975 688,697,706,714,737,746,784,785,78 8,800,807 VAD_TNDX1 Local 133 R(4) 4 1 45 330,334,626,630,769,773 VAD_TNDX2 Local 133 R(4) 4 1 45 784,788 VBM2 Scalar 3 R(4) 4 2 17955 COM 227,237,681,730 VBM3 Scalar 3 R(4) 4 2 17955 COM VM Local 105 R(4) 4 2 17955 VMK Local 681 R(4) 4 scalar 681,683,687,730,732,736 VRBLS Common 2 13143060 VSTL Local 120 R(4) 4 3 807975 706,714,736 VTADV Subr 2 VTM Scalar 7 R(4) 4 3 807975 COM 681,730 VZ0 Scalar 4 R(4) 4 2 17955 COM WATER Local 128 R(4) 4 3 1 ALC 359,374,383,394,423,424,445,454,46 9,475,486,498,501,505,506,512,513, 520,521,522,523,530,538,547,555 WFA Local 95 R(4) 4 1 44 149,433,434 Page 26 Source Listing VTADV 2014-12-17 20:47 Symbol Table VTADV.F Name Object Declared Type Bytes Dimen Elements Attributes References WFB Local 95 R(4) 4 1 44 150,433,434 WPDAR Scalar 8 R(4) 4 2 17955 COM Z0 Scalar 3 R(4) 4 2 17955 COM Page 27 Source Listing VTADV 2014-12-17 20:47 Subprograms/Common Blocks VTADV.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CLDWTR Common 2 3447720 CONTIN Common 2 12999420 CTLBLK Common 4 112 DYNAM Common 2 725404 GLB_TABLE Common 37 128 INDX Common 2 6564 INDXG Common 8 8084 MAPPINGS Common 56 5024 MASKS Common 2 6894720 MPPCOM Common 1 1464 PVRBLS Common 2 5530140 TEMPCOM Common 42 6603768 TOPO Common 50 17478548 VRBLS Common 2 13143060 VTADV Subr 2 COMPILER OPTIONS BEING USED -align nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores no -auto -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ Page 28 Source Listing VTADV 2014-12-17 20:47 VTADV.F -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest no -ftz -fp_model precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits nono_precise_libm_functions -heap_arrays 0 no -threadprivate_compat -g2 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, Page 29 Source Listing VTADV 2014-12-17 20:47 VTADV.F /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : VTADV.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100