Page 1 Source Listing RADTN 2014-12-17 20:47 /tmpnwprd/ifortG0lQe6.i 1 # 1 "RADTN.F" 2 C&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 3 SUBROUTINE RADTN 4 C ****************************************************************** 5 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 6 C . . . 7 C SUBPROGRAM: RADTN THE OUTER RADIATION DRIVER 8 C PRGRMMR: BLACK ORG: W/NP22 DATE: 93-12-?? 9 C 10 C ABSTRACT: 11 C RADTN PRIMARILY SERVES TO SET UP THE ARRAYS NEEDED AS INPUT 12 C FOR RADFS (THE INNER RADIATION DRIVER). GROUPS OF MODEL COLUMNS 13 C ARE SENT TO RADFS BUT FIRST THEY ARE "LIFTED" SO THAT THE LOWEST 14 C LAYER ABOVE THE GROUND HAS A VERTICAL INDEX VALUE OF LM NOT LMH. 15 C THIS ROUTINE IS CALLED AS OFTEN AS DESIRED (EVERY 1 TO 2 HOURS) 16 C FOR BOTH THE SHORT AND LONGWAVE EFFECTS. THE RESULTING TEMPER- 17 C ATURE TENDENCIES, TOTAL DOWNWARD AND SHORTWAVE UPWARD FLUXES ARE 18 C COLLECTED. 19 C THE INITIAL GROUND POTENTIAL TEMPERATURE IS ALSO COMPUTED HERE 20 C AND IS SIMPLY AN ADIABATIC EXTRAPOLATION FROM THE LOWEST MID- 21 C LAYER VALUE ABOVE THE GROUND. 22 C 23 C PROGRAM HISTORY LOG: 24 C 87-09-?? BLACK - ORIGINATOR 25 C 92-10-?? BALDWIN - VARIOUS CLOUD EFFECTS WERE INCLUDED 26 C WHICH WERE ALREADY IN THE MRF 27 C 93-11-?? ZHAO - TIED TO UPDATED GFDL RADIATION SCHEME 28 C USING MODEL-PREDICTED CLOUD 29 C 95-03-25 BLACK - CONVERSION FROM 1-D TO 2-D IN HORIZONTAL 30 C 95-04-13 BLACK - PARALLELIZED THE LARGE LOOP STEPPING 31 C THROUGH THE DOMAIN THAT CALLS RADFS 32 C 95-10-10 ZHAO - i) THE CALCULATION OF CLOUD FRACTION WAS 33 C CHANGED TO USE BOTH CLOUD WATER/ICE 34 C MIXING RATIO AND RELATIVE HUMIDITY 35 C (RANDALL, 1994); 36 C ii) THE CLOUD INPUTS WERE CHANGED TO USE 37 C CLOUD FRACTION IN EACH MODEL LAYER 38 C AFTER Y.T. HOU (1995). 39 C 96-06-03 ZHAO - SNOW ALBEDO IS CHANGED ACCORDING TO 40 C SUGGESTIONS FROM KEN MITCHELL AND FEI CHEN 41 C 96-07-23 ZHAO - ADD CALL TO SOLARD TO CALCULATE THE NON- 42 C DIMENSIONAL SUN-EARTH DISTANCE R1 WHICH 43 C WILL BE USED IN RADFS TO COMPUTE SOLAR 44 C CONSTANT SOLC ON EACH DAY 45 C 96-07-26 BLACK - ADDED OZONE COMPUTATIONS 46 C 97-05-19 ZHAO - DIAGNOSTIC CLOUDS (LOW, MIDDLE, AND HIGH) 47 C ARE MODIFIED TO USE THE MAXIMUM OF 48 C CONVECTIVE AND STRATIFORM. THIS WILL REPLACE 49 C THE PREVIOUS SCHEME WHICH USES ONLY CONVECTIVE 50 C CLOUDS AT CONVECTIVE POINTS. THIS WILL 51 C AFFECT CFRACL, CFRACM, CFRACH, AND WILL 52 C AFFECT THE TOTAL CLOUD FRACTION CALCULATION 53 C IN THE POST PROCESSORS. 54 C 98-??-?? TUCCILLO - ADDED PARALLELISM FOR CLASS VIII 55 C 98-10-27 BLACK - PARALLELISM INTO NEWEST VERSION 56 C 57 C Page 2 Source Listing RADTN 2014-12-17 20:47 RADTN.F 58 C 59 C USAGE: CALL RADTN FROM MAIN PROGRAM EBU 60 C INPUT ARGUMENT LIST: 61 C NONE 62 C 63 C OUTPUT ARGUMENT LIST: 64 C NONE 65 C 66 C OUTPUT FILES: 67 C NONE 68 C 69 C SUBPROGRAMS CALLED: 70 C 71 C UNIQUE: 72 C ZENITH 73 C RADFS 74 C GRADFS 75 C SOLARD 76 C O3CLIM 77 C OZON2D 78 C 79 C LIBRARY: 80 C NONE 81 C 82 C COMMON BLOCKS: CTLBLK 83 C LOOPS 84 C MASKS 85 C DYNAMD 86 C PHYS 87 C VRBLS 88 C PVRBLS 89 C CLDWTR 90 C CNVCLD 91 C INDX 92 C ACMCLD 93 C ACMRDS 94 C ACMRDL 95 C 96 C ATTRIBUTES: 97 C LANGUAGE: FORTRAN 90 98 C MACHINE : IBM SP 99 C$$$ 100 C ****************************************************************** 101 C * Note: Convective clouds are added in this subroutine * 102 C * for use in the eta model in which model-predicted * 103 C * clouds are not used in the convective precipitation * 104 C * processes. * 105 C * For use with the version of the eta model in which * 106 C * the model-predicted clouds are linked into the model's * 107 C * convective precipitation processes, just set: * 108 C * CNCLD=.FALSE. * 109 C * Qingyun Zhao 12-9-94 * 110 C ****************************************************************** 111 C----------------------------------------------------------------------- 112 INCLUDE "parmeta" 135 INCLUDE "parm.tbl" 137 INCLUDE "parmsoil" Page 3 Source Listing RADTN 2014-12-17 20:47 RADTN.F 139 INCLUDE "mpp.h" 140 # 1 "./sp.h" 1 141 # 4 142 143 # 116 "RADTN.F" 2 144 C----------------------------------------------------------------------- 203 P A R A M E T E R 204 & (CAPA=0.28589641,RTD=57.2957795 205 &, WA=.10,WG=1.-WA,KSMUD=0) 206 C-------------------------CLOUD---------------------------------------- 207 P A R A M E T E R 208 & (A1=610.78,A2=17.2693882,A3=273.16,A4=35.86 209 c &, PQ0=379.90516,SNOALB=0.55) 210 &, PQ0=379.90516) 211 C-------------------------CLOUD---------------------------------------- 212 P A R A M E T E R 213 & (IMJM=IM*JM-JM/2,JAM=6+2*(JM-10),LM1=LM-1,LP1=LM+1) 214 C 215 P A R A M E T E R 216 & (SLPM=1.01325E5,EPSQ1=1.E-5,EPSQ=2.E-12,EPSO3=1.E-10,HPINC=1.E1 217 &, CLDRH0=0.80,TRESH=1.00,RNRM=1./(TRESH-CLDRH0) 218 &, CLDRH2=0.90,TRESH2=1.00,RNRM2=1./(TRESH2-CLDRH2) 219 &, CLAPSE=-0.0005,CLPSE=-0.0006,DCLPS=-0.0001 220 &, CM1=2937.4,CM2=4.9283,CM3=23.5518,EPS=0.622,PBOT=10000.0 221 &, STBOL=5.67E-8,PI2=2.*3.14159265,RLAG=14.8125) 222 C 223 P A R A M E T E R 224 & (NB=12) 225 C----------------------------------------------------------------------- 226 P A R A M E T E R 227 & (K15=SELECTED_REAL_KIND(15)) 228 C 229 R E A L 230 & (KIND=K15) PROD,DDX,EEX 231 C----------------------------------------------------------------------- 232 L O G I C A L 233 & RUN,FIRST,RESTRT,SIGMA,CALL1,SHORT,LONG 234 &,BCLD(IDIM1:IDIM2),BTEMP1(IDIM1:IDIM2) 235 &,BITX,BITY,BITZ,BITW,BIT1,BIT2,BITC,BITS,BITCP1,BITSP1 236 C-------------------------CONVECTION------------------------------------ 237 L O G I C A L 238 & CNCLD 239 C-------------------------CONVECTION------------------------------------ 240 INCLUDE "CTLBLK.comm" 241 C----------------------------------------------------------------------- 251 INCLUDE "LOOPS.comm" 252 C----------------------------------------------------------------------- 259 INCLUDE "MASKS.comm" 260 C----------------------------------------------------------------------- 271 INCLUDE "DYNAMD.comm" 272 C----------------------------------------------------------------------- 288 INCLUDE "PHYS.comm" 289 C----------------------------------------------------------------------- 317 INCLUDE "VRBLS.comm" 318 C----------------------------------------------------------------------- 328 INCLUDE "PVRBLS.comm" 329 C----------------------------------------------------------------------- Page 4 Source Listing RADTN 2014-12-17 20:47 RADTN.F 359 INCLUDE "SOIL.comm" 360 C----------------------------------------------------------------------- 373 INCLUDE "CLDWTR.comm" 374 C----------------------------------------------------------------------- 382 INCLUDE "CNVCLD.comm" 383 C----------------------------------------------------------------------- 390 INCLUDE "INDX.comm" 391 C----------------------------------------------------------------------- 403 INCLUDE "ACMCLD.comm" 404 C---------------------------------------------------------------------- 411 INCLUDE "ACMRDL.comm" 412 C---------------------------------------------------------------------- 421 INCLUDE "ACMRDS.comm" 422 C----------------------------------------------------------------------- 431 COMMON /CUINIT/ CURAD 432 LOGICAL CURAD 433 C----------------------------------------------------------------------- 434 C O M M O N 435 & /SWRSAV/ABCFF(NB),PWTS(NB),CFCO2,CFO3,REFLO3,RRAYAV 436 C O M M O N 437 & /RD1TIM/K400,CTHK(3),LTOP(3),PTOPC(4),TAUCV(3),R1 438 &, LVL(IDIM1:IDIM2,JDIM1:JDIM2) 439 C----------------------------------------------------------------------- 440 D I M E N S I O N 441 & TENDK (LM),CLDAMT(0:LM) 442 &, PSFC (IDIM1:IDIM2),TSKN (IDIM1:IDIM2) 443 c &, ALBEDO(IDIM1:IDIM2),XLAT(IDIM1:IDIM2),COSZ (IDIM1:IDIM2) 444 &, ALBDO (IDIM1:IDIM2),XLAT(IDIM1:IDIM2),COSZ (IDIM1:IDIM2) 445 &, CLDCFR(IDIM1:IDIM2,3),MBOT(IDIM1:IDIM2,3) 446 &, CLDF (IDIM1:IDIM2,LP1),SLMSK (IDIM1:IDIM2) 447 &, TENDS (IDIM1:IDIM2,LM),TENDL (IDIM1:IDIM2,LM) 448 C 449 &, PMID (IDIM1:IDIM2,LM),TMID (IDIM1:IDIM2,LM) 450 &, QMID (IDIM1:IDIM2,LM),THMID(IDIM1:IDIM2,LM) 451 &, OZN (IDIM1:IDIM2,LM),POZN (IDIM1:IDIM2,LM) 452 &, MTOP(IDIM1:IDIM2,3),ICVB(IDIM1:IDIM2), ICVT(IDIM1:IDIM2) 453 &, CV(IDIM1:IDIM2),SV(IDIM1:IDIM2) 454 C 455 &, FLWUP (IDIM1:IDIM2),FSWDN (IDIM1:IDIM2),FSWUP (IDIM1:IDIM2) 456 &, FSWDNS(IDIM1:IDIM2),FSWUPS(IDIM1:IDIM2) 457 &, FLWDNS(IDIM1:IDIM2),FLWUPS(IDIM1:IDIM2) 458 &, PDSL (IDIM1:IDIM2,JDIM1:JDIM2) 459 &, FNE(IDIM1:IDIM2,JDIM1:JDIM2),FSE(IDIM1:IDIM2,JDIM1:JDIM2) 460 &, TL (IDIM1:IDIM2,JDIM1:JDIM2) 461 D I M E N S I O N 462 & PBOTL(IDIM1:IDIM2,JDIM1:JDIM2), PTOPL(IDIM1:IDIM2,JDIM1:JDIM2) 463 &,PBOTM(IDIM1:IDIM2,JDIM1:JDIM2), PTOPM(IDIM1:IDIM2,JDIM1:JDIM2) 464 &,PBOTH(IDIM1:IDIM2,JDIM1:JDIM2), PTOPH(IDIM1:IDIM2,JDIM1:JDIM2) 465 &,TOT (IDIM1:IDIM2,JDIM1:JDIM2) 466 D I M E N S I O N 467 & CC(9),PPT(9) 468 D I M E N S I O N 469 & PINT(IDIM1:IDIM2,LP1),PHALF(LP1),CSTR(IDIM1:IDIM2) 470 &, EMIS(IDIM1:IDIM2,LP1), TAUC(IDIM1:IDIM2) 471 &, CVB(IDIM1:IDIM2),CVT(IDIM1:IDIM2),TAUDAR(IDIM1:IDIM2) 472 D I M E N S I O N 473 & CAMT(IDIM1:IDIM2,LP1),NCLDS(IDIM1:IDIM2) Page 5 Source Listing RADTN 2014-12-17 20:47 RADTN.F 474 &,ITYP(IDIM1:IDIM2,LP1),KTOP(IDIM1:IDIM2,LP1) 475 &,KBTM(IDIM1:IDIM2,LP1),RRCL(IDIM1:IDIM2,NB,LP1) 476 &,TTCL(IDIM1:IDIM2,NB,LP1),KCLD(IDIM1:IDIM2) 477 C--------------------CLOUD---------------------------------------------- 478 D I M E N S I O N 479 & CCR(IDIM1:IDIM2,LM),IW(IDIM1:IDIM2,LM),CSMID(IDIM1:IDIM2,LM) 480 &,WMID(IDIM1:IDIM2,LM),HMID(IDIM1:IDIM2,LM) 481 &,BMID(IDIM1:IDIM2),UMID(IDIM1:IDIM2) 482 &,CCMID(IDIM1:IDIM2,LM) 483 C--------------------CLOUD---------------------------------------------- 484 D A T A 485 1 PLOMD/64200./,PMDHI/35000./,PHITP/15000./,P400/40000./ 486 2, PLBTM/105000./ 487 D A T A 488 1 NFILE/14/ 489 DATA CC/0.,0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8/ 490 DATA PPT/.14,.31,.70,1.6,3.4,7.7,17.,38.,85./ 491 C---------------------------------------------------------------------- 492 UTIM=1. 493 CNCLD=.TRUE. 494 C*** 495 C*** ASSIGN THE PRESSURES FOR CLOUD DOMAIN BOUNDARIES 496 C*** 497 PTOPC(1)=PLBTM 498 PTOPC(2)=PLOMD 499 PTOPC(3)=PMDHI 500 PTOPC(4)=PHITP 501 C*** 502 C*** FIND THE 'SEA LEVEL PRESSURE'. 503 C*** 504 DO J=MYJS,MYJE 505 DO I=MYIS,MYIE 506 PDSL(I,J)=RES(I,J)*PD(I,J) 507 ENDDO 508 ENDDO 509 C********************************************************************** 510 C*** THE FOLLOWING CODE IS EXECUTED EACH TIME THE RADIATION IS CALLED. 511 C********************************************************************** 512 C----------------------CONVECTION-------------------------------------- 513 C NRADPP IS THE NUMBER OF TIME STEPS TO ACCUMULATE CONVECTIVE PRECIP 514 C FOR RADIATION 515 C NOTE: THIS WILL NOT WORK IF NRADS AND NRADL ARE DIFFERENT UNLESS 516 C THEY ARE INTEGER MULTIPLES OF EACH OTHER 517 C CLSTP IS THE NUMBER OF HOURS OF THE ACCUMULATION PERIOD 518 C 519 NTSPH=NINT(3600./DT) 520 NRADPP=MIN(NRADS,NRADL) 521 CLSTP=1.0*NRADPP/NTSPH 522 C----------------------CONVECTION-------------------------------------- 523 C*** 524 C*** STATE WHETHER THE SHORT OR LONGWAVE COMPUTATIONS ARE TO BE DONE. 525 C*** 526 SHORT=.FALSE. 527 LONG=.FALSE. 528 IF(MOD(NTSD,NRADS).EQ.1.OR.RESTRT)SHORT=.TRUE. 529 IF(MOD(NTSD,NRADL).EQ.1.OR.RESTRT)LONG=.TRUE. 530 ITIMSW=0 Page 6 Source Listing RADTN 2014-12-17 20:47 RADTN.F 531 ITIMLW=0 532 IF(SHORT)ITIMSW=1 533 IF(LONG) ITIMLW=1 534 C----------------------------------------------------------------------- 535 C*** 536 C*** FLAG FOR RESETTING CUPPT,HTOP,HBOT IN CHKOUT 537 C*** 538 IF (MOD(NTSD,NRADPP).EQ.1) CURAD=.TRUE. 539 C*** 540 C*** FIND THE MEAN COSINE OF THE SOLAR ZENITH ANGLE 541 C*** BETWEEN THE CURRENT TIME AND THE NEXT TIME RADIATION IS 542 C*** CALLED. ONLY AVERAGE IF THE SUN IS ABOVE THE HORIZON. 543 C*** 544 TIME=(NTSD-1)*DT 545 CALL ZENITH(TIME,DAYI,HOUR) 546 JD=INT(DAYI+0.50) 547 ADDL=0. 548 IF(MOD(IDAT(3),4).EQ.0)ADDL=1. 549 RANG=PI2*(DAYI-RLAG)/(365.25+ADDL) 550 RSIN1=SIN(RANG) 551 RCOS1=COS(RANG) 552 RCOS2=COS(2.*RANG) 553 IF(SHORT)THEN 554 !$omp parallel do private(i,j) 555 DO J=MYJS,MYJE 556 DO I=MYIS,MYIE 557 CZMEAN(I,J)=0. 558 TOT(I,J)=0. 559 ENDDO 560 ENDDO 561 C 562 DO II=0,NRADS,NPHS 563 TIMES=(NTSD-1)*DT+II*DT 564 CALL ZENITH(TIMES,DAYI,HOUR) 565 !$omp parallel do private(i,j) 566 DO J=MYJS,MYJE 567 DO I=MYIS,MYIE 568 IF(CZEN(I,J).GT.0.)THEN 569 CZMEAN(I,J)=CZMEAN(I,J)+CZEN(I,J) 570 TOT(I,J)=TOT(I,J)+1. 571 ENDIF 572 ENDDO 573 ENDDO 574 ENDDO 575 !$omp parallel do private(i,j) 576 DO J=MYJS,MYJE 577 DO I=MYIS,MYIE 578 IF(TOT(I,J).GT.0.)CZMEAN(I,J)=CZMEAN(I,J)/TOT(I,J) 579 ENDDO 580 ENDDO 581 ENDIF 582 C 583 C2345678901234567890123456789012345678901234567890123456789012345678901 584 !$omp parallel do 585 c!$omp& private (aa,albedo,apes,bb,bcld,bit1,bit2,bitc,bitcp1,bits, 586 !$omp& private (aa,albdo,apes,bb,bcld,bit1,bit2,bitc,bitcp1,bits, 587 !$omp& bitsp1,bitw,bitx,bity,bitz,bmid,btemp1,camt, Page 7 Source Listing RADTN 2014-12-17 20:47 RADTN.F 588 !$omp& cc1,cc2,ccmid,ccr,cfravg,cl1,cl2, 589 !$omp& cldamt,cldcfr,cldmax,clpfil,cosz,cr1, 590 !$omp& csmid,cstr,cv,cwmkl,dd,ddp,delp,denom, 591 !$omp& dpcl,dthdp,ee,emis,exner,fctra, 592 !$omp& fctrb,ff,fiq,fiw,flwdns,flwup,flwups,fswdn,fswdns, 593 !$omp& fswup,fswups,gg,hh,hmid,i,icvb,icvt,ir,ityp, 594 !$omp& iw,iwkl,j,kbt1,kbt2,kbtm,kcld,kntlyr, 595 !$omp& kth1,kth2,ktop,ktop1,l,l400,lbase,lin,ll, 596 !$omp& llbot,lltop,lml,ltrop,lvlij,malvl,mbot,mtop, 597 !$omp& n,nband,nbtm,nc,ncld,nclds,nktp,nlvl, 598 !$omp& nmod,ozn,p1,p2,pdslij,pint, 599 !$omp& pmid,pmod,pozn,pp,prs1,prs2,psfc,qc,qi,qint,qkl, 600 !$omp& qmid,qsum,qw,rqkl,rrcl,slmsk,snofac,sv, 601 !$omp& tauc,taudar,tcld,tendl,tends,thmid, 602 !$omp& tkl,tmid,tmt0,tmt15,tskn,ttcl, 603 !$omp& u00kl,umid,us,wmid,xlat) 604 c 605 C******************************************************************** 606 C*** THIS IS THE BEGINNING OF THE PRIMARY LOOP THROUGH THE DOMAIN 607 C******************************************************************** 608 C ********************* 609 DO 700 J = MYJS, MYJE 610 C ********************* 611 C 612 DO 125 L=1,LM 613 DO I=MYIS,MYIE 614 IR=IRAD(I) 615 TMID(I,L)=T(I,J,1) 616 QMID(I,L)=EPSQ 617 CSMID(I,L)=0. 618 WMID(I,L)=0. 619 CCMID(I,L)=0. 620 IW(I,L)=0. 621 CCR(I,L)=0. 622 HMID(I,L)=0. 623 OZN(I,L)=EPSO3 624 TENDS(I,L)=0. 625 TENDL(I,L)=0. 626 ENDDO 627 125 CONTINUE 628 C 629 DO 140 N=1,3 630 DO I=MYIS,MYIE 631 CLDCFR(I,N)=0. 632 MTOP(I,N)=0 633 MBOT(I,N)=0 634 ENDDO 635 140 CONTINUE 636 C*** 637 C*** FILL IN WORKING ARRAYS WHERE VALUES AT L=LM ARE THOSE THAT 638 C*** ARE ACTUALLY AT ETA LEVEL L=LMH. 639 C*** 640 DO 200 I=MYIS,MYIE 641 IR=IRAD(I) 642 LML=LMH(I,J) 643 LVLIJ=LVL(I,J) 644 ctlb BMID(I)=HBM2(IR,J) Page 8 Source Listing RADTN 2014-12-17 20:47 RADTN.F 645 BMID(I)=HBM2(I,J) 646 UMID(I)=U00(I,J) 647 C 648 DO L=1,LML 649 PMID(I,L+LVLIJ)=AETA(L)*PDSL(I,J)+PT 650 PINT(I,L+LVLIJ+1)=ETAD(L+1)*PDSL(I,J)+PT 651 EXNER=(1.E5/PMID(I,L+LVLIJ))**CAPA 652 TMID(I,L+LVLIJ)=T(I,J,L) 653 THMID(I,L+LVLIJ)=T(I,J,L)*EXNER 654 QMID(I,L+LVLIJ)=Q(I,J,L) 655 WMID(I,L+LVLIJ)=CWM(I,J,L) 656 HMID(I,L+LVLIJ)=HTM(I,J,L) 657 ENDDO 658 C*** 659 C*** FILL IN ARTIFICIAL VALUES ABOVE THE TOP OF THE DOMAIN. 660 C*** PRESSURE DEPTHS OF THESE LAYERS IS 1 HPA. 661 C*** TEMPERATURES ABOVE ARE ALREADY ISOTHERMAL WITH (TRUE) LAYER 1. 662 C*** 663 IF(LVLIJ.GT.0)THEN 664 KNTLYR=0 665 C 666 DO L=LVLIJ,1,-1 667 KNTLYR=KNTLYR+1 668 PMID(I,L)=PT-REAL(2*KNTLYR-1)*0.5*HPINC 669 PINT(I,L+1)=PMID(I,L)+0.5*HPINC 670 EXNER=(1.E5/PMID(I,L))**CAPA 671 THMID(I,L)=TMID(I,L)*EXNER 672 ENDDO 673 ENDIF 674 C 675 IF(LVLIJ.EQ.0) THEN 676 PINT(I,1)=PT 677 ELSE 678 PINT(I,1)=PMID(I,1)-0.5*HPINC 679 ENDIF 680 200 CONTINUE 681 C*** 682 C*** FILL IN THE SURFACE PRESSURE, SKIN TEMPERATURE, GEODETIC LATITUDE, 683 C*** ZENITH ANGLE, SEA MASK, AND ALBEDO. THE SKIN TEMPERATURE IS 684 C*** NEGATIVE OVER WATER. 685 C*** 686 DO 250 I=MYIS,MYIE 687 PSFC(I)=PD(I,J)+PT 688 APES=(PSFC(I)*1.E-5)**CAPA 689 TSKN(I)=THS(I,J)*APES*(1.-2.*SM(I,J)) 690 if(mype.eq.13.and.i.eq.40.and.j.eq.77) then 691 print*,'apes,sm(i,j),ths(i,j)=',apes,sm(i,j),ths(i,j) 692 endif 693 SLMSK(I)=SM(I,J) 694 C 695 C ---------------------------------------------------------------------- 696 c turn off snow albedo calculation since it is now calculated in SFLX. 697 c SNO(I,J)=AMAX1(SNO(I,J),0.) 698 c SNOFAC=AMIN1(SNO(I,J)/0.02, 1.0) 699 c ALBEDO(I)=ALB(I,J)+(1.0-VEGFRC(I,J))*SNOFAC*(SNOALB-ALB(I,J)) 700 ALBDO(I)=ALBEDO(I,J) 701 C Page 9 Source Listing RADTN 2014-12-17 20:47 RADTN.F 702 XLAT(I)=GLAT(I,J)*RTD 703 COSZ(I)=CZMEAN(I,J) 704 250 CONTINUE 705 C----------------------------------------------------------------------- 706 C*******************STRATIFORM CLOUD SECTION*************************** 707 C----------------------------------------------------------------------- 708 C CALCULATE STRATIFORM CLOUD COVERAGE AT EACH MODEL GRID POINT WHICH 709 C WILL BE USED IN THE MODEL RADIATION PARAMETERIZATION SCHEME. 710 C----------------------------------------------------------------------- 711 US=1. 712 CCLIMIT=1.0E-3 713 CLIMIT =1.0E-20 714 C------------------QW, QI AND QINT-------------------------------------- 715 DO 280 I=MYIS,MYIE 716 LML=LMH(I,J) 717 LVLIJ=LVL(I,J) 718 C 719 DO 275 L=1,LML 720 LL=L+LVLIJ 721 HH=HMID(I,LL)*BMID(I) 722 TKL=TMID(I,LL) 723 QKL=QMID(I,LL) 724 CWMKL=WMID(I,LL) 725 TMT0=(TKL-273.16)*HH 726 TMT15=AMIN1(TMT0,-15.)*HH 727 AI=0.008855 728 BI=1. 729 C 730 IF(TMT0.LT.-20.)THEN 731 AI=0.007225 732 BI=0.9674 733 ENDIF 734 C 735 PP=PMID(I,LL) 736 QW=HH*PQ0/PP*EXP(HH*A2*(TKL-A3)/(TKL-A4)) 737 QI=QW*(BI+AI*AMIN1(TMT0,0.)) 738 QINT=QW*(1.-0.00032*TMT15*(TMT15+15.)) 739 IF(TMT0.LE.-40.) QINT=QI 740 C-------------------ICE-WATER ID NUMBER IW------------------------------ 741 U00KL=UMID(I)+UL(L)*(0.95-UMID(I))*UTIM 742 IF(TMT0.LT.-15.0)THEN 743 FIQ=QKL-U00KL*QI 744 IF(FIQ.GT.0..OR.CWMKL.GT.CLIMIT)THEN 745 IW(I,LL)=1 746 ELSE 747 IW(I,LL)=0 748 ENDIF 749 ENDIF 750 C 751 IF(TMT0.GE.0.)THEN 752 IW(I,LL)=0 753 ENDIF 754 C 755 IF(TMT0.LT.0.0.AND.TMT0.GE.-15.0)THEN 756 IW(I,LL)=0 757 IF(IW(I,LL-1).EQ.1.AND.CWMKL.GT.CLIMIT) IW(I,LL)=1 758 ENDIF Page 10 Source Listing RADTN 2014-12-17 20:47 RADTN.F 759 C 760 IWKL=IW(I,LL) 761 C 762 C----------------THE SATURATION SPECIFIC HUMIDITY------------------------ 763 C 764 FIW=FLOAT(IWKL) 765 QC=(1.-FIW)*QINT+FIW*QI 766 C 767 C----------------THE RELATIVE HUMIDITY---------------------------------- 768 C 769 IF(QC.LE.EPSQ1.OR.QKL.LE.EPSQ1)THEN 770 RQKL=0. 771 ELSE 772 RQKL=QKL/QC 773 ENDIF 774 C 775 C----------------CLOUD COVER RATIO CCR---------------------------------- 776 C 777 IF(RQKL.GE.0.9999)THEN 778 CCR(I,LL)=AMIN1(US,RQKL) 779 ELSE 780 ARG=-1000.*CWMKL/(US-RQKL) 781 ARG=AMAX1(ARG,-25.) 782 CCR(I,LL)= RQKL*(1.-EXP(ARG)) 783 ENDIF 784 CSMID(I,LL)=AMIN1(US,CCR(I,LL)) 785 C---------------------------------------------------------------------- 786 275 CONTINUE 787 280 CONTINUE 788 C---------------------------------------------------------------------- 789 C********************************************************************** 790 C NOW CHECK THE CLOUDS PRODUCED ABOVE TO MAKE SURE THEY ARE GOOD 791 C ENOUGH FOR RADIATION CALCULATIONS 792 C********************************************************************** 793 C*** 794 C*** NO STRATIFORM CLOUDS FOR THIS TYPE 795 C*** 796 DO 350 I=MYIS,MYIE 797 C 798 LML=LMH(I,J) 799 LVLIJ=LVL(I,J) 800 C*** 801 C*** ZERO OUT CLDAMT IF LAND AND BELOW PBOT ABOVE GROUND 802 C*** 803 IF(SM(I,J).LT.0.5)THEN 804 DO L=1,LML 805 LL=LML-L+1+LVLIJ 806 DDP=PSFC(I)-PMID(I,LL) 807 IF(DDP.GE.PBOT) GO TO 290 808 CSMID(I,LL)=0. 809 ENDDO 810 290 CONTINUE 811 ENDIF 812 C*** 813 C*** CHECK FOR OCEAN STRATUS (LOW CLOUD) 814 C*** LOOK ONLY OVER OCEAN AND ONLY IF AN INVERSION (DTHDP.LE.-0.05) 815 C*** IS PRESENT WITH AT LEAST 2 CLOUD FREE LAYERS ABOVE IT Page 11 Source Listing RADTN 2014-12-17 20:47 RADTN.F 816 C*** 817 IF(SM(I,J).GT.0.5)THEN 818 C 819 C*** FIND BASE OF INVERSION 820 C 821 LBASE=LM 822 DO L=1,LML-1 823 LL=LML-L+1+LVLIJ 824 DTHDP=(THMID(I,LL-1)-THMID(I,LL)) 825 & /(PMID(I,LL-1)-PMID(I,LL)) 826 IF(DTHDP.LE.CLAPSE)THEN 827 LBASE=LL 828 GO TO 300 829 ENDIF 830 ENDDO 831 300 CONTINUE 832 C 833 C*** CHECK 2 LAYERS ABOVE LBASE FOR DRYNESS 834 C 835 IF(CSMID(I,LBASE-1).LE.0..AND.CSMID(I,LBASE-2).LE.0. 836 1 .AND.LBASE.LT.LM)THEN 837 IF(DTHDP.GT.CLPSE)THEN 838 CLPFIL=1.-((CLPSE-DTHDP)/DCLPS) 839 ELSE 840 CLPFIL=1. 841 ENDIF 842 C 843 DO L=1,LML 844 LL=LML-L+1+LVLIJ 845 DDP=PSFC(I)-PMID(I,LL) 846 IF(DDP.GE.PBOT) GO TO 310 847 CSMID(I,LL)=CSMID(I,LL)*CLPFIL 848 ENDDO 849 310 CONTINUE 850 C 851 C*** IF NO INVERSION OR IF CLDS EXIST IN EITHER OF THE 2 LAYERS ABOVE 852 C*** INVERSION, ZERO OUT CLOUD BELOW PBOT 853 C 854 ELSE 855 DO L=1,LML 856 LL=LML-L+1+LVLIJ 857 DDP=PSFC(I)-PMID(I,LL) 858 IF(DDP.GE.PBOT) GO TO 320 859 CSMID(I,LL)=0. 860 ENDDO 861 320 CONTINUE 862 C 863 ENDIF 864 C------------ 865 ENDIF 866 C------------ 867 C*** 868 C*** REMOVE HIGH CLOUDS ABOVE THE TROPOPAUSE 869 C*** 870 L400=LM 871 DO L=1,LML 872 LL=LML-L+1+LVLIJ Page 12 Source Listing RADTN 2014-12-17 20:47 RADTN.F 873 IF(PMID(I,LL).LE.40000.0)THEN 874 L400=LL 875 GO TO 330 876 ENDIF 877 ENDDO 878 330 CONTINUE 879 C 880 LTROP=LM 881 DO LL=L400,2,-1 882 DTHDP=(THMID(I,LL-1)-THMID(I,LL)) 883 1 /(PMID(I,LL-1)-PMID(I,LL)) 884 C 885 IF(DTHDP.LT.-0.0025.OR.QMID(I,LL).LE.EPSQ1)THEN 886 LTROP=LL 887 GOTO 340 888 ENDIF 889 ENDDO 890 340 IF(LTROP.LT.LM)THEN 891 DO LL=LTROP,1,-1 892 CSMID(I,LL)=0. 893 ENDDO 894 ENDIF 895 350 CONTINUE 896 C 897 C********************************************************************* 898 C*****************END OF STRATIFORM CLOUD SECTION***************** 899 C------------------------CONVECTION-------------------------------- 900 C*** 901 C*** CONVECTIVE CLOUD SECTION 902 C*** 903 C*** THIS PART WAS MODIFIED TO COMPUTE CONVECTIVE CLOUDS AT EACH 904 C*** MODEL LAYER BASED ON CONVECTIVE PRECIPITATION RATES. CURRENTLY, 905 C*** CLOUDS ARE SET TO 0.75*CV(I) BELOW 400MB 906 C*** AND 0.90*CV(I) ABOVE 400MB TO ACCOUNT FOR CIRRUS CAP 907 C*** Q.ZHAO 95-3-22 908 C 909 C*** 910 C*** NON-PRECIPITATING CLOUD FRACTION OF 20 PERCENT IS ADDED AT 911 C*** AT POINTS WHERE THE SHALLOW AND DEEP CONVECTIONS ACCUR. 912 C*** Q. ZHAO 97-5-2 913 C 914 C COMPUTE THE CONVECTIVE CLOUD COVER FOR RADIATION 915 C 916 C----------------------------------------------------------------- 917 IF(CNCLD)THEN 918 C----------------------------------------------------------------- 919 DO 375 I=MYIS,MYIE 920 IF(HBOT(I,J)-HTOP(I,J).GT.1.0)THEN 921 C IF(HTOP(I,J).LT.HBOT(I,J))THEN 922 SV(I)=0.0 923 ELSE 924 SV(I)=0.0 925 ENDIF 926 C 927 PMOD=CUPPT(I,J)*24.0*1000.0/CLSTP 928 NMOD=0 929 C Page 13 Source Listing RADTN 2014-12-17 20:47 RADTN.F 930 DO NC=1,9 931 IF(PMOD.GT.PPT(NC)) NMOD=NC 932 ENDDO 933 C 934 C*** CLOUD TOPS AND BOTTOMS COME FROM CUCNVC 935 C*** ADD LVL TO BE CONSISTENT WITH OTHER WORKING ARRAYS 936 C 937 IF(NMOD.EQ.0)THEN 938 CV(I)=0. 939 ELSEIF(NMOD.EQ.9)THEN 940 CV(I)=CC(9) 941 ELSE 942 CC1=CC(NMOD) 943 CC2=CC(NMOD+1) 944 P1=PPT(NMOD) 945 P2=PPT(NMOD+1) 946 CV(I)=CC1+(CC2-CC1)*(PMOD-P1)/(P2-P1) 947 ENDIF 948 C 949 CV(I)=AMAX1(SV(I),CV(I)) 950 CV(I)=AMIN1(1.0,CV(I)) 951 C 952 IF(CV(I).EQ.0.0)THEN 953 ICVT(I)=0 954 ICVB(I)=0 955 ELSE 956 ICVT(I)=INT(HTOP(I,J)+0.50)+LVL(I,J) 957 ICVB(I)=INT(HBOT(I,J)+0.50)+LVL(I,J) 958 ENDIF 959 375 CONTINUE 960 C*** 961 C*** MAKE SURE CLOUDS ARE DEEP ENOUGH 962 C*** 963 DO I=MYIS,MYIE 964 BCLD(I)=CV(I).GT.0..AND. 965 1 (ICVB(I)-ICVT(I)).GE.1 966 BTEMP1(I)=BCLD(I) 967 ENDDO 968 C*** 969 C*** COMPUTE CONVECTIVE CLOUD FRACTION 970 C*** 971 DO 390 I=MYIS,MYIE 972 IF(BCLD(I)) THEN 973 LML=LMH(I,J) 974 LVLIJ=LVL(I,J) 975 C 976 DO L=1,LML 977 LL=L+LVLIJ 978 IF(LL.GT.ICVB(I).OR.LL.LT.ICVT(I))THEN 979 CCMID(I,LL)=0. 980 ELSE 981 CCMID(I,LL)=CV(I) 982 ENDIF 983 CCMID(I,LL)=AMIN1(1.0,CCMID(I,LL)) 984 ENDDO 985 ENDIF 986 390 CONTINUE Page 14 Source Listing RADTN 2014-12-17 20:47 RADTN.F 987 C*** 988 C*** REMOVE HIGH CLOUDS ABOVE THE TROPOPAUSE 989 C*** 990 L400=LM 991 DO 425 I=MYIS,MYIE 992 LML=LMH(I,J) 993 LVLIJ=LVL(I,J) 994 C 995 DO L = 1, LML 996 LL=LML-L+1+LVLIJ 997 IF(PMID(I,LL).LE.40000.0)THEN 998 L400=LL 999 GO TO 400 1000 ENDIF 1001 ENDDO 1002 400 CONTINUE 1003 C 1004 LTROP=LM 1005 DO LL=L400,2,-1 1006 DTHDP=(THMID(I,LL-1)-THMID(I,LL)) 1007 1 /(PMID(I,LL-1)-PMID(I,LL)) 1008 IF(DTHDP.LT.-0.0025.OR.QMID(I,LL).LE.EPSQ1)THEN 1009 LTROP=LL 1010 GOTO 410 1011 ENDIF 1012 ENDDO 1013 410 IF(LTROP.LT.LM)THEN 1014 DO LL=LTROP,1,-1 1015 CCMID(I,LL)=0. 1016 ENDDO 1017 ENDIF 1018 425 CONTINUE 1019 C*** 1020 C----------------------------------------------------------------- 1021 ENDIF 1022 C------------------------CONVECTION-------------------------------- 1023 C*****************END OF CONVECTIVE CLOUD SECTION***************** 1024 C********************************************************************* 1025 C*** 1026 C*** DETERMINE THE FRACTIONAL CLOUD COVERAGE FOR HIGH, MID 1027 C*** AND LOW OF CLOUDS FROM THE CLOUD COVERAGE AT EACH LEVEL 1028 C*** 1029 C*** NOTE: THIS IS FOR DIAGNOSTICS ONLY!!! 1030 C*** 1031 C*** 1032 DO 500 I=MYIS,MYIE 1033 C 1034 CSTR(I)=0.0 1035 C 1036 DO L=0,LM 1037 CLDAMT(L)=0. 1038 ENDDO 1039 C 1040 C*** NOW GOES LOW, MIDDLE, HIGH 1041 C 1042 DO 480 NLVL=1,3 1043 CLDMAX=0. Page 15 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1044 MALVL=LM 1045 LLTOP=LTOP(NLVL)+LVL(I,J) 1046 C*** 1047 C*** GO TO THE NEXT CLOUD LAYER IF THE TOP OF THE CLOUD-TYPE IN 1048 C*** QUESTION IS BELOW GROUND OR IS IN THE LOWEST LAYER ABOVE GROUND. 1049 C*** 1050 IF(LLTOP.GE.LM)GO TO 480 1051 C 1052 IF(NLVL.GT.1)THEN 1053 LLBOT=LTOP(NLVL-1)-1+LVL(I,J) 1054 LLBOT=MIN(LLBOT,LM1) 1055 ELSE 1056 LLBOT=LM1 1057 ENDIF 1058 C 1059 DO 435 L=LLTOP,LLBOT 1060 CLDAMT(L)=AMAX1(CSMID(I,L),CCMID(I,L)) 1061 IF(CLDAMT(L).GT.CLDMAX)THEN 1062 MALVL=L 1063 CLDMAX=CLDAMT(L) 1064 ENDIF 1065 435 CONTINUE 1066 C********************************************************************* 1067 C NOW, CALCULATE THE TOTAL CLOUD FRACTION IN THIS PRESSURE DOMAIN 1068 C USING THE METHOD DEVELOPED BY Y.H., K.A.C. AND A.K. (NOV., 1992). 1069 C IN THIS METHOD, IT IS ASSUMED THAT SEPERATED CLOUD LAYERS ARE 1070 C RADOMLY OVERLAPPED AND ADJACENT CLOUD LAYERS ARE MAXIMUM OVERLAPPED. 1071 C VERTICAL LOCATION OF EACH TYPE OF CLOUD IS DETERMINED BY THE THICKEST 1072 C CONTINUING CLOUD LAYERS IN THE DOMAIN. 1073 C********************************************************************* 1074 CL1=0.0 1075 CL2=0.0 1076 KBT1=LLBOT 1077 KBT2=LLBOT 1078 KTH1=0 1079 KTH2=0 1080 C 1081 DO 450 LL=LLTOP,LLBOT 1082 L=LLBOT-LL+LLTOP 1083 BIT1=.FALSE. 1084 CR1=CLDAMT(L) 1085 BITX=(PINT(I,L).GE.PTOPC(NLVL+1)).AND. 1086 1 (PINT(I,L).LT.PTOPC(NLVL)).AND. 1087 2 (CLDAMT(L).GT.0.0) 1088 BIT1=BIT1.OR.BITX 1089 IF(.NOT.BIT1)GO TO 450 1090 C*** 1091 C*** BITY=T: FIRST CLOUD LAYER; BITZ=T:CONSECUTIVE CLOUD LAYER 1092 C*** NOTE: WE ASSUME THAT THE THICKNESS OF EACH CLOUD LAYER IN THE 1093 C*** DOMAIN IS LESS THAN 200 MB TO AVOID TOO MUCH COOLING OR 1094 C*** HEATING. SO WE SET CTHK(NLVL)=200*E2. BUT THIS LIMIT MAY 1095 C*** WORK WELL FOR CONVECTIVE CLOUDS. MODIFICATION MAY BE 1096 C*** NEEDED IN THE FUTURE. 1097 C*** 1098 BITY=BITX.AND.(KTH2.LE.0) 1099 BITZ=BITX.AND.(KTH2.GT.0) 1100 C Page 16 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1101 IF(BITY)THEN 1102 KBT2=L 1103 KTH2=1 1104 ENDIF 1105 C 1106 IF(BITZ)THEN 1107 KTOP1=KBT2-KTH2+1 1108 DPCL=PMID(I,KBT2)-PMID(I,KTOP1) 1109 IF(DPCL.LT.CTHK(NLVL))THEN 1110 KTH2=KTH2+1 1111 ELSE 1112 KBT2=KBT2-1 1113 ENDIF 1114 ENDIF 1115 IF(BITX)CL2=AMAX1(CL2,CR1) 1116 C*** 1117 C*** AT THE DOMAIN BOUNDARY OR SEPARATED CLD LAYERS, RANDOM OVERLAP. 1118 C*** CHOOSE THE THICKEST OR THE LARGEST FRACTION AMT AS THE CLD 1119 C*** LAYER IN THAT DOMAIN. 1120 C*** 1121 BIT2=.FALSE. 1122 BITY=BITX.AND.(CLDAMT(L-1).LE.0.0.OR. 1123 1 PINT(I,L-1).LT.PTOPC(NLVL+1)) 1124 BITZ=BITY.AND.CL1.GT.0.0 1125 BITW=BITY.AND.CL1.LE.0.0 1126 BIT2=BIT2.OR.BITY 1127 IF(.NOT.BIT2)GO TO 450 1128 C 1129 IF(BITZ)THEN 1130 KBT1=INT((CL1*KBT1+CL2*KBT2)/(CL1+CL2)) 1131 KTH1=INT((CL1*KTH1+CL2*KTH2)/(CL1+CL2))+1 1132 CL1=CL1+CL2-CL1*CL2 1133 ENDIF 1134 C 1135 IF(BITW)THEN 1136 KBT1=KBT2 1137 KTH1=KTH2 1138 CL1=CL2 1139 ENDIF 1140 C 1141 IF(BITY)THEN 1142 KBT2=LLBOT 1143 KTH2=0 1144 CL2=0.0 1145 ENDIF 1146 450 CONTINUE 1147 C*** 1148 CLDCFR(I,NLVL)=AMIN1(1.0,CL1) 1149 MTOP(I,NLVL)=MIN(KBT1,KBT1-KTH1+1) 1150 MBOT(I,NLVL)=KBT1 1151 480 CONTINUE 1152 500 CONTINUE 1153 C*** 1154 C*** SET THE UN-NEEDED TAUDAR TO ONE 1155 C*** 1156 DO I=MYIS,MYIE 1157 TAUDAR(I)=1.0 Page 17 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1158 ENDDO 1159 C---------------------------------------------------------------------- 1160 C NOW, CALCULATE THE CLOUD RADIATIVE PROPERTIES AFTER DAVIS (1982), 1161 C HARSHVARDHAN ET AL (1987) AND Y.H., K.A.C. AND A.K. (1993). 1162 C 1163 C UPDATE: THE FOLLOWING PARTS ARE MODIFIED, AFTER Y.T.H. (1994), TO 1164 C CALCULATE THE RADIATIVE PROPERTIES OF CLOUDS ON EACH MODEL 1165 C LAYER. BOTH CONVECTIVE AND STRATIFORM CLOUDS ARE USED 1166 C IN THIS CALCULATIONS. 1167 C 1168 C QINGYUN ZHAO 95-3-22 1169 C 1170 C---------------------------------------------------------------------- 1171 C 1172 C*** 1173 C*** INITIALIZE ARRAYS FOR USES LATER 1174 C*** 1175 1176 DO 600 I=MYIS,MYIE 1177 LML=LMH(I,J) 1178 LVLIJ=LVL(I,J) 1179 C 1180 C*** 1181 C*** NOTE: LAYER=1 IS THE SURFACE, AND LAYER=2 IS THE FIRST CLOUD 1182 C*** LAYER ABOVE THE SURFACE AND SO ON. 1183 C*** 1184 EMIS(I,1)=1.0 1185 KTOP(I,1)=LP1 1186 KBTM(I,1)=LP1 1187 CAMT(I,1)=1.0 1188 ITYP(I,1)=0 1189 KCLD(I)=2 1190 C 1191 DO NBAND=1,NB 1192 RRCL(I,NBAND,1)=0.0 1193 TTCL(I,NBAND,1)=1.0 1194 ENDDO 1195 C 1196 DO 510 L=2,LP1 1197 ITYP(I,L)=0 1198 CAMT(I,L)=0.0 1199 KTOP(I,L)=1 1200 KBTM(I,L)=1 1201 EMIS(I,L)=0.0 1202 C 1203 DO NBAND=1,NB 1204 RRCL(I,NBAND,L)=0.0 1205 TTCL(I,NBAND,L)=1.0 1206 ENDDO 1207 510 CONTINUE 1208 C*** 1209 C*** NOW CALCULATE THE AMOUNT, TOP, BOTTOM AND TYPE OF EACH CLOUD LAYER 1210 C*** CLOUD TYPE=1: STRATIFORM CLOUD 1211 C*** TYPE=2: CONVECTIVE CLOUD 1212 C*** WHEN BOTH CONVECTIVE AND STRATIFORM CLOUDS EXIST AT THE SAME POINT, 1213 C*** SELECT CONVECTIVE CLOUD (TYPE=2),IN OTHER WORDS, CONVECTIVE CLOUDS 1214 C*** HAVE THE HIGHER PRIORITY THAN STRATIFORM CLOUDS. Page 18 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1215 C*** CLOUD LAYERS ARE SEPARATED BY: 1216 C*** 1. NO-CLOUD LAYER 1217 C*** 2. DIFFERENT CLOUD TYPE 1218 C*** NOTE: THERE IS ONLY ONE CONVECTIVE CLOUD LAYER IN ONE COLUMN. 1219 C*** KTOP AND KBTM ARE THE TOP AND BOTTOM OF EACH CLOUD LAYER IN TERMS O 1220 C*** ETA MODEL LEVEL. 1221 C*** 1222 DO 540 L=2,LML 1223 LL=LML-L+1+LVLIJ 1224 BITC=CCMID(I,LL).GT.0.1 1225 BITS=CSMID(I,LL).GT.0.1 1226 BITCP1=CCMID(I,LL+1).GT.0.1 1227 BITSP1=CSMID(I,LL+1).GT.0.1 1228 BIT1=BITS.OR.BITC 1229 C------------------- 1230 IF(BIT1)THEN 1231 C------------------- 1232 IF(ITYP(I,KCLD(I)).EQ.0)THEN 1233 CAMT(I,KCLD(I))=CSMID(I,LL) 1234 ITYP(I,KCLD(I))=1 1235 KBTM(I,KCLD(I))=LL 1236 C 1237 IF(BITC)THEN 1238 CAMT(I,KCLD(I))=CCMID(I,LL) 1239 ITYP(I,KCLD(I))=2 1240 ENDIF 1241 ELSE 1242 IF(BITC)THEN 1243 IF(BITCP1)THEN 1244 CAMT(I,KCLD(I))=AMAX1(CAMT(I,KCLD(I)),CCMID(I,LL)) 1245 ELSE 1246 KCLD(I)=KCLD(I)+1 1247 CAMT(I,KCLD(I))=CCMID(I,LL) 1248 ITYP(I,KCLD(I))=2 1249 KTOP(I,KCLD(I)-1)=LL+1 1250 KBTM(I,KCLD(I))=LL 1251 ENDIF 1252 ELSE 1253 IF(BITCP1)THEN 1254 KCLD(I)=KCLD(I)+1 1255 CAMT(I,KCLD(I))=CSMID(I,LL) 1256 ITYP(I,KCLD(I))=1 1257 KTOP(I,KCLD(I)-1)=LL+1 1258 KBTM(I,KCLD(I))=LL 1259 ELSE 1260 CAMT(I,KCLD(I))=AMAX1(CAMT(I,KCLD(I)),CSMID(I,LL)) 1261 ENDIF 1262 ENDIF 1263 ENDIF 1264 C------------------- 1265 ELSE 1266 C------------------- 1267 IF(BITCP1.OR.BITSP1)THEN 1268 KCLD(I)=KCLD(I)+1 1269 KTOP(I,KCLD(I)-1)=LL+1 1270 ITYP(I,KCLD(I))=0 1271 CAMT(I,KCLD(I))=0.0 Page 19 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1272 ENDIF 1273 C------------------- 1274 ENDIF 1275 C------------------- 1276 540 CONTINUE 1277 C*** 1278 C*** THE REAL NUMBER OF CLOUD LAYERS IS (THE FIRST IS THE GROUNG; 1279 C*** THE LAST IS THE SKY): 1280 C*** 1281 NCLDS(I)=KCLD(I)-2 1282 NCLD=NCLDS(I) 1283 C*** 1284 C*** NOW CALCULATE CLOUD RADIATIVE PROPERTIES 1285 C*** 1286 IF(NCLD.GE.1)THEN 1287 C*** 1288 C*** NOTE: THE FOLLOWING CALCULATIONS, THE UNIT FOR PRESSURE IS MB!!! 1289 C*** 1290 DO 580 NC=2,NCLD+1 1291 C 1292 TAUC(I)=0.0 1293 QSUM=0.0 1294 NKTP=LP1 1295 NBTM=0 1296 BITX=CAMT(I,NC).GT.0.1 1297 NKTP=MIN(NKTP,KTOP(I,NC)) 1298 NBTM=MAX(NBTM,KBTM(I,NC)) 1299 C 1300 DO 560 LL=NKTP,NBTM 1301 IF(LL.GE.KTOP(I,NC).AND.LL.LE.KBTM(I,NC).AND.BITX)THEN 1302 PRS1=PINT(I,LL)*0.01 1303 PRS2=PINT(I,LL+1)*0.01 1304 DELP=PRS2-PRS1 1305 TCLD=TMID(I,LL)-273.16 1306 QSUM=QSUM+QMID(I,LL)*DELP*(PRS1+PRS2) 1307 1 /(120.1612*SQRT(TMID(I,LL))) 1308 C*** 1309 C*** FOR CONVECTIVE CLOUD OR STARTIFORM CLOUD WITH TOP ABOVE 500MB 1310 C*** 1311 IF(ITYP(I,NC).EQ.2 1312 1 .OR.PINT(I,KTOP(I,NC)).LE.PTOPC(3))THEN 1313 IF(TCLD.LE.-10.0)THEN 1314 TAUC(I)=TAUC(I)+DELP*AMAX1(0.1E-3, 1315 1 2.0E-6*(TCLD+82.5)**2) 1316 ELSE 1317 TAUC(I)=TAUC(I)+DELP*AMIN1(0.08,6.949E-3*TCLD+0.1) 1318 ENDIF 1319 ELSE 1320 C*** 1321 C*** FOR LOW AND MID STRATIFORM CLOUDS 1322 C*** 1323 IF(TCLD.LE.-20.0)THEN 1324 TAUC(I)=TAUC(I)+DELP*AMAX1(0.1E-3,2.56E-5* 1325 1 (TCLD+82.5)**2) 1326 ELSE 1327 TAUC(I)=TAUC(I)+DELP*0.1 1328 ENDIF Page 20 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1329 ENDIF 1330 ENDIF 1331 560 CONTINUE 1332 C 1333 IF(BITX)EMIS(I,NC)=1.0-EXP(-0.75*TAUC(I)) 1334 IF(QSUM.GE.EPSQ1)THEN 1335 C 1336 DO 570 NBAND=1,NB 1337 IF(BITX)THEN 1338 PROD=ABCFF(NBAND)*QSUM 1339 DDX=TAUC(I)/(TAUC(I)+PROD) 1340 EEX=1.0-DDX 1341 IF(ABS(EEX).GE.1.E-8)THEN 1342 DD=DDX 1343 EE=EEX 1344 FF=1.0-DD*0.85 1345 AA=MIN(50.0,SQRT(3.0*EE*FF)*TAUC(I)) 1346 AA=EXP(-AA) 1347 BB=FF/EE 1348 GG=SQRT(BB) 1349 DD=(GG+1.0)*(GG+1.0)-(GG-1.0)*(GG-1.0)*AA*AA 1350 RRCL(I,NBAND,NC)=MAX(0.1E-5,(BB-1.0)*(1.0-AA*AA)/DD) 1351 TTCL(I,NBAND,NC)=AMAX1(0.1E-5,4.0*GG*AA/DD) 1352 ENDIF 1353 ENDIF 1354 570 CONTINUE 1355 ENDIF 1356 580 CONTINUE 1357 C 1358 ENDIF 1359 C 1360 600 CONTINUE 1361 C********************************************************************* 1362 C****************** COMPUTE OZONE AT MIDLAYERS ********************* 1363 C********************************************************************* 1364 C 1365 C*** MODIFY PRESSURES SO THAT THE ENTIRE COLUMN OF OZONE (TO 0 MB) 1366 C*** IS INCLUDED IN THE MODEL COLUMN EVEN WHEN PT > 0 MB 1367 C*** 1368 DO L=1,LM 1369 DO I=MYIS,MYIE 1370 DENOM=1./(PINT(I,LP1)-PINT(I,1)) 1371 FCTRA=PINT(I,LP1)*DENOM 1372 FCTRB=-PINT(I,1)*PINT(I,LP1)*DENOM 1373 POZN(I,L)=PMID(I,L)*FCTRA+FCTRB 1374 ENDDO 1375 ENDDO 1376 C 1377 CALL OZON2D(LM,POZN,XLAT,RSIN1,RCOS1,RCOS2,OZN) 1378 C 1379 C*** 1380 C*** NOW THE VARIABLES REQUIRED BY RADFS HAVE BEEN CALCULATED. 1381 C*** 1382 C---------------------------------------------------------------------- 1383 C*** 1384 C*** CALL THE GFDL RADIATION DRIVER 1385 C*** Page 21 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1386 C*** 1387 CALL RADFS 1388 c 1 (PSFC,PMID,PINT,QMID,TMID,OZN,TSKN,SLMSK,ALBEDO,XLAT 1389 1 (PSFC,PMID,PINT,QMID,TMID,OZN,TSKN,SLMSK,ALBDO,XLAT 1390 2, CAMT,ITYP,KTOP,KBTM,NCLDS,EMIS,RRCL,TTCL 1391 3, COSZ,TAUDAR,1 1392 4, 1,0 1393 5, ETAD,AETA,ITIMSW,ITIMLW,JD,R1,HOUR,TENDS,TENDL 1394 6, FLWUP,FSWUP,FSWDN,FSWDNS,FSWUPS,FLWDNS,FLWUPS) 1395 C---------------------------------------------------------------------- 1396 DO 650 I=MYIS,MYIE 1397 PDSLIJ=PDSL(I,J) 1398 PMOD=CUPPT(I,J)*24.0*1000.0/CLSTP 1399 CFRACL(I,J)=CLDCFR(I,1) 1400 CFRACM(I,J)=CLDCFR(I,2) 1401 CFRACH(I,J)=CLDCFR(I,3) 1402 C 1403 C*** ARRAYS ACFRST AND ACFRCV ACCUMULATE AVERAGE STRATIFORM AND 1404 C*** CONVECTIVE CLOUD FRACTIONS, RESPECTIVELY. THIS INFORMATION 1405 C*** IS PASSED TO THE POST PROCESSOR VIA COMMON BLOCK ACMCLD. 1406 C 1407 CFRAVG=AMAX1(CFRACL(I,J),AMAX1(CFRACM(I,J),CFRACH(I,J))) 1408 IF(CNCLD)THEN 1409 IF(PMOD.LE.PPT(1))THEN 1410 ACFRST(I,J)=ACFRST(I,J)+CFRAVG 1411 NCFRST(I,J)=NCFRST(I,J)+1 1412 ELSE 1413 ACFRCV(I,J)=ACFRCV(I,J)+CFRAVG 1414 NCFRCV(I,J)=NCFRCV(I,J)+1 1415 ENDIF 1416 ELSE 1417 ACFRST(I,J)=ACFRST(I,J)+CFRAVG 1418 NCFRST(I,J)=NCFRST(I,J)+1 1419 ENDIF 1420 650 CONTINUE 1421 C*** 1422 C*** COLLECT ATMOSPHERIC TEMPERATURE TENDENCIES DUE TO RADIATION. 1423 C*** ALSO COLLECT THE TOTAL SW AND INCOMING LW RADIATION (W/M**2) 1424 C*** AND CONVERT TO FORM NEEDED FOR PREDICTION OF THS IN SURFCE. 1425 C*** 1426 DO 660 I=MYIS,MYIE 1427 DO L=1,LM 1428 LL=LVL(I,J)+L 1429 IF(SHORT)RSWTT(I,J,L)=TENDS(I,LL) 1430 IF(LONG) RLWTT(I,J,L)=TENDL(I,LL) 1431 c if(mype.eq.13.and.i.eq.40.and.j.eq.77) then 1432 c print*,'l,ll,tendl(i,ll),rlwtt(i,j,l)=', 1433 c * l,ll,tendl(i,ll),rlwtt(i,j,l) 1434 c endif 1435 IF(LL.EQ.LM)GO TO 660 1436 ENDDO 1437 660 CONTINUE 1438 C*** 1439 C*** SUM THE LW INCOMING AND SW RADIATION (W/M**2) FOR RADIN. 1440 C*** 1441 DO 675 I=MYIS,MYIE 1442 IF(LONG)THEN Page 22 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1443 SIGT4(I,J)=STBOL*TMID(I,LM)*TMID(I,LM)* 1444 1 TMID(I,LM)*TMID(I,LM) 1445 ENDIF 1446 C 1447 C*** ACCUMULATE VARIOUS LW AND SW RADIATIVE FLUXES FOR POST 1448 C*** PROCESSOR. PASSED VIA COMMON ACMRDL AND ACMRDS. 1449 C 1450 IF(LONG)THEN 1451 RLWIN(I,J) =FLWDNS(I) 1452 RLWOUT(I,J)=FLWUPS(I) 1453 RLWTOA(I,J)=FLWUP(I) 1454 ENDIF 1455 IF(SHORT)THEN 1456 RSWIN(I,J) =FSWDNS(I) 1457 RSWOUT(I,J)=FSWUPS(I) 1458 RSWTOA(I,J)=FSWUP(I) 1459 ENDIF 1460 675 CONTINUE 1461 C*** 1462 C*** THIS ROW IS FINISHED. GO TO NEXT 1463 C*** 1464 C ********************* 1465 700 CONTINUE 1466 C ********************* 1467 C---------------------------------------------------------------------- 1468 C*** 1469 C*** CALLS TO RADIATION THIS TIME STEP ARE COMPLETE. 1470 C*** 1471 C---------------------------------------------------------------------- 1472 C---------------------------------------------------------------------- 1473 C*** 1474 C*** HORIZONTAL SMOOTHING OF TEMPERATURE TENDENCIES 1475 C*** 1476 C---------------------------------------------------------------------- 1477 IF(SHORT) THEN 1478 DO 800 L=1,LM 1479 CALL ZERO2(TL) 1480 CALL ZERO2(FNE) 1481 CALL ZERO2(FSE) 1482 C 1483 IF(KSMUD.GE.1)THEN 1484 DO 750 KS=1,KSMUD 1485 C 1486 DO J=MYJS,MYJE 1487 DO I=MYIS,MYIE 1488 TL(I,J)=RSWTT(I,J,L)*HTM(I,J,L) 1489 ENDDO 1490 ENDDO 1491 C 1492 DO J=MYJS,MYJE 1493 DO I=MYIS,MYIE 1494 FNE(I,J)=(TL(I+IHE(J),J+1)-TL(I,J)) 1495 1 *HTM(I,J,L)*HTM(I+IHE(J),J+1,L) 1496 ENDDO 1497 ENDDO 1498 C 1499 DO J=MYJS1,MYJE Page 23 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1500 DO I=MYIS,MYIE 1501 FSE(I,J)=(TL(I+IHE(J),J-1)-TL(I,J)) 1502 1 *HTM(I+IHE(J),J-1,L)*HTM(I,J,L) 1503 ENDDO 1504 ENDDO 1505 C 1506 DO J=MYJS2,MYJE2 1507 DO I=MYIS,MYIE 1508 TL(I,J)=(FNE(I,J)-FNE(I+IHW(J),J-1) 1509 1 +FSE(I,J)-FSE(I+IHW(J),J+1)) 1510 2 *HBM2(I,J)*0.125+TL(I,J) 1511 ENDDO 1512 ENDDO 1513 C 1514 DO J=MYJS,MYJE 1515 DO I=MYIS,MYIE 1516 RSWTT(I,J,L)=TL(I,J) 1517 ENDDO 1518 ENDDO 1519 C 1520 750 CONTINUE 1521 ENDIF 1522 C 1523 800 CONTINUE 1524 ENDIF 1525 C---------------------------------------------------------------------- 1526 C 1527 IF(LONG)THEN 1528 C 1529 DO 900 L=1,LM 1530 CALL ZERO2(TL) 1531 CALL ZERO2(FNE) 1532 CALL ZERO2(FSE) 1533 C 1534 IF(KSMUD.GE.1)THEN 1535 DO 850 KS=1,KSMUD 1536 C 1537 DO J=MYJS,MYJE 1538 DO I=MYIS,MYIE 1539 TL(I,J)=RLWTT(I,J,L)*HTM(I,J,L) 1540 ENDDO 1541 ENDDO 1542 C 1543 DO J=MYJS,MYJE1 1544 DO I=MYIS,MYIE 1545 FNE(I,J)=(TL(I+IHE(J),J+1)-TL(I,J)) 1546 1 *HTM(I,J,L)*HTM(I+IHE(J),J+1,L) 1547 ENDDO 1548 ENDDO 1549 C 1550 DO J=MYJS1,MYJE 1551 DO I=MYIS,MYIE 1552 FSE(I,J)=(TL(I+IHE(J),J-1)-TL(I,J)) 1553 1 *HTM(I+IHE(J),J-1,L)*HTM(I,J,L) 1554 ENDDO 1555 ENDDO 1556 C Page 24 Source Listing RADTN 2014-12-17 20:47 RADTN.F 1557 DO J=MYJS2,MYJE2 1558 DO I=MYIS,MYIE 1559 TL(I,J)=(FNE(I,J)-FNE(I+IHW(J),J-1) 1560 1 +FSE(I,J)-FSE(I+IHW(J),J+1)) 1561 2 *HBM2(I,J)*0.125+TL(I,J) 1562 ENDDO 1563 ENDDO 1564 C 1565 DO J=MYJS,MYJE 1566 DO I=MYIS,MYIE 1567 RLWTT(I,J,L)=TL(I,J) 1568 ENDDO 1569 ENDDO 1570 C 1571 850 CONTINUE 1572 ENDIF 1573 900 CONTINUE 1574 ENDIF 1575 C----------------------------------------------------------------------- 1576 C----------------------------------------------------------------------- 1577 RETURN 1578 END ENTRY POINTS Name radtn_ Page 25 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 125 Label 378 363 140 Label 386 380 200 Label 431 391 250 Label 455 437 275 Label 537 470 280 Label 538 466 290 Label 561 558 300 Label 582 579 310 Label 600 597 320 Label 612 609 330 Label 629 626 340 Label 641 638 350 Label 646 547 375 Label 710 670 390 Label 737 722 400 Label 753 750 410 Label 764 761 425 Label 769 742 435 Label 816 810 450 Label 897 832,840,878 480 Label 902 793,801 500 Label 903 783 510 Label 958 947 540 Label 1027 973 560 Label 1082 1051 570 Label 1105 1087 580 Label 1107 1041 600 Label 1111 927 650 Label 1171 1147 660 Label 1188 1177,1186 675 Label 1211 1192 700 Label 1216 360 750 Label 1271 1235 800 Label 1274 1229 850 Label 1322 1286 900 Label 1324 1280 A1 Param 122 R(4) 4 scalar A2 Param 122 R(4) 4 scalar 487 A3 Param 122 R(4) 4 scalar 487 A4 Param 122 R(4) 4 scalar 487 AA Local 1096 R(4) 4 scalar 1096,1097,1100,1101,1102 ABCFF Scalar 186 R(4) 4 1 12 COM 1089 ABS Func 1092 scalar 1092 ACCLIQ Scalar 13 R(4) 4 2 17955 COM ACFRCV Scalar 4 R(4) 4 2 17955 COM 1164 ACFRST Scalar 5 R(4) 4 2 17955 COM 1161,1168 ACMCLD Common 2 287288 ACMRDL Common 2 3662832 ACMRDS Common 2 3662832 ACPREC Scalar 12 R(4) 4 2 17955 COM ADDL Local 298 R(4) 4 scalar 298,299,300 Page 26 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References AETA Scalar 9 R(4) 4 1 45 COM 400,1144 AETAD Scalar 4 R(4) 4 1 45 COM AFSI Scalar 26 R(4) 4 2 17955 COM AI Local 478 R(4) 4 scalar 478,482,488 AKHS Scalar 7 R(4) 4 2 17955 COM AKMS Scalar 7 R(4) 4 2 17955 COM ALBASE Scalar 20 R(4) 4 2 17955 COM ALBDO Local 195 R(4) 4 1 45 451,1140 ALBEDO Scalar 21 R(4) 4 2 17955 COM 451 ALWIN Scalar 5 R(4) 4 2 17955 COM ALWOUT Scalar 6 R(4) 4 2 17955 COM ALWTOA Scalar 6 R(4) 4 2 17955 COM AMAX1 Func 532 scalar 532,700,811,866,995,1011,1065,1075 ,1102,1158 AMIN1 Func 477 scalar 477,488,529,535,701,734,899,1068 APES Local 439 R(4) 4 scalar 439,440,442 ARDLW Scalar 3 R(4) 4 scalar COM ARDSW Scalar 3 R(4) 4 scalar COM ARG Local 531 R(4) 4 scalar 531,532,533 ASWIN Scalar 5 R(4) 4 2 17955 COM ASWOUT Scalar 6 R(4) 4 2 17955 COM ASWTOA Scalar 6 R(4) 4 2 17955 COM BB Local 1098 R(4) 4 scalar 1098,1099,1101 BCLD Local 148 L(4) 4 1 45 715,717,723 BI Local 479 R(4) 4 scalar 479,483,488 BIT1 Local 149 L(4) 4 scalar 834,839,840,979,981 BIT2 Local 149 L(4) 4 scalar 872,877,878 BITC Local 149 L(4) 4 scalar 975,979,988,993 BITCP1 Local 149 L(4) 4 scalar 977,994,1004,1018 BITS Local 149 L(4) 4 scalar 976,979 BITSP1 Local 149 L(4) 4 scalar 978,1018 BITW Local 149 L(4) 4 scalar 876,886 BITX Local 149 L(4) 4 scalar 836,839,849,850,866,873,1047,1052, 1084,1088 BITY Local 149 L(4) 4 scalar 849,852,873,875,876,877,892 BITZ Local 149 L(4) 4 scalar 850,857,875,880 BMID Local 232 R(4) 4 1 45 396,472 BTEMP1 Local 148 L(4) 4 1 45 717 CALL1 Local 147 L(4) 4 scalar CAMT Local 224 R(4) 4 2 2070 938,949,984,989,995,998,1006,1011, 1022,1047,1141 CAPA Param 118 R(4) 4 scalar 402,421,439 CC Local 218 R(4) 4 1 9 240,691,693,694 CC1 Local 693 R(4) 4 scalar 693,697 CC2 Local 694 R(4) 4 scalar 694,697 CCLIMIT Local 463 R(4) 4 scalar 463 CCMID Local 233 R(4) 4 2 2025 370,730,732,734,766,811,975,977,98 9,995,998 CCR Local 230 R(4) 4 2 2025 372,529,533,535 CFCO2 Scalar 186 R(4) 4 scalar COM CFO3 Scalar 186 R(4) 4 scalar COM CFRACH Scalar 4 R(4) 4 2 17955 COM 1152,1158 CFRACL Scalar 3 R(4) 4 2 17955 COM 1150,1158 CFRACM Scalar 4 R(4) 4 2 17955 COM 1151,1158 CFRAVG Local 1158 R(4) 4 scalar 1158,1161,1164,1168 Page 27 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References CI Scalar 6 R(4) 4 scalar COM CL1 Local 825 R(4) 4 scalar 825,875,876,881,882,883,889,899 CL2 Local 826 R(4) 4 scalar 826,866,881,882,883,889,895 CLAPSE Param 133 R(4) 4 scalar 577 CLDAMT Local 192 R(4) 4 1 46 788,811,812,814,835,838,873 CLDCFR Local 196 R(4) 4 2 135 382,899,1150,1151,1152 CLDEFI Scalar 11 R(4) 4 2 17955 COM CLDF Local 197 R(4) 4 2 2070 CLDMAX Local 794 R(4) 4 scalar 794,812,814 CLDRH0 Param 131 R(4) 4 scalar 131 CLDRH2 Param 132 R(4) 4 scalar 132 CLDWTR Common 2 3447720 CLIMIT Local 464 R(4) 4 scalar 464,495,508 CLPFIL Local 589 R(4) 4 scalar 589,591,598 CLPSE Param 133 R(4) 4 scalar 588,589 CLSTP Local 272 R(4) 4 scalar 272,678,1149 CM1 Param 134 R(4) 4 scalar CM2 Param 134 R(4) 4 scalar CM3 Param 134 R(4) 4 scalar CMC Scalar 8 R(4) 4 2 17955 COM CNCLD Local 152 L(4) 4 scalar 244,668,1159 CNVBOT Scalar 18 R(4) 4 2 17955 COM CNVCLD Common 2 287280 CNVTOP Scalar 18 R(4) 4 2 17955 COM COS Func 302 scalar 302,303 COSZ Local 195 R(4) 4 1 45 454,1142 CPGFUX Scalar 9 R(4) 4 2 17955 COM CPGFV Scalar 3 R(4) 4 scalar COM CR1 Local 835 R(4) 4 scalar 835,866 CS Scalar 6 R(4) 4 scalar COM CSMID Local 230 R(4) 4 2 2025 368,535,559,586,598,610,643,811,97 6,978,984,1006,1011 CSTR Local 220 R(4) 4 1 45 785 CTHK Scalar 188 R(4) 4 1 3 COM 860 CTLBLK Common 4 112 CUINIT Common 182 4 CUPPT Scalar 3 R(4) 4 2 17955 COM 678,1149 CUPREC Scalar 13 R(4) 4 2 17955 COM CURAD Scalar 182 L(4) 4 scalar COM 289 CURVX Scalar 9 R(4) 4 2 17955 COM CV Local 204 R(4) 4 1 45 689,691,697,700,701,703,715,732 CVB Local 222 R(4) 4 1 45 CVT Local 222 R(4) 4 1 45 CWM Scalar 3 R(4) 4 3 807975 COM 406 CWMKL Local 475 R(4) 4 scalar 475,495,508,531 CZEN Scalar 16 R(4) 4 2 17955 COM 319,320 CZMEAN Scalar 23 R(4) 4 2 17955 COM 308,320,329,454 DAYI Local 296 R(4) 4 scalar 296,297,300,315 DCLPS Param 133 R(4) 4 scalar 589 DD Local 1093 R(4) 4 scalar 1093,1095,1100,1101,1102 DDMPUX Scalar 12 R(4) 4 2 17955 COM DDMPVX Scalar 12 R(4) 4 2 17955 COM DDP Local 557 R(4) 4 scalar 557,558,596,597,608,609 DDX Local 144 R(8) 8 scalar 1090,1091,1093 DELP Local 1055 R(4) 4 scalar 1055,1057,1065,1068,1075,1078 Page 28 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References DENOM Local 1121 R(4) 4 scalar 1121,1122,1123 DETA Scalar 9 R(4) 4 1 45 COM DETAD Scalar 4 R(4) 4 1 45 COM DFL Scalar 4 R(4) 4 1 46 COM DFRLG Scalar 9 R(4) 4 1 46 COM DI Scalar 6 R(4) 4 scalar COM DPCL Local 859 R(4) 4 scalar 859,860 DS Scalar 6 R(4) 4 scalar COM DT Scalar 6 R(4) 4 scalar COM 270,295,314 DTD Scalar 5 R(4) 4 scalar COM DTHDP Local 575 R(4) 4 scalar 575,577,588,589,633,636,757,759 DTQ2 Scalar 5 R(4) 4 scalar COM DXX Scalar 8 R(4) 4 2 17955 COM DY Scalar 3 R(4) 4 scalar COM DYNAM Common 2 725404 EE Local 1094 R(4) 4 scalar 1094,1096,1098 EEX Local 144 R(8) 8 scalar 1091,1092,1094 EF4T Scalar 3 R(4) 4 scalar COM EM Scalar 6 R(4) 4 1 760 COM EMIS Local 221 R(4) 4 2 2070 935,952,1084,1141 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 134 R(4) 4 scalar EPSO3 Param 130 R(4) 4 scalar 374 EPSQ Param 130 R(4) 4 scalar 367 EPSQ1 Param 130 R(4) 4 scalar 520,636,759,1085 EPSR Scalar 13 R(4) 4 2 17955 COM ETAD Scalar 4 R(4) 4 1 46 COM 401,1144 EXNER Local 402 R(4) 4 scalar 402,404,421,422 EXP Func 487 scalar 487,533,1084,1097 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 FADX Scalar 13 R(4) 4 2 17955 COM FCP Scalar 10 R(4) 4 2 17955 COM FCTRA Local 1122 R(4) 4 scalar 1122,1124 FCTRB Local 1123 R(4) 4 scalar 1123,1124 FDIVX Scalar 10 R(4) 4 2 17955 COM FF Local 1095 R(4) 4 scalar 1095,1096,1098 FIQ Local 494 R(4) 4 scalar 494,495 FIRST Scalar 147 L(4) 4 scalar COM FIS Scalar 7 R(4) 4 2 17955 COM FIW Local 515 R(4) 4 scalar 515,516 FLOAT Func 515 scalar 515 FLWDNS Local 208 R(4) 4 1 45 1145,1202 FLWUP Local 206 R(4) 4 1 45 1145,1204 FLWUPS Local 208 R(4) 4 1 45 1145,1203 FNE Local 210 R(4) 4 2 17955 1231,1245,1259,1282,1296,1310 FSE Local 210 R(4) 4 2 17955 1232,1252,1260,1283,1303,1311 FSWDN Local 206 R(4) 4 1 45 1145 FSWDNS Local 207 R(4) 4 1 45 1145,1207 FSWUP Local 206 R(4) 4 1 45 1145,1209 FSWUPS Local 207 R(4) 4 1 45 1145,1208 FX Scalar 11 R(4) 4 2 17955 COM Page 29 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GFFC Scalar 19 R(4) 4 2 17955 COM GG Local 1099 R(4) 4 scalar 1099,1100,1102 GLAT Scalar 15 R(4) 4 2 17955 COM 453 GLB_TABLE Common 37 128 GLON Scalar 15 R(4) 4 2 17955 COM GRNFLX Scalar 5 R(4) 4 2 17955 COM HBM2 Scalar 8 R(4) 4 2 17955 COM 396,1261,1312 HBM3 Scalar 8 R(4) 4 2 17955 COM HBOT Scalar 17 R(4) 4 2 17955 COM 671,708 HDAC Scalar 22 R(4) 4 2 17955 COM HDACV Scalar 22 R(4) 4 2 17955 COM HH Local 472 R(4) 4 scalar 472,476,477,487 HMID Local 231 R(4) 4 2 2025 373,407,472 HOUR Local 296 R(4) 4 scalar 296,315,1144 HPINC Param 130 R(4) 4 scalar 419,420,429 HTM Scalar 6 R(4) 4 3 807975 COM 407,1239,1246,1253,1290,1297,1304 HTMG Scalar 52 R(4) 4 3 4183695 COM HTOP Scalar 17 R(4) 4 2 17955 COM 671,707 I Local 256 I(4) 4 scalar 256,257,307,308,309,318,319,320,32 1,328,329,364,365,366,367,368,369, 370,371,372,373,374,375,376,381,38 2,383,384,391,392,393,394,396,397, 400,401,402,403,404,405,406,407,41 9,420,421,422,427,429,437,438,439, 440,441,442,444,451,453,454,466,46 7,468,472,473,474,475,486,492,496, 498,503,507,508,511,529,533,535,54 7,549,550,554,557,559,568,575,576, 586,596,598,608,610,624,633,634,63 6,643,670,671,673,675,678,689,691, 697,700,701,703,704,705,707,708,71 4,715,716,717,722,723,724,725,729, 730,732,734,742,743,744,748,757,75 8,759,766,783,785,796,804,811,836, 837,859,874,899,900,901,907,908,92 7,928,929,935,936,937,938,939,940, 943,944,948,949,950,951,952,955,95 6,975,976,977,978,983,984,985,986, 989,990,995,997,998,999,1000,1001, 1005,1006,1007,1008,1009,1011,1019 ,1020,1021,1022,1032,1033,1043,104 7,1048,1049,1052,1053,1054,1056,10 57,1058,1062,1063,1065,1068,1075,1 078,1084,1090,1096,1101,1102,1120, 1121,1122,1123,1124,1147,1148,1149 ,1150,1151,1152,1158,1161,1162,116 4,1165,1168,1169,1177,1179,1180,11 81,1192,1194,1195,1202,1203,1204,1 207,1208,1209,1238,1239,1244,1245, 1246,1251,1252,1253,1258,1259,1260 ,1261,1266,1267,1289,1290,1295,129 6,1297,1302,1303,1304,1309,1310,13 11,1312,1317,1318 Page 30 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References IBROW Scalar 22 I(4) 4 scalar COM ICHUNKTAB Scalar 29 I(4) 4 1 8 COM ICVB Local 203 I(4) 4 1 45 705,708,716,729 ICVT Local 203 I(4) 4 1 45 704,707,716,729 IDAT Scalar 5 I(4) 4 1 3 COM 299 IDIM1 Param 19 I(4) 4 scalar 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,20,21,22,23,26,148,189,1 93,195,196,197,198,200,201,202,203 ,204,206,207,208,209,210,211,213,2 14,215,216,220,221,222,224,225,226 ,227,230,231,232,233 IDIM2 Param 19 I(4) 4 scalar 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,20,21,22,23,26,148,189,1 93,195,196,197,198,200,201,202,203 ,204,206,207,208,209,210,211,213,2 14,215,216,220,221,222,224,225,226 ,227,230,231,232,233 IDTAD Scalar 6 I(4) 4 scalar COM IE_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 IHE Scalar 3 I(4) 4 1 399 COM 1245,1246,1252,1253,1296,1297,1303 ,1304 IHEG Scalar 9 I(4) 4 1 387 COM IHHA Scalar 3 I(4) 4 1 760 COM IHLA Scalar 3 I(4) 4 1 760 COM IHRST Scalar 5 I(4) 4 scalar COM IHW Scalar 3 I(4) 4 1 399 COM 1259,1260,1310,1311 IHWG Scalar 9 I(4) 4 1 387 COM II Local 313 I(4) 4 scalar 313,314 ILCOL Scalar 22 I(4) 4 scalar COM ILPAD1 Scalar 23 I(4) 4 scalar COM ILPAD2 Scalar 23 I(4) 4 scalar COM ILPAD3 Scalar 23 I(4) 4 scalar COM ILPAD4 Scalar 23 I(4) 4 scalar COM ILPAD5 Scalar 23 I(4) 4 scalar COM IM Param 3 I(4) 4 scalar 9,16,19,43,44,45,46,47,51,52,57,12 7 IMJM Param 127 I(4) 4 scalar INDX Common 2 6564 INDXG Common 8 8084 INPES Param 11 I(4) 4 scalar 16,19,27,28,29,38,39 INT Func 297 scalar 297,707,708,881,882 INUMQ Scalar 34 I(4) 4 1 100 COM IOUT Scalar 5 I(4) 4 scalar COM IQUILT_GROUP Scalar 33 I(4) 4 scalar COM IR Local 365 I(4) 4 scalar 365,392 IRAD Scalar 5 I(4) 4 1 45 COM 365,392 IRADG Scalar 9 I(4) 4 1 473 COM IRCOL Scalar 22 I(4) 4 scalar COM IRPAD1 Scalar 24 I(4) 4 scalar COM IRPAD2 Scalar 24 I(4) 4 scalar COM IRPAD3 Scalar 24 I(4) 4 scalar COM IRPAD4 Scalar 24 I(4) 4 scalar COM Page 31 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References IRPAD5 Scalar 24 I(4) 4 scalar COM ISLOPE Scalar 7 I(4) 4 2 17955 COM ISLTYP Scalar 6 I(4) 4 2 17955 COM IS_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM ITAIL Param 16 I(4) 4 scalar ITB Param 1 I(4) 4 scalar 10,25 ITBQ Param 1 I(4) 4 scalar 11,25 ITEMP Scalar 47 I(4) 4 2 91719 COM ITEMP2 Scalar 47 I(4) 4 2 91719 COM ITIMLW Local 282 I(4) 4 scalar 282,284,1144 ITIMSW Local 281 I(4) 4 scalar 281,283,1144 ITROW Scalar 22 I(4) 4 scalar COM ITYP Local 225 I(4) 4 2 2070 939,948,983,985,990,999,1007,1021, 1062,1141 IVE Scalar 4 I(4) 4 1 399 COM IVEG Scalar 9 I(4) 4 1 387 COM IVGTYP Scalar 6 I(4) 4 2 17955 COM IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IVW Scalar 4 I(4) 4 1 399 COM IVWG Scalar 9 I(4) 4 1 387 COM IW Local 230 I(4) 4 2 2025 371,496,498,503,507,508,511 IWKL Local 511 I(4) 4 scalar 511,515 J Local 255 I(4) 4 scalar 255,257,306,308,309,317,319,320,32 1,327,329,360,366,393,394,396,397, 400,401,403,404,405,406,407,438,44 0,441,442,444,451,453,454,467,468, 549,550,554,568,671,678,707,708,72 4,725,743,744,796,804,928,929,1148 ,1149,1150,1151,1152,1158,1161,116 2,1164,1165,1168,1169,1179,1180,11 81,1194,1202,1203,1204,1207,1208,1 209,1237,1239,1243,1245,1246,1250, 1252,1253,1257,1259,1260,1261,1265 ,1267,1288,1290,1294,1296,1297,130 1,1303,1304,1308,1310,1311,1312,13 16,1318 JAM Param 127 I(4) 4 scalar 3,6 JBPAD1 Scalar 25 I(4) 4 scalar COM JBPAD2 Scalar 25 I(4) 4 scalar COM JBPAD3 Scalar 25 I(4) 4 scalar COM JBPAD4 Scalar 25 I(4) 4 scalar COM JBPAD5 Scalar 25 I(4) 4 scalar COM JD Local 297 I(4) 4 scalar 297,1144 JDIM1 Param 20 I(4) 4 scalar 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,20,21,22,23,26,189,209,2 10,211,213,214,215,216 JDIM2 Param 20 I(4) 4 scalar 3,4,5,6,7,8,9,10,11,12,13,14,15,16 ,17,18,19,20,21,22,23,26,189,209,2 10,211,213,214,215,216 JE_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 Page 32 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References JM Param 3 I(4) 4 scalar 9,17,20,43,44,45,46,47,51,52,57,12 7 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 JRA Scalar 3 I(4) 4 1 760 COM JS_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM JTAIL Param 17 I(4) 4 scalar JTB Param 1 I(4) 4 scalar 10,25 JTBQ Param 1 I(4) 4 scalar 25 JTPAD1 Scalar 26 I(4) 4 scalar COM JTPAD2 Scalar 26 I(4) 4 scalar COM JTPAD3 Scalar 26 I(4) 4 scalar COM JTPAD4 Scalar 26 I(4) 4 scalar COM JTPAD5 Scalar 26 I(4) 4 scalar COM K15 Param 141 I(4) 4 scalar 144 K400 Scalar 188 I(4) 4 scalar COM KBT1 Local 827 I(4) 4 scalar 827,881,887,900,901 KBT2 Local 828 I(4) 4 scalar 828,853,858,859,863,881,887,893 KBTM Local 226 I(4) 4 2 2070 937,951,986,1001,1009,1049,1052,11 41 KCLD Local 227 I(4) 4 1 45 940,983,984,985,986,989,990,995,99 7,998,999,1000,1001,1005,1006,1007 ,1008,1009,1011,1019,1020,1021,102 2,1032 KNTLYR Local 415 I(4) 4 scalar 415,418,419 KS Local 1235 I(4) 4 scalar 1235,1286 KSMUD Param 119 I(4) 4 scalar 1234,1235,1285,1286 KTH1 Local 829 I(4) 4 scalar 829,882,888,900 KTH2 Local 830 I(4) 4 scalar 830,849,850,854,858,861,882,888,89 4 KTM Scalar 5 I(4) 4 scalar COM KTOP Local 225 I(4) 4 2 2070 936,950,1000,1008,1020,1048,1052,1 063,1141 KTOP1 Local 858 I(4) 4 scalar 858,859 L Local 363 I(4) 4 scalar 363,366,367,368,369,370,371,372,37 3,374,375,376,399,400,401,402,403, 404,405,406,407,417,419,420,421,42 2,470,471,492,555,556,573,574,594, 595,606,607,622,623,727,728,746,74 7,787,788,810,811,812,813,814,833, 835,836,837,838,853,873,874,947,94 8,949,950,951,952,955,956,973,974, 1119,1124,1178,1179,1180,1181,1229 ,1239,1246,1253,1267,1280,1290,129 7,1304,1318 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM L400 Local 621 I(4) 4 scalar 621,625,632,741,749,756 LBASE Local 572 I(4) 4 scalar 572,578,586,587 LC Scalar 5 I(4) 4 2 17955 COM LIST Scalar 5 I(4) 4 scalar COM LL Local 471 I(4) 4 scalar 471,472,473,474,475,486,496,498,50 3,507,508,511,529,533,535,556,557, 559,574,575,576,578,595,596,598,60 7,608,610,623,624,625,632,633,634, Page 33 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References 636,637,642,643,728,729,730,732,73 4,747,748,749,756,757,758,759,760, 765,766,832,833,974,975,976,977,97 8,984,986,989,995,998,1000,1001,10 06,1008,1009,1011,1020,1051,1052,1 053,1054,1056,1057,1058,1179,1180, 1181,1186 LLBOT Local 804 I(4) 4 scalar 804,805,807,810,827,828,832,833,89 3 LLTOP Local 796 I(4) 4 scalar 796,801,810,832,833 LM Param 3 I(4) 4 scalar 3,4,5,6,7,9,21,52,127,192,198,200, 201,202,230,231,233,363,572,587,62 1,631,641,741,755,764,787,795,801, 1119,1128,1178,1186,1194,1195,1229 ,1280 LM1 Param 127 I(4) 4 scalar 805,807 LMH Scalar 4 I(4) 4 2 17955 COM 393,467,549,724,743,928 LML Local 393 I(4) 4 scalar 393,399,467,470,549,555,556,573,57 4,594,595,606,607,622,623,724,727, 743,746,747,928,973,974 LMV Scalar 4 I(4) 4 2 17955 COM LONG Local 147 L(4) 4 scalar 278,280,284,1181,1193,1201,1278 LOOPS Common 2 158840 LP1 Param 127 I(4) 4 scalar 4,9,197,220,221,224,225,226,227,93 6,937,947,1045,1121,1122,1123 LSM Param 3 I(4) 4 scalar LTOP Scalar 188 I(4) 4 1 3 COM 796,804 LTROP Local 631 I(4) 4 scalar 631,637,641,642,755,760,764,765 LVL Scalar 189 I(4) 4 2 17955 COM 394,468,550,707,708,725,744,796,80 4,929,1179 LVLIJ Local 394 I(4) 4 scalar 394,400,401,402,403,404,405,406,40 7,414,417,426,468,471,550,556,574, 595,607,623,725,728,744,747,929,97 4 MALVL Local 795 I(4) 4 scalar 795,813 MAPPINGS Common 56 5024 MASKS Common 2 6894720 MAX Func 1049 scalar 1049,1101 MBOT Local 196 I(4) 4 2 135 384,901 MIN Func 271 scalar 271,805,900,1048,1096 MOD Func 279 scalar 279,280,289,299 MPI_COMM_COMP Scalar 30 I(4) 4 scalar COM MPI_COMM_INTER Scalar 31 I(4) 4 scalar COM MPI_COMM_INTER_ARRAY Scalar 32 I(4) 4 1 100 COM MPPCOM Common 1 1464 MTOP Local 203 I(4) 4 2 135 383,900 MXSNAL Scalar 2 R(4) 4 2 17955 COM MYIE Scalar 4 I(4) 4 scalar COM 256,307,318,328,364,381,391,437,46 6,547,670,714,722,742,783,907,927, 1120,1147,1177,1192,1238,1244,1251 ,1258,1266,1289,1295,1302,1309,131 7 MYIE1 Scalar 4 I(4) 4 scalar COM MYIE1_P1 Scalar 8 I(4) 4 scalar COM MYIE1_P2 Scalar 8 I(4) 4 scalar COM Page 34 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References MYIE1_P3 Scalar 8 I(4) 4 scalar COM MYIE1_P4 Scalar 8 I(4) 4 scalar COM MYIE2 Scalar 4 I(4) 4 scalar COM MYIE2_P1 Scalar 9 I(4) 4 scalar COM MYIE_P1 Scalar 7 I(4) 4 scalar COM MYIE_P2 Scalar 7 I(4) 4 scalar COM MYIE_P3 Scalar 7 I(4) 4 scalar COM MYIE_P4 Scalar 7 I(4) 4 scalar COM MYIE_P5 Scalar 7 I(4) 4 scalar COM MYIS Scalar 4 I(4) 4 scalar COM 256,307,318,328,364,381,391,437,46 6,547,670,714,722,742,783,907,927, 1120,1147,1177,1192,1238,1244,1251 ,1258,1266,1289,1295,1302,1309,131 7 MYIS1 Scalar 4 I(4) 4 scalar COM MYIS1_P1 Scalar 6 I(4) 4 scalar COM MYIS1_P2 Scalar 6 I(4) 4 scalar COM MYIS1_P3 Scalar 6 I(4) 4 scalar COM MYIS1_P4 Scalar 6 I(4) 4 scalar COM MYIS2 Scalar 4 I(4) 4 scalar COM MYIS_P1 Scalar 5 I(4) 4 scalar COM MYIS_P2 Scalar 5 I(4) 4 scalar COM MYIS_P3 Scalar 5 I(4) 4 scalar COM MYIS_P4 Scalar 5 I(4) 4 scalar COM MYIS_P5 Scalar 5 I(4) 4 scalar COM MYJE Scalar 15 I(4) 4 scalar COM 255,306,317,327,360,1237,1243,1250 ,1265,1288,1301,1316 MYJE1 Scalar 15 I(4) 4 scalar COM 1294 MYJE1_P1 Scalar 17 I(4) 4 scalar COM MYJE1_P2 Scalar 17 I(4) 4 scalar COM MYJE1_P3 Scalar 17 I(4) 4 scalar COM MYJE1_P4 Scalar 17 I(4) 4 scalar COM MYJE2 Scalar 15 I(4) 4 scalar COM 1257,1308 MYJE2_P1 Scalar 18 I(4) 4 scalar COM MYJE2_P2 Scalar 18 I(4) 4 scalar COM MYJE2_P3 Scalar 18 I(4) 4 scalar COM MYJE2_P4 Scalar 18 I(4) 4 scalar COM MYJE3 Scalar 15 I(4) 4 scalar COM MYJE3_P4 Scalar 18 I(4) 4 scalar COM MYJE4 Scalar 15 I(4) 4 scalar COM MYJE4_P1 Scalar 19 I(4) 4 scalar COM MYJE4_P4 Scalar 19 I(4) 4 scalar COM MYJE5 Scalar 15 I(4) 4 scalar COM MYJE5_P1 Scalar 19 I(4) 4 scalar COM MYJE5_P2 Scalar 19 I(4) 4 scalar COM MYJE_P1 Scalar 16 I(4) 4 scalar COM MYJE_P2 Scalar 16 I(4) 4 scalar COM MYJE_P3 Scalar 16 I(4) 4 scalar COM MYJE_P4 Scalar 16 I(4) 4 scalar COM MYJE_P5 Scalar 16 I(4) 4 scalar COM MYJS Scalar 10 I(4) 4 scalar COM 255,306,317,327,360,1237,1243,1265 ,1288,1294,1316 MYJS1 Scalar 10 I(4) 4 scalar COM 1250,1301 MYJS1_P1 Scalar 12 I(4) 4 scalar COM MYJS1_P2 Scalar 12 I(4) 4 scalar COM Page 35 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References MYJS1_P3 Scalar 12 I(4) 4 scalar COM MYJS1_P4 Scalar 12 I(4) 4 scalar COM MYJS2 Scalar 10 I(4) 4 scalar COM 1257,1308 MYJS2_P1 Scalar 13 I(4) 4 scalar COM MYJS2_P2 Scalar 13 I(4) 4 scalar COM MYJS2_P3 Scalar 13 I(4) 4 scalar COM MYJS2_P4 Scalar 13 I(4) 4 scalar COM MYJS3 Scalar 10 I(4) 4 scalar COM MYJS3_P4 Scalar 13 I(4) 4 scalar COM MYJS4 Scalar 10 I(4) 4 scalar COM MYJS4_P1 Scalar 14 I(4) 4 scalar COM MYJS4_P4 Scalar 14 I(4) 4 scalar COM MYJS5 Scalar 10 I(4) 4 scalar COM MYJS5_P1 Scalar 14 I(4) 4 scalar COM MYJS5_P2 Scalar 14 I(4) 4 scalar COM MYJS_P1 Scalar 11 I(4) 4 scalar COM MYJS_P2 Scalar 11 I(4) 4 scalar COM MYJS_P3 Scalar 11 I(4) 4 scalar COM MYJS_P4 Scalar 11 I(4) 4 scalar COM MYJS_P5 Scalar 11 I(4) 4 scalar COM MYPE Scalar 2 I(4) 4 scalar COM 441 MY_E Scalar 20 I(4) 4 scalar COM MY_IE_GLB Scalar 2 I(4) 4 scalar COM MY_IE_LOC Scalar 3 I(4) 4 scalar COM MY_IS_GLB Scalar 2 I(4) 4 scalar COM MY_IS_LOC Scalar 3 I(4) 4 scalar COM MY_JE_GLB Scalar 2 I(4) 4 scalar COM MY_JE_LOC Scalar 3 I(4) 4 scalar COM MY_JS_GLB Scalar 2 I(4) 4 scalar COM MY_JS_LOC Scalar 3 I(4) 4 scalar COM MY_N Scalar 20 I(4) 4 scalar COM MY_NE Scalar 21 I(4) 4 scalar COM MY_NEB Scalar 21 I(4) 4 1 8 COM MY_NW Scalar 21 I(4) 4 scalar COM MY_S Scalar 20 I(4) 4 scalar COM MY_SE Scalar 21 I(4) 4 scalar COM MY_SW Scalar 21 I(4) 4 scalar COM MY_W Scalar 20 I(4) 4 scalar COM N Local 380 I(4) 4 scalar 380,382,383,384 NB Param 138 I(4) 4 scalar 186,226,227,942,954,1087 NBAND Local 942 I(4) 4 scalar 942,943,944,954,955,956,1087,1089, 1101,1102 NBC Scalar 5 I(4) 4 scalar COM NBOCO Scalar 6 I(4) 4 scalar COM NBTM Local 1046 I(4) 4 scalar 1046,1049,1051 NC Local 681 I(4) 4 scalar 681,682,1041,1047,1048,1049,1052,1 062,1063,1084,1101,1102 NCFRCV Scalar 4 I(4) 4 2 17955 COM 1165 NCFRST Scalar 5 I(4) 4 2 17955 COM 1162,1169 NCLD Local 1033 I(4) 4 scalar 1033,1037,1041 NCLDS Local 224 I(4) 4 1 45 1032,1033,1141 NCLOD Scalar 3 I(4) 4 scalar COM NCNVC Scalar 7 I(4) 4 scalar COM NCP Scalar 6 I(4) 4 scalar COM NEST Scalar 7 I(4) 4 scalar COM Page 36 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References NFCST Scalar 5 I(4) 4 scalar COM NFILE Local 239 I(4) 4 scalar 239 NINT Func 270 scalar 270 NKTP Local 1045 I(4) 4 scalar 1045,1048,1051 NLVL Local 793 I(4) 4 scalar 793,796,803,804,836,837,860,874,89 9,900,901 NMOD Local 679 I(4) 4 scalar 679,682,688,690,693,694,695,696 NPES Scalar 2 I(4) 4 scalar COM NPHS Scalar 7 I(4) 4 scalar COM 313 NPREC Scalar 6 I(4) 4 scalar COM NRADL Scalar 7 I(4) 4 scalar COM 271,280 NRADPP Local 271 I(4) 4 scalar 271,272,289 NRADS Scalar 7 I(4) 4 scalar COM 271,279,313 NRDLW Scalar 3 I(4) 4 scalar COM NRDSW Scalar 3 I(4) 4 scalar COM NROOT Param 1 I(4) 4 scalar NSHDE Scalar 6 I(4) 4 scalar COM NSOIL Param 1 I(4) 4 scalar 8,9,10 NSTART Scalar 6 I(4) 4 scalar COM NTDDMP Scalar 6 I(4) 4 scalar COM NTSD Scalar 6 I(4) 4 scalar COM 279,280,289,295,314 NTSPH Local 270 I(4) 4 scalar 270,272 NTSTM Scalar 6 I(4) 4 scalar COM OZN Local 202 R(4) 4 2 2025 374,1128,1140 OZON2D Subr 1128 1128 P1 Local 695 R(4) 4 scalar 695,697 P2 Local 696 R(4) 4 scalar 696,697 P400 Local 236 R(4) 4 scalar 236 PBOT Param 134 R(4) 4 scalar 558,597,609 PBOTH Local 215 R(4) 4 2 17955 PBOTL Local 213 R(4) 4 2 17955 PBOTM Local 214 R(4) 4 2 17955 PCTSNO Scalar 5 R(4) 4 2 17955 COM PD Scalar 3 R(4) 4 2 17955 COM 257,438 PDSL Local 209 R(4) 4 2 17955 257,400,401,1148 PDSLIJ Local 1148 R(4) 4 scalar 1148 PHALF Local 220 R(4) 4 1 46 PHITP Local 236 R(4) 4 scalar 236,251 PHYS Common 4 1788916 PI2 Param 135 R(4) 4 scalar 300 PINT Local 220 R(4) 4 2 2070 401,420,427,429,836,837,874,1053,1 054,1063,1121,1122,1123,1140 PL Scalar 7 R(4) 4 scalar COM PLBTM Local 237 R(4) 4 scalar 237,248 PLOMD Local 236 R(4) 4 scalar 236,249 PLQ Scalar 8 R(4) 4 scalar COM PMDHI Local 236 R(4) 4 scalar 236,250 PMID Local 200 R(4) 4 2 2025 400,402,419,420,421,429,486,557,57 6,596,608,624,634,748,758,859,1124 ,1140 PMOD Local 678 R(4) 4 scalar 678,682,697,1149,1160 POZN Local 202 R(4) 4 2 2025 1124,1128 PP Local 486 R(4) 4 scalar 486,487 PPT Local 218 R(4) 4 1 9 241,682,695,696,1160 PQ0 Param 124 R(4) 4 scalar 487 Page 37 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References PREC Scalar 12 R(4) 4 2 17955 COM PROD Local 144 R(8) 8 scalar 1089,1090 PRS1 Local 1053 R(4) 4 scalar 1053,1055,1057 PRS2 Local 1054 R(4) 4 scalar 1054,1055,1057 PSFC Local 193 R(4) 4 1 45 438,439,557,596,608,1140 PSHLTR Scalar 19 R(4) 4 2 17955 COM PT Scalar 5 R(4) 4 scalar COM 400,401,419,427,438 PTBL Scalar 25 R(4) 4 2 10184 COM PTDYN Scalar 3 R(4) 4 scalar COM PTOPC Scalar 188 R(4) 4 1 4 COM 248,249,250,251,836,837,874,1063 PTOPH Local 215 R(4) 4 2 17955 PTOPL Local 213 R(4) 4 2 17955 PTOPM Local 214 R(4) 4 2 17955 PVRBLS Common 2 5530140 PWTS Scalar 186 R(4) 4 1 12 COM Q Scalar 6 R(4) 4 3 807975 COM 405 Q10 Scalar 16 R(4) 4 2 17955 COM Q2 Scalar 21 R(4) 4 3 807975 COM Q30 Scalar 14 R(4) 4 2 17955 COM QC Local 516 R(4) 4 scalar 516,520,523 QI Local 488 R(4) 4 scalar 488,490,494,516 QINT Local 489 R(4) 4 scalar 489,490,516 QKL Local 474 R(4) 4 scalar 474,494,520,523 QMID Local 201 R(4) 4 2 2025 367,405,474,636,759,1057,1140 QS Scalar 6 R(4) 4 2 17955 COM QS0 Scalar 10 R(4) 4 1 134 COM QSHLTR Scalar 18 R(4) 4 2 17955 COM QSUM Local 1044 R(4) 4 scalar 1044,1057,1085,1089 QW Local 487 R(4) 4 scalar 487,488,489 QWBS Scalar 9 R(4) 4 2 17955 COM QZ0 Scalar 5 R(4) 4 2 17955 COM R1 Scalar 188 R(4) 4 scalar COM 1144 RADFS Subr 1138 1138 RADIN Scalar 14 R(4) 4 2 17955 COM RADOT Scalar 14 R(4) 4 2 17955 COM RADTN Subr 2 RANG Local 300 R(4) 4 scalar 300,301,302,303 RCOS1 Local 302 R(4) 4 scalar 302,1128 RCOS2 Local 303 R(4) 4 scalar 303,1128 RD1TIM Common 188 71880 RDETA Scalar 4 R(4) 4 1 45 COM RDP Scalar 7 R(4) 4 scalar COM RDPQ Scalar 8 R(4) 4 scalar COM RDQ Scalar 7 R(4) 4 scalar COM RDTH Scalar 7 R(4) 4 scalar COM RDTHE Scalar 7 R(4) 4 scalar COM RDTHEQ Scalar 8 R(4) 4 scalar COM RDYN Scalar 3 R(4) 4 scalar COM REAL Func 419 scalar 419 REFLO3 Scalar 186 R(4) 4 scalar COM RES Scalar 7 R(4) 4 2 17955 COM 257 RESTRT Scalar 147 L(4) 4 scalar COM 279,280 RF Scalar 8 R(4) 4 2 17955 COM RLAG Param 135 R(4) 4 scalar 300 RLWIN Scalar 4 R(4) 4 2 17955 COM 1202 Page 38 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References RLWOUT Scalar 4 R(4) 4 2 17955 COM 1203 RLWTOA Scalar 5 R(4) 4 2 17955 COM 1204 RLWTT Scalar 7 R(4) 4 3 807975 COM 1181,1290,1318 RNRM Param 131 R(4) 4 scalar RNRM2 Param 132 R(4) 4 scalar ROI Scalar 6 R(4) 4 scalar COM ROS Scalar 6 R(4) 4 scalar COM RQKL Local 521 R(4) 4 scalar 521,523,528,529,531,533 RRAYAV Scalar 186 R(4) 4 scalar COM RRCL Local 226 R(4) 4 3 24840 943,955,1101,1141 RSIN1 Local 301 R(4) 4 scalar 301,1128 RSWIN Scalar 4 R(4) 4 2 17955 COM 1207 RSWOUT Scalar 4 R(4) 4 2 17955 COM 1208 RSWTOA Scalar 5 R(4) 4 2 17955 COM 1209 RSWTT Scalar 7 R(4) 4 3 807975 COM 1180,1239,1267 RTD Param 118 R(4) 4 scalar 453 RTDPTH Scalar 10 R(4) 4 1 4 COM RUN Scalar 147 L(4) 4 scalar COM SELECTED_REAL_KIND Func 141 scalar 141 SFCEXC Scalar 3 R(4) 4 2 17955 COM SH2O Scalar 10 R(4) 4 3 71820 COM SHORT Local 147 L(4) 4 scalar 277,279,283,304,1180,1206,1228 SI Scalar 10 R(4) 4 2 17955 COM SICE Scalar 4 R(4) 4 2 17955 COM SIGMA Scalar 147 L(4) 4 scalar COM SIGT4 Scalar 23 R(4) 4 2 17955 COM 1194 SIN Func 301 scalar 301 SINGLRST Local 2 L(4) 4 scalar SLDPTH Scalar 10 R(4) 4 1 4 COM SLMSK Local 197 R(4) 4 1 45 444,1140 SLPM Param 130 R(4) 4 scalar SM Scalar 4 R(4) 4 2 17955 COM 440,442,444,554,568 SMC Scalar 8 R(4) 4 3 71820 COM SMSTAV Scalar 4 R(4) 4 2 17955 COM SMSTOT Scalar 4 R(4) 4 2 17955 COM SNO Scalar 10 R(4) 4 2 17955 COM SOIL Common 2 1651892 SOILTB Scalar 3 R(4) 4 2 17955 COM SQRT Func 1058 scalar 1058,1096,1099 SQS Scalar 10 R(4) 4 1 134 COM SR Scalar 5 R(4) 4 2 17955 COM SST Scalar 20 R(4) 4 2 17955 COM STBOL Param 135 R(4) 4 scalar 1194 STC Scalar 9 R(4) 4 3 71820 COM STHE Scalar 10 R(4) 4 1 76 COM STHEQ Scalar 11 R(4) 4 1 152 COM SUBPOST Scalar 2 L(4) 4 scalar COM SV Local 204 R(4) 4 1 45 673,675,700 SWRSAV Common 186 112 T Scalar 5 R(4) 4 3 807975 COM 366,403,404 TAUC Local 221 R(4) 4 1 45 1043,1065,1068,1075,1078,1084,1090 ,1096 TAUCV Scalar 188 R(4) 4 1 3 COM TAUDAR Local 222 R(4) 4 1 45 908,1142 TCLD Local 1056 R(4) 4 scalar 1056,1064,1066,1068,1074,1076 Page 39 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References TCLOD Scalar 3 R(4) 4 scalar COM TDTD Scalar 5 R(4) 4 scalar COM TDTQ2 Scalar 5 R(4) 4 scalar COM TEMP1 Scalar 43 R(4) 4 2 91719 COM TEMP10 Scalar 45 R(4) 4 2 91719 COM TEMP11 Scalar 45 R(4) 4 2 91719 COM TEMP12 Scalar 45 R(4) 4 2 91719 COM TEMP13 Scalar 46 R(4) 4 2 91719 COM TEMP14 Scalar 46 R(4) 4 2 91719 COM TEMP15 Scalar 46 R(4) 4 2 91719 COM TEMP16 Scalar 46 R(4) 4 2 91719 COM TEMP2 Scalar 43 R(4) 4 2 91719 COM TEMP2X Scalar 51 R(4) 4 2 92971 COM TEMP3 Scalar 43 R(4) 4 2 91719 COM TEMP4 Scalar 43 R(4) 4 2 91719 COM TEMP5 Scalar 44 R(4) 4 2 91719 COM TEMP6 Scalar 44 R(4) 4 2 91719 COM TEMP7 Scalar 44 R(4) 4 2 91719 COM TEMP8 Scalar 44 R(4) 4 2 91719 COM TEMP9 Scalar 45 R(4) 4 2 91719 COM TEMPCOM Common 42 6603768 TENDK Local 192 R(4) 4 1 45 TENDL Local 198 R(4) 4 2 2025 376,1144,1181 TENDS Local 198 R(4) 4 2 2025 375,1144,1180 TG Scalar 19 R(4) 4 2 17955 COM TH10 Scalar 16 R(4) 4 2 17955 COM TH30 Scalar 14 R(4) 4 2 17955 COM THE0 Scalar 10 R(4) 4 1 76 COM THE0Q Scalar 11 R(4) 4 1 152 COM THL Scalar 7 R(4) 4 scalar COM THMID Local 201 R(4) 4 2 2025 404,422,575,633,757 THS Scalar 6 R(4) 4 2 17955 COM 440,442 THZ0 Scalar 5 R(4) 4 2 17955 COM TIME Local 295 R(4) 4 scalar 295,296 TIMES Local 314 R(4) 4 scalar 314,315 TKL Local 473 R(4) 4 scalar 473,476,487 TL Local 211 R(4) 4 2 17955 1230,1239,1245,1252,1259,1261,1267 ,1281,1290,1296,1303,1310,1312,131 8 TMID Local 200 R(4) 4 2 2025 366,403,422,473,1056,1058,1140,119 4,1195 TMT0 Local 476 R(4) 4 scalar 476,477,481,488,490,493,502,506 TMT15 Local 477 R(4) 4 scalar 477,489 TOPO Common 50 17478548 TOT Local 216 R(4) 4 2 17955 309,321,329 TRDLW Scalar 3 R(4) 4 scalar COM TRDSW Scalar 3 R(4) 4 scalar COM TRESH Param 131 R(4) 4 scalar 131 TRESH2 Param 132 R(4) 4 scalar 132 TSHLTR Scalar 18 R(4) 4 2 17955 COM TSKN Local 193 R(4) 4 1 45 440,1140 TTBL Scalar 25 R(4) 4 2 10184 COM TTBLQ Scalar 25 R(4) 4 2 66880 COM TTCL Local 227 R(4) 4 3 24840 944,956,1102,1141 TTVG Scalar 51 R(4) 4 2 92971 COM Page 40 Source Listing RADTN 2014-12-17 20:47 Symbol Table RADTN.F Name Object Declared Type Bytes Dimen Elements Attributes References TWBS Scalar 9 R(4) 4 2 17955 COM U Scalar 5 R(4) 4 3 807975 COM U00 Scalar 3 R(4) 4 2 17955 COM 397 U00KL Local 492 R(4) 4 scalar 492,494 U10 Scalar 17 R(4) 4 2 17955 COM U30 Scalar 15 R(4) 4 2 17955 COM UL Scalar 4 R(4) 4 1 90 COM 492 UMID Local 232 R(4) 4 1 45 397,492 US Local 462 R(4) 4 scalar 462,529,531,535 USTAR Scalar 3 R(4) 4 2 17955 COM UTIM Local 243 R(4) 4 scalar 243,492 UZ0 Scalar 4 R(4) 4 2 17955 COM V Scalar 6 R(4) 4 3 807975 COM V10 Scalar 17 R(4) 4 2 17955 COM V30 Scalar 15 R(4) 4 2 17955 COM VBM2 Scalar 3 R(4) 4 2 17955 COM VBM3 Scalar 3 R(4) 4 2 17955 COM VEGFRC Scalar 7 R(4) 4 2 17955 COM VRBLS Common 2 13143060 VTM Scalar 7 R(4) 4 3 807975 COM VZ0 Scalar 4 R(4) 4 2 17955 COM WA Param 119 R(4) 4 scalar 119 WG Param 119 R(4) 4 scalar WMID Local 231 R(4) 4 2 2025 369,406,475 WPDARX Scalar 8 R(4) 4 2 17955 COM XLAT Local 195 R(4) 4 1 45 453,1128,1140 Z0 Scalar 3 R(4) 4 2 17955 COM ZENITH Subr 296 296,315 ZERO2 Subr 1230 1230,1231,1232,1281,1282,1283 Page 41 Source Listing RADTN 2014-12-17 20:47 Subprograms/Common Blocks RADTN.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References ACMCLD Common 2 287288 ACMRDL Common 2 3662832 ACMRDS Common 2 3662832 CLDWTR Common 2 3447720 CNVCLD Common 2 287280 CTLBLK Common 4 112 CUINIT Common 182 4 DYNAM Common 2 725404 GLB_TABLE Common 37 128 INDX Common 2 6564 INDXG Common 8 8084 LOOPS Common 2 158840 MAPPINGS Common 56 5024 MASKS Common 2 6894720 MPPCOM Common 1 1464 PHYS Common 4 1788916 PVRBLS Common 2 5530140 RADTN Subr 2 RD1TIM Common 188 71880 SOIL Common 2 1651892 SWRSAV Common 186 112 TEMPCOM Common 42 6603768 TOPO Common 50 17478548 VRBLS Common 2 13143060 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 Page 42 Source Listing RADTN 2014-12-17 20:47 RADTN.F -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__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant -fpe3 -fprm nearest no -ftz -fp_model precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits nono_precise_libm_functions -heap_arrays 0 no -threadprivate_compat -g2 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations Page 43 Source Listing RADTN 2014-12-17 20:47 RADTN.F -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : RADTN.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100