Page 1 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 1 SUBROUTINE MISCLN(IMOUT,JMOUT) 2 C 3 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 4 C . . . 5 C SUBPROGRAM: MISCLN POSTS MISCELLANEOUS FIELDS 6 C PRGRMMR: TREADON ORG: W/NP2 DATE: 92-12-20 7 C 8 C ABSTRACT: 9 C THIS ROUTINE HAS BECOME THE CATCH-ALL FOR MISCELLANEOUS 10 C OUTPUT FIELDS POSTED BY THE ETA POST PROCESSOR. 11 C CURRENTLY THIS ROUTINE POSTS THE FOLLOWING FIELDS: 12 C (1) TROPOPAUSE LEVEL Z,P, T, U, V, AND VERTICAL WIND SHEAR, 13 C (2) MAX WIND LEVEL Z, P, U, AND V, 14 C (3) FD LEVEL T, U, AND V, 15 C (4) FREEZING LEVEL Z AND RH, 16 C (5) CONSTANT MASS (BOUNDARY) FIELDS, 17 C (6) LFM LOOK-ALIKE FIELDS, AND 18 C (7) NGM LOOK-ALIKE FIELDS. 19 C 20 C . 21 C 22 C PROGRAM HISTORY LOG: 23 C 92-12-20 RUSS TREADON 24 C 93-06-19 RUSS TREADON - ADDED TYPE 2 CAPE POSTING. 25 C 94-11-07 MIKE BALDWIN - ADDED HELICITY POSTING. 26 C 96-03-26 MIKE BALDWIN - CHANGE ETA BOUNDARY LAYER LABELS FOR GRIB 27 C 96-11-19 MIKE BALDWIN - BACK OUT PREVIOUS CHANGE 28 C 97-04-25 MIKE BALDWIN - CHANGE ETA BOUNDARY LAYER LABELS FOR GRIB 29 C 97-04-29 GEOFF MANIKIN - ADDED TROPOPAUSE HEIGHT AND 30 C MAX WIND LEVEL FIELDS 31 C 98-06-15 T BLACK - CONVERSION FROM 1-D TO 2-D 32 C 98-07-17 MIKE BALDWIN - REMOVED LABL84 33 C 00-01-04 JIM TUCCILLO - MPI VERSION 34 C 35 C USAGE: CALL MISCLN(IMOUT,JMOUT) 36 C INPUT ARGUMENT LIST: 37 C IMOUT - FIRST DIMENSION OF OUTPUT GRID. 38 C JMOUT - SECOND DIMENSION OF OUTPUT GRID. 39 C 40 C OUTPUT ARGUMENT LIST: 41 C NONE 42 C 43 C OUTPUT FILES: 44 C STDOUT - RUN-TIME MESSAGES (FOR DEBUGGING, NOW TURNED OFF) 45 C 46 C SUBPROGRAMS CALLED: 47 C UTILITIES: 48 C TRPAUS - COMPUTE TROPOPAUSE LEVEL FIELDS. 49 C CALMXW - COMPUTE MAX WIND LEVEL FIELDS. 50 C E2OUT - INTERPOLATE STAGGERED E-GRID TO OUTPUT GRID. 51 C SCLFLD - SCALE ARRAY ELEMENTS BY CONSTANT. 52 C OUTPUT - DRIVER FOR TYPE OF FILE OUTPUT. 53 C CALPOT2 - CALCULATE POTENTIAL TEMPERATURE. 54 C FDLVL - COMPUTE FD LEVEL DATA (AGL OR MSL). 55 C FRZLVL - COMPUTE FREEZING LEVEL DATA. 56 C BOUND - BOUND ARRAY ELEMENTS BETWEEN MINIMUM AND MAXIMUM VALUES. 57 C BNDLYR2 - COMPUTE BOUNDARY LAYER FIELDS. Page 2 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 58 C CALDWP2 - CALCULATE DEWPOINT TEMPERATURE. 59 C CALMCVG - CALCULATE MOISTURE CONVERGENCE. 60 C OTLFT2 - COMPUTE LIFTED INDEX AT 500MB. 61 C CALLCL - COMPUTE LCL DATA. 62 C LFMFLD - COMPUTE LFM LOOK-ALIKE FIELDS. 63 C NGMFLD2 - COMPUTE NGM LOOK-ALIKE FIELDS. 64 C CALTHTE - COMPUTE THETA-E. 65 C CALHEL - COMPUTE HELICITY AND STORM MOTION. 66 C 67 C LIBRARY: 68 C COMMON - RQSTFLD 69 C MASKS 70 C IOUNIT 71 C OPTIONS 72 C 73 C ATTRIBUTES: 74 C LANGUAGE: FORTRAN 75 C MACHINE : CRAY C-90 76 C$$$ 77 C 78 C 79 C INCLUDE GRID DIMENSIONS. DERIVE DEPENDENT PARAMETERS. 80 C 81 INCLUDE "parmeta" 99 INCLUDE "parmout" 103 INCLUDE "params" 104 C 105 C SET LOCAL PARAMETERS. MAKE SURE NFD AND NBND AGREE 106 C WITH THE VALUES SET IN SUBROUTINES FDLVL AND BNDLYR2, 107 C RESPECTIVELY. 149 PARAMETER (NFD=6,NBND=6,SMALL=1.E-12) 150 PARAMETER (C2K=273.15) 151 C 152 C DECLARE VARIABLES. 153 C 154 LOGICAL NORTH 155 LOGICAL RUN,FIRST,RESTRT,SIGMA,OLDRD,STDRD 156 LOGICAL FIELD1,FIELD2 157 INTEGER LVLBND(IM,JM,NBND),LB2(IM,JM) 158 REAL P1D(IM,JM),T1D(IM,JM),Q1D(IM,JM),U1D(IM,JM),V1D(IM,JM) 159 REAL SHR1D(IM,JM),Z1D(IM,JM),RH1D(IM,JM) 160 REAL OMGBND(IM,JM,NBND),PWTBND(IM,JM,NBND) 161 REAL PBND(IM,JM,NBND),TBND(IM,JM,NBND),QBND(IM,JM,NBND) 162 REAL UBND(IM,JM,NBND),VBND(IM,JM,NBND),RHBND(IM,JM,NBND) 163 REAL T78483(IM,JM),T89671(IM,JM),P78483(IM,JM),P89671(IM,JM) 164 REAL QM8510(IM,JM),RH4710(IM,JM),RH8498(IM,JM) 165 REAL RH4796(IM,JM),RH1847(IM,JM),UST(IM,JM),VST(IM,JM) 166 REAL RH3310(IM,JM),RH6610(IM,JM),RH3366(IM,JM),PW3310(IM,JM) 167 REAL HTFD(NFD),T6D(IM,JM,NFD),U6D(IM,JM,NFD),V6D(IM,JM,NFD) 168 REAL PETABND(NBND),SIGBND(NBND),HELI(IM,JM) 169 REAL EGRID1(IM,JM),EGRID2(IM,JM) 170 REAL GRID1(IMOUT,JMOUT),GRID2(IMOUT,JMOUT) 171 REAL MAXWP(IM,JM),MAXWZ(IM,JM),MAXWU(IM,JM), MAXWV(IM,JM) 172 173 C 174 C INCLUDE COMMON BLOCKS. 175 INCLUDE "RQSTFLD.comm" Page 3 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 188 INCLUDE "MASKS.comm" 196 INCLUDE "IOUNIT.comm" 204 INCLUDE "OPTIONS.comm" 208 INCLUDE "VRBLS.comm" 215 INCLUDE "EXTRA.comm" 227 INCLUDE "DYNAM.comm" 239 INCLUDE "CTLBLK.comm" 240 C 241 C 242 C EQUIVALENCES FOR SUBROUTINE MISCLN. 257 EQUIVALENCE (TBND(1,1,1),T6D(1,1,1)) 258 EQUIVALENCE (TBND(1,1,2),T6D(1,1,2)) 259 EQUIVALENCE (TBND(1,1,3),T6D(1,1,3)) 260 EQUIVALENCE (TBND(1,1,4),T6D(1,1,4)) 261 EQUIVALENCE (TBND(1,1,5),T6D(1,1,5)) 262 EQUIVALENCE (TBND(1,1,6),T6D(1,1,6)) 263 EQUIVALENCE (UBND(1,1,1),U6D(1,1,1)) 264 EQUIVALENCE (UBND(1,1,2),U6D(1,1,2)) 265 EQUIVALENCE (UBND(1,1,3),U6D(1,1,3)) 266 EQUIVALENCE (UBND(1,1,4),U6D(1,1,4)) 267 EQUIVALENCE (UBND(1,1,5),U6D(1,1,5)) 268 EQUIVALENCE (UBND(1,1,6),U6D(1,1,6)) 269 EQUIVALENCE (VBND(1,1,1),V6D(1,1,1)) 270 EQUIVALENCE (VBND(1,1,2),V6D(1,1,2)) 271 EQUIVALENCE (VBND(1,1,3),V6D(1,1,3)) 272 EQUIVALENCE (VBND(1,1,4),V6D(1,1,4)) 273 EQUIVALENCE (VBND(1,1,5),V6D(1,1,5)) 274 EQUIVALENCE (VBND(1,1,6),V6D(1,1,6)) 275 C 276 C SET FD LEVEL HEIGHTS IN GEOPOTENTAL METERS. 277 DATA HTFD / 914.E0,1524.E0,1829.E0, 278 X 2134.E0,2743.E0,3658.E0/ 279 C 280 C SET MIDPOINT "SIGMA" VALUES FOR ETA BOUNDARY LAYERS. 281 DATA SIGBND / 0.985,0.955,0.925,0.895,0.865,0.835 / 282 DATA PETABND / 15.,45.,75.,105.,135.,165. / 283 C 284 C**************************************************************************** 285 C START MISCLN HERE. 286 C 287 C HELICITY AND STORM MOTION. 288 IF (IGET(162).GT.0.OR.IGET(163).GT.0.OR.IGET(164).GT.0) THEN 289 CALL CALHEL(UST,VST,HELI) 290 IF (IGET(162).GT.0) THEN 291 CALL E2OUT(162,000,HELI,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 292 ID(1:25) = 0 293 ID(10) = 30 294 ID(11) = 0 295 CALL OUTPUT(IOUTYP,IGET(162),LVLS(1,IGET(162)), 296 X GRID1,IMOUT,JMOUT) 297 ENDIF 298 IF ((IGET(163).GT.0).OR.(IGET(164).GT.0)) THEN 299 CALL E2OUT(163,164,UST,VST,GRID1,GRID2,IMOUT,JMOUT) 300 ID(1:25) = 0 301 ID(10) = 60 302 ID(11) = 0 303 IF (IGET(163).GT.0) CALL OUTPUT(IOUTYP,IGET(163), Page 4 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 304 X LVLS(1,IGET(163)),GRID1,IMOUT,JMOUT) 305 ID(1:25) = 0 306 ID(10) = 60 307 ID(11) = 0 308 IF (IGET(164).GT.0) CALL OUTPUT(IOUTYP,IGET(164), 309 X LVLS(1,IGET(164)),GRID2,IMOUT,JMOUT) 310 ENDIF 311 ENDIF 312 C 313 C 314 C 315 C ***BLOCK 1: TROPOPAUSE P, Z, T, U, V, AND WIND SHEAR. 316 C 317 IF ( (IGET(054).GT.0).OR.(IGET(055).GT.0).OR. 318 X (IGET(056).GT.0).OR.(IGET(057).GT.0).OR. 319 X (IGET(177).GT.0).OR. 320 X (IGET(058).GT.0).OR.(IGET(108).GT.0) ) THEN 321 CALL TRPAUS(P1D,T1D,Z1D,U1D,V1D,SHR1D) 322 C 323 C TROPOPAUSE PRESSURE. 324 IF (IGET(054).GT.0) THEN 325 CALL E2OUT(054,000,P1D,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 326 ID(1:25) = 0 327 CALL OUTPUT(IOUTYP,IGET(054),LVLS(1,IGET(054)), 328 X GRID1,IMOUT,JMOUT) 329 ENDIF 330 331 C TROPOPAUSE HEIGHT. 332 IF (IGET(177).GT.0) THEN 333 CALL E2OUT(177,000,Z1D,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 334 ID(1:25) = 0 335 CALL OUTPUT(IOUTYP,IGET(177),LVLS(1,IGET(177)), 336 X GRID1,IMOUT,JMOUT) 337 ENDIF 338 C 339 C TROPOPAUSE TEMPERATURE. 340 IF (IGET(055).GT.0) THEN 341 CALL E2OUT(055,000,T1D,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 342 ID(1:25) = 0 343 CALL OUTPUT(IOUTYP,IGET(055),LVLS(1,IGET(055)), 344 X GRID1,IMOUT,JMOUT) 345 ENDIF 346 C 347 C TROPOPAUSE POTENTIAL TEMPERATURE. 348 IF (IGET(108).GT.0) THEN 349 CALL CALPOT2(P1D,T1D,EGRID1,IM,JM) 350 CALL E2OUT(108,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 351 ID(1:25) = 0 352 CALL OUTPUT(IOUTYP,IGET(108),LVLS(1,IGET(108)), 353 X GRID1,IMOUT,JMOUT) 354 ENDIF 355 C 356 C TROPOPAUSE U WIND AND/OR V WIND. 357 IF ((IGET(056).GT.0).OR.(IGET(057).GT.0)) THEN 358 CALL E2OUT(056,057,U1D,V1D,GRID1,GRID2,IMOUT,JMOUT) 359 ID(1:25) = 0 360 IF (IGET(056).GT.0) CALL OUTPUT(IOUTYP,IGET(056), Page 5 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 361 X LVLS(1,IGET(056)),GRID1,IMOUT,JMOUT) 362 ID(1:25) = 0 363 IF (IGET(057).GT.0) CALL OUTPUT(IOUTYP,IGET(057), 364 X LVLS(1,IGET(057)),GRID2,IMOUT,JMOUT) 365 ENDIF 366 C 367 C TROPOPAUSE WIND SHEAR. 368 IF (IGET(058).GT.0) THEN 369 CALL E2OUT(058,000,SHR1D,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 370 ID(1:25) = 0 371 CALL OUTPUT(IOUTYP,IGET(058),LVLS(1,IGET(058)), 372 X GRID1,IMOUT,JMOUT) 373 ENDIF 374 ENDIF 375 C 376 C 377 C 378 C ***BLOCK 2: MAX WIND LEVEL P, Z, U, AND V 379 C 380 C MAX WIND LEVEL CALCULATIONS 381 IF ((IGET(173).GT.0) .OR. (IGET(174).GT.0) .OR. 382 X (IGET(175).GT.0) .OR. (IGET(176).GT.0)) THEN 383 CALL CALMXW(MAXWP,MAXWZ,MAXWU,MAXWV) 384 ENDIF 385 C PRESSURE OF MAX WIND LEVEL 386 IF (IGET(173).GT.0) THEN 387 CALL E2OUT(173,000,MAXWP,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 388 ID(1:25) = 0 389 CALL OUTPUT(IOUTYP,IGET(173),LVLS(1,IGET(173)), 390 X GRID1,IMOUT,JMOUT) 391 ENDIF 392 C HEIGHT OF MAX WIND LEVEL 393 IF (IGET(174).GT.0) THEN 394 CALL E2OUT(174,000,MAXWZ,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 395 ID(1:25) = 0 396 CALL OUTPUT(IOUTYP,IGET(174),LVLS(1,IGET(174)), 397 X GRID1,IMOUT,JMOUT) 398 ENDIF 399 400 C MAX WIND LEVEL U WIND AND/OR V WIND. 401 IF ((IGET(175).GT.0).OR.(IGET(176).GT.0)) THEN 402 CALL E2OUT(175,176,MAXWU,MAXWV,GRID1,GRID2,IMOUT,JMOUT) 403 ID(1:25) = 0 404 IF (IGET(175).GT.0) CALL OUTPUT(IOUTYP,IGET(175), 405 X LVLS(1,IGET(175)),GRID1,IMOUT,JMOUT) 406 ID(1:25) = 0 407 IF (IGET(176).GT.0) CALL OUTPUT(IOUTYP,IGET(176), 408 X LVLS(1,IGET(176)),GRID2,IMOUT,JMOUT) 409 ENDIF 410 C 411 C 412 C 413 C ***BLOCK 3: FD LEVEL T, U, AND V. 414 C 415 IF ( (IGET(059).GT.0).OR.(IGET(060).GT.0).OR. 416 X (IGET(061).GT.0) ) THEN 417 C Page 6 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 418 C DETERMINE WHETHER TO DO MSL OR AGL FD LEVELS 419 C 420 ITYPE=1 421 DO IFD = 1,NFD 422 IF (IGET(059).GT.0) THEN 423 IF (LVLS(IFD,IGET(059)).GT.1) ITYPE=2 424 ENDIF 425 IF (IGET(060).GT.0) THEN 426 IF (LVLS(IFD,IGET(060)).GT.1) ITYPE=2 427 ENDIF 428 IF (IGET(061).GT.0) THEN 429 IF (LVLS(IFD,IGET(061)).GT.1) ITYPE=2 430 ENDIF 431 ENDDO 432 CALL FDLVL(ITYPE,T6D,U6D,V6D) 433 C 434 DO 10 IFD = 1,NFD 435 ID(1:25) = 0 436 ISVALUE = NINT(HTFD(IFD)) 437 ID(11) = ISVALUE 438 C 439 C FD LEVEL TEMPERATURE. 440 IF (IGET(059).GT.0) THEN 441 IF (LVLS(IFD,IGET(059)).GT.0) THEN 442 CALL E2OUT(059,000,T6D(1,1,IFD),EGRID2, 443 X GRID1,GRID2,IMOUT,JMOUT) 444 CALL OUTPUT(IOUTYP,IGET(059), 445 X LVLS(IFD,IGET(059)),GRID1,IMOUT,JMOUT) 446 ENDIF 447 ENDIF 448 C 449 C FD LEVEL U WIND AND/OR V WIND. 450 IF ((IGET(060).GT.0).OR.(IGET(061).GT.0)) THEN 451 CALL E2OUT(060,061,U6D(1,1,IFD),V6D(1,1,IFD), 452 X GRID1,GRID2,IMOUT,JMOUT) 453 IF (IGET(060).GT.0) THEN 454 IF (LVLS(IFD,IGET(060)).GT.0) CALL OUTPUT(IOUTYP, 455 X IGET(060),LVLS(IFD,IGET(060)),GRID1,IMOUT,JMOUT) 456 ENDIF 457 IF (IGET(061).GT.0) THEN 458 IF (LVLS(IFD,IGET(061)).GT.0) CALL OUTPUT(IOUTYP, 459 X IGET(061),LVLS(IFD,IGET(061)),GRID2,IMOUT,JMOUT) 460 ENDIF 461 ENDIF 462 10 CONTINUE 463 ENDIF 464 C 465 C 466 C 467 C ***BLOCK 4: FREEZING LEVEL Z AND RH. 468 C 469 IF ( (IGET(062).GT.0).OR.(IGET(063).GT.0) ) THEN 470 CALL FRZLVL(Z1D,RH1D) 471 C 472 C FREEZING LEVEL HEIGHT. 473 IF (IGET(062).GT.0) THEN 474 CALL E2OUT(062,000,Z1D,EGRID2,GRID1,GRID2,IMOUT,JMOUT) Page 7 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 475 ID(1:25) = 0 476 CALL BOUND (GRID1,D00,H99999,IMOUT,JMOUT) 477 CALL OUTPUT(IOUTYP,IGET(062),LVLS(1,IGET(062)), 478 X GRID1,IMOUT,JMOUT) 479 ENDIF 480 C 481 C FREEZING LEVEL RELATIVE HUMIDITY. 482 IF (IGET(063).GT.0) THEN 483 CALL E2OUT(063,000,RH1D,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 484 ID(1:25) = 0 485 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 486 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 487 CALL OUTPUT(IOUTYP,IGET(063),LVLS(1,IGET(063)), 488 X GRID1,IMOUT,JMOUT) 489 ENDIF 490 ENDIF 491 IF (IGET(165).GT.0) THEN 492 CALL FRZLVL2(Z1D,RH1D) 493 C 494 C HIGHEST FREEZING LEVEL HEIGHT. 495 CALL E2OUT(165,000,Z1D,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 496 ID(1:25) = 0 497 CALL BOUND (GRID1,D00,H99999,IMOUT,JMOUT) 498 CALL OUTPUT(IOUTYP,IGET(165),LVLS(1,IGET(165)), 499 X GRID1,IMOUT,JMOUT) 500 ENDIF 501 C 502 C 503 C 504 C ***BLOCK 5: BOUNDARY LAYER FIELDS. 505 C 506 IF ( (IGET(067).GT.0).OR.(IGET(068).GT.0).OR. 507 X (IGET(069).GT.0).OR.(IGET(070).GT.0).OR. 508 X (IGET(071).GT.0).OR.(IGET(072).GT.0).OR. 509 X (IGET(073).GT.0).OR.(IGET(074).GT.0).OR. 510 X (IGET(088).GT.0).OR.(IGET(089).GT.0).OR. 511 X (IGET(090).GT.0).OR.(IGET(075).GT.0).OR. 512 X (IGET(109).GT.0).OR.(IGET(110).GT.0).OR. 513 X (IGET(031).GT.0).OR.(IGET(032).GT.0).OR. 514 X (IGET(107).GT.0).OR.(IGET(091).GT.0).OR. 515 X (IGET(092).GT.0).OR.(IGET(093).GT.0).OR. 516 X (IGET(094).GT.0).OR.(IGET(095).GT.0).OR. 517 X (IGET(096).GT.0).OR.(IGET(097).GT.0).OR. 518 X (IGET(098).GT.0) ) THEN 519 C 520 C COMPUTE ETA BOUNDARY LAYER FIELDS. 521 CALL BNDLYR2(PBND,TBND,QBND,RHBND,UBND,VBND, 522 X OMGBND,PWTBND,LVLBND) 523 C 524 C LOOP OVER NBND BOUNDARY LAYERS. 525 DO 20 LBND = 1,NBND 526 ID(1:25) = 0 527 ID(10) = NINT(PETABND(LBND)+15.) 528 ID(11) = NINT(PETABND(LBND)-15.) 529 C 530 C BOUNDARY LAYER PRESSURE. 531 IF (IGET(067).GT.0) THEN Page 8 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 532 IF (LVLS(LBND,IGET(067)).GT.0) THEN 533 CALL E2OUT(067,000,PBND(1,1,LBND),EGRID2, 534 X GRID1,GRID2,IMOUT,JMOUT) 535 CALL OUTPUT(IOUTYP,IGET(067), 536 X LVLS(LBND,IGET(067)),GRID1,IMOUT,JMOUT) 537 ENDIF 538 ENDIF 539 C 540 C BOUNDARY LAYER TEMPERATURE. 541 IF (IGET(068).GT.0) THEN 542 IF (LVLS(LBND,IGET(068)).GT.0) THEN 543 CALL E2OUT(068,000,TBND(1,1,LBND),EGRID2, 544 X GRID1,GRID2,IMOUT,JMOUT) 545 CALL OUTPUT(IOUTYP,IGET(068),LVLS(LBND,IGET(068)), 546 X GRID1,IMOUT,JMOUT) 547 ENDIF 548 ENDIF 549 C 550 C BOUNDARY LAYER POTENTIAL TEMPERATURE. 551 IF (IGET(069).GT.0) THEN 552 IF (LVLS(LBND,IGET(069)).GT.0) THEN 553 CALL CALPOT2(PBND(1,1,LBND),TBND(1,1,LBND),EGRID1,IM,JM) 554 CALL E2OUT(069,000,EGRID1,EGRID2,GRID1,GRID2, 555 X IMOUT,JMOUT) 556 CALL OUTPUT(IOUTYP,IGET(069),LVLS(LBND,IGET(069)), 557 X GRID1,IMOUT,JMOUT) 558 ENDIF 559 ENDIF 560 C 561 C BOUNDARY LAYER RELATIVE HUMIDITY. 562 IF (IGET(072).GT.0) THEN 563 IF (LVLS(LBND,IGET(072)).GT.0) THEN 564 CALL E2OUT(072,000,RHBND(1,1,LBND),EGRID2, 565 X GRID1,GRID2,IMOUT,JMOUT) 566 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 567 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 568 CALL OUTPUT(IOUTYP,IGET(072),LVLS(LBND,IGET(072)), 569 X GRID1,IMOUT,JMOUT) 570 ENDIF 571 ENDIF 572 C 573 C BOUNDARY LAYER DEWPOINT TEMPERATURE. 574 IF (IGET(070).GT.0) THEN 575 IF (LVLS(LBND,IGET(070)).GT.0) THEN 576 CALL CALDWP2(PBND(1,1,LBND),QBND(1,1,LBND),EGRID1, 577 X TBND(1,1,LBND)) 578 CALL E2OUT(070,000,EGRID1,EGRID2,GRID1,GRID2, 579 X IMOUT,JMOUT) 580 CALL OUTPUT(IOUTYP,IGET(070),LVLS(LBND,IGET(070)), 581 X GRID1,IMOUT,JMOUT) 582 ENDIF 583 ENDIF 584 C 585 C BOUNDARY LAYER SPECIFIC HUMIDITY. 586 IF (IGET(071).GT.0) THEN 587 IF (LVLS(LBND,IGET(071)).GT.0) THEN 588 CALL E2OUT(071,000,QBND(1,1,LBND),EGRID2, Page 9 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 589 X GRID1,GRID2,IMOUT,JMOUT) 590 CALL BOUND(GRID1,H1M12,H99999,IMOUT,JMOUT) 591 CALL OUTPUT(IOUTYP,IGET(071),LVLS(LBND,IGET(071)), 592 X GRID1,IMOUT,JMOUT) 593 ENDIF 594 ENDIF 595 C 596 C BOUNDARY LAYER MOISTURE CONVERGENCE. 597 IF (IGET(088).GT.0) THEN 598 IF (LVLS(LBND,IGET(088)).GT.0) THEN 599 CALL CALMCVG(QBND(1,1,LBND),UBND(1,1,LBND), 600 X VBND(1,1,LBND),-1,EGRID1) 601 C CONVERT TO DIVERGENCE FOR GRIB 602 DO J=JSTA,JEND 603 DO I=1,IM 604 EGRID1(I,J)=-1.0*EGRID1(I,J) 605 ENDDO 606 ENDDO 607 CALL E2OUT(088,000,EGRID1,EGRID2,GRID1,GRID2, 608 X IMOUT,JMOUT) 609 CALL OUTPUT(IOUTYP,IGET(088),LVLS(LBND,IGET(088)), 610 X GRID1,IMOUT,JMOUT) 611 ENDIF 612 ENDIF 613 C 614 C BOUNDARY LAYER U WIND AND/OR V WIND. 615 C 616 FIELD1=.FALSE. 617 FIELD2=.FALSE. 618 C 619 IF(IGET(073).GT.0)THEN 620 IF(LVLS(LBND,IGET(073)).GT.0)FIELD1=.TRUE. 621 ENDIF 622 IF(IGET(074).GT.0)THEN 623 IF(LVLS(LBND,IGET(074)).GT.0)FIELD2=.TRUE. 624 ENDIF 625 C 626 IF(FIELD1.OR.FIELD2)THEN 627 CALL E2OUT(073,074,UBND(1,1,LBND),VBND(1,1,LBND), 628 X GRID1,GRID2,IMOUT,JMOUT) 629 C 630 C TWO ROUTINES WERE WRITTEN DURING POST PROCESSOR 631 C DEVELOPMENT TO MIMIC LFM BOUNDARY LAYER WINDS. 632 C BOTH APPLY AN EKMAN SPIRAL ROTATION TO AN ARTIFICIALLY 633 C REDUCED GEOSTROPHIC WIND FIELD. SUBROUTINE CALEKM 634 C BASES THE GEOSTROPHIC WIND ON 1000MB HEIGHTS WHILE 635 C CALEKM2 USES THE SEA LEVEL PRESSURE FIELD. IF 636 C LBND EQUALS ONE AND LVLS IS 8 POST THE CALEKM WINDS 637 C AS BOUNDARY LAYER ONE WINDS. IF LBND EQUALS ONE AND 638 C LVLS IS 9 POST CALEKM2 WINDS AS BOUNDARY LAYER ONE 639 C WINDS. OTHERWISE, POST THE "TRUE" BOUNDARY LAYER 640 C ONE WINDS. 641 C - 19 OCTOBER 1993, RUSS TREADON 642 C 643 IF (IGET(073).GT.0) THEN 644 IF ((LBND.EQ.1).AND.(LVLS(LBND,IGET(073)).EQ.8)) THEN 645 WRITE(STDOUT,*)'MISCLN: BNDLYR ',LBND, Page 10 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 646 X ' UV USING EKMAN-GEO Z1000 UV' 647 CALL CALEKM(U1D,V1D) 648 CALL E2OUT(073,074,U1D,V1D,GRID1,GRID2,IMOUT,JMOUT) 649 ENDIF 650 ENDIF 651 C 652 IF (IGET(073).GT.0) THEN 653 IF ((LBND.EQ.1).AND.(LVLS(LBND,IGET(073)).EQ.9)) THEN 654 WRITE(STDOUT,*)'MISCLN: BNDLYR ',LBND, 655 X ' UV USING EKMAN-GEO PMSL UV' 656 CALL CALEKM2(U1D,V1D) 657 CALL E2OUT(073,074,U1D,V1D,GRID1,GRID2,IMOUT,JMOUT) 658 ENDIF 659 ENDIF 660 C 661 IF (IGET(073).GT.0) THEN 662 IF (LVLS(LBND,IGET(073)).GT.0) 663 X CALL OUTPUT(IOUTYP,IGET(073), 664 X LVLS(LBND,IGET(073)),GRID1,IMOUT,JMOUT) 665 ENDIF 666 IF (IGET(074).GT.0) THEN 667 IF (LVLS(LBND,IGET(074)).GT.0) 668 X CALL OUTPUT(IOUTYP,IGET(074), 669 X LVLS(LBND,IGET(074)),GRID2,IMOUT,JMOUT) 670 ENDIF 671 ENDIF 672 C 673 C BOUNDARY LAYER OMEGA. 674 IF (IGET(090).GT.0) THEN 675 IF (LVLS(LBND,IGET(090)).GT.0) THEN 676 CALL E2OUT(090,000,OMGBND(1,1,LBND),EGRID2, 677 X GRID1,GRID2,IMOUT,JMOUT) 678 CALL OUTPUT(IOUTYP,IGET(090),LVLS(LBND,IGET(090)), 679 X GRID1,IMOUT,JMOUT) 680 ENDIF 681 ENDIF 682 C 683 C BOUNDARY LAYER PRECIPITBLE WATER. 684 IF (IGET(089).GT.0) THEN 685 IF (LVLS(LBND,IGET(089)).GT.0) THEN 686 CALL E2OUT(089,000,PWTBND(1,1,LBND),EGRID2, 687 X GRID1,GRID2,IMOUT,JMOUT) 688 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) 689 CALL OUTPUT(IOUTYP,IGET(089),LVLS(LBND,IGET(089)), 690 X GRID1,IMOUT,JMOUT) 691 ENDIF 692 ENDIF 693 C 694 C BOUNDARY LAYER LIFTED INDEX. 695 IF (IGET(075).GT.0) THEN 696 IF (LVLS(LBND,IGET(075)).GT.0) THEN 697 CALL OTLFT2(PBND(1,1,LBND),TBND(1,1,LBND), 698 X QBND(1,1,LBND),EGRID1) 699 CALL E2OUT(075,000,EGRID1,EGRID2, 700 X GRID1,GRID2,IMOUT,JMOUT) 701 C 702 C 23 SEPTEMBER 1993, RUSS TREADON. Page 11 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 703 C ON LFM FORECAST GRID 026 WE POST THE FIRST ETA 704 C BOUNDARY LAYER LI AS THE SURFACE TO 500MB LI. 705 C THE CALL TO E2OUT INTERPOLATES THE LI FROM 706 C THE E-GRID TO THE OUTPUT GRID. IF WE ARE POSTING 707 C DATA THE THE LFM GRID WE NEED TO SAVE THIS 708 C INTERPOLATED FIELD FOR USE BELOW. WHY? WHEN 709 C PACKING THE DATA IN GRIB1 THE VALUES IN GRID1 710 C ARE MULTIPLIED BY 10. IF WE SEND GRID1 TO THE 711 C GRIB PACKER A SECOND TIME THE LI'S ARE 10 TIMES 712 C TOO BIG. 713 C 714 CALL E2OUT(075,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 715 IF (KGTYPE.EQ.026) THEN 716 DO J=JSTA,JEND 717 DO I=1,IMOUT 718 GRID2(I,J) = GRID1(I,J) 719 ENDDO 720 ENDDO 721 ENDIF 722 723 724 CALL OUTPUT(IOUTYP,IGET(075),LVLS(LBND,IGET(075)), 725 X GRID1,IMOUT,JMOUT) 726 C 727 C 25 JUNE 1993, RUSS TREADON. 728 C ON THE LFM FORECAST GRID (026) WE POST THE FIRST 729 C BOUNDARY LAYER LIFTED INDEX AS THE SURFACE LI. 730 C 731 IF (KGTYPE.EQ.026) THEN 732 DO J=JSTA,JEND 733 DO I=1,IMOUT 734 GRID2(I,J) = GRID2(I,J) + C2K 735 ENDDO 736 ENDDO 737 C 738 ID(1:25) = 0 739 ID(9) = 101 740 ID(10) = 50 741 ID(11) = 0 742 CALL OUTPUT(IOUTYP,IGET(030),LVLS(1,IGET(030)), 743 X GRID2,IMOUT,JMOUT) 744 ENDIF 745 ENDIF 746 ENDIF 747 C 748 C END OF ETA BOUNDARY LAYER LOOP. 749 20 CONTINUE 750 C 751 C BEST LIFTED INDEX FROM BOUNDARY LAYER FIELDS. 752 C 753 IF (IGET(031).GT.0) THEN 754 DO J=JSTA,JEND 755 DO I=1,IM 756 EGRID1(I,J) = H99999 757 EGRID2(I,J) = H99999 758 ENDDO 759 ENDDO Page 12 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 760 C 761 DO 50 LBND = 1,NBND 762 CALL OTLFT2(PBND(1,1,LBND),TBND(1,1,LBND), 763 X QBND(1,1,LBND),EGRID2) 764 DO J=JSTA,JEND 765 DO I=1,IM 766 EGRID1(I,J)=AMIN1(EGRID1(I,J),EGRID2(I,J)) 767 ENDDO 768 ENDDO 769 50 CONTINUE 770 CALL E2OUT(031,000,EGRID1,EGRID2, 771 X GRID1,GRID2,IMOUT,JMOUT) 772 ID(1:25) = 0 773 ID(10) = PETABND(NBND)+15. 774 ID(11) = PETABND(1)-15. 775 CALL OUTPUT(IOUTYP,IGET(031),LVLS(1,IGET(031)), 776 X GRID1,IMOUT,JMOUT) 777 ENDIF 778 C 779 C BEST BOUNDARY LAYER CAPE AND CINS. 780 C 781 FIELD1=.FALSE. 782 FIELD2=.FALSE. 783 C 784 IF(IGET(032).GT.0)THEN 785 IF(LVLS(2,IGET(032)).GT.0)FIELD1=.TRUE. 786 ENDIF 787 IF(IGET(107).GT.0)THEN 788 IF(LVLS(2,IGET(107)).GT.0)FIELD2=.TRUE. 789 ENDIF 790 C 791 IF(FIELD1.OR.FIELD2)THEN 792 ITYPE = 2 793 C 794 DO J=JSTA,JEND 795 DO I=1,IM 796 EGRID1(I,J) = -H99999 797 EGRID2(I,J) = -H99999 798 ENDDO 799 ENDDO 800 C 801 DO 80 LBND = 1,NBND 802 CALL CALTHTE(PBND(1,1,LBND),TBND(1,1,LBND), 803 X QBND(1,1,LBND),EGRID1) 804 DO J=JSTA,JEND 805 DO I=1,IM 806 IF (EGRID1(I,J).GT.EGRID2(I,J)) THEN 807 EGRID2(I,J) = EGRID1(I,J) 808 LB2(I,J) = LVLBND(I,J,LBND) 809 P1D(I,J) = PBND(I,J,LBND) 810 T1D(I,J) = TBND(I,J,LBND) 811 Q1D(I,J) = QBND(I,J,LBND) 812 ENDIF 813 ENDDO 814 ENDDO 815 80 CONTINUE 816 C Page 13 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 817 CALL CALCAPE(ITYPE,P1D,T1D,Q1D,LB2,EGRID1,EGRID2) 818 C 819 IF (IGET(032).GT.0) THEN 820 CALL E2OUT(032,000,EGRID1,EGRID2, 821 X GRID1,GRID2,IMOUT,JMOUT) 822 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) 823 ID(1:25) = 0 824 ID(09) = 116 825 ID(10) = PETABND(NBND)+15. 826 ID(11) = PETABND(1)-15. 827 CALL OUTPUT(IOUTYP,IGET(032),LVLS(1,IGET(032)), 828 X GRID1,IMOUT,JMOUT) 829 ENDIF 830 C 831 IF (IGET(107).GT.0) THEN 832 CALL E2OUT(107,000,EGRID2,EGRID1, 833 X GRID1,GRID2,IMOUT,JMOUT) 834 DO J=JSTA,JEND 835 DO I=1,IMOUT 836 GRID1(I,J) = -1.*GRID1(I,J) 837 ENDDO 838 ENDDO 839 C 840 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) 841 C 842 DO J=JSTA,JEND 843 DO I=1,IMOUT 844 GRID1(I,J) = -1.*GRID1(I,J) 845 ENDDO 846 ENDDO 847 C 848 ID(1:25) = 0 849 ID(09) = 116 850 ID(10) = PETABND(NBND)+15. 851 ID(11) = PETABND(1)-15. 852 CALL OUTPUT(IOUTYP,IGET(107),LVLS(1,IGET(107)), 853 X GRID1,IMOUT,JMOUT) 854 ENDIF 855 ENDIF 856 C 857 C BOUNDARY LAYER LIFTING CONDENSATION PRESSURE AND HEIGHT. 858 C EGRID1 IS LCL PRESSURE. EGRID2 IS LCL HEIGHT. 859 C 860 IF ( (IGET(109).GT.0).OR.(IGET(110).GT.0) ) THEN 861 CALL CALLCL(PBND(1,1,1),TBND(1,1,1), 862 X QBND(1,1,1),EGRID1,EGRID2) 863 IF (IGET(109).GT.0) THEN 864 CALL E2OUT(109,000,EGRID2,EGRID1, 865 X GRID1,GRID2,IMOUT,JMOUT) 866 ID(1:25) = 0 867 CALL OUTPUT(IOUTYP,IGET(109),ILVL, 868 X GRID1,IMOUT,JMOUT) 869 ENDIF 870 IF (IGET(110).GT.0) THEN 871 CALL E2OUT(110,000,EGRID1,EGRID2, 872 X GRID1,GRID2,IMOUT,JMOUT) 873 ID(1:25) = 0 Page 14 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 874 CALL OUTPUT(IOUTYP,IGET(110),ILVL, 875 X GRID1,IMOUT,JMOUT) 876 ENDIF 877 ENDIF 878 C 879 C NGM BOUNDARY LAYER FIELDS. 880 C 881 IF ( (IGET(091).GT.0).OR.(IGET(092).GT.0).OR. 882 X (IGET(093).GT.0).OR.(IGET(094).GT.0).OR. 883 X (IGET(095).GT.0).OR.(IGET(095).GT.0).OR. 884 X (IGET(096).GT.0).OR.(IGET(097).GT.0).OR. 885 X (IGET(098).GT.0) ) THEN 886 C 887 C COMPUTE SIGMA 0.89671 AND 0.78483 TEMPERATURES 888 C INTERPOLATE LINEAR IN LOG P 889 IF (IGET(097).GT.0.OR.IGET(098).GT.0) THEN 890 DO J=JSTA,JEND 891 DO I=1,IM 892 P78483(I,J)=ALOG((PD(I,J)+PT)*0.78483) 893 P89671(I,J)=ALOG((PD(I,J)+PT)*0.89671) 894 ENDDO 895 ENDDO 896 !$omp parallel do 897 !$omp& private(fac1,fac2,pkl1,pku1,t78483,t89671) 898 DO L=2,LM 899 DO J=JSTA,JEND 900 DO I=1,IM 901 PKL1=0.5*(ALPINT(I,J,L)+ALPINT(I,J,L+1)) 902 PKU1=0.5*(ALPINT(I,J,L)+ALPINT(I,J,L-1)) 903 IF(P78483(I,J).LT.PKL1.AND.P78483(I,J).GT.PKU1)THEN 904 FAC1=(PKL1-P78483(I,J))/(PKL1-PKU1) 905 FAC2=(P78483(I,J)-PKU1)/(PKL1-PKU1) 906 T78483(I,J)=T(I,J,L)*FAC2+T(I,J,L-1)*FAC1 907 ENDIF 908 IF(P89671(I,J).LT.PKL1.AND.P89671(I,J).GT.PKU1)THEN 909 FAC1=(PKL1-P89671(I,J))/(PKL1-PKU1) 910 FAC2=(P89671(I,J)-PKU1)/(PKL1-PKU1) 911 T89671(I,J)=T(I,J,L)*FAC2+T(I,J,L-1)*FAC1 912 ENDIF 913 ENDDO 914 ENDDO 915 ENDDO 916 C 917 C SIGMA 0.89671 TEMPERATURE 918 IF (IGET(097).GT.0) THEN 919 ID(1:25) = 0 920 ISVALUE = 8967 921 ID(11) = ISVALUE 922 CALL E2OUT(097,000,T89671,EGRID2,GRID1,GRID2, 923 X IMOUT,JMOUT) 924 CALL OUTPUT(IOUTYP,IGET(097),LVLS(1,IGET(097)), 925 X GRID1,IMOUT,JMOUT) 926 ENDIF 927 C 928 C SIGMA 0.78483 TEMPERATURE 929 IF (IGET(098).GT.0) THEN 930 ID(1:25) = 0 Page 15 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 931 ISVALUE = 7848 932 ID(11) = ISVALUE 933 CALL E2OUT(098,000,T78483,EGRID2,GRID1,GRID2, 934 X IMOUT,JMOUT) 935 CALL OUTPUT(IOUTYP,IGET(098),LVLS(1,IGET(098)), 936 X GRID1,IMOUT,JMOUT) 937 ENDIF 938 ENDIF 939 C 940 C NGM SIGMA LAYER 0.98230 FIELDS. THESE FIELDS ARE 941 C THE FIRST ETA LAYER BOUNDARY LAYER FIELDS. AN 942 C EXCEPTION IS MADE FOR RELATIVE HUMIDITY ON THE LFM 943 C FORECAST GRID (026). ON THIS OUTPUT GRID THE SIGMA 944 C 0.98230 RELATIVE HUMIDITY IS THE MEAN RELATIVE HUMIDITY 945 C COMPUTED BY AVERAGING ETA BOUNDARY LAYER 1, 2, AND 3 946 C RELATIVE HUMIDITIES. THIS AD HOC APPROACH WAS TAKEN 947 C TO PRODUCE A MORE "LFM-LIKE" BOUNDARY LAYER RELATIVE 948 C HUMIDITY FIELD FOR THE AVIATION BRANCH OF MOD. 949 C 950 C THIS IS NOT A SATISFACTORY SOLUTION. THE SIGMA 951 C 0.98230 RELATIVE HUMIDITY SHOULD BE CONSISTENT 952 C WITH THE OTHER SIGMA 0.98230 FIELDS. 953 C 954 C 955 IF ( (IGET(091).GT.0).OR.(IGET(092).GT.0).OR. 956 X (IGET(093).GT.0).OR.(IGET(094).GT.0).OR. 957 X (IGET(095).GT.0).OR.(IGET(095).GT.0).OR. 958 X (IGET(096).GT.0) ) THEN 959 C 960 ID(1:25) = 0 961 ISVALUE = 9823 962 ID(11) = ISVALUE 963 C 964 C PRESSURE. 965 IF (IGET(091).GT.0) THEN 966 CALL E2OUT(091,000,PBND(1,1,1),EGRID2, 967 X GRID1,GRID2,IMOUT,JMOUT) 968 CALL OUTPUT(IOUTYP,IGET(091),LVLS(1,IGET(091)), 969 X GRID1,IMOUT,JMOUT) 970 ENDIF 971 C 972 C TEMPERATURE. 973 IF (IGET(092).GT.0) THEN 974 CALL E2OUT(092,000,TBND(1,1,1),EGRID2, 975 X GRID1,GRID2,IMOUT,JMOUT) 976 CALL OUTPUT(IOUTYP,IGET(092),LVLS(1,IGET(092)), 977 X GRID1,IMOUT,JMOUT) 978 ENDIF 979 C 980 C SPECIFIC HUMIDITY. 981 IF (IGET(093).GT.0) THEN 982 CALL E2OUT(093,000,QBND(1,1,1),EGRID2, 983 X GRID1,GRID2,IMOUT,JMOUT) 984 CALL BOUND(GRID1,H1M12,H99999,IMOUT,JMOUT) 985 CALL OUTPUT(IOUTYP,IGET(093),LVLS(1,IGET(093)), 986 X GRID1,IMOUT,JMOUT) 987 ENDIF Page 16 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 988 C 989 C RELATIVE HUMIDITY. 990 IF (IGET(094).GT.0) THEN 991 IF (KGTYPE.EQ.026) THEN 992 DO J=JSTA,JEND 993 DO I=1,IM 994 RH1D(I,J)=(RHBND(I,J,1)+RHBND(I,J,2)+ 995 & RHBND(I,J,3))*0.33333333 996 ENDDO 997 ENDDO 998 ELSE 999 DO J=JSTA,JEND 1000 DO I=1,IM 1001 RH1D(I,J)=RHBND(I,J,1) 1002 ENDDO 1003 ENDDO 1004 ENDIF 1005 CALL E2OUT(094,000,RH1D,EGRID2,GRID1,GRID2, 1006 X IMOUT,JMOUT) 1007 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 1008 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 1009 CALL OUTPUT(IOUTYP,IGET(094),LVLS(1,IGET(094)), 1010 X GRID1,IMOUT,JMOUT) 1011 ENDIF 1012 C 1013 C U AND/OR V WIND. 1014 IF ((IGET(095).GT.0).OR.(IGET(096).GT.0)) THEN 1015 CALL E2OUT(095,096,UBND(1,1,1),VBND(1,1,1), 1016 X GRID1,GRID2,IMOUT,JMOUT) 1017 IF (IGET(095).GT.0) CALL OUTPUT(IOUTYP,IGET(095), 1018 X LVLS(1,IGET(095)),GRID1,IMOUT,JMOUT) 1019 IF (IGET(096).GT.0) CALL OUTPUT(IOUTYP,IGET(096), 1020 X LVLS(1,IGET(096)),GRID2,IMOUT,JMOUT) 1021 ENDIF 1022 ENDIF 1023 ENDIF 1024 C 1025 C ENDIF FOR BOUNDARY LAYER BLOCK. 1026 C 1027 ENDIF 1028 C 1029 C 1030 C 1031 C ***BLOCK 6: MISCELLANEOUS LAYER MEAN LFM AND NGM FIELDS. 1032 C 1033 IF ( (IGET(066).GT.0).OR.(IGET(081).GT.0).OR. 1034 X (IGET(082).GT.0).OR.(IGET(104).GT.0).OR. 1035 X (IGET(099).GT.0).OR.(IGET(100).GT.0).OR. 1036 X (IGET(101).GT.0).OR.(IGET(102).GT.0).OR. 1037 X (IGET(103).GT.0) ) THEN 1038 C 1039 C LFM "MEAN" RELATIVE HUMIDITIES AND PRECIPITABLE WATER. 1040 C 1041 IF ( (IGET(066).GT.0).OR.(IGET(081).GT.0).OR. 1042 X (IGET(082).GT.0).OR.(IGET(104).GT.0) ) THEN 1043 CALL LFMFLD(RH3310,RH6610,RH3366,PW3310) 1044 ID(1:25) = 0 Page 17 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 1045 C 1046 C SIGMA 0.33-1.00 MEAN RELATIVE HUMIIDITY. 1047 IF (IGET(066).GT.0) THEN 1048 CALL E2OUT(066,000,RH3310,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1049 ID(1:25) = 0 1050 ID(10) = 33 1051 ID(11) = 100 1052 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 1053 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 1054 CALL OUTPUT(IOUTYP,IGET(066),LVLS(1,IGET(066)), 1055 X GRID1,IMOUT,JMOUT) 1056 ENDIF 1057 C 1058 C SIGMA 0.66-1.00 MEAN RELATIVE HUMIIDITY. 1059 IF (IGET(081).GT.0) THEN 1060 CALL E2OUT(081,000,RH6610,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1061 ID(1:25) = 0 1062 ID(10) = 67 1063 ID(11) = 100 1064 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 1065 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 1066 CALL OUTPUT(IOUTYP,IGET(081),LVLS(1,IGET(081)), 1067 X GRID1,IMOUT,JMOUT) 1068 ENDIF 1069 C 1070 C SIGMA 0.33-0.66 MEAN RELATIVE HUMIIDITY. 1071 IF (IGET(082).GT.0) THEN 1072 CALL E2OUT(082,000,RH3366,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1073 ID(1:25) = 0 1074 ID(10) = 33 1075 ID(11) = 67 1076 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 1077 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 1078 CALL OUTPUT(IOUTYP,IGET(082),LVLS(1,IGET(082)), 1079 X GRID1,IMOUT,JMOUT) 1080 ENDIF 1081 C 1082 C SIGMA 0.33-1.00 PRECIPITABLE WATER. 1083 IF (IGET(104).GT.0) THEN 1084 CALL E2OUT(104,000,PW3310,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 1085 ID(1:25) = 0 1086 ID(10) = 33 1087 ID(11) = 100 1088 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) 1089 CALL OUTPUT(IOUTYP,IGET(104),LVLS(1,IGET(104)), 1090 X GRID1,IMOUT,JMOUT) 1091 ENDIF 1092 ENDIF 1093 C 1094 C VARIOUS LAYER MEAN NGM SIGMA FIELDS. 1095 C 1096 IF ( (IGET(099).GT.0).OR.(IGET(100).GT.0).OR. 1097 X (IGET(101).GT.0).OR.(IGET(102).GT.0).OR. 1098 X (IGET(103).GT.0) ) THEN 1099 CALL NGMFLD2(RH4710,RH4796,RH1847,RH8498,QM8510) 1100 C 1101 C SIGMA 0.47191-1.00000 RELATIVE HUMIDITY. Page 18 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 1102 IF (IGET(099).GT.0) THEN 1103 CALL E2OUT(099,000,RH4710,EGRID2,GRID1,GRID2, 1104 X IMOUT,JMOUT) 1105 ID(1:25) = 0 1106 ID(10) = 47 1107 ID(11) = 100 1108 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 1109 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 1110 CALL OUTPUT(IOUTYP,IGET(099),LVLS(1,IGET(099)), 1111 X GRID1,IMOUT,JMOUT) 1112 ENDIF 1113 C 1114 C SIGMA 0.47191-0.96470 RELATIVE HUMIDITY. 1115 IF (IGET(100).GT.0) THEN 1116 CALL E2OUT(100,000,RH4796,EGRID2,GRID1,GRID2, 1117 X IMOUT,JMOUT) 1118 ID(1:25) = 0 1119 ID(10) = 47 1120 ID(11) = 96 1121 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 1122 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 1123 CALL OUTPUT(IOUTYP,IGET(100),LVLS(1,IGET(100)), 1124 X GRID1,IMOUT,JMOUT) 1125 ENDIF 1126 C 1127 C SIGMA 0.18019-0.47191 RELATIVE HUMIDITY. 1128 IF (IGET(101).GT.0) THEN 1129 CALL E2OUT(101,000,RH1847,EGRID2,GRID1,GRID2, 1130 X IMOUT,JMOUT) 1131 ID(1:25) = 0 1132 ID(10) = 18 1133 ID(11) = 47 1134 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 1135 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 1136 CALL OUTPUT(IOUTYP,IGET(101),LVLS(1,IGET(101)), 1137 X GRID1,IMOUT,JMOUT) 1138 ENDIF 1139 C 1140 C SIGMA 0.84368-0.98230 RELATIVE HUMIDITY. 1141 IF (IGET(102).GT.0) THEN 1142 CALL E2OUT(102,000,RH8498,EGRID2,GRID1,GRID2, 1143 X IMOUT,JMOUT) 1144 ID(1:25) = 0 1145 ID(10) = 84 1146 ID(11) = 98 1147 CALL SCLFLD(GRID1,H100,IMOUT,JMOUT) 1148 CALL BOUND(GRID1,H1,H100,IMOUT,JMOUT) 1149 CALL OUTPUT(IOUTYP,IGET(102),LVLS(1,IGET(102)), 1150 X GRID1,IMOUT,JMOUT) 1151 ENDIF 1152 C 1153 C SIGMA 0.85000-1.00000 MOISTURE CONVERGENCE. 1154 IF (IGET(103).GT.0) THEN 1155 CALL E2OUT(103,000,QM8510,EGRID2,GRID1,GRID2, 1156 X IMOUT,JMOUT) 1157 C CONVERT TO DIVERGENCE FOR GRIB 1158 DO J = JSTA,JEND Page 19 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f 1159 DO I = 1,IMOUT 1160 GRID1(I,J) = -1.0 * GRID1(I,J) 1161 ENDDO 1162 ENDDO 1163 ID(1:25) = 0 1164 ID(10) = 85 1165 ID(11) = 100 1166 CALL OUTPUT(IOUTYP,IGET(103),LVLS(1,IGET(103)), 1167 X GRID1,IMOUT,JMOUT) 1168 ENDIF 1169 ENDIF 1170 ENDIF 1171 C 1172 C END OF ROUTINE. 1173 C 1174 RETURN 1175 END ENTRY POINTS Name miscln_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 330 302 20 Label 617 393 50 Label 637 629 80 Label 683 669 A2 Param 41 R(4) 4 scalar A3 Param 41 R(4) 4 scalar A4 Param 41 R(4) 4 scalar AD05 Param 33 R(4) 4 scalar AETA Scalar 4 R(4) 4 1 45 COM ALOG Func 760 scalar 760,761 ALPINT Scalar 4 R(4) 4 3 4219074 COM 769,770 AMIN1 Func 634 scalar 634 AVBL Scalar 2 CHAR 20 1 250 COM BETA Param 37 R(4) 4 scalar 37 BNDLYR2 Subr 389 389 BOUND Subr 344 344,354,365,435,458,556,690,708,85 2,876,921,933,945,956,977,990,1003 ,1016 BTG Param 37 R(4) 4 scalar C2K Param 89 R(4) 4 scalar 602 CALCAPE Subr 685 685 CALDWP2 Subr 444 444 CALEKM Subr 515 515 CALEKM2 Subr 524 524 CALHEL Subr 157 157 CALLCL Subr 729 729 CALMCVG Subr 467 467 Page 20 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References CALMXW Subr 251 251 CALPOT2 Subr 217 217,421 CALTHTE Subr 670 670 CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CM1 Param 23 R(4) 4 scalar CM2 Param 23 R(4) 4 scalar CM3 Param 23 R(4) 4 scalar CP Param 34 R(4) 4 scalar 36,40 CPDR Param 36 R(4) 4 scalar CPGFU Scalar 8 R(4) 4 2 91719 COM CPGFV Scalar 3 R(4) 4 scalar COM CTLBLK Common 2 96 CURV Scalar 8 R(4) 4 2 91719 COM D00 Param 25 R(4) 4 scalar 344,365,556,690,708,956 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 D608 Param 29 R(4) 4 scalar DATSET Scalar 3 CHAR 6 scalar COM DDMPU Scalar 9 R(4) 4 2 91719 COM DDMPV Scalar 9 R(4) 4 2 91719 COM DEC Scalar 8 R(4) 4 1 250 COM DETA Scalar 4 R(4) 4 1 45 COM DFL Scalar 4 R(4) 4 1 46 COM DP10M Param 35 R(4) 4 scalar DT Scalar 4 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 DX Scalar 8 R(4) 4 2 91719 COM DY Scalar 3 R(4) 4 scalar COM DYNAM Common 2 3675964 E2OUT Subr 159 159,167,193,201,209,218,226,237,25 5,262,270,310,319,342,351,363,401, 411,422,432,446,456,475,495,516,52 5,544,554,567,582,638,688,700,732, 739,790,801,834,842,850,873,883,91 6,928,940,952,971,984,997,1010,102 3 EF4T Scalar 3 R(4) 4 scalar COM EGRID1 Local 108 R(4) 4 2 91719 217,218,421,422,444,446,468,472,47 5,566,567,582,624,634,638,664,671, 674,675,685,688,700,730,732,739 EGRID2 Local 108 R(4) 4 2 91719 159,193,201,209,218,237,255,262,31 0,342,351,363,401,411,422,432,446, 456,475,544,554,567,582,625,631,63 4,638,665,674,675,685,688,700,730, 732,739,790,801,834,842,850,873,91 6,928,940,952,971,984,997,1010,102 3 Page 21 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References ELDR Param 36 R(4) 4 scalar ELWV Param 36 R(4) 4 scalar 36 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 EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar ETA Scalar 4 R(4) 4 1 46 COM EXTRA Common 3 71506268 F Scalar 9 R(4) 4 2 91719 COM F4D Scalar 3 R(4) 4 scalar COM F4Q Scalar 3 R(4) 4 scalar COM F4Q2 Scalar 4 R(4) 4 1 45 COM FAC1 Local 772 R(4) 4 scalar 772,774,777,779 FAC2 Local 773 R(4) 4 scalar 773,774,778,779 FAD Scalar 9 R(4) 4 2 91719 COM FCP Scalar 8 R(4) 4 2 91719 COM FDIV Scalar 9 R(4) 4 2 91719 COM FDLVL Subr 300 300 FIELD Scalar 2 CHAR 20 1 250 COM FIELD1 Local 95 L(4) 4 scalar 484,488,494,649,653,659 FIELD2 Local 95 L(4) 4 scalar 485,491,494,650,656,659 FIRST Scalar 94 L(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 FRZLVL Subr 338 338 FRZLVL2 Subr 360 360 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 GRID1 Local 109 R(4) 4 2 0 159,164,167,172,193,196,201,204,20 9,212,218,221,226,229,237,240,255, 258,262,265,270,273,311,313,320,32 3,342,344,346,351,353,354,356,363, 365,367,402,404,412,414,422,425,43 3,434,435,437,446,449,457,458,460, 475,478,496,516,525,532,545,547,55 5,556,558,568,582,586,593,639,644, 689,690,696,701,704,708,712,721,73 3,736,740,743,790,793,801,804,835, 837,843,845,851,852,854,873,875,87 6,878,884,886,916,920,921,923,928, 932,933,935,940,944,945,947,952,95 6,958,971,976,977,979,984,989,990, 992,997,1002,1003,1005,1010,1015,1 016,1018,1023,1028,1035 GRID2 Local 109 R(4) 4 2 0 159,167,177,193,201,209,218,226,23 2,237,255,262,270,276,311,320,327, 342,351,363,402,412,422,433,446,45 7,475,496,516,525,537,545,555,568, 582,586,602,611,639,689,701,733,74 0,790,801,835,843,851,873,884,888, Page 22 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References 916,928,940,952,971,984,997,1010,1 023 H1 Param 26 R(4) 4 scalar 33,354,435,876,921,933,945,977,990 ,1003,1016 H100 Param 27 R(4) 4 scalar 353,354,434,435,875,876,920,921,93 2,933,944,945,976,977,989,990,1002 ,1003,1015,1016 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 458,852 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 344,365,458,556,624,625,664,665,69 0,708,852,956 HBM2 Scalar 3 R(4) 4 2 91719 COM HELI Local 107 R(4) 4 2 91719 157,159 HTFD Local 106 R(4) 4 1 6 145,304 HTM Scalar 5 R(4) 4 3 4127355 COM I Local 471 I(4) 4 scalar 471,472,585,586,601,602,623,624,62 5,633,634,663,664,665,673,674,675, 676,677,678,679,703,704,711,712,75 9,760,761,768,769,770,771,772,773, 774,776,777,778,779,861,862,863,86 8,869,1027,1028 IBESSL Scalar 2 I(4) 4 scalar COM ICNT Scalar 7 I(4) 4 1 1024 COM ID Scalar 11 I(4) 4 1 25 COM 160,161,162,168,169,170,173,174,17 5,194,202,210,219,227,230,238,256, 263,271,274,303,305,343,352,364,39 4,395,396,606,607,608,609,640,641, 642,691,692,693,694,716,717,718,71 9,734,741,787,789,798,800,828,830, 912,917,918,919,929,930,931,941,94 2,943,953,954,955,973,974,975,986, 987,988,999,1000,1001,1012,1013,10 14,1031,1032,1033 IDAT Scalar 3 I(4) 4 1 3 COM IDENT Scalar 10 I(4) 4 1 250 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 IFD Local 289 I(4) 4 scalar 289,291,294,297,302,304,309,310,31 3,319,322,323,326,327 IFILV Scalar 11 I(4) 4 1 250 COM IFLAG Scalar 2 I(4) 4 scalar COM IGET Scalar 8 I(4) 4 1 250 COM 156,158,163,166,171,172,176,177,18 5,186,187,188,192,195,200,203,208, Page 23 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References 211,216,220,225,228,229,231,232,23 6,239,249,250,254,257,261,264,269, 272,273,275,276,283,284,290,291,29 3,294,296,297,308,309,312,313,318, 321,322,323,325,326,327,337,341,34 5,350,355,359,366,374,375,376,377, 378,379,380,381,382,383,384,385,38 6,399,400,403,404,409,410,413,419, 420,424,430,431,436,442,443,448,45 4,455,459,465,466,477,487,488,490, 491,511,512,520,521,529,530,531,53 2,534,535,536,537,542,543,546,552, 553,557,563,564,592,610,621,643,65 2,653,655,656,687,695,699,720,728, 731,735,738,742,749,750,751,752,75 3,757,786,792,797,803,823,824,825, 826,833,836,841,844,849,853,858,87 7,882,885,886,887,888,901,902,903, 904,905,909,910,915,922,927,934,93 9,946,951,957,964,965,966,970,978, 983,991,996,1004,1009,1017,1022,10 34 IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 IHRST Scalar 3 I(4) 4 scalar COM ILVL Local 735 I(4) 4 scalar 735,742 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,14,96,97,98,99, 100,101,102,103,104,105,106,107,10 8,110,217,421,471,623,633,663,673, 759,768,861,868 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 109,159,164,167,172,177,193,196,20 1,204,209,212,218,221,226,229,232, 237,240,255,258,262,265,270,273,27 6,311,313,320,323,327,342,344,346, 351,353,354,356,363,365,367,402,40 4,412,414,423,425,433,434,435,437, 447,449,457,458,460,476,478,496,51 6,525,532,537,545,547,555,556,558, 568,582,585,593,601,611,639,644,68 9,690,696,701,703,708,711,721,733, 736,740,743,791,793,802,804,835,83 7,843,845,851,852,854,874,875,876, 878,884,886,888,916,920,921,923,92 8,932,933,935,940,944,945,947,952, 956,958,972,976,977,979,985,989,99 0,992,998,1002,1003,1005,1011,1015 ,1016,1018,1024,1027,1035 IMT Param 4 I(4) 4 scalar IMX Param 3 I(4) 4 scalar INPES Param 6 I(4) 4 scalar 11,14 IOFFS Scalar 2 I(4) 4 scalar COM IOUNIT Common 4 20 IOUT Scalar 3 I(4) 4 scalar COM Page 24 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References IOUTYP Scalar 8 I(4) 4 scalar COM 163,171,176,195,203,211,220,228,23 1,239,257,264,272,275,312,322,326, 345,355,366,403,413,424,436,448,45 9,477,531,536,546,557,592,610,643, 695,720,735,742,792,803,836,844,85 3,877,885,887,922,934,946,957,978, 991,1004,1017,1034 IQ Scalar 9 I(4) 4 1 250 COM IS Scalar 9 I(4) 4 1 250 COM ISMFUL Scalar 9 I(4) 4 1 250 COM ISMOUT Scalar 10 I(4) 4 1 250 COM ISMSTG Scalar 9 I(4) 4 1 250 COM ISVALUE Local 304 I(4) 4 scalar 304,305,788,789,799,800,829,830 ITAIL Param 11 I(4) 4 scalar ITYPE Local 288 I(4) 4 scalar 288,291,294,297,300,660,685 IUP Scalar 7 I(4) 4 scalar COM J Local 470 I(4) 4 scalar 470,472,584,586,600,602,622,624,62 5,632,634,662,664,665,672,674,675, 676,677,678,679,702,704,710,712,75 8,760,761,767,769,770,771,772,773, 774,776,777,778,779,860,862,863,86 7,869,1026,1028 JAM Param 6 I(4) 4 scalar 6 JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 470,584,600,622,632,662,672,702,71 0,758,767,860,867,1026 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 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,96,97,98,99, 100,101,102,103,104,105,106,107,10 8,110,217,421 JM2 Param 6 I(4) 4 scalar JMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 109,159,164,167,172,177,193,196,20 1,204,209,212,218,221,226,229,232, 237,240,255,258,262,265,270,273,27 6,311,313,320,323,327,342,344,346, 351,353,354,356,363,365,367,402,40 4,412,414,423,425,433,434,435,437, 447,449,457,458,460,476,478,496,51 6,525,532,537,545,547,555,556,558, 568,582,593,611,639,644,689,690,69 6,701,708,721,733,736,740,743,791, 793,802,804,835,837,843,845,851,85 2,854,874,875,876,878,884,886,888, 916,920,921,923,928,932,933,935,94 0,944,945,947,952,956,958,972,976, 977,979,985,989,990,992,998,1002,1 003,1005,1011,1015,1016,1018,1024, 1035 JMT Param 4 I(4) 4 scalar JMX Param 3 I(4) 4 scalar Page 25 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References JNPES Param 6 I(4) 4 scalar 12,15 JSTA Scalar 6 I(4) 4 scalar COM 470,584,600,622,632,662,672,702,71 0,758,767,860,867,1026 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM JSTA_M2 Scalar 7 I(4) 4 scalar COM JTAIL Param 12 I(4) 4 scalar KGTYPE Scalar 8 I(4) 4 scalar COM 583,599,859 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 L Local 766 I(4) 4 scalar 766,769,770,774,779 LB Param 9 I(4) 4 scalar LB2 Local 96 I(4) 4 2 91719 676,685 LBND Local 393 I(4) 4 scalar 393,395,396,400,401,404,410,411,41 3,420,421,424,431,432,436,443,444, 445,448,455,456,459,466,467,468,47 7,488,491,495,512,513,521,522,530, 532,535,537,543,544,546,553,554,55 7,564,565,566,592,629,630,631,669, 670,671,676,677,678,679 LCNTRL Scalar 5 I(4) 4 scalar COM LFMFLD Subr 911 911 LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LIST Scalar 3 I(4) 4 scalar COM LM Param 3 I(4) 4 scalar 3,4,5,8,19,766 LM1 Param 8 I(4) 4 scalar LP1 Param 8 I(4) 4 scalar 4 LSM Param 3 I(4) 4 scalar LUNCO Scalar 5 I(4) 4 scalar COM LUNLL Scalar 5 I(4) 4 scalar COM LUNOUT Scalar 5 I(4) 4 scalar COM LVLBND Local 96 I(4) 4 3 550314 390,676 LVLS Scalar 10 I(4) 4 2 11250 COM 163,172,177,195,203,211,220,229,23 2,239,257,264,273,276,291,294,297, 309,313,322,323,326,327,345,355,36 6,400,404,410,413,420,424,431,436, 443,448,455,459,466,477,488,491,51 2,521,530,532,535,537,543,546,553, 557,564,592,610,643,653,656,695,72 0,792,803,836,844,853,877,886,888, 922,934,946,957,978,991,1004,1017, 1034 MASKS Common 2 34853220 Page 26 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References MAXWP Local 110 R(4) 4 2 91719 251,255 MAXWU Local 110 R(4) 4 2 91719 251,270 MAXWV Local 110 R(4) 4 2 91719 251,270 MAXWZ Local 110 R(4) 4 2 91719 251,262 ME Scalar 6 I(4) 4 scalar COM MISCLN Subr 1 MXFLD Param 3 I(4) 4 scalar 2,8,9,10,11 MXLVL Param 3 I(4) 4 scalar 10 NBC Scalar 3 I(4) 4 scalar COM NBND Param 88 I(4) 4 scalar 96,99,100,101,107,393,629,641,669, 693,718 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 NFD Param 88 I(4) 4 scalar 106,289,302 NFLD Scalar 8 I(4) 4 scalar COM NGMFLD2 Subr 967 967 NINC Param 7 I(4) 4 scalar NINT Func 304 scalar 304,395,396 NORTH Local 93 L(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 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 NWORD Param 19 I(4) 4 scalar 19 NWORD3 Param 19 I(4) 4 scalar NWORD6 Param 19 I(4) 4 scalar OLDRD Scalar 94 L(4) 4 scalar COM OMGBND Local 99 R(4) 4 3 550314 390,544 ONEPS Param 39 R(4) 4 scalar OPTIONS Common 2 24 OTLFT2 Subr 565 565,630 OUTPUT Subr 163 163,171,176,195,203,211,220,228,23 1,239,257,264,272,275,312,322,326, 345,355,366,403,413,424,436,448,45 9,477,531,536,546,557,592,610,643, 695,720,735,742,792,803,836,844,85 3,877,885,887,922,934,946,957,978, 991,1004,1017,1034 OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar P1D Local 97 R(4) 4 2 91719 189,193,217,677,685 P78483 Local 102 R(4) 4 2 91719 760,771,772,773 P89671 Local 102 R(4) 4 2 91719 761,776,777,778 PARALLEL Common 9 8240 PBND Local 100 R(4) 4 3 550314 389,401,421,444,565,630,670,677,72 9,834 PD Scalar 3 R(4) 4 2 91719 COM 760,761 Page 27 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References PDSL Scalar 6 R(4) 4 2 91719 COM PDVP1 Scalar 6 R(4) 4 2 91719 COM PETABND Local 107 R(4) 4 1 6 150,395,396,641,642,693,694,718,71 9 PFSL Scalar 7 R(4) 4 2 91719 COM PINT Scalar 4 R(4) 4 3 4219074 COM PKL1 Local 769 R(4) 4 scalar 769,771,772,773,776,777,778 PKU1 Local 770 R(4) 4 scalar 770,771,772,773,776,777,778 PQ0 Param 41 R(4) 4 scalar PSLP Scalar 6 R(4) 4 2 91719 COM PT Scalar 3 R(4) 4 scalar COM 760,761 PTSL Scalar 7 R(4) 4 2 91719 COM PW3310 Local 105 R(4) 4 2 91719 911,952 PWTBND Local 99 R(4) 4 3 550314 390,554 Q Scalar 5 R(4) 4 3 4127355 COM Q1D Local 97 R(4) 4 2 91719 679,685 QBND Local 100 R(4) 4 3 550314 389,444,456,467,566,631,671,679,73 0,850 QM8510 Local 103 R(4) 4 2 91719 967,1023 QSL Scalar 9 R(4) 4 2 91719 COM R Scalar 3 R(4) 4 scalar COM RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDETA Scalar 4 R(4) 4 1 45 COM RDOCP Param 40 R(4) 4 scalar RES Scalar 3 R(4) 4 2 91719 COM RESTRT Scalar 94 L(4) 4 scalar COM RG Param 36 R(4) 4 scalar RH1847 Local 104 R(4) 4 2 91719 967,997 RH1D Local 98 R(4) 4 2 91719 338,351,360,862,869,873 RH3310 Local 105 R(4) 4 2 91719 911,916 RH3366 Local 105 R(4) 4 2 91719 911,940 RH4710 Local 103 R(4) 4 2 91719 967,971 RH4796 Local 104 R(4) 4 2 91719 967,984 RH6610 Local 105 R(4) 4 2 91719 911,928 RH8498 Local 103 R(4) 4 2 91719 967,1010 RHBND Local 101 R(4) 4 3 550314 389,432,862,863,869 RITE2 Scalar 5 L(4) 4 scalar COM RITEHD Scalar 5 L(4) 4 scalar COM ROG Param 32 R(4) 4 scalar RQSTFLD Common 7 64130 RTD Param 34 R(4) 4 scalar RUN Scalar 94 L(4) 4 scalar COM SATDEL Scalar 2 R(4) 4 scalar COM SCLFLD Subr 353 353,434,875,920,932,944,976,989,10 02,1015 SHR1D Local 98 R(4) 4 2 91719 189,237 SICE Scalar 3 R(4) 4 2 91719 COM SIGBND Local 107 R(4) 4 1 6 149 SIGMA Scalar 94 L(4) 4 scalar COM SLP Scalar 5 R(4) 4 2 91719 COM SM Scalar 3 R(4) 4 2 91719 COM SMALL Param 88 R(4) 4 scalar SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar Page 28 Source Listing MISCLN 2014-12-17 20:48 Symbol Table MISCLN.f Name Object Declared Type Bytes Dimen Elements Attributes References SPVAL Scalar 2 R(4) 4 scalar COM STDOUT Scalar 2 I(4) 4 scalar COM 513,522 STDRD Scalar 94 L(4) 4 scalar COM SVALUE Scalar 8 R(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM 774,779 T1D Local 97 R(4) 4 2 91719 189,209,217,678,685 T500 Scalar 6 R(4) 4 2 91719 COM T6D Local 106 R(4) 4 3 550314 300,310 T78483 Local 102 R(4) 4 2 91719 774,801 T89671 Local 102 R(4) 4 2 91719 779,790 TBND Local 100 R(4) 4 3 550314 389,411,421,445,565,630,670,678,72 9,842 TFRZ Param 39 R(4) 4 scalar TRPAUS Subr 189 189 TSL Scalar 8 R(4) 4 2 91719 COM TTV Scalar 8 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM U1D Local 97 R(4) 4 2 91719 189,226,515,516,524,525 U6D Local 106 R(4) 4 3 550314 300,319 UBND Local 101 R(4) 4 3 550314 389,467,495,883 UST Local 104 R(4) 4 2 91719 157,167 V Scalar 5 R(4) 4 3 4127355 COM V1D Local 97 R(4) 4 2 91719 189,226,515,516,524,525 V6D Local 106 R(4) 4 3 550314 300,319 VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VBND Local 101 R(4) 4 3 550314 389,468,495,883 VRBLS Common 2 67138308 VST Local 104 R(4) 4 2 91719 157,167 VTM Scalar 5 R(4) 4 3 4127355 COM WPDAR Scalar 8 R(4) 4 2 91719 COM Z1000 Scalar 9 R(4) 4 2 91719 COM Z1D Local 98 R(4) 4 2 91719 189,201,338,342,360,363 ZINT Scalar 4 R(4) 4 3 4219074 COM Page 29 Source Listing MISCLN 2014-12-17 20:48 Subprograms/Common Blocks MISCLN.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CTLBLK Common 2 96 DYNAM Common 2 3675964 EXTRA Common 3 71506268 IOUNIT Common 4 20 MASKS Common 2 34853220 MISCLN Subr 1 OPTIONS Common 2 24 PARALLEL Common 9 8240 RQSTFLD Common 7 64130 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__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ Page 30 Source Listing MISCLN 2014-12-17 20:48 MISCLN.f -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 : MISCLN.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100