Page 1 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f 1 SUBROUTINE BNDLYR2(PBND,TBND,QBND,RHBND,UBND,VBND, 2 X OMGBND,PWTBND,LVLBND) 3 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 4 C . . . 5 C SUBPROGRAM: BNDLYR2 COMPUTES CONSTANT MASS MEAN FIELDS 6 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-01-29 7 C 8 C ABSTRACT: THIS ROUTINE COMPUTES CONSTANT MASS (BOUNDARY LAYER) 9 C FIELDS. THE FIELDS ARE A MEAN OVER LAYERS PARAMETER DPBND 10 C (PASCALS) THICK. THERE ARE NBND CONSTANT MASS LAYERS, EACH 11 C DPBND THICK STARTING FROM THE SURFACE UP. COMPUTED BOUNDARY 12 C LAYER FIELDS ARE PRESSURE, TEMPERATURE, SPECIFIC HUMIDITY, 13 C RELATIVE HUMIDITY, U AND V WINDS, VERTICAL VELOCITY, 14 C AND PRECIPITABLE WATER. GIVEN THESE FUNDAMENTAL VARIABLES 15 C OTHER FIELDS MAY BE COMPUTED. 16 C 17 C ***WARNING*** IF YOU CHANGE PARAMETER NBND IN THIS ROUTINE 18 C DON'T FOREGET TO CHANGE IT ALSO IN THE CALLING 19 C SUBPROGRAM, MISCLN. 20 C . 21 C 22 C PROGRAM HISTORY LOG: 23 C 93-01-29 RUSS TREADON 24 C 93-05-07 RUSS TREADON - ADDED DOC BLOCK AND MORE COMMENTS. 25 C 93-06-19 RUSS TREADON - ADDED LVLBND TO PARAMETER LIST. 26 C 96-03-07 MIKE BALDWIN - CHANGE PWTR CALC TO INCLUDE CLD WTR 27 C SPEED UP CODE 28 C 98-06-16 T BLACK - CONVERSION FROM 1-D TO 2-D 29 C 98-08-18 MIKE BALDWIN - CHANGE QSBND TO RHBND IN CALL, 30 C COMPUTE RH OVER ICE 31 C 98-12-22 MIKE BALDWIN - BACK OUT RH OVER ICE 32 C 00-01-04 JIM TUCCILLO - MPI VERSION 33 C 34 C USAGE: CALL BNDLYR2(PBND,TBND,QBND,RHBND,UBND,VBND, 35 C OMGBND,PWTBND,LVLBND) 36 C 37 C INPUT ARGUMENT LIST: 38 C NONE 39 C 40 C OUTPUT ARGUMENT LIST: 41 C PBND - LAYER MEAN PRESSURE IN NBND BOUNDARY LAYERS (NBL). 42 C TBND - LAYER MEAN TEMPERATURE IN NBL. 43 C QBND - LAYER MEAN SPECIFIC HUMIDITY IN NBL. 44 C RHBND - LAYER MEAN RELATIVE HUM. (QBND/QSBND) IN NBL. 45 C UBND - LAYER MEAN U WIND COMPONENT IN NBL. 46 C VBND - LAYER MEAN V WIND COMPONENT IN NBL. 47 C OMGBND - LAYER MEAN VERTICAL VELOCITY IN NBL. 48 C PWTBND - LAYER PRECIPITABLE WATER IN NBL. 49 C LVLBND - ETA LAYER AT MIDPOINT OF NBL. 50 C 51 C OUTPUT FILES: 52 C NONE 53 C 54 C SUBPROGRAMS CALLED: 55 C UTILITIES: 56 C H2V 57 C Page 2 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f 58 C LIBRARY: 59 C COMMON - LOOPS 60 C VRBLS 61 C MAPOT 62 C EXTRA 63 C MASKS 64 C PVRBLS 65 C OMGAOT 66 C OPTIONS 67 C INDX 68 C 69 C ATTRIBUTES: 70 C LANGUAGE: FORTRAN 90 71 C MACHINE : CRAY C-90 72 C$$$ 73 C 74 C 75 C INCLUDE GLOBAL PARAMETERS. SET LOCAL PARAMETERS. 76 C DPBND: DEPTH OF BOUNDARY LAYER (IN PASCALS). 77 C 78 INCLUDE "parmeta" 96 INCLUDE "params" 97 C 98 C 99 C 141 PARAMETER (DPBND=30.E2,RHOWAT=1.E3,rrhowt=1./rhowat) 142 PARAMETER (NBND=6) 143 C 144 C DECLARE VARIABLES. 145 C 146 INTEGER LVLBND(IM,JM,NBND) 147 REAL PBND(IM,JM,NBND),TBND(IM,JM,NBND),QBND(IM,JM,NBND) 148 REAL UBND(IM,JM,NBND),VBND(IM,JM,NBND),OMGBND(IM,JM,NBND) 149 REAL PWTBND(IM,JM,NBND),PBINT(IM,JM,NBND+1) 150 REAL QSBND(IM,JM,NBND),RHBND(IM,JM,NBND),IWL(IM,JM),IWM1 151 REAL PSUM(IM,JM,NBND),PVSUM(IM,JM,NBND),NSUM(IM,JM,NBND) 152 C 153 C INCLUDE COMMONS. 154 INCLUDE "LOOPS.comm" 162 INCLUDE "VRBLS.comm" 169 INCLUDE "MAPOT.comm" 182 INCLUDE "EXTRA.comm" 194 INCLUDE "MASKS.comm" 202 INCLUDE "PVRBLS.comm" 214 INCLUDE "CLDWTR.comm" 220 INCLUDE "OMGAOT.comm" 224 INCLUDE "OPTIONS.comm" 228 INCLUDE "INDX.comm" 234 INCLUDE "CTLBLK.comm" 235 C 236 C***************************************************************************** 237 C START BNDLYR HERE 252 CLIMIT =1.0E-20 253 C 254 C LOOP OVER HORIZONTAL GRID. AT EACH MASS POINT COMPUTE 255 C PRESSURE AT THE INTERFACE OF EACH BOUNDARY LAYER. 256 C Page 3 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f 257 !$omp parallel do 258 DO J=JSTA,JEND 259 DO I=1,IM 260 PBINT(I,J,1)=PD(I,J)+PT 261 ENDDO 262 ENDDO 263 C 264 DO LBND=2,NBND+1 265 !$omp parallel do 266 DO J=JSTA,JEND 267 DO I=1,IM 268 PBINT(I,J,LBND)=PBINT(I,J,LBND-1)-DPBND 269 ENDDO 270 ENDDO 271 ENDDO 272 C 273 C LOOP OVER HORIZONTAL. AT EACH MASS POINT COMPUTE 274 C MASS WEIGHTED LAYER MEAN P, T, Q, U, V, OMEGA, 275 C WAND PRECIPITABLE WATER IN EACH BOUNDARY LAYER FROM THE SURFACE UP. 276 C 277 !$omp parallel do 278 !$omp& private(ai,bi,dp,iwl,iwm1,pm,qi,qint,qsat,qw,tmt0,tmt15) 279 DO LBND=1,NBND 280 DO J=JSTA,JEND 281 DO I=1,IM 282 PBND(I,J,LBND) = D00 283 TBND(I,J,LBND) = D00 284 QBND(I,J,LBND) = D00 285 QSBND(I,J,LBND) = D00 286 RHBND(I,J,LBND) = D00 287 UBND(I,J,LBND) = D00 288 VBND(I,J,LBND) = D00 289 OMGBND(I,J,LBND) = D00 290 LVLBND(I,J,LBND) = 0 291 NSUM(I,J,LBND) = 0 292 PSUM(I,J,LBND) = D00 293 PVSUM(I,J,LBND) = D00 294 PWTBND(I,J,LBND) = D00 295 ENDDO 296 ENDDO 297 C 298 IWL=0. 299 DO L=1,LM 300 DO J=JSTA,JEND 301 DO I=1,IM 302 C 303 C COMPUTE IW 304 C 305 IF (L.LE.LMH(I,J)) THEN 306 IWM1=IWL(I,J) 307 IF(CWM(I,J,L).GT.CLIMIT) THEN 308 IF(T(I,J,L).LT.258.15)THEN 309 IWL(I,J)=1. 310 ELSEIF(T(I,J,L).GE.273.15)THEN 311 IWL(I,J)=0. 312 ELSE 313 IF(IWM1.EQ.1.0)IWL(I,J)=1. Page 4 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f 314 ENDIF 315 ELSE 316 IWL(I,J)=0. 317 ENDIF 318 ENDIF 319 C 320 PM=D50*(PINT(I,J,L)+PINT(I,J,L+1)) 321 IF((PBINT(I,J,LBND).GE.PM).AND. 322 1 (PBINT(I,J,LBND+1).LE.PM)) THEN 323 DP = PINT(I,J,L+1)-PINT(I,J,L) 324 PSUM(I,J,LBND) =PSUM(I,J,LBND)+DP 325 NSUM(I,J,LBND) =NSUM(I,J,LBND)+1 326 LVLBND(I,J,LBND)=LVLBND(I,J,LBND)+L 327 TBND(I,J,LBND) =TBND(I,J,LBND)+T(I,J,L)*DP 328 QBND(I,J,LBND) =QBND(I,J,LBND)+Q(I,J,L)*DP 329 OMGBND(I,J,LBND)=OMGBND(I,J,LBND)+OMGA(I,J,L)*DP 330 PWTBND(I,J,LBND)=PWTBND(I,J,LBND) 331 1 +(Q(I,J,L)+CWM(I,J,L))*DP*GI 332 C 333 TMT0=T(I,J,L)-273.16 334 TMT15=AMIN1(TMT0,-15.) 335 AI=0.008855 336 BI=1. 337 IF(TMT0.LT.-20.)THEN 338 AI=0.007225 339 BI=0.9674 340 ENDIF 341 QW=PQ0/PM 342 1 *EXP(A2*(T(I,J,L)-A3)/(T(I,J,L)-A4)) 343 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 344 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 345 IF(TMT0.LT.-15.)THEN 346 QSAT=QI 347 ELSEIF(TMT0.GE.0.)THEN 348 QSAT=QINT 349 ELSE 350 IF(IWL(I,J).GT.0.0) THEN 351 QSAT=QI 352 ELSE 353 QSAT=QINT 354 ENDIF 355 ENDIF 356 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 357 C DELETE THIS LINE TO SWITCH BACK TO RH VS ICE 358 QSAT=QW 359 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 360 QSBND(I,J,LBND)=QSBND(I,J,LBND)+QSAT*DP 361 ENDIF 362 ENDDO 363 ENDDO 364 ENDDO 365 C 366 C NEED TO MAKE SURE THAT PINT HALOS EXIST 367 C 368 DO L = 1, LP1 369 CALL EXCH(PINT(1,1,L)) 370 end do Page 5 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f 371 C 372 DO L=1,LM 373 C 374 DO J=JSTA_M,JEND_M 375 DO I=2,IM-1 376 IE=I+IVE(J) 377 IW=I+IVW(J) 378 PV1=0.25*(PINT(IW,J,L)+PINT(IE,J,L) 379 1 +PINT(I,J+1,L)+PINT(I,J-1,L)) 380 PV2=0.25*(PINT(IW,J,L+1)+PINT(IE,J,L+1) 381 1 +PINT(I,J+1,L+1)+PINT(I,J-1,L+1)) 382 DP=PV2-PV1 383 PMV=0.5*(PV1+PV2) 384 IF((PBINT(IW,J,LBND).GE.PMV).AND. 385 X (PBINT(IW,J,LBND+1).LE.PMV)) THEN 386 PVSUM(I,J,LBND)=PVSUM(I,J,LBND)+DP 387 UBND(I,J,LBND)=UBND(I,J,LBND)+U(I,J,L)*DP 388 VBND(I,J,LBND)=VBND(I,J,LBND)+V(I,J,L)*DP 389 ENDIF 390 C 391 ENDDO 392 ENDDO 393 ENDDO 394 ENDDO 395 C 396 !$omp parallel do 397 !$omp& private(rpsum) 398 DO LBND=1,NBND 399 DO J=JSTA,JEND 400 DO I=1,IM 401 IF(PSUM(I,J,LBND).NE.0.)THEN 402 RPSUM = 1./PSUM(I,J,LBND) 403 LVLBND(I,J,LBND)= LVLBND(I,J,LBND)/NSUM(I,J,LBND) 404 PBND(I,J,LBND) = (PBINT(I,J,LBND)+PBINT(I,J,LBND+1))*0.5 405 TBND(I,J,LBND) = TBND(I,J,LBND)*RPSUM 406 QBND(I,J,LBND) = QBND(I,J,LBND)*RPSUM 407 QSBND(I,J,LBND) = QSBND(I,J,LBND)*RPSUM 408 OMGBND(I,J,LBND)= OMGBND(I,J,LBND)*RPSUM 409 ENDIF 410 ENDDO 411 ENDDO 412 C 413 DO J=JSTA_M,JEND_M 414 DO I=2,IM-1 415 IF(PVSUM(I,J,LBND).NE.0.)THEN 416 RPVSUM =1./PVSUM(I,J,LBND) 417 UBND(I,J,LBND)=UBND(I,J,LBND)*RPVSUM 418 VBND(I,J,LBND)=VBND(I,J,LBND)*RPVSUM 419 ENDIF 420 ENDDO 421 ENDDO 422 ENDDO 423 C 424 C IF NO ETA MID LAYER PRESSURES FELL WITHIN A BND LYR, 425 C FIND THE CLOSEST LAYER TO THE BND LYR AND ASSIGN THE VALUES THERE 426 C 427 !$omp parallel do Page 6 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f 428 !$omp& private(ai,bi,delp,delpv,dp,ie,iw,l,lv,pm,pmin,pminv,pmv, 429 !$omp& qi,qint,qsat,qw,tmt0,tmt15) 430 DO LBND=1,NBND 431 DO J=JSTA,JEND 432 DO I=1,IM 433 IF(PSUM(I,J,LBND).EQ.0.)THEN 434 L=LM 435 PMIN=9999999. 436 PBND(I,J,LBND)=(PBINT(I,J,LBND)+PBINT(I,J,LBND+1))*0.5 437 C 438 DO LL=1,LM 439 PM=D50*(PINT(I,J,LL)+PINT(I,J,LL+1)) 440 DELP=ABS(PM-PBND(I,J,LBND)) 441 IF(DELP.LT.PMIN)THEN 442 PMIN=DELP 443 L=LL 444 ENDIF 445 ENDDO 446 C 447 DP=PINT(I,J,L+1)-PINT(I,J,L) 448 PM=D50*(PINT(I,J,L)+PINT(I,J,L+1)) 449 LVLBND(I,J,LBND)=L 450 TBND(I,J,LBND)=T(I,J,L) 451 QBND(I,J,LBND)=Q(I,J,L) 452 C 453 TMT0=T(I,J,L)-273.16 454 TMT15=AMIN1(TMT0,-15.) 455 AI=0.008855 456 BI=1. 457 IF(TMT0.LT.-20.)THEN 458 AI=0.007225 459 BI=0.9674 460 ENDIF 461 QW=PQ0/PM 462 1 *EXP(A2*(T(I,J,L)-A3)/(T(I,J,L)-A4)) 463 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 464 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 465 IF(TMT0.LT.-15.)THEN 466 QSAT=QI 467 ELSEIF(TMT0.GE.0.)THEN 468 QSAT=QINT 469 ELSE 470 IF(IWL(I,J).GT.0.0) THEN 471 QSAT=QI 472 ELSE 473 QSAT=QINT 474 ENDIF 475 ENDIF 476 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 477 C DELETE THIS LINE TO SWITCH BACK TO RH VS ICE 478 QSAT=QW 479 CMEB 12/22/98 SWITCH TO RH VS WATER NO MATTER WHAT 480 QSBND(I,J,LBND)=QSAT 481 OMGBND(I,J,LBND)=OMGA(I,J,L) 482 PWTBND(I,J,LBND)=(Q(I,J,L)+CWM(I,J,L))*DP*GI 483 ENDIF 484 C Page 7 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f 485 C RH, BOUNDS CHECK 486 C 487 RHBND(I,J,LBND)=QBND(I,J,LBND)/QSBND(I,J,LBND) 488 IF (RHBND(I,J,LBND).GT.1.0) THEN 489 RHBND(I,J,LBND)=1.0 490 QBND(I,J,LBND)=RHBND(I,J,LBND)*QSBND(I,J,LBND) 491 ENDIF 492 IF (RHBND(I,J,LBND).LT.0.01) THEN 493 RHBND(I,J,LBND)=0.01 494 QBND(I,J,LBND)=RHBND(I,J,LBND)*QSBND(I,J,LBND) 495 ENDIF 496 ENDDO 497 ENDDO 498 C 499 DO J=JSTA_M,JEND_M 500 DO I=2,IM-1 501 IF(PVSUM(I,J,LBND).EQ.0.)THEN 502 LV=LM 503 PMINV=9999999. 504 IE=I+IVE(J) 505 IW=I+IVW(J) 506 C 507 C PINT HALOS UPDATED ALREADY 508 C 509 DO LL=1,LM 510 PMV=0.125*(PINT(IW,J,LL)+PINT(IE,J,LL)+ 511 1 PINT(I,J+1,LL)+PINT(I,J-1,LL)+ 512 2 PINT(IW,J,LL+1)+PINT(IE,J,LL+1)+ 513 3 PINT(I,J+1,LL+1)+PINT(I,J-1,LL+1)) 514 DELPV=ABS(PMV-PBND(I,J,LBND)) 515 IF(DELPV.LT.PMINV)THEN 516 PMINV=DELPV 517 LV=LL 518 ENDIF 519 ENDDO 520 C 521 UBND(I,J,LBND)=U(I,J,LV) 522 VBND(I,J,LBND)=V(I,J,LV) 523 ENDIF 524 ENDDO 525 ENDDO 526 ENDDO 527 C 528 C END OF ROUTINE 529 C 530 RETURN 531 END Page 8 Source Listing BNDLYR2 2014-12-17 20:47 Entry Points BNDLYR2.f ENTRY POINTS Name bndlyr2_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References A2 Param 41 R(4) 4 scalar 200,320 A3 Param 41 R(4) 4 scalar 200,320 A4 Param 41 R(4) 4 scalar 200,320 ABS Func 298 scalar 298,372 ACCLIQ Scalar 7 R(4) 4 2 91719 COM ACPREC Scalar 7 R(4) 4 2 91719 COM AD05 Param 33 R(4) 4 scalar AETA Scalar 8 R(4) 4 1 45 COM AI Local 193 R(4) 4 scalar 193,196,201,313,316,321 AKHS Scalar 3 R(4) 4 2 91719 COM AKMS Scalar 3 R(4) 4 2 91719 COM ALPINT Scalar 4 R(4) 4 3 4219074 COM ALSL Scalar 7 R(4) 4 1 39 COM AMIN1 Func 192 scalar 192,201,312,321 BETA Param 37 R(4) 4 scalar 37 BI Local 194 R(4) 4 scalar 194,197,201,314,317,321 BNDLYR2 Subr 1 BTG Param 37 R(4) 4 scalar CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CLDEFI Scalar 6 R(4) 4 2 91719 COM CLDWTR Common 2 17610408 CLIMIT Local 110 R(4) 4 scalar 110,165 CM1 Param 23 R(4) 4 scalar CM2 Param 23 R(4) 4 scalar CM3 Param 23 R(4) 4 scalar CMLD Scalar 6 R(4) 4 scalar COM CP Param 34 R(4) 4 scalar 36,40 CPDR Param 36 R(4) 4 scalar CPHI0 Scalar 9 R(4) 4 scalar COM CTLBLK Common 2 96 CUPREC Scalar 7 R(4) 4 2 91719 COM CWM Scalar 3 R(4) 4 3 4127355 COM 165,189,340 D00 Param 25 R(4) 4 scalar 140,141,142,143,144,145,146,147,15 0,151,152 D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar D05 Param 29 R(4) 4 scalar 33 D1 Param 25 R(4) 4 scalar D115 Param 29 R(4) 4 scalar D125 Param 27 R(4) 4 scalar D25 Param 27 R(4) 4 scalar D50 Param 25 R(4) 4 scalar 178,297,306 Page 9 Source Listing BNDLYR2 2014-12-17 20:47 Symbol Table BNDLYR2.f Name Object Declared Type Bytes Dimen Elements Attributes References D608 Param 29 R(4) 4 scalar DELP Local 298 R(4) 4 scalar 298,299,300 DELPV Local 372 R(4) 4 scalar 372,373,374 DETA Scalar 8 R(4) 4 1 45 COM DISLP Scalar 7 R(4) 4 scalar COM DLMD Scalar 6 R(4) 4 scalar COM DP Local 181 R(4) 4 scalar 181,182,185,186,187,189,218,240,24 4,245,246,305,340 DP10M Param 35 R(4) 4 scalar DP30 Scalar 6 R(4) 4 scalar COM DPBND Param 83 R(4) 4 scalar 126 DPHD Scalar 6 R(4) 4 scalar COM DT Scalar 4 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 ELDR Param 36 R(4) 4 scalar ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar ERLAM0 Scalar 9 R(4) 4 scalar COM ETA Scalar 8 R(4) 4 1 46 COM EXCH Subr 227 227 EXP Func 200 scalar 200,320 EXTRA Common 3 71506268 FIRST Scalar 3 R(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM FMD Param 39 R(4) 4 scalar 39 FMW Param 39 R(4) 4 scalar 39 FSL Scalar 8 R(4) 4 2 91719 COM G Param 32 R(4) 4 scalar 32 GI Param 32 R(4) 4 scalar 36,189,340 H1 Param 26 R(4) 4 scalar 33 H100 Param 27 R(4) 4 scalar H1000 Param 27 R(4) 4 scalar H10000 Param 28 R(4) 4 scalar H10E5 Param 28 R(4) 4 scalar H1E2 Param 30 R(4) 4 scalar H1E4 Param 30 R(4) 4 scalar H1M12 Param 26 R(4) 4 scalar H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar H2 Param 26 R(4) 4 scalar H4 Param 26 R(4) 4 scalar H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar HBM2 Scalar 3 R(4) 4 2 91719 COM HTM Scalar 5 R(4) 4 3 4127355 COM I Local 117 I(4) 4 scalar 117,118,125,126,139,140,141,142,14 3,144,145,146,147,148,149,150,151, 152,159,163,164,165,166,167,168,16 9,171,174,178,179,180,181,182,183, 184,185,186,187,188,189,191,200,20 8,218,233,234,235,237,239,244,245, 246,258,259,260,261,262,263,264,26 5,266,272,273,274,275,276,290,291, Page 10 Source Listing BNDLYR2 2014-12-17 20:47 Symbol Table BNDLYR2.f Name Object Declared Type Bytes Dimen Elements Attributes References 294,297,298,305,306,307,308,309,31 1,320,328,338,339,340,345,346,347, 348,350,351,352,358,359,362,363,36 9,371,372,379,380 IBESSL Scalar 2 I(4) 4 scalar COM ICNT Scalar 7 I(4) 4 1 1024 COM IDAT Scalar 3 I(4) 4 1 3 COM IDIM1 Param 14 I(4) 4 scalar IDIM2 Param 14 I(4) 4 scalar IDN Scalar 7 I(4) 4 scalar COM IDSP Scalar 7 I(4) 4 1 1024 COM IDTAD Scalar 4 I(4) 4 scalar COM IE Local 234 I(4) 4 scalar 234,236,238,362,368,370 IFLAG Scalar 2 I(4) 4 scalar COM IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 IHE Scalar 3 I(4) 4 1 387 COM IHHA Scalar 3 I(4) 4 1 760 COM IHLA Scalar 3 I(4) 4 1 760 COM IHRST Scalar 3 I(4) 4 scalar COM IHW Scalar 3 I(4) 4 1 387 COM IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,14,88,89,90,91, 92,93,117,125,139,159,233,258,272, 290,358 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMT Param 4 I(4) 4 scalar INDX Common 2 6192 INPES Param 6 I(4) 4 scalar 11,14 IOFFS Scalar 2 I(4) 4 scalar COM IOUT Scalar 3 I(4) 4 scalar COM ISHDE Scalar 4 I(4) 4 1 61 COM ITAIL Param 11 I(4) 4 scalar IUNIT Scalar 10 I(4) 4 scalar COM IUP Scalar 7 I(4) 4 scalar COM IVE Scalar 3 I(4) 4 1 387 COM 234,362 IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IVW Scalar 3 I(4) 4 1 387 COM 235,363 IW Local 235 I(4) 4 scalar 235,236,238,242,243,363,368,370 IWL Local 92 R(4) 4 2 91719 156,164,167,169,171,174,208,328 IWM1 Local 92 R(4) 4 scalar 164,171 IXM Scalar 3 I(4) 4 scalar COM IYM Scalar 3 I(4) 4 scalar COM J Local 116 I(4) 4 scalar 116,118,124,126,138,140,141,142,14 3,144,145,146,147,148,149,150,151, 152,158,163,164,165,166,167,168,16 9,171,174,178,179,180,181,182,183, 184,185,186,187,188,189,191,200,20 8,218,232,234,235,236,237,238,239, 242,243,244,245,246,257,259,260,26 1,262,263,264,265,266,271,273,274, 275,276,289,291,294,297,298,305,30 6,307,308,309,311,320,328,338,339, 340,345,346,347,348,350,351,352,35 Page 11 Source Listing BNDLYR2 2014-12-17 20:47 Symbol Table BNDLYR2.f Name Object Declared Type Bytes Dimen Elements Attributes References 7,359,362,363,368,369,370,371,372, 379,380 JAM Param 6 I(4) 4 scalar 3 JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 116,124,138,158,257,289 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 232,271,357 JEND_M2 Scalar 7 I(4) 4 scalar COM JGSTL Param 10 I(4) 4 scalar 15 JGSTR Param 10 I(4) 4 scalar 15 JM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,12,15,88,89,90,91, 92,93 JM2 Param 6 I(4) 4 scalar JMT Param 4 I(4) 4 scalar JNPES Param 6 I(4) 4 scalar 12,15 JRA Scalar 3 I(4) 4 1 760 COM JSTA Scalar 6 I(4) 4 scalar COM 116,124,138,158,257,289 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 232,271,357 JSTA_M2 Scalar 7 I(4) 4 scalar COM JTAIL Param 12 I(4) 4 scalar KMNT Scalar 5 I(4) 4 2 4118670 COM KMNTM Scalar 5 I(4) 4 1 45 COM KSAV Param 14 I(4) 4 scalar KSB Scalar 2 I(4) 4 scalar COM KSFI Param 13 I(4) 4 scalar KSLI Param 14 I(4) 4 scalar KSLP Param 14 I(4) 4 scalar KSLPD Param 13 I(4) 4 scalar KSMUD Param 13 I(4) 4 scalar KSRH Param 14 I(4) 4 scalar KSSLP Param 13 I(4) 4 scalar KST Param 13 I(4) 4 scalar KSUV Param 14 I(4) 4 scalar KUNIT Scalar 10 I(4) 4 scalar COM L Local 157 I(4) 4 scalar 157,163,165,166,168,178,181,184,18 5,186,187,189,191,200,226,227,230, 236,237,238,239,245,246,292,301,30 5,306,307,308,309,311,320,339,340 LB Param 9 I(4) 4 scalar LBND Local 122 I(4) 4 scalar 122,126,137,140,141,142,143,144,14 5,146,147,148,149,150,151,152,179, 180,182,183,184,185,186,187,188,21 8,242,243,244,245,246,256,259,260, 261,262,263,264,265,266,273,274,27 5,276,288,291,294,298,307,308,309, 338,339,340,345,346,347,348,350,35 1,352,359,372,379,380 LC Scalar 3 I(4) 4 2 91719 COM LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LIST Scalar 3 I(4) 4 scalar COM LL Local 296 I(4) 4 scalar 296,297,301,367,368,369,370,371,37 5 Page 12 Source Listing BNDLYR2 2014-12-17 20:47 Symbol Table BNDLYR2.f Name Object Declared Type Bytes Dimen Elements Attributes References LM Param 3 I(4) 4 scalar 2,3,5,8,9,19,157,230,292,296,360,3 67 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 163 LMV Scalar 5 I(4) 4 2 91719 COM LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4,8,226 LSL Scalar 3 I(4) 4 scalar COM LSM Param 3 I(4) 4 scalar 7 LUNIT Scalar 10 I(4) 4 scalar COM LV Local 360 I(4) 4 scalar 360,375,379,380 LVLBND Dummy 2 I(4) 4 3 550314 ARG,INOUT 148,184,261,307 MAPOT Common 2 1448 MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM MUNIT Scalar 10 I(4) 4 scalar COM NBC Scalar 3 I(4) 4 scalar COM NBND Param 84 I(4) 4 scalar 88,89,90,91,92,93,122,137,256,288 NBOCO Scalar 4 I(4) 4 scalar COM NCP Scalar 4 I(4) 4 scalar COM NDDAMP Scalar 4 I(4) 4 scalar COM NFCST Scalar 3 I(4) 4 scalar COM NINC Param 7 I(4) 4 scalar NPHS Scalar 5 I(4) 4 scalar COM NPREC Scalar 4 I(4) 4 scalar COM NRADL Scalar 5 I(4) 4 scalar COM NRADS Scalar 5 I(4) 4 scalar COM NSHDE Scalar 4 I(4) 4 scalar COM NSTAT Param 19 I(4) 4 scalar NSUM Local 93 R(4) 4 3 550314 149,183,261 NSUMD Param 15 I(4) 4 scalar NTSD Scalar 4 I(4) 4 scalar COM NTSTM Scalar 4 I(4) 4 scalar COM NUM_PROCS Scalar 6 I(4) 4 scalar COM NUNIT Scalar 10 I(4) 4 scalar COM NWORD Param 19 I(4) 4 scalar 19 NWORD3 Param 19 I(4) 4 scalar NWORD6 Param 19 I(4) 4 scalar OLDRD Scalar 10 R(4) 4 scalar COM OMGA Scalar 2 R(4) 4 3 4127355 COM 187,339 OMGAOT Common 2 16509420 OMGBND Dummy 2 R(4) 4 3 550314 ARG,INOUT 147,187,266,339 ONEPS Param 39 R(4) 4 scalar OPTIONS Common 2 24 OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar PARALLEL Common 9 8240 PBINT Local 91 R(4) 4 3 642033 118,126,179,180,242,243,262,294 PBND Dummy 1 R(4) 4 3 550314 ARG,INOUT 140,262,294,298,372 PD Scalar 3 R(4) 4 2 91719 COM 118 PDSL Scalar 6 R(4) 4 2 91719 COM PDVP1 Scalar 6 R(4) 4 2 91719 COM PFSL Scalar 7 R(4) 4 2 91719 COM PINT Scalar 4 R(4) 4 3 4219074 COM 178,181,227,236,237,238,239,297,30 5,306,368,369,370,371 Page 13 Source Listing BNDLYR2 2014-12-17 20:47 Symbol Table BNDLYR2.f Name Object Declared Type Bytes Dimen Elements Attributes References PM Local 178 R(4) 4 scalar 178,179,180,199,297,298,306,319 PMIN Local 293 R(4) 4 scalar 293,299,300 PMINV Local 361 R(4) 4 scalar 361,373,374 PMV Local 241 R(4) 4 scalar 241,242,243,368,372 PQ0 Param 41 R(4) 4 scalar 199,319 PREC Scalar 7 R(4) 4 2 91719 COM PSLP Scalar 6 R(4) 4 2 91719 COM PSUM Local 93 R(4) 4 3 550314 150,182,259,260,291 PT Scalar 5 R(4) 4 scalar COM 118 PTSL Scalar 7 R(4) 4 2 91719 COM PV1 Local 236 R(4) 4 scalar 236,240,241 PV2 Local 238 R(4) 4 scalar 238,240,241 PVRBLS Common 2 23480064 PVSUM Local 93 R(4) 4 3 550314 151,244,273,274,359 PWTBND Dummy 2 R(4) 4 3 550314 ARG,INOUT 152,188,340 Q Scalar 5 R(4) 4 3 4127355 COM 186,189,309,340 Q2 Scalar 9 R(4) 4 3 4127355 COM QBND Dummy 1 R(4) 4 3 550314 ARG,INOUT 142,186,264,309,345,348,352 QI Local 201 R(4) 4 scalar 201,204,209,321,324,329 QINT Local 202 R(4) 4 scalar 202,206,211,322,326,331 QS Scalar 3 R(4) 4 2 91719 COM QSAT Local 204 R(4) 4 scalar 204,206,209,211,216,218,324,326,32 9,331,336,338 QSBND Local 92 R(4) 4 3 550314 143,218,265,338,345,348,352 QSL Scalar 9 R(4) 4 2 91719 COM QW Local 199 R(4) 4 scalar 199,201,202,216,319,321,322,336 QWBS Scalar 5 R(4) 4 2 91719 COM QZ0 Scalar 4 R(4) 4 2 91719 COM R Scalar 5 R(4) 4 scalar COM RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDOCP Param 40 R(4) 4 scalar RES Scalar 3 R(4) 4 2 91719 COM RESTRT Scalar 3 R(4) 4 scalar COM RF Scalar 5 R(4) 4 2 91719 COM RG Param 36 R(4) 4 scalar RHBND Dummy 1 R(4) 4 3 550314 ARG,INOUT 144,345,346,347,348,350,351,352 RHOWAT Param 83 R(4) 4 scalar 83 ROG Param 32 R(4) 4 scalar RPSUM Local 260 R(4) 4 scalar 260,263,264,265,266 RPVSUM Local 274 R(4) 4 scalar 274,275,276 RRHOWT Param 83 R(4) 4 scalar RTD Param 34 R(4) 4 scalar RUN Scalar 3 R(4) 4 scalar COM SATDEL Scalar 2 R(4) 4 scalar COM SBD Scalar 6 R(4) 4 scalar COM SI Scalar 6 R(4) 4 2 91719 COM SICE Scalar 3 R(4) 4 2 91719 COM SIGMA Scalar 3 R(4) 4 scalar COM SLP Scalar 5 R(4) 4 2 91719 COM SM Scalar 3 R(4) 4 2 91719 COM SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar SNO Scalar 6 R(4) 4 2 91719 COM SPHI0 Scalar 9 R(4) 4 scalar COM Page 14 Source Listing BNDLYR2 2014-12-17 20:47 Symbol Table BNDLYR2.f Name Object Declared Type Bytes Dimen Elements Attributes References SPL Scalar 7 R(4) 4 1 39 COM SPVAL Scalar 2 R(4) 4 scalar COM SR Scalar 3 R(4) 4 2 91719 COM STDRD Scalar 10 R(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM 166,168,185,191,200,308,311,320 T500 Scalar 6 R(4) 4 2 91719 COM TBND Dummy 1 R(4) 4 3 550314 ARG,INOUT 141,185,263,308 TFRZ Param 39 R(4) 4 scalar THS Scalar 3 R(4) 4 2 91719 COM THZ0 Scalar 4 R(4) 4 2 91719 COM TLM0D Scalar 6 R(4) 4 scalar COM TMT0 Local 191 R(4) 4 scalar 191,192,195,201,203,205,311,312,31 5,321,323,325 TMT15 Local 192 R(4) 4 scalar 192,202,312,322 TPH0D Scalar 6 R(4) 4 scalar COM TSHDE Scalar 8 R(4) 4 1 61 COM TSL Scalar 8 R(4) 4 2 91719 COM TSPH Scalar 5 R(4) 4 scalar COM TTV Scalar 8 R(4) 4 2 91719 COM TWBS Scalar 5 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM 245,379 U00 Scalar 3 R(4) 4 2 91719 COM UBND Dummy 1 R(4) 4 3 550314 ARG,INOUT 145,245,275,379 UL Scalar 3 R(4) 4 1 90 COM UZ0 Scalar 4 R(4) 4 2 91719 COM V Scalar 5 R(4) 4 3 4127355 COM 246,380 VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VBND Dummy 1 R(4) 4 3 550314 ARG,INOUT 146,246,276,380 VRBLS Common 2 67138308 VTM Scalar 5 R(4) 4 3 4127355 COM VZ0 Scalar 4 R(4) 4 2 91719 COM WBD Scalar 6 R(4) 4 scalar COM X1P Scalar 6 R(4) 4 scalar COM Y1P Scalar 6 R(4) 4 scalar COM Z0 Scalar 3 R(4) 4 2 91719 COM Z0SLP Scalar 7 R(4) 4 scalar COM Z1000 Scalar 9 R(4) 4 2 91719 COM ZINT Scalar 4 R(4) 4 3 4219074 COM Page 15 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f 532 C Page 16 Source Listing BNDLYR2 2014-12-17 20:47 Subprograms/Common Blocks BNDLYR2.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References BNDLYR2 Subr 1 CLDWTR Common 2 17610408 CTLBLK Common 2 96 EXTRA Common 3 71506268 INDX Common 2 6192 LOOPS Common 2 748952 MAPOT Common 2 1448 MASKS Common 2 34853220 OMGAOT Common 2 16509420 OPTIONS Common 2 24 PARALLEL Common 9 8240 PVRBLS Common 2 23480064 VRBLS Common 2 67138308 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__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ Page 17 Source Listing BNDLYR2 2014-12-17 20:47 BNDLYR2.f -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 noprecise -fp_model fast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -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 : /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/,.f,./.f,/usrx/local/intel/composerxe/mkl/include/.f, /usrx/local/intel/composerxe/tbb/include/.f,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/.f, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/.f,/usr/local/include/.f,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/.f, /usr/include/.f,/usr/include/.f -list filename : BNDLYR2.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100