Page 1 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 1 SUBROUTINE CLDRAD(IMOUT,JMOUT) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CLDRAD POST SNDING/CLOUD/RADTN FIELDS 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-08-30 6 C 7 C ABSTRACT: THIS ROUTINE COMPUTES/POSTS SOUNDING, CLOUD 8 C RELATED, AND RADIATION FIELDS. UNDER THE HEADING OF 9 C SOUNDING FIELDS FALL THE THREE ETA MODEL LIFTED INDICES, 10 C CAPE, CIN, AND TOTAL COLUMN PRECIPITABLE WATER. 11 C 12 C THE THREE ETA MODEL LIFTED INDICES DIFFER ONLY IN THE 13 C DEFINITION OF THE PARCEL TO LIFT. ONE LIFTS PARCELS FROM 14 C THE LOWEST ABOVE GROUND ETA LAYER. ANOTHER LIFTS MEAN 15 C PARCELS FROM ANY OF NBND BOUNDARY LAYERS (SEE SUBROUTINE 16 C BNDLYR2). THE FINAL TYPE OF LIFTED INDEX IS A BEST LIFTED 17 C INDEX BASED ON THE NBND BOUNDARY LAYER LIFTED INDICES. 18 C 19 C TWO TYPES OF CAPE/CIN ARE AVAILABLE. ONE IS BASED ON PARCELS 20 C IN THE LOWEST ETA LAYER ABOVE GROUND. THE OTHER IS BASED 21 C ON A LAYER MEAN PARCEL IN THE N-TH BOUNDARY LAYER ABOVE 22 C THE GROUND. SEE SUBROUTINE CALCAPE FOR DETAILS. 23 C 24 C THE CLOUD FRACTION AND LIQUID CLOUD WATER FIELDS ARE DIRECTLY 25 C FROM THE MODEL WITH MINIMAL POST PROCESSING. THE LIQUID 26 C CLOUD WATER, 3-D CLOUD FRACTION, AND TEMPERATURE TENDENCIES 27 C DUE TO PRECIPITATION ARE NOT POSTED IN THIS ROUTINE. SEE 28 C SUBROUTINE ETAFLD2 FOR THESE FIELDS. LIFTING CONDENSATION 29 C LEVEL HEIGHT AND PRESSURE ARE COMPUTED AND POSTED IN 30 C SUBROUTINE MISCLN. 31 C 32 C THE RADIATION FIELDS POSTED BY THIS ROUTINE ARE THOSE COMPUTED 33 C DIRECTLY IN THE MODEL. 34 C 35 C PROGRAM HISTORY LOG: 36 C 93-08-30 RUSS TREADON 37 C 94-08-04 MICHAEL BALDWIN - ADDED OUTPUT OF INSTANTANEOUS SFC 38 C FLUXES OF NET SW AND LW DOWN RADIATION 39 C 97-04-25 MICHAEL BALDWIN - FIX PDS FOR PRECIPITABLE WATER 40 C 97-04-29 GEOFF MANIKIN - MOVED CLOUD TOP TEMPS CALCULATION 41 C TO THIS SUBROUTINE. CHANGED METHOD 42 C OF DETERMINING WHERE CLOUD BASE AND 43 C TOP ARE FOUND AND ADDED HEIGHT OPTION 44 C FOR TOP AND BASE. 45 C 98-04-29 GEOFF MANIKIN - CHANGED VALUE FOR CLOUD BASE/TOP PRESSURES 46 C AND HEIGHTS FROM SPVAL TO -500 47 C 98-06-15 T BLACK - CONVERSION FROM 1-D TO 2-D 48 C 98-07-17 MIKE BALDWIN - REMOVED LABL84 49 C 00-01-04 JIM TUCCILLO - MPI VERSION 50 C 00-02-22 GEOFF MANIKIN - CHANGED VALUE FOR CLOUD BASE/TOP PRESSURES 51 C AND HEIGHTS FROM SPVAL TO -500 (WAS NOT IN 52 C PREVIOUS IBM VERSION) 53 C 54 C 55 C 56 C USAGE: CALL CLDRAD(IMOUT,JMOUT) 57 C INPUT ARGUMENT LIST: Page 2 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 58 C IMOUT - FIRST DIMENSION OF OUTPUT GRID. 59 C JMOUT - SECOND DIMENSION OF OUTPUT GRID. 60 C 61 C OUTPUT ARGUMENT LIST: 62 C NONE 63 C 64 C OUTPUT FILES: 65 C NONE 66 C 67 C SUBPROGRAMS CALLED: 68 C UTILITIES: 69 C NONE 70 C LIBRARY: 71 C COMMON - RQSTFLD 72 C EXTRA 73 C MAPOT 74 C PHYS 75 C LOOPS 76 C MASKS 77 C VRBLS 78 C CNVCLD 79 C CTLBLK 80 C ACMCLD 81 C ACMRDS 82 C ACMRDL 83 C OPTIONS 84 C CLDWTR 85 C 86 C ATTRIBUTES: 87 C LANGUAGE: FORTRAN 88 C MACHINE : CRAY C-90 89 C$$$ 90 C 91 C 92 C INCLUDE GRID DIMENSIONS. SET/DERIVE OTHER PARAMETERS. 93 C 94 INCLUDE "parmeta" 112 INCLUDE "parmout" 116 INCLUDE "params" 158 INCLUDE "parm.tbl" 159 C 160 C SET CELSIUS TO KELVIN CONVERSION. 162 PARAMETER (C2K=273.15) 163 C 164 C DECLARE VARIABLES. 165 C 166 LOGICAL RUN,FIRST,RESTRT,SIGMA,OLDRD,STDRD 167 LOGICAL NEED(IM,JM) 168 INTEGER L1D(IM,JM) 169 REAL EGRID1(IM,JM),EGRID2(IM,JM) 170 REAL GRID1(IMOUT,JMOUT),GRID2(IMOUT,JMOUT),CLDTP(IM,JM), 171 & CLDTZ(IM,JM),CLDBP(IM,JM),CLDBZ(IM,JM),CLDTT(IM,JM) 172 C 173 C INCLUDE COMMON BLOCKS. 174 INCLUDE "RQSTFLD.comm" 187 INCLUDE "EXTRA.comm" 199 INCLUDE "MAPOT.comm" Page 3 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 212 INCLUDE "MASKS.comm" 220 INCLUDE "PHYS1.comm" 241 INCLUDE "CNVCLD.comm" 247 INCLUDE "LOOPS.comm" 255 INCLUDE "VRBLS.comm" 262 INCLUDE "CTLBLK.comm" 277 INCLUDE "ACMCLD.comm" 284 INCLUDE "ACMRDS.comm" 290 INCLUDE "ACMRDL.comm" 296 INCLUDE "OPTIONS.comm" 300 INCLUDE "CLDWTR.comm" 301 C 302 C 303 C************************************************************************* 304 C START CLDRAD HERE. 305 C 306 C*** BLOCK 1. SOUNDING DERIVED FIELDS. 307 C 308 C ETA SURFACE TO 500MB LIFTED INDEX. TO BE CONSISTENT WITH THE 309 310 C LFM AND NGM POSTING WE ADD 273.15 TO THE LIFTED INDEX 311 C NOTE: 25 JUNE 1993, RUSS TREADON. 312 C ON THE LFM FORECAST GRID (026) WE POST 313 C THE FIRST BOUNDARY LAYER LIFTED INDEX. 314 C SEE SUBROUTINE MISCLN. 315 C 316 C THE BEST (SIX LAYER) AND BOUNDARY LAYER LIFTED INDICES ARE 317 C COMPUTED AND POSTED IN SUBROUTINE MISCLN. 323 print*,'Start of CLDRAD' 324 C 325 IF ( (IGET(030).GT.0).AND.(KGTYPE.NE.026) ) THEN 326 DO J=JSTA,JEND 327 DO I=1,IM 328 EGRID1(I,J) = SPVAL 329 ENDDO 330 ENDDO 331 C 332 CALL OTLIFT2(T500,EGRID1) 333 C 334 DO J=JSTA,JEND 335 DO I=1,IM 336 EGRID1(I,J) = EGRID1(I,J) + C2K 337 ENDDO 338 ENDDO 339 C 340 CALL E2OUT(030,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 341 ID(1:25)=0 342 ID(10) =50 343 ID(11) =100 344 CALL OUTPUT(IOUTYP,IGET(030),LVLS(1,IGET(030)), 345 X GRID1,IMOUT,JMOUT) 346 ENDIF 347 C 348 C SOUNDING DERIVED AREA INTEGRATED ENERGIES - CAPE AND CIN. 349 IF ((IGET(032).GT.0).OR.(IGET(107).GT.0)) THEN 350 IF ( (LVLS(1,IGET(032)).GT.0) .OR. 351 X (LVLS(1,IGET(107)).GT.0) ) THEN Page 4 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 352 ITYPE = 1 353 CALL CALCAPE(ITYPE,P1D,T1D,Q1D,L1D,EGRID1,EGRID2) 354 C 355 C CONVECTIVE AVAILABLE POTENTIAL ENERGY. 356 IF (IGET(032).GT.0) THEN 357 CALL E2OUT(032,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 358 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) 359 ID(1:25)=0 360 CALL OUTPUT(IOUTYP,IGET(032),LVLS(1,IGET(032)), 361 X GRID1,IMOUT,JMOUT) 362 ENDIF 363 C 364 C CONVECTIVE INHIBITION. 365 IF (IGET(107).GT.0) THEN 366 CALL E2OUT(107,000,EGRID2,EGRID1,GRID1,GRID2,IMOUT,JMOUT) 367 DO J=JSTA,JEND 368 DO I=1,IMOUT 369 GRID1(I,J) = -1.*GRID1(I,J) 370 ENDDO 371 ENDDO 372 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) 373 DO J=JSTA,JEND 374 DO I=1,IMOUT 375 GRID1(I,J) = -1.*GRID1(I,J) 376 ENDDO 377 ENDDO 378 ID(1:25)=0 379 CALL OUTPUT(IOUTYP,IGET(107),LVLS(1,IGET(107)), 380 X GRID1,IMOUT,JMOUT) 381 ENDIF 382 ENDIF 383 ENDIF 384 C 385 C TOTAL COLUMN PRECIPITABLE WATER. 386 IF (IGET(080).GT.0) THEN 387 CALL CALPW(EGRID1) 388 CALL E2OUT(080,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 389 ID(1:25)=0 390 CALL BOUND(GRID1,D00,H99999,IMOUT,JMOUT) 391 CALL OUTPUT(IOUTYP,IGET(080),LVLS(1,IGET(080)), 392 X GRID1,IMOUT,JMOUT) 393 ENDIF 394 C 395 C 396 C 397 C 398 C*** BLOCK 2. 2-D CLOUD FIELDS. 399 C 400 C LOW CLOUD FRACTION. 401 IF (IGET(037).GT.0) THEN 402 CALL E2OUT(037,000,CFRACL,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 403 ID(1:25)=0 404 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 405 CALL OUTPUT(IOUTYP,IGET(037),LVLS(1,IGET(037)), 406 X GRID1,IMOUT,JMOUT) 407 ENDIF 408 C Page 5 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 409 C MIDDLE CLOUD FRACTION. 410 IF (IGET(038).GT.0) THEN 411 CALL E2OUT(038,000,CFRACM,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 412 ID(1:25)=0 413 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 414 CALL OUTPUT(IOUTYP,IGET(038),LVLS(1,IGET(038)), 415 X GRID1,IMOUT,JMOUT) 416 ENDIF 417 C 418 C HIGH CLOUD FRACTION. 419 IF (IGET(039).GT.0) THEN 420 CALL E2OUT(039,000,CFRACH,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 421 ID(1:25)=0 422 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 423 CALL OUTPUT(IOUTYP,IGET(039),LVLS(1,IGET(039)), 424 X GRID1,IMOUT,JMOUT) 425 ENDIF 426 C 427 C TOTAL CLOUD FRACTION (INSTANTANEOUS). 428 IF (IGET(161).GT.0) THEN 429 DO J=JSTA,JEND 430 DO I=1,IM 431 EGRID1(I,J)=AMAX1(CFRACL(I,J), 432 1 AMAX1(CFRACM(I,J),CFRACH(I,J))) 433 ENDDO 434 ENDDO 435 CALL E2OUT(161,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 436 ID(1:25)=0 437 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 438 CALL OUTPUT(IOUTYP,IGET(161),LVLS(1,IGET(161)), 439 X GRID1,IMOUT,JMOUT) 440 ENDIF 441 C 442 C TIME AVERAGED TOTAL CLOUD FRACTION. 443 IF (IGET(144).GT.0) THEN 444 DO 40 J=JSTA,JEND 445 DO 40 I=1,IM 446 ISUM = NCFRST(I,J)+NCFRCV(I,J) 447 IF (ISUM.GT.0) THEN 448 EGRID1(I,J)=(ACFRST(I,J)+ACFRCV(I,J))/ISUM 449 ELSE 450 EGRID1(I,J) = D00 451 ENDIF 452 40 CONTINUE 453 C 454 CALL E2OUT(144,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 455 ID(1:25)= 0 456 IFHR = NTSD/TSPH+0.5 457 ITCLOD = INT(TCLOD) 458 IFINCR = MOD(IFHR,ITCLOD) 459 ID(19) = IFHR 460 ID(20) = 3 461 IF (IFINCR.EQ.0) THEN 462 ID(18) = IFHR-ITCLOD 463 ELSE 464 ID(18) = IFHR-IFINCR 465 ENDIF Page 6 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 466 IF (ID(18).LT.0) ID(18) = 0 467 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 468 CALL OUTPUT(IOUTYP,IGET(144),LVLS(1,IGET(144)), 469 X GRID1,IMOUT,JMOUT) 470 ENDIF 471 C 472 C TIME AVERAGED STRATIFORM CLOUD FRACTION. 473 IF (IGET(139).GT.0) THEN 474 DO 50 J=JSTA,JEND 475 DO 50 I=1,IM 476 IF (NCFRST(I,J).GT.0) THEN 477 EGRID1(I,J) = ACFRST(I,J)/NCFRST(I,J) 478 ELSE 479 EGRID1(I,J) = D00 480 ENDIF 481 50 CONTINUE 482 C 483 CALL E2OUT(139,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 484 ID(1:25)=0 485 IFHR = NTSD/TSPH+0.5 486 ITCLOD = INT(TCLOD) 487 IFINCR = MOD(IFHR,ITCLOD) 488 ID(19) = IFHR 489 ID(20) = 3 490 IF (IFINCR.EQ.0) THEN 491 ID(18) = IFHR-ITCLOD 492 ELSE 493 ID(18) = IFHR-IFINCR 494 ENDIF 495 IF (ID(18).LT.0) ID(18) = 0 496 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 497 CALL OUTPUT(IOUTYP,IGET(139),LVLS(1,IGET(139)), 498 X GRID1,IMOUT,JMOUT) 499 ENDIF 500 C 501 C TIME AVERAGED CONVECTIVE CLOUD FRACTION. 502 IF (IGET(143).GT.0) THEN 503 DO 60 J=JSTA,JEND 504 DO 60 I=1,IM 505 IF (NCFRCV(I,J).GT.0) THEN 506 EGRID1(I,J) = ACFRCV(I,J)/NCFRCV(I,J) 507 ELSE 508 EGRID1(I,J) = D00 509 ENDIF 510 60 CONTINUE 511 C 512 CALL E2OUT(143,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 513 ID(1:25)=0 514 IFHR = NTSD/TSPH+0.5 515 ITCLOD = INT(TCLOD) 516 IFINCR = MOD(IFHR,ITCLOD) 517 ID(19) = IFHR 518 ID(20) = 3 519 IF (IFINCR.EQ.0) THEN 520 ID(18) = IFHR-ITCLOD 521 ELSE 522 ID(18) = IFHR-IFINCR Page 7 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 523 ENDIF 524 IF (ID(18).LT.0) ID(18) = 0 525 CALL SCLFLD(GRID1,100.,IMOUT,JMOUT) 526 CALL OUTPUT(IOUTYP,IGET(143),LVLS(1,IGET(143)), 527 X GRID1,IMOUT,JMOUT) 528 ENDIF 529 C 530 C CLOUD BASE FIELDS 531 IF ((IGET(148).GT.0) .OR. (IGET(178).GT.0)) THEN 532 CLIMIT =1.0E-06 533 DO J=JSTA,JEND 534 DO I=1,IM 535 NEED(I,J)=.TRUE. 536 CLDBP(I,J) = -50000. 537 CLDBZ(I,J) = -5000. 538 ENDDO 539 ENDDO 540 541 !$omp parallel do 542 !$omp& private(cbot,lev,llmh) 543 DO 70 J=JSTA,JEND 544 DO 70 I=1,IM 545 LLMH=LMH(I,J) 546 CBOT = 500 547 DO L=LLMH,1,-1 548 C GSM 549 C START AT THE FIRST LAYER ABOVE GROUND, AND FIND THE 550 C FIRST LAYER WITH A VALUE OF CLOUD WATER GREATER THAN 551 C THE SIGNIFICANT LIMIT (VALUE DESIGNATED BY Q. ZHAO). 552 C THIS LAYER WILL BE THE CLOUD BOTTOM UNLESS THE BOTTOM 553 C OF THE CONVECTIVE CLOUD (HBOT) IS FOUND BELOW IN WHICH 554 C CASE HBOT BECOMES THE CLOUD BASE LAYER. 555 IF (CWM(I,J,L).GT.CLIMIT.AND.NEED(I,J)) THEN 556 CBOT=L 557 IF (HBOT(I,J).GT.CBOT) THEN 558 CBOT = HBOT(I,J) 559 ENDIF 560 NEED(I,J)=.FALSE. 561 ENDIF 562 ENDDO 563 C 564 IF (CBOT.EQ.500.) THEN 565 CLDBP(I,J) = -50000. 566 CLDBZ(I,J) = -5000. 567 ELSE IF (CBOT.EQ.LM) THEN 568 CLDBP(I,J) = AETA1(INT(CBOT))*PDSL(I,J)+PT1 569 CLDBZ(I,J) = ZINT(I,J,LM) 570 ELSE 571 LEV = CBOT 572 CLDBP(I,J) = AETA1(INT(LEV))*PDSL(I,J)+PT1 573 CLDBZ(I,J) = HTM(I,J,LEV+1)*T(I,J,LEV+1) 574 1 *(Q(I,J,LEV+1)*D608+H1)*ROG* 575 2 (LOG(PINT(I,J,LEV+1))-LOG(CLDBP(I,J))) 576 3 +ZINT(I,J,LEV+1) 577 ENDIF 578 70 CONTINUE 579 C CLOUD BOTTOM PRESSURE Page 8 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 580 IF (IGET(148).GT.0) THEN 581 CALL E2OUT(148,000,CLDBP,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 582 ID(1:25)=0 583 CALL OUTPUT(IOUTYP,IGET(148),LVLS(1,IGET(148)), 584 X GRID1,IMOUT,JMOUT) 585 ENDIF 586 587 C CLOUD BOTTOM HEIGHT 588 IF (IGET(178).GT.0) THEN 589 CALL E2OUT(148,000,CLDBZ,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 590 ID(1:25)=0 591 CALL OUTPUT(IOUTYP,IGET(178),LVLS(1,IGET(178)), 592 X GRID1,IMOUT,JMOUT) 593 ENDIF 594 ENDIF 595 C 596 C CLOUD TOP FIELDS 597 IF ((IGET(149).GT.0) .OR. (IGET(179).GT.0) .OR. 598 X (IGET(168).GT.0)) THEN 599 CLIMIT =1.0E-06 600 DO J=JSTA,JEND 601 DO I=1,IM 602 NEED(I,J)=.TRUE. 603 ENDDO 604 ENDDO 605 C GSM 606 C START AT THE TOP OF THE ATMOSPHERE. FIND THE 607 C FIRST LAYER WITH A VALUE OF CLOUD WATER GREATER THAN 608 C THE SIGNIFICANT LIMIT (VALUE DESIGNATED BY Q. ZHAO). 609 C THIS LAYER WILL BE THE CLOUD TOP UNLESS THE TOP 610 C OF THE CONVECTIVE CLOUD (HTOP) IS FOUND ABOVE IN WHICH 611 C CASE HTOP BECOMES THE CLOUD TOP LAYER. 612 613 !$omp parallel do 614 !$omp& private(ctop,lev,lmhk) 615 DO 80 J=JSTA,JEND 616 DO 80 I=1,IM 617 CTOP = 0. 618 LMHK=LMH(I,J) 619 DO L=1,LMHK 620 ! print *, i,j,l,CWM(I,J,L),CLIMIT,"-----" 621 IF (CWM(I,J,L).GT.CLIMIT.AND.NEED(I,J)) THEN 622 CTOP=L 623 IF (HTOP(I,J).LT.CTOP) THEN 624 CTOP = HTOP(I,J) 625 ENDIF 626 NEED(I,J)=.FALSE. 627 ENDIF 628 ENDDO 629 C 630 IF (CTOP.EQ.0.)THEN 631 CLDTP(I,J) = -50000. 632 CLDTZ(I,J) = -5000. 633 LMHK=LMH(I,J) 634 CLDTT(I,J) = T(I,J,LMHK) 635 ELSE IF (CTOP.EQ.LM) THEN 636 CLDTP(I,J) = AETA1(INT(LM))*PDSL(I,J)+PT1 Page 9 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 637 CLDTZ(I,J) = ZINT(I,J,LM) 638 CLDTT(I,J) = T(I,J,LM) 639 ELSE 640 LEV = CTOP 641 CLDTP(I,J) = AETA1(INT(LEV))*PDSL(I,J)+PT1 642 CLDTZ(I,J) = HTM(I,J,LEV+1)*T(I,J,LEV+1) 643 1 *(Q(I,J,LEV+1)*D608+H1)*ROG* 644 2 (LOG(PINT(I,J,LEV+1))-LOG(CLDTP(I,J))) 645 3 +ZINT(I,J,LEV+1) 646 CLDTT(I,J) = T(I,J,LEV) 647 ENDIF 648 ! print *, i,j,CLDTP(I,J), CTOP, NEED(I,J) 649 80 CONTINUE 650 C CLOUD TOP PRESSURE 651 IF (IGET(149).GT.0) THEN 652 CALL E2OUT(149,000,CLDTP,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 653 ID(1:25)=0 654 CALL OUTPUT(IOUTYP,IGET(149),LVLS(1,IGET(149)), 655 X GRID1,IMOUT,JMOUT) 656 ENDIF 657 658 C CLOUD TOP HEIGHT 659 IF (IGET(179).GT.0) THEN 660 CALL E2OUT(179,000,CLDTZ,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 661 ID(1:25)=0 662 CALL OUTPUT(IOUTYP,IGET(179),LVLS(1,IGET(179)), 663 X GRID1,IMOUT,JMOUT) 664 ENDIF 665 666 C CLOUD TOP TEMPS 667 IF (IGET(168).GT.0) THEN 668 CALL E2OUT(168,000,CLDTT,EGRID2, 669 X GRID1,GRID2,IMOUT,JMOUT) 670 ID(1:25)=0 671 CALL OUTPUT(IOUTYP,IGET(168),LVLS(1,IGET(168)), 672 X GRID1,IMOUT,JMOUT) 673 ENDIF 674 ENDIF 675 676 C*** BLOCK 3. RADIATION FIELDS. 677 C 678 C 679 C TIME AVERAGED SURFACE SHORT WAVE INCOMING RADIATION. 680 IF (IGET(126).GT.0) THEN 681 IF(ARDSW.GT.0.)THEN 682 RRNUM=1./ARDSW 683 ELSE 684 RRNUM=0. 685 ENDIF 686 DO J=JSTA,JEND 687 DO I=1,IM 688 EGRID1(I,J) = ASWIN(I,J)*RRNUM 689 ENDDO 690 ENDDO 691 CALL E2OUT(126,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 692 ID(1:25)=0 693 IFHR = NTSD/TSPH+0.5 Page 10 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 694 ITRDSW = INT(TRDSW) 695 IFINCR = MOD(IFHR,ITRDSW) 696 ID(19) = IFHR 697 ID(20) = 3 698 IF (IFINCR.EQ.0) THEN 699 ID(18) = IFHR-ITRDSW 700 ELSE 701 ID(18) = IFHR-IFINCR 702 ENDIF 703 IF (ID(18).LT.0) ID(18) = 0 704 CALL OUTPUT(IOUTYP,IGET(126),LVLS(1,IGET(126)), 705 X GRID1,IMOUT,JMOUT) 706 ENDIF 707 C 708 C TIME AVERAGED SURFACE LONG WAVE INCOMING RADIATION. 709 IF (IGET(127).GT.0) THEN 710 IF(ARDLW.GT.0.)THEN 711 RRNUM=1./ARDLW 712 ELSE 713 RRNUM=0. 714 ENDIF 715 DO J=JSTA,JEND 716 DO I=1,IM 717 EGRID1(I,J) = ALWIN(I,J)*RRNUM 718 ENDDO 719 ENDDO 720 CALL E2OUT(127,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 721 ID(1:25)=0 722 IFHR = NTSD/TSPH+0.5 723 ITRDLW = INT(TRDLW) 724 IFINCR = MOD(IFHR,ITRDLW) 725 ID(19) = IFHR 726 ID(20) = 3 727 IF (IFINCR.EQ.0) THEN 728 ID(18) = IFHR-ITRDLW 729 ELSE 730 ID(18) = IFHR-IFINCR 731 ENDIF 732 IF (ID(18).LT.0) ID(18) = 0 733 CALL OUTPUT(IOUTYP,IGET(127),LVLS(1,IGET(127)), 734 X GRID1,IMOUT,JMOUT) 735 ENDIF 736 C 737 C TIME AVERAGED SURFACE SHORT WAVE OUTGOING RADIATION. 738 IF (IGET(128).GT.0) THEN 739 IF(ARDSW.GT.0.)THEN 740 RRNUM=1./ARDSW 741 ELSE 742 RRNUM=0. 743 ENDIF 744 DO J=JSTA,JEND 745 DO I=1,IM 746 EGRID1(I,J) = -1.0*ASWOUT(I,J)*RRNUM 747 ENDDO 748 ENDDO 749 CALL E2OUT(128,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 750 ID(1:25)=0 Page 11 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 751 IFHR = NTSD/TSPH+0.5 752 ITRDSW = INT(TRDSW) 753 IFINCR = MOD(IFHR,ITRDSW) 754 ID(19) = IFHR 755 ID(20) = 3 756 IF (IFINCR.EQ.0) THEN 757 ID(18) = IFHR-ITRDSW 758 ELSE 759 ID(18) = IFHR-IFINCR 760 ENDIF 761 IF (ID(18).LT.0) ID(18) = 0 762 CALL OUTPUT(IOUTYP,IGET(128),LVLS(1,IGET(128)), 763 X GRID1,IMOUT,JMOUT) 764 ENDIF 765 C 766 C TIME AVERAGED SURFACE LONG WAVE OUTGOING RADIATION. 767 IF (IGET(129).GT.0) THEN 768 IF(ARDLW.GT.0.)THEN 769 RRNUM=1./ARDLW 770 ELSE 771 RRNUM=0. 772 ENDIF 773 DO J=JSTA,JEND 774 DO I=1,IM 775 EGRID1(I,J) = -1.0*ALWOUT(I,J)*RRNUM 776 ENDDO 777 ENDDO 778 CALL E2OUT(129,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 779 ID(1:25)=0 780 IFHR = NTSD/TSPH+0.5 781 ITRDLW = INT(TRDLW) 782 IFINCR = MOD(IFHR,ITRDLW) 783 ID(19) = IFHR 784 ID(20) = 3 785 IF (IFINCR.EQ.0) THEN 786 ID(18) = IFHR-ITRDLW 787 ELSE 788 ID(18) = IFHR-IFINCR 789 ENDIF 790 IF (ID(18).LT.0) ID(18) = 0 791 CALL OUTPUT(IOUTYP,IGET(129),LVLS(1,IGET(129)), 792 X GRID1,IMOUT,JMOUT) 793 ENDIF 794 C 795 C TIME AVERAGED TOP OF THE ATMOSPHERE SHORT WAVE RADIATION. 796 IF (IGET(130).GT.0) THEN 797 IF(ARDSW.GT.0.)THEN 798 RRNUM=1./ARDSW 799 ELSE 800 RRNUM=0. 801 ENDIF 802 DO J=JSTA,JEND 803 DO I=1,IM 804 EGRID1(I,J) = ASWTOA(I,J)*RRNUM 805 ENDDO 806 ENDDO 807 CALL E2OUT(130,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) Page 12 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 808 ID(1:25)=0 809 IFHR = NTSD/TSPH+0.5 810 ITRDSW = INT(TRDSW) 811 IFINCR = MOD(IFHR,ITRDSW) 812 ID(19) = IFHR 813 ID(20) = 3 814 IF (IFINCR.EQ.0) THEN 815 ID(18) = IFHR-ITRDSW 816 ELSE 817 ID(18) = IFHR-IFINCR 818 ENDIF 819 IF (ID(18).LT.0) ID(18) = 0 820 CALL OUTPUT(IOUTYP,IGET(130),LVLS(1,IGET(130)), 821 X GRID1,IMOUT,JMOUT) 822 ENDIF 823 C 824 C TIME AVERAGED TOP OF THE ATMOSPHERE LONG WAVE RADIATION. 825 IF (IGET(131).GT.0) THEN 826 IF(ARDLW.GT.0.)THEN 827 RRNUM=1./ARDLW 828 ELSE 829 RRNUM=0. 830 ENDIF 831 DO J=JSTA,JEND 832 DO I=1,IM 833 EGRID1(I,J) = ALWTOA(I,J)*RRNUM 834 ENDDO 835 ENDDO 836 CALL E2OUT(131,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 837 ID(1:25)=0 838 IFHR = NTSD/TSPH+0.5 839 ITRDLW = INT(TRDLW) 840 IFINCR = MOD(IFHR,ITRDLW) 841 ID(19) = IFHR 842 ID(20) = 3 843 IF (IFINCR.EQ.0) THEN 844 ID(18) = IFHR-ITRDLW 845 ELSE 846 ID(18) = IFHR-IFINCR 847 ENDIF 848 IF (ID(18).LT.0) ID(18) = 0 849 CALL OUTPUT(IOUTYP,IGET(131),LVLS(1,IGET(131)), 850 X GRID1,IMOUT,JMOUT) 851 ENDIF 852 C 853 C CURRENT INCOMING SW RADIATION AT THE SURFACE. 854 IF (IGET(156).GT.0) THEN 855 DO J=JSTA,JEND 856 DO I=1,IM 857 IF(CZMEAN(I,J).GT.1.E-6) THEN 858 FACTRS=CZEN(I,J)/CZMEAN(I,J) 859 ELSE 860 FACTRS=0.0 861 ENDIF 862 EGRID1(I,J)=HBM2(I,J)*RSWIN(I,J)*FACTRS 863 ENDDO 864 ENDDO Page 13 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 865 C 866 CALL E2OUT(156,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 867 ID(1:25)=0 868 CALL OUTPUT(IOUTYP,IGET(156),LVLS(1,IGET(156)), 869 X GRID1,IMOUT,JMOUT) 870 ENDIF 871 C 872 C CURRENT INCOMING LW RADIATION AT THE SURFACE. 873 IF (IGET(157).GT.0) THEN 874 DO J=JSTA,JEND 875 DO I=1,IM 876 IF(SIGT4(I,J).GT.0.0) THEN 877 LLMH=LMH(I,J) 878 TLMH=T(I,J,LLMH) 879 FACTRL=5.67E-8*TLMH*TLMH*TLMH*TLMH/SIGT4(I,J) 880 ELSE 881 FACTRL=0.0 882 ENDIF 883 EGRID1(I,J)=HBM2(I,J)*RLWIN(I,J)*FACTRL 884 ENDDO 885 ENDDO 886 C 887 CALL E2OUT(157,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 888 ID(1:25)=0 889 CALL OUTPUT(IOUTYP,IGET(157),LVLS(1,IGET(157)), 890 X GRID1,IMOUT,JMOUT) 891 ENDIF 892 C 893 C CURRENT OUTGOING SW RADIATION AT THE SURFACE. 894 IF (IGET(141).GT.0) THEN 895 DO J=JSTA,JEND 896 DO I=1,IM 897 IF(CZMEAN(I,J).GT.1.E-6) THEN 898 FACTRS=CZEN(I,J)/CZMEAN(I,J) 899 ELSE 900 FACTRS=0.0 901 ENDIF 902 EGRID1(I,J)=HBM2(I,J)*RSWOUT(I,J)*FACTRS 903 ENDDO 904 ENDDO 905 C 906 CALL E2OUT(141,000,EGRID1,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 907 ID(1:25)=0 908 CALL OUTPUT(IOUTYP,IGET(141),LVLS(1,IGET(141)), 909 X GRID1,IMOUT,JMOUT) 910 ENDIF 911 C 912 C CURRENT OUTGOING LW RADIATION AT THE SURFACE. 913 IF (IGET(142).GT.0) THEN 914 CALL E2OUT(142,000,RADOT,EGRID2,GRID1,GRID2,IMOUT,JMOUT) 915 ID(1:25)=0 916 CALL OUTPUT(IOUTYP,IGET(142),LVLS(1,IGET(142)), 917 X GRID1,IMOUT,JMOUT) 918 ENDIF 919 C 920 C END OF ROUTINE. 921 print*,'End of CLDRAD' Page 14 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f 922 C 923 RETURN 924 END ENTRY POINTS Name cldrad_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 40 Label 272 264,265 50 Label 301 294,295 60 Label 330 323,324 70 Label 398 363,364 80 Label 469 435,436 A2 Param 41 R(4) 4 scalar A3 Param 41 R(4) 4 scalar A4 Param 41 R(4) 4 scalar ACFRCV Scalar 3 R(4) 4 2 91719 COM 268,326 ACFRST Scalar 4 R(4) 4 2 91719 COM 268,297 ACMCLD Common 2 1467512 ACMRDL Common 2 2201268 ACMRDS Common 2 2201268 AD05 Param 33 R(4) 4 scalar AETA Scalar 8 R(4) 4 1 45 COM AETA1 Scalar 7 R(4) 4 1 45 COM 388,392,456,461 ALBASE Scalar 14 R(4) 4 2 91719 COM ALBEDO Scalar 15 R(4) 4 2 91719 COM ALPINT Scalar 4 R(4) 4 3 4219074 COM ALSL Scalar 7 R(4) 4 1 39 COM ALWIN Scalar 4 R(4) 4 2 91719 COM 537 ALWOUT Scalar 4 R(4) 4 2 91719 COM 595 ALWTOA Scalar 4 R(4) 4 2 91719 COM 653 AMAX1 Func 251 scalar 251,252 ARDLW Scalar 3 R(4) 4 scalar COM 530,531,588,589,646,647 ARDSW Scalar 3 R(4) 4 scalar COM 501,502,559,560,617,618 ASWIN Scalar 4 R(4) 4 2 91719 COM 508 ASWOUT Scalar 4 R(4) 4 2 91719 COM 566 ASWTOA Scalar 4 R(4) 4 2 91719 COM 624 AVBL Scalar 2 CHAR 20 1 250 COM BETA Param 37 R(4) 4 scalar 37 BOUND Subr 178 178,192,210 BTG Param 37 R(4) 4 scalar C2K Param 100 R(4) 4 scalar 156 CALCAPE Subr 173 173 CALPW Subr 207 207 CAPA Param 28 R(4) 4 scalar 35 CBOT Local 366 R(4) 4 scalar 366,376,377,378,384,387,388,391 CFRACH Scalar 3 R(4) 4 2 91719 COM 240,252 CFRACL Scalar 3 R(4) 4 2 91719 COM 222,251 CFRACM Scalar 3 R(4) 4 2 91719 COM 231,252 Page 15 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References CFT0 Param 33 R(4) 4 scalar CI Scalar 4 R(4) 4 scalar COM CLDBP Local 109 R(4) 4 2 91719 356,385,388,392,395,401 CLDBZ Local 109 R(4) 4 2 91719 357,386,389,393,409 CLDRAD Subr 1 CLDTP Local 108 R(4) 4 2 91719 451,456,461,464,472 CLDTT Local 109 R(4) 4 2 91719 454,458,466,488 CLDTZ Local 109 R(4) 4 2 91719 452,457,462,480 CLDWTR Common 2 17610408 CLIMIT Local 352 R(4) 4 scalar 352,375,419,441 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 CNVCLD Common 2 1467504 CP Param 34 R(4) 4 scalar 36,40 CPDR Param 36 R(4) 4 scalar CPHI0 Scalar 9 R(4) 4 scalar COM CS Scalar 4 R(4) 4 scalar COM CTLBLK Common 2 96 CTOP Local 437 R(4) 4 scalar 437,442,443,444,450,455,460 CUPPT Scalar 3 R(4) 4 2 91719 COM CWM Scalar 3 R(4) 4 3 4127355 COM 375,441 CZEN Scalar 12 R(4) 4 2 91719 COM 678,718 CZMEAN Scalar 16 R(4) 4 2 91719 COM 677,678,717,718 D00 Param 25 R(4) 4 scalar 178,192,210,270,299,328 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 394,463 DATSET Scalar 3 CHAR 6 scalar COM DEC Scalar 8 R(4) 4 1 250 COM DETA Scalar 8 R(4) 4 1 45 COM DETA1 Scalar 7 R(4) 4 1 45 COM DFRLG Scalar 7 R(4) 4 1 46 COM DI Scalar 4 R(4) 4 scalar COM DISLP Scalar 7 R(4) 4 scalar COM DLMD Scalar 6 R(4) 4 scalar COM DP10M Param 35 R(4) 4 scalar DP30 Scalar 6 R(4) 4 scalar COM DPHD Scalar 6 R(4) 4 scalar COM DS Scalar 4 R(4) 4 scalar COM DT Scalar 4 R(4) 4 scalar COM DTD Scalar 3 R(4) 4 scalar COM DTQ2 Scalar 3 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 E2OUT Subr 160 160,177,186,208,222,231,240,255,27 4,303,332,401,409,472,480,488,511, 540,569,598,627,656,686,707,726,73 Page 16 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References 4 EGRID1 Local 107 R(4) 4 2 91719 148,152,156,160,173,177,186,207,20 8,251,255,268,270,274,297,299,303, 326,328,332,508,511,537,540,566,56 9,595,598,624,627,653,656,682,686, 703,707,722,726 EGRID2 Local 107 R(4) 4 2 91719 160,173,177,186,208,222,231,240,25 5,274,303,332,401,409,472,480,488, 511,540,569,598,627,656,686,707,72 6,734 ELDR Param 36 R(4) 4 scalar ELWDWN Scalar 13 R(4) 4 2 91719 COM ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar EPSR Scalar 11 R(4) 4 2 91719 COM ERAD Param 34 R(4) 4 scalar ERLAM0 Scalar 9 R(4) 4 scalar COM ETA Scalar 8 R(4) 4 1 46 COM EXTRA Common 3 71506268 FACTRL Local 699 R(4) 4 scalar 699,701,703 FACTRS Local 678 R(4) 4 scalar 678,680,682,718,720,722 FIELD Scalar 2 CHAR 20 1 250 COM FIRST Scalar 104 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 FSL Scalar 8 R(4) 4 2 91719 COM G Param 32 R(4) 4 scalar 32 GFFC Scalar 14 R(4) 4 2 91719 COM GI Param 32 R(4) 4 scalar 36 GLAT Scalar 12 R(4) 4 2 91719 COM GLON Scalar 12 R(4) 4 2 91719 COM GRID1 Local 108 R(4) 4 2 0 160,165,177,178,181,186,189,192,19 5,200,208,210,212,222,224,226,231, 233,235,240,242,244,255,257,259,27 4,287,289,303,316,318,332,345,347, 401,404,409,412,472,475,480,483,48 9,492,511,525,540,554,569,583,598, 612,627,641,656,670,686,689,707,71 0,726,729,734,737 GRID2 Local 108 R(4) 4 2 0 160,177,186,208,222,231,240,255,27 4,303,332,401,409,472,480,489,511, 540,569,598,627,656,686,707,726,73 4 H1 Param 26 R(4) 4 scalar 33,394,463 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 Page 17 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 178,192,210 HBM2 Scalar 3 R(4) 4 2 91719 COM 682,703,722 HBOT Scalar 13 R(4) 4 2 91719 COM 377,378 HDAC Scalar 16 R(4) 4 2 91719 COM HDACV Scalar 16 R(4) 4 2 91719 COM HTM Scalar 5 R(4) 4 3 4127355 COM 393,462 HTOP Scalar 13 R(4) 4 2 91719 COM 443,444 I Local 147 I(4) 4 scalar 147,148,155,156,188,189,194,195,25 0,251,252,265,266,268,270,295,296, 297,299,324,325,326,328,354,355,35 6,357,364,365,375,377,378,380,385, 386,388,389,392,393,394,395,396,42 1,422,436,438,441,443,444,446,451, 452,453,454,456,457,458,461,462,46 3,464,465,466,507,508,536,537,565, 566,594,595,623,624,652,653,676,67 7,678,682,695,696,697,698,699,703, 716,717,718,722 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 161,162,163,179,198,209,223,232,24 1,256,275,279,280,282,284,286,304, 308,309,311,313,315,333,337,338,34 0,342,344,402,410,473,481,490,512, 516,517,519,521,523,541,545,546,54 8,550,552,570,574,575,577,579,581, 599,603,604,606,608,610,628,632,63 3,635,637,639,657,661,662,664,666, 668,687,708,727,735 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 IFHR Local 276 I(4) 4 scalar 276,278,279,282,284,305,307,308,31 1,313,334,336,337,340,342,513,515, 516,519,521,542,544,545,548,550,57 1,573,574,577,579,600,602,603,606, 608,629,631,632,635,637,658,660,66 1,664,666 IFILV Scalar 11 I(4) 4 1 250 COM IFINCR Local 278 I(4) 4 scalar 278,281,284,307,310,313,336,339,34 2,515,518,521,544,547,550,573,576, 579,602,605,608,631,634,637,660,66 3,666 IFLAG Scalar 2 I(4) 4 scalar COM IGET Scalar 8 I(4) 4 1 250 COM 145,164,169,170,171,176,180,185,19 9,206,211,221,225,230,234,239,243, 248,258,263,288,293,317,322,346,35 1,400,403,408,411,417,418,471,474, Page 18 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References 479,482,487,491,500,524,529,553,55 8,582,587,611,616,640,645,669,674, 688,693,709,714,728,733,736 IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 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 IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,12,13,14,15,16, 105,106,107,108,109,147,155,250,26 5,295,324,354,364,421,436,507,536, 565,594,623,652,676,695,716 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 108,160,165,177,178,181,186,188,19 2,194,200,208,210,212,222,224,226, 231,233,235,240,242,244,255,257,25 9,274,287,289,303,316,318,332,345, 347,401,404,409,412,472,475,480,48 3,489,492,511,525,540,554,569,583, 598,612,627,641,656,670,686,689,70 7,710,726,729,734,737 IMT Param 4 I(4) 4 scalar IMX Param 3 I(4) 4 scalar INPES Param 6 I(4) 4 scalar 11,14 INT Func 277 scalar 277,306,335,388,392,456,461,514,54 3,572,601,630,659 IOFFS Scalar 2 I(4) 4 scalar COM IOUT Scalar 3 I(4) 4 scalar COM IOUTYP Scalar 8 I(4) 4 scalar COM 164,180,199,211,225,234,243,258,28 8,317,346,403,411,474,482,491,524, 553,582,611,640,669,688,709,728,73 6 IQ Scalar 9 I(4) 4 1 250 COM IS Scalar 9 I(4) 4 1 250 COM ISHDE Scalar 4 I(4) 4 1 61 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 ISUM Local 266 I(4) 4 scalar 266,267,268 ITAIL Param 11 I(4) 4 scalar ITB Param 1 I(4) 4 scalar 8,18 ITBQ Param 1 I(4) 4 scalar 9,18 ITCLOD Local 277 I(4) 4 scalar 277,278,282,306,307,311,335,336,34 0 ITRDLW Local 543 I(4) 4 scalar 543,544,548,601,602,606,659,660,66 4 ITRDSW Local 514 I(4) 4 scalar 514,515,519,572,573,577,630,631,63 5 ITYPE Local 172 I(4) 4 scalar 172,173 IUNIT Scalar 10 I(4) 4 scalar COM IUP Scalar 7 I(4) 4 scalar COM IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IXM Scalar 3 I(4) 4 scalar COM Page 19 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References IYM Scalar 3 I(4) 4 scalar COM J Local 146 I(4) 4 scalar 146,148,154,156,187,189,193,195,24 9,251,252,264,266,268,270,294,296, 297,299,323,325,326,328,353,355,35 6,357,363,365,375,377,378,380,385, 386,388,389,392,393,394,395,396,42 0,422,435,438,441,443,444,446,451, 452,453,454,456,457,458,461,462,46 3,464,465,466,506,508,535,537,564, 566,593,595,622,624,651,653,675,67 7,678,682,694,696,697,698,699,703, 715,717,718,722 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 146,154,187,193,249,264,294,323,35 3,363,420,435,506,535,564,593,622, 651,675,694,715 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,11,12,13,14,15,16, 105,106,107,108,109 JM2 Param 6 I(4) 4 scalar JMOUT Dummy 1 I(4) 4 scalar ARG,INOUT 108,160,165,177,178,181,186,192,20 0,208,210,212,222,224,226,231,233, 235,240,242,244,255,257,259,274,28 7,289,303,316,318,332,345,347,401, 404,409,412,472,475,480,483,489,49 2,511,525,540,554,569,583,598,612, 627,641,656,670,686,689,707,710,72 6,729,734,737 JMT Param 4 I(4) 4 scalar JMX Param 3 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 146,154,187,193,249,264,294,323,35 3,363,420,435,506,535,564,593,622, 651,675,694,715 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 JTB Param 1 I(4) 4 scalar 8,18 JTBQ Param 1 I(4) 4 scalar 18 KGTYPE Scalar 8 I(4) 4 scalar COM 145 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 Page 20 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 KTM Scalar 3 I(4) 4 scalar COM KUNIT Scalar 10 I(4) 4 scalar COM L Local 367 I(4) 4 scalar 367,375,376,439,441,442 L1D Local 106 I(4) 4 2 91719 173 LB Param 9 I(4) 4 scalar LC Scalar 3 I(4) 4 2 91719 COM LEV Local 391 I(4) 4 scalar 391,392,393,394,395,396,460,461,46 2,463,464,465,466 LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LIST Scalar 3 I(4) 4 scalar COM LLMH Local 365 I(4) 4 scalar 365,367,697,698 LM Param 3 I(4) 4 scalar 3,5,7,8,19,387,389,455,456,457,458 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 365,438,453,697 LMHK Local 438 I(4) 4 scalar 438,439,453,454 LMV Scalar 5 I(4) 4 2 91719 COM LOG Func 395 scalar 395,464 LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4,7,8 LSL Scalar 3 I(4) 4 scalar COM LSM Param 3 I(4) 4 scalar 7 LUNIT Scalar 10 I(4) 4 scalar COM LVLS Scalar 10 I(4) 4 2 11250 COM 164,170,171,180,199,211,225,234,24 3,258,288,317,346,403,411,474,482, 491,524,553,582,611,640,669,688,70 9,728,736 MAPOT Common 2 1448 MASKS Common 2 34853220 ME Scalar 6 I(4) 4 scalar COM MOD Func 278 scalar 278,307,336,515,544,573,602,631,66 0 MUNIT Scalar 10 I(4) 4 scalar COM 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 NBOCO Scalar 4 I(4) 4 scalar COM NCFRCV Scalar 3 I(4) 4 2 91719 COM 266,325,326 NCFRST Scalar 4 I(4) 4 2 91719 COM 266,296,297 NCLOD Scalar 3 I(4) 4 scalar COM NCP Scalar 4 I(4) 4 scalar COM NDDAMP Scalar 4 I(4) 4 scalar COM NEED Local 105 L(4) 4 2 91719 355,375,380,422,441,446 NFCST Scalar 3 I(4) 4 scalar COM NFLD Scalar 8 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 Page 21 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References NRADS Scalar 5 I(4) 4 scalar COM NRDLW Scalar 3 I(4) 4 scalar COM NRDSW Scalar 3 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 276,305,334,513,542,571,600,629,65 8 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 104 L(4) 4 scalar COM ONEPS Param 39 R(4) 4 scalar OPTIONS Common 2 24 OTLIFT2 Subr 152 152 OUTPUT Subr 164 164,180,199,211,225,234,243,258,28 8,317,346,403,411,474,482,491,524, 553,582,611,640,669,688,709,728,73 6 OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar P1D Local 173 R(4) 4 scalar 173 PARALLEL Common 9 8240 PD Scalar 3 R(4) 4 2 91719 COM PDSL Scalar 6 R(4) 4 2 91719 COM 388,392,456,461 PDVP1 Scalar 6 R(4) 4 2 91719 COM PFSL Scalar 7 R(4) 4 2 91719 COM PHYS Common 2 7690036 PINT Scalar 4 R(4) 4 3 4219074 COM 395,464 PL Scalar 5 R(4) 4 scalar COM PLQ Scalar 6 R(4) 4 scalar COM PQ0 Param 41 R(4) 4 scalar PSLP Scalar 6 R(4) 4 2 91719 COM PT Scalar 5 R(4) 4 scalar COM PT1 Scalar 3 R(4) 4 scalar COM 388,392,456,461 PTBL Scalar 18 R(4) 4 2 10184 COM PTSL Scalar 7 R(4) 4 2 91719 COM Q Scalar 5 R(4) 4 3 4127355 COM 394,463 Q1D Local 173 R(4) 4 scalar 173 QS0 Scalar 8 R(4) 4 1 134 COM QSL Scalar 9 R(4) 4 2 91719 COM R Scalar 5 R(4) 4 scalar COM RADIN Scalar 11 R(4) 4 2 91719 COM RADOT Scalar 11 R(4) 4 2 91719 COM 734 RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDOCP Param 40 R(4) 4 scalar RDP Scalar 5 R(4) 4 scalar COM RDPQ Scalar 6 R(4) 4 scalar COM RDQ Scalar 5 R(4) 4 scalar COM RDTH Scalar 5 R(4) 4 scalar COM RDTHE Scalar 5 R(4) 4 scalar COM Page 22 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References RDTHEQ Scalar 6 R(4) 4 scalar COM RES Scalar 3 R(4) 4 2 91719 COM RESTRT Scalar 104 L(4) 4 scalar COM RG Param 36 R(4) 4 scalar RITE2 Scalar 5 L(4) 4 scalar COM RITEHD Scalar 5 L(4) 4 scalar COM RLWIN Scalar 3 R(4) 4 2 91719 COM 703 RLWOUT Scalar 3 R(4) 4 2 91719 COM RLWTOA Scalar 3 R(4) 4 2 91719 COM RMXSNAL Scalar 11 R(4) 4 2 91719 COM ROG Param 32 R(4) 4 scalar 394,463 ROI Scalar 4 R(4) 4 scalar COM ROS Scalar 4 R(4) 4 scalar COM RQSTFLD Common 7 64130 RRNUM Local 502 R(4) 4 scalar 502,504,508,531,533,537,560,562,56 6,589,591,595,618,620,624,647,649, 653 RSWIN Scalar 3 R(4) 4 2 91719 COM 682 RSWOUT Scalar 3 R(4) 4 2 91719 COM 722 RSWTOA Scalar 3 R(4) 4 2 91719 COM RTD Param 34 R(4) 4 scalar RUN Scalar 104 L(4) 4 scalar COM SATDEL Scalar 2 R(4) 4 scalar COM SBD Scalar 6 R(4) 4 scalar COM SCLFLD Subr 224 224,233,242,257,287,316,345 SICE Scalar 3 R(4) 4 2 91719 COM SIGMA Scalar 104 L(4) 4 scalar COM SIGT4 Scalar 16 R(4) 4 2 91719 COM 696,699 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 SPHI0 Scalar 9 R(4) 4 scalar COM SPL Scalar 7 R(4) 4 1 39 COM SPVAL Scalar 2 R(4) 4 scalar COM 148 SQS Scalar 8 R(4) 4 1 134 COM SR Scalar 3 R(4) 4 2 91719 COM SST Scalar 14 R(4) 4 2 91719 COM STDRD Scalar 104 L(4) 4 scalar COM STHE Scalar 8 R(4) 4 1 76 COM STHEQ Scalar 9 R(4) 4 1 152 COM SVALUE Scalar 8 R(4) 4 scalar COM SWNET Scalar 12 R(4) 4 2 91719 COM T Scalar 5 R(4) 4 3 4127355 COM 393,454,458,462,466,698 T1D Local 173 R(4) 4 scalar 173 T500 Scalar 6 R(4) 4 2 91719 COM 152 TCLOD Scalar 3 R(4) 4 scalar COM 277,306,335 TDTD Scalar 3 R(4) 4 scalar COM TDTQ2 Scalar 3 R(4) 4 scalar COM TFRZ Param 39 R(4) 4 scalar TG Scalar 14 R(4) 4 2 91719 COM THE0 Scalar 8 R(4) 4 1 76 COM THE0Q Scalar 9 R(4) 4 1 152 COM THL Scalar 5 R(4) 4 scalar COM TLM0D Scalar 6 R(4) 4 scalar COM Page 23 Source Listing CLDRAD 2014-12-17 20:47 Symbol Table CLDRAD.f Name Object Declared Type Bytes Dimen Elements Attributes References TLMH Local 698 R(4) 4 scalar 698,699 TPH0D Scalar 6 R(4) 4 scalar COM TRDLW Scalar 3 R(4) 4 scalar COM 543,601,659 TRDSW Scalar 3 R(4) 4 scalar COM 514,572,630 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 276,305,334,513,542,571,600,629,65 8 TTBL Scalar 18 R(4) 4 2 10184 COM TTBLQ Scalar 18 R(4) 4 2 66880 COM TTV Scalar 8 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM U00 Scalar 3 R(4) 4 2 91719 COM UL Scalar 3 R(4) 4 1 90 COM V Scalar 5 R(4) 4 3 4127355 COM VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VRBLS Common 2 67138308 VTM Scalar 5 R(4) 4 3 4127355 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 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 389,396,457,465 Page 24 Source Listing CLDRAD 2014-12-17 20:47 Subprograms/Common Blocks CLDRAD.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References ACMCLD Common 2 1467512 ACMRDL Common 2 2201268 ACMRDS Common 2 2201268 CLDRAD Subr 1 CLDWTR Common 2 17610408 CNVCLD Common 2 1467504 CTLBLK Common 2 96 EXTRA Common 3 71506268 LOOPS Common 2 748952 MAPOT Common 2 1448 MASKS Common 2 34853220 OPTIONS Common 2 24 PARALLEL Common 9 8240 PHYS Common 2 7690036 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__ Page 25 Source Listing CLDRAD 2014-12-17 20:47 CLDRAD.f -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest no -ftz -fp_model 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 : CLDRAD.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100