ACALC.f90#52 INTEGER :: I , J AERMOD.f90#615 INTEGER :: I , IDSTAT #673 WRITE (DBGUNT,200) (MODOPS(I),I=1,18) ANYAVG.f90#50 INTEGER I , NLVLS , NDXABV , NDXBLW #84 DO I = NDXABV + 1 , NDXBLW !153850K #85 SUM = SUM + (HTS(I)-HTS(I-1))*0.5*(PARRAY(I)+PARRAY(I-1)) !544507K APARM.f90#34 INTEGER :: I , J , ISDX #155 DO I = 1 , NVERTS(ISDX) #156 AXCNTR(ISDX) = AXCNTR(ISDX) + AXVERT(I,ISDX) ! 36 #157 AYCNTR(ISDX) = AYCNTR(ISDX) + AYVERT(I,ISDX) ARDIST.f90#37 INTEGER :: I , INDX #57 DO I = 1 , NVERT + 1 #58 XSRC = XVERT(I) !6575040 #59 YSRC = YVERT(I) #60 SPA(I,1) = -((XR-XSRC)*WDSIN+(YR-YSRC)*WDCOS) #61 SPA(I,2) = (XR-XSRC)*WDCOS - (YR-YSRC)*WDSIN #62 XMINREC = MIN(XMINREC,SPA(I,1)) #63 XMAXREC = MAX(XMAXREC,SPA(I,1)) #64 YMINREC = MIN(YMINREC,SPA(I,2)) #65 YMAXREC = MAX(YMAXREC,SPA(I,2)) AVETIM.f90#29 INTEGER :: I , J , K #41 DO I = 3 , IFC #42 KOPT = FIELD(I) ! 9 #65 DO I = 3 , IFC #66 KOPT = FIELD(I) ! 9 CAV_SRC.f90#96 DO I = 2 , 3 #97 QC(I) = 0.0 ! 99136 #98 YRC(I) = YR #99 ZRC(I) = ZR #100 SZC(I) = SZCAV0 #101 SYC(I) = SYCAV0 #145 DO I = 2 , 3 #146 QC(I) = 0.0 ! 642172 #147 YRC(I) = YR #148 ZRC(I) = ZR #149 SZC(I) = SZCAV0 #150 SYC(I) = SYCAV0 #186 DO I = N1 , N2 ! 3572 #187 ZRC(I) = 0.0 ! 3834 #188 SZC(I) = SZCAV0 #189 SYC(I) = SYCAV0 CHK_ENDYR.f90#26 INTEGER :: IEND_DAY , I , J , K , L , M #90 DO I = 1 , NHIVAL ! 0 #91 HIVALU(J,I,K,L,M) = 0.0 ! 0 #92 NHIDAT(J,I,K,L,M) = 0 #93 HCLMSG(J,I,K,L,M) = ' ' DAYRNG.f90#31 INTEGER :: I , K , IMN , IDY , IMN1 , IDY1 , IMN2 , IDY2 , JDAYB ,& #45 DO I = 3 , IFC ! 0 #54 CALL FSPLIT(PATH,KEYWRD,FIELD(I),ILEN_FLD,'-',RMARK,BEGRNG, & DEBOPT.f90#44 DO I = 3 , IFC ! 0 #45 KOPT = FIELD(I) ! 0 DEFINE.f90#59 DO I = LOCB(1) , LOCE(2) + 1 #60 IF ( RUNST(I).NE.' ' ) BLINE = .FALSE. ! 9792 #64 DO I = LOCB(1) + 12 , ISTRG ! 816 #66 IF ( .NOT.INFLD .AND. RUNST(I).NE.' ' ) THEN ! 97920 #76 ELSEIF ( INFLD .AND. RUNST(I).EQ.' ' ) THEN #81 LOCE(IFC) = I - 1 #86 IF ( INFLD .AND. I.EQ.ISTRG ) LOCE(IFC) = ISTRG ! 97920 DPUNIT.f90#55 DO I = 1 , NTYP ! 0 #56 EMIFAC(I) = FNUM ! 0 #57 EMILBL(I) = FIELD(4) #58 OUTLBL(I) = FIELD(5) #60 PERLBL(I) = RUNST1(LOCB(5):LOCE(5))//'/YR' ! 0 #62 PERLBL(I) = FIELD(5) ! 0 #66 DO I = 2 , NTYP ! 0 #67 EMIFAC(I) = FNUM ! 0 #68 EMILBL(I) = FIELD(4) #69 OUTLBL(I) = FIELD(5) #71 PERLBL(I) = RUNST1(LOCB(5):LOCE(5))//'/YR' ! 0 #73 PERLBL(I) = FIELD(5) ! 0 DSBLDG.f90#26 INTEGER :: I , IH , IL , ISDX #76 DO I = 1 , NUMSRC ! 0 #78 CALL ASNGRP(SRCID(I),LID,HID,INGRP) ! 0 #79 IF ( INGRP .AND. SRCTYP(I).EQ.'POINT' ) THEN #80 ISDX = I ! 0 EMVARY.f90#32 INTEGER :: I , IH , IL , ISDX , IQMAX #105 DO I = 1 , NUMSRC ! 0 #107 CALL ASNGRP(SRCID(I),LID,HID,INGRP) ! 0 #109 ISDX = I ! 0 ERRHDL.f90#32 I = 1 #54 DO WHILE ( .NOT.FIND .AND. I.LE.IERRN ) ! 1282 #55 IF ( INERCD.EQ.ERRCOD(I) ) THEN ! 189994 #56 ERRMG1 = ERRMSG(I) ! 1282 #59 I = I + 1 ! 189994 EVALCK.f90#28 INTEGER :: I , INDEX EVALFL.f90#31 INTEGER :: I , INDEX #37 DO I = 1 , NUMARC #52 CWIC = SRT2PI*SYMAX(I)*ARCCL(I) ! 0 #55 CWICN = CWIC*UMAX(I)*AMAX1(HEMAX(I),ZI) #61 CWICW = SRT2PI*SYMAX(I)*(CHIDMW(I)+CHINMW(I)) & #62 & + SRT2PI*SY3MX(I)*CHI3MW(I) #63 CWICL = SRT2PI*SYMAX(I)*(CHIDML(I)+CHINML(I)) & #64 & + SRT2PI*SY3MX(I)*CHI3ML(I) #71 CWICW = SRT2PI*SYMAX(I)*UMAX(I)*ZI*(CHIDMW(I)+CHINMW(I)) & #72 & + SRT2PI*SY3MX(I)*U3MAX(I)*HPEN*CHI3MW(I) #73 CWICL = SRT2PI*SYMAX(I)*UMAX(I)*ZI*(CHIDML(I)+CHINML(I)) & #74 & + SRT2PI*SY3MX(I)*U3MAX(I)*HPEN*CHI3ML(I) #82 UOUST = UMAX(I)/USTAR ! 0 #88 IF ( UMAX(I).GE.1.0E-10 ) THEN ! 0 #89 SVOU = SVMAX(I)/UMAX(I) ! 0 #96 HEOZI = HEMAX(I)/ZI ! 0 #123 UOWST = UMAX(I)/WSTAR ! 0 #129 IF ( UMAX(I).GE.1.0E-10 .AND. ZI.GE.1.0E-10 ) THEN ! 0 #130 XNDIM = DXMAX(I)*WSTAR/(UMAX(I)*ZI) ! 0 #143 PWSTAR = SZMAX(I) #163 UOUT = U3MAX(I) ! 0 #164 SYOUT = SY3MX(I) #166 UOUT = UMAX(I) ! 0 #167 SYOUT = SYMAX(I) #179 WRITE (IELUNT(ISRC),9000) SRCID(ISRC) , KURDAT , ARCID(I) , & #180 & ARCMAX(I) , QMAX(I) , CWIC , CWICN , & #181 & DXMAX(I) , UOUT , SVMAX(I) , SWMAX(I)& #182 & , SYOUT , HEMAX(I) , OBUOUT , ZI , & #184 & PPF , CHIDML(I) , CHINML(I) , & #185 & CHI3ML(I) , XNDIM , HEOZI , FSTAR , & #186 & AHS(ISRC) , ARCCL(I) , AFV , & #187 & HSBLMX(I) EVCARD.f90#37 INTEGER :: I , ILSAVE , ITEMPDATE , ITEMPYEAR #54 DO I = 1 , 366 ! 0 #55 IPROC(I) = 0 ! 0 #109 DO I = 1 , NUMEVE ! 0 #110 ITEMPDATE = EVDATE(I) ! 0 EVLINI.f90#34 DO I = 1 , NUMARC #35 ARCMAX(I) = 0.0 ! 0 #40 ARCCL(I) = 0.0 #41 QMAX(I) = 0.0 #42 DXMAX(I) = 0.0 #43 UMAX(I) = 0.0 #44 U3MAX(I) = 0.0 #45 SVMAX(I) = 0.0 #46 SWMAX(I) = 0.0 #47 SYMAX(I) = 0.0 #48 SY3MX(I) = 0.0 #49 HEMAX(I) = 0.0 #50 CHIDMW(I) = 0.0 #51 CHINMW(I) = 0.0 #52 CHI3MW(I) = 0.0 #53 CHIDML(I) = 0.0 #54 CHINML(I) = 0.0 #55 CHI3ML(I) = 0.0 EV_FLUSH.f90#31 DO I = 1 , NUMTYP #32 HRVAL(I) = 0.0 ! 0 EV_METEXT.f90#97 DO I = 1 , NPLVLS #99 IFLAG(I) = AIFLAG(IHOUR,I) ! 0 #100 PFLHT(I) = APFLHT(IHOUR,I) #101 PFLWD(I) = APFLWD(IHOUR,I) #102 PFLWS(I) = APFLWS(IHOUR,I) #103 PFLTA(I) = APFLTA(IHOUR,I) #104 PFLSA(I) = APFLSA(IHOUR,I) #105 PFLSV(I) = APFLSV(IHOUR,I) #106 PFLSW(I) = APFLSW(IHOUR,I) #111 DO I = 1 , NTGLVL #113 PFLTG(I) = APFLTG(IHOUR,I) ! 0 #114 PFLTGZ(I) = APFLTGZ(IHOUR,I) EV_SETUP.f90#38 INTEGER :: I , IFSTAT #74 READ (INUNIT,RDFRM,END=999) RUNST1 , (RUNST(I),I=1,ISTRG) EXKEY.f90#39 DO I = 1 , IKN #41 IF ( INPFLD.EQ.KEYWD(I) ) THEN ! 20598 EXPATH.f90#40 DO I = 1 , IPN #42 IF ( INPFLD.EQ.PATHWY(I) ) THEN ! 1212 FLAGDF.f90#42 DO I = 1 , NREC ! 0 #43 AZFLAG(I) = ZFLG ! 0 FLGHGT.f90#32 INTEGER :: I , J #39 DO I = 1 , IFC #40 IF ( FIELD(I).EQ.'FLAG' ) ISC = I + 1 ! 0 #66 DO I = ISC + 1 , IFC #67 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 FSPLIT.f90#31 INTEGER :: I , LENGTH , IDELM #38 I = LENGTH #47 DO WHILE ( .NOT.MEND .AND. I.GE.1 ) #48 CHK = INPFLD(I:I) ! 18480 #55 ENDFLD = INPFLD(I+1:LENGTH) #56 IF ( I.EQ.1 ) THEN #73 I = I - 1 ! 18480 GASDEP.f90#28 INTEGER :: I , IH , IL , ISDX #126 DO I = 1 , NUMSRC ! 0 #128 CALL ASNGRP(SRCID(I),LID,HID,INGRP) ! 0 #130 ISDX = I ! 0 GDLAND.f90#24 INTEGER I , J , ILAND_NDX #47 DO I = 3 , IFC #49 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 GDSEAS.f90#23 INTEGER I , J , ISEA_NDX #46 DO I = 3 , IFC #48 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 GENCAR.f90#27 INTEGER :: I , J , K #36 DO I = 1 , IFC #37 IF ( FIELD(I).EQ.'XYINC' ) ISC = I + 1 ! 0 #84 DO I = 1 , ICOUNT ! 0 #85 XCOORD(I,INNET) = XINT + XDELTA*FLOAT(I-1) ! 0 GENCIR.f90#27 INTEGER :: I , ISDX , NSIDES , NEWRAD #44 DO I = 1 , NSIDES #45 IF ( I.NE.1 ) ANG = ANG + ANGINC ! 0 #47 AXVERT(I,ISDX) = (NEWRAD*SIN(ANG/RTODEG)) + AXS(ISDX) #48 AYVERT(I,ISDX) = (NEWRAD*COS(ANG/RTODEG)) + AYS(ISDX) GENPOL.f90#26 INTEGER :: I , J , K #35 DO I = 1 , IFC #36 IF ( FIELD(I).EQ.'GDIR' ) ISC = I + 1 ! 18 GETFLD.f90#25 INTEGER :: I , J #36 DO I = 1 , IFC #37 IF ( LOCE(I)-LOCB(I).LE.(ILEN_FLD-1) ) THEN ! 5646 #39 WRITE (FIELD(I),WRTFRM) (RUNST(J),J=LOCB(I),LOCE(I)) ! 5646 #43 WRITE (FIELD(I),WRTFRM) (RUNST(J),J=LOCB(I),LOCB(I)+ILEN_FLD& GRDDEN.f90#31 INTEGER :: I , NBELOW #52 DO I = 1 , MXGLVL #56 TAMB = GRIDPT(I) - GOVRCP*(GRIDHT(I)+ZBASE) ! 264828 #58 GRIDRHO(I) = RAMB0*(TAMB0/TAMB)*EXP(-G*GRIDHT(I)/(RGASM*TBAR)) GRDEPS.f90#36 DO I = 1 , MXGLVL #38 TSUBLR = AR1*ZI/GRIDSW(I) ! 0 #39 GRIDEPS(I) = 0.78*GRIDSW(I)*GRIDSW(I)/TSUBLR GRDURBAN.f90#39 DO I = 1 , MXGLVL #41 GRDSVR(I) = GRIDSV(I) ! 0 #42 GRDSWR(I) = GRIDSW(I) #43 GRDTGR(I) = GRIDTG(I) #44 GRDPTR(I) = GRIDPT(I) #48 IF ( GRIDHT(I).LE.ZIURB ) THEN #51 ELSEIF ( GRIDHT(I).GT.ZIURB .AND. GRIDHT(I).LE.ZDCRS ) THEN #55 CALL GINTRP(ZIURB,SV2,ZDCRS,SV2DCR,GRIDHT(I),VAL2) #65 IF ( GRIDHT(I).LE.0.1*ZIURB ) THEN ! 0 #66 SW2 = 1.6*(GRIDHT(I)/ZIURB)**0.667*URBWSTR**2 ! 0 #69 ELSEIF ( GRIDHT(I).GT.0.1*ZIURB .AND. GRIDHT(I).LE.ZIURB ) THEN #72 ELSEIF ( GRIDHT(I).GT.ZIURB ) THEN #76 SW2 = 0.35*URBWSTR**2*EXP(-(6.*(GRIDHT(I)-ZIURB)/ZIURB)) ! 0 #81 GRDSVU(I) = SQRT(GRIDSV(I)**2+SVURB**2) ! 0 #82 GRDSWU(I) = SQRT(GRIDSW(I)**2+SWURB**2) #84 IF ( GRIDHT(I).LE.ZIURB ) THEN #85 GRDTGU(I) = 1.0E-5 ! 0 #87 GRDTGU(I) = GRIDTG(I) ! 0 HEADER.f90#66 WRITE (IOUNIT,9040) (MODOPS(I),I=1,18) HILHGT.f90#26 INTEGER :: I , J , IZH1 #34 DO I = 1 , IFC #35 IF ( FIELD(I).EQ.'HILL' ) ISC = I + 1 ! 0 #61 DO I = ISC + 1 , IFC #62 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 #81 DO I = IZH1 , IZH ! 0 #82 ZHTMP2(I) = 0.3048*ZHTMP2(I) ! 0 HPSORT.f90#30 INTEGER :: I , J , IDIM , NVAR , ILMID , IR #55 I = ILMID ! 134720 #63 UVAR(I) = UVAR(J) ! 118584 #64 I = J #71 UVAR(I) = RU ! 134720 HQREAD.f90#28 INTEGER :: I , IS , IHR , IHYEAR , IHMON , IHDAY , IHHOUR #47 & (RUNST(I),I=1,ISTRG) HREMIS.f90#25 INTEGER :: I , K , IH , IL #68 DO I = 1 , NUMSRC ! 0 #69 QFLAG(I) = 'HOURLY' ! 0 #73 DO I = 4 , IFC ! 0 #74 CALL FSPLIT(PATH,KEYWRD,FIELD(I),ILEN_FLD,'-',RMARK,LOWID, & HRLOOP.f90#54 INTEGER :: IEND_DAY , I , J , K , L , M , ILSAVE #414 DO I = IPAIR + 1 , NPAIR ! 0 #415 IDCONC(IAVE,I) = IDUM ! 0 #416 TXCONC(IAVE,I) = RDUM #418 WRITE (ITXUNT(IAVE)) (IDCONC(IAVE,I),I=1,NPAIR) ! 0 #419 WRITE (ITXUNT(IAVE)) (TXCONC(IAVE,I),I=1,NPAIR) HRQEXT.f90#28 INTEGER :: I , IS , IHYEAR , IHMON , IHDAY , IHHOUR #42 READ (IHREMI,RDFRM,ERR=99,END=999) RUNST1 , (RUNST(I),I=1,ISTRG) INCLUD.f90#30 INTEGER :: I , ILREAL #88 READ (INCUNT,RDFRM,END=999) RUNST1 , (RUNST(I),I=1,ISTRG) INPPDM.f90#28 INTEGER :: I , J , K , IH , IL , ISDX , IING #81 SOURCE_LOOP:DO I = 1 , NUMSRC ! 0 #83 CALL ASNGRP(SRCID(I),LID,HID,INGRP) ! 0 #85 IING = I ! 0 #98 APDIAM(ISET,I) = FNUM ! 0 INPPDN.f90#28 INTEGER :: I , J , K , IH , IL , ISDX , IING #84 SOURCE_LOOP:DO I = 1 , NUMSRC ! 0 #86 CALL ASNGRP(SRCID(I),LID,HID,INGRP) ! 0 #88 IING = I ! 0 #104 APDENS(ISET,I) = FNUM ! 0 INPPHI.f90#26 INTEGER :: I , J , K , IH , IL , ISDX , IING #82 SOURCE_LOOP:DO I = 1 , NUMSRC ! 0 #84 CALL ASNGRP(SRCID(I),LID,HID,INGRP) ! 0 #86 IING = I ! 0 #102 APHI(ISET,I) = FNUM ! 0 LWRUPR.f90#30 INTEGER :: I , INDCHK #39 DO I = 1 , ISTRG #40 IF ( RUNST(I).NE.' ' ) THEN ! 107712 #41 INDCHK = INDEX(LWCASE,RUNST(I)) ! 26946 #42 IF ( INDCHK.NE.0 ) RUNST(I) = UPCASE(INDCHK:INDCHK) MAIN1.f90#907 DATA (KEYWD(I),I=1,IKN)/'STARTING' , 'FINISHED' , 'TITLEONE' , & MARCHING.f90#23 DO I = 1 , 7 !6441956 #24 FNEW(I) = FOLD(I) + RHS(I)*DS ! 45093K MECARD.f90#31 INTEGER :: I , ND , NDYS #112 DO I = 1 , 366 ! 0 #113 IPROC(I) = 0 ! 0 #225 DO I = 1 , ND #226 IF ( IPROC(I).EQ.1 ) NDYS = NDYS + 1 ! 1095 MEREAD.f90#41 INTEGER :: I , J , IHR , IJDAY , IDATCHK , JUYI , JUSI , JSYI , & #251 DO I = 1 , NTGLVL #252 APFLTG(IHR,I) = PFLTG(I) ! 0 #253 APFLTGZ(IHR,I) = PFLTGZ(I) #310 DO I = 1 , NTGLVL #311 APFLTG(IHR,I) = PFLTG(I) ! 0 #312 APFLTGZ(IHR,I) = PFLTGZ(I) METDAT.f90#36 INTEGER :: I , IHR , ILMAX , IJDAY #79 DO I = 1 , 24 ! 2 #84 IF ( I.EQ.1 ) THEN ! 48 #180 DO I = 1 , NPLVLS #181 WRITE (IOUNIT,9036) KYEAR , KMONTH , KDAY , KHOUR , PFLHT(I) , & #182 & IFLAG(I) , PFLWD(I) , PFLWS(I) , PFLTA(I) ,& #183 & PFLSA(I) , PFLSW(I) , PFLSV(I) METDEB.f90#64 DO I = MXGLVL , 1 , -1 #65 WRITE (DBMUNT,6120) I , GRIDHT(I) , GRIDWD(I) , GRIDWS(I) , & #66 & GRIDSV(I) , GRIDSW(I) , GRIDPT(I) , & #67 & GRIDTG(I) METH_2.f90#25 INTEGER :: I , IH , IL , ISDX #101 DO I = 1 , NUMSRC ! 0 #103 CALL ASNGRP(SRCID(I),LID,HID,INGRP) ! 0 #105 ISDX = I ! 0 METSUM.f90#31 INTEGER :: I , ILMAX #82 DO I = 1 , NPLVLS #83 IF ( PFLTA(I).EQ.-999.0 ) THEN ! 0 #84 PFLTEMP = PFLTA(I) ! 0 #86 PFLTEMP = PFLTA(I) - DCTODK ! 0 #88 WRITE (IPUNIT,99026) KYEAR , KMONTH , KDAY , KHOUR , PFLHT(I) ,& #89 & IFLAG(I) , PFLWD(I) , PFLWS(I) , PFLTEMP ,& #90 & PFLSW(I) , PFLSV(I) MODOPT.f90#68 DO I = 3 , IFC ! 3 #69 KOPT = FIELD(I) ! 6 #85 DO I = 3 , IFC ! 3 #86 KOPT = FIELD(I) ! 6 #93 DO I = 3 , IFC #94 KOPT = FIELD(I) ! 6 NO2RAT.f90#26 INTEGER :: I , IH , IL , ISDX , IING , IPINDEX #84 DO I = 1 , NUMSRC ! 0 #86 CALL ASNGRP(SRCID(I),LID,HID,INGRP) ! 0 #88 IING = I ! 0 NUMGRAD.f90#35 DO I = NTRM1 , 1 , -1 #36 IF ( X.GE.XTR(I) ) THEN ! 14580K #37 IP1 = I + 1 ! 669780 #38 ZEFF = ZTR(IP1) - (ZTR(IP1)-ZTR(I))*(XTR(IP1)-X) & #39 & /(XTR(IP1)-XTR(I)) NUMMET.f90#51 DO I = 1 , NZA #52 UAMB(I) = WS*(ZGPTA(I)/ZANEM)**P ! 0 #59 DO I = 1 , NZA #60 TAMB(I) = TAMB0 + (GRADACT-ADIA)*ZGPTA(I) ! 0 #65 DO I = 2 , NZA #66 DEDZ(I) = ADIA + (TAMB(I)-TAMB(I-1))/(ZGPTA(I)-ZGPTA(I-1)) ! 0 #72 DO I = 1 , NZA #74 TBAR = 0.5*(TAMB(I)+TAMB0) ! 0 #76 RAMB(I) = RAMB0*(TAMB0/TAMB(I)) & #77 & *EXP(-GRAVI*ZGPTA(I)/(RGAS*TBAR)) #102 DO I = 1 , NZA #103 WRITE (IO6,10) I , ZGPTA(I) , UAMB(I) , TAMB(I) ! 0 #110 DO I = 1 , NZA #111 WRITE (IO6,22) I , ZFACEA(I) , DEDZ(I) , ZGPTA(I) , RAMB(I) ! 0 NUMPR1.f90#101 DO I = 2 , 20 #107 DO I = 21 , 26 #113 DO I = 27 , 41 #119 DO I = 42 , 45 #126 DO I = 2 , NZA #127 ZFACEA(I) = 0.5*(ZGPTA(I)+ZGPTA(I-1)) ! 132 NUMRISE.f90#263 DO I = 1 , 7 #264 RHS(I) = 0.5*(RHSTEMP(I)-RHS(I)) ! 22546K #566 DO I = 1 , NTR #567 WRITE (IO6,32) I , XTR(I) , YTR(I) , ZTR(I) , RTR(I) , & #568 & (RTR(I)*0.8) OCALC.f90#43 INTEGER :: I , II , J , ICAT , INOUT , NDXR #211 DO I = 1 , NVERT + 1 ! 0 #212 XVM(I) = AXVERT(I,ISRC) ! 0 #213 YVM(I) = AYVERT(I,ISRC) OLMGRP.f90#26 INTEGER :: I , K , IH , IL #47 DO I = 1 , NUMOLM #48 IF ( TEMPID.EQ.OLMID(I) ) CONT = .TRUE. ! 0 #67 DO I = 1 , NUMSRC ! 0 #68 IGRP_OLM(I,IOLM) = 1 ! 0 #69 L_OLMGRP(I) = .TRUE. #73 DO I = 4 , IFC ! 0 #74 CALL FSPLIT(PATH,KEYWRD,FIELD(I),ILEN_FLD,'-',RMARK,LOWID, & OPARM.f90#30 INTEGER :: I , ISDX #145 DO I = 1 , NVERT #146 AXCNTR(ISDX) = AXCNTR(ISDX) + AXVERT(I,ISDX) ! 0 #147 AYCNTR(ISDX) = AYCNTR(ISDX) + AYVERT(I,ISDX) OUDALY.f90#26 INTEGER :: I , J , IPRDT #46 DO I = 1 , NUMAVE ! 0 #47 IDYTAB(I) = 1 ! 0 #52 DO I = 3 , IFC ! 0 #53 IF ( FIELD(I).EQ.'MONTH' .AND. MONTH ) THEN ! 0 #58 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) OUEVAL.f90#25 INTEGER :: I , J , K , INDSRC #103 DO I = 1 , NUMSRC #104 IF ( I.NE.INDSRC ) THEN ! 0 #105 IF ( EVLFIL(INDSRC).EQ.EVLFIL(I) .AND. IELUNT(INDSRC) & #106 & .EQ.IELUNT(I) ) THEN #108 ELSEIF ( EVLFIL(INDSRC).EQ.EVLFIL(I) .AND. IELUNT(INDSRC) & #109 & .NE.IELUNT(I) ) THEN #113 ELSEIF ( EVLFIL(INDSRC).NE.EVLFIL(I) .AND. IELUNT(INDSRC) & #114 & .EQ.IELUNT(I) ) THEN OUHIGH.f90#26 INTEGER :: I , J , ILOCH(NAVE) , IPRDT , IPRDT1 , IPRDT2 , ISPRD ,& #32 DATA (NCHR1(I),I=1,10)/'FIRST' , 'SECOND' , 'THIRD' , 'FOURTH' , & #35 DATA (NCHR2(I),I=1,10)/'1ST' , '2ND' , '3RD' , '4TH' , '5TH' , & #40 DO I = 1 , NVAL #41 HIGHST(I) = 0 ! 6 #44 DO I = 1 , NAVE ! 3 #45 ILOCH(I) = 0 ! 6 #66 DO I = 1 , NUMAVE ! 3 #67 INHI(I) = 1 ! 6 #68 ILOCH(I) = 1 #75 DO I = 1 , NUMAVE #76 IF ( IPRDT.EQ.KAVE(I) ) THEN ! 0 #78 INHI(I) = 1 #79 ILOCH(I) = 1 #94 DO I = 1 , NUMAVE #95 IF ( IPRDT1.EQ.KAVE(I) ) THEN ! 0 #97 INHI(I) = 1 #98 ILOCH(I) = 1 #119 DO I = 1 , NUMAVE #120 IF ( KAVE(I).GE.IPRDT1 .AND. KAVE(I).LE.IPRDT2 ) THEN ! 0 #122 INHI(I) = 1 #123 ILOCH(I) = 1 #142 DO I = 4 , IFC ! 3 #144 CALL FSPLIT(PATH,KEYWRD,FIELD(I),ILEN_FLD,'-',RMARK,LPRD,HPRD) ! 3 #176 DO I = 1 , NUMAVE ! 3 #178 IF ( HIGHST(J).EQ.1 .AND. ILOCH(I).EQ.1 ) THEN ! 12 #179 NHIAVE(J,I) = 1 ! 12 OUMXFL.f90#30 INTEGER :: I , J , IPRDT , IDAT #161 DO I = 1 , NUMGRP ! 0 #162 IF ( I.NE.INDGRP .OR. J.NE.INDAVE ) THEN ! 0 #163 IF ( THRFIL(INDGRP,INDAVE).EQ.THRFIL(I,J) .AND. & #164 & IMXUNT(INDGRP,INDAVE).EQ.IMXUNT(I,J) ) THEN #166 ELSEIF ( THRFIL(INDGRP,INDAVE).EQ.THRFIL(I,J) .AND. & #167 & IMXUNT(INDGRP,INDAVE).NE.IMXUNT(I,J) ) THEN #171 ELSEIF ( THRFIL(INDGRP,INDAVE).NE.THRFIL(I,J) .AND. & #172 & IMXUNT(INDGRP,INDAVE).EQ.IMXUNT(I,J) ) THEN #222 WRITE (IMXUNT(INDGRP,INDAVE),9007) (MODOPS(I),I=1,18) OUMXVL.f90#26 INTEGER :: I , J , IPRDT #51 DO I = 1 , NUMAVE ! 3 #52 MAXAVE(I) = 1 ! 6 #106 DO I = 1 , NUMAVE ! 3 #107 IMXVAL(I) = INUM ! 6 OUPLOT.f90#27 INTEGER :: I , J , K , IPRDT #32 DATA (NCHR1(I),I=1,10)/'FIRST' , 'SECOND' , 'THIRD' , 'FOURTH' , & #35 DATA (NCHR2(I),I=1,10)/'1ST' , '2ND' , '3RD' , '4TH' , '5TH' , & #119 DO I = 1 , 10 #120 IF ( FIELD(5).EQ.NCHR1(I) .OR. FIELD(5).EQ.NCHR2(I) ) THEN ! 0 #191 DO I = 1 , NHIVAL ! 0 #192 IF ( I.NE.INDVAL .OR. J.NE.INDGRP .OR. K.NE.INDAVE ) THEN! 0 #193 IF ( PLTFIL(INDVAL,INDGRP,INDAVE).EQ.PLTFIL(I,J,K) & #195 & .EQ.IPLUNT(I,J,K) ) THEN #197 ELSEIF ( PLTFIL(INDVAL,INDGRP,INDAVE).EQ.PLTFIL(I,J,K)& #199 & .NE.IPLUNT(I,J,K) ) THEN #203 ELSEIF ( PLTFIL(INDVAL,INDGRP,INDAVE).NE.PLTFIL(I,J,K)& #205 & .EQ.IPLUNT(I,J,K) ) THEN #216 DO I = 1 , NUMGRP ! 0 #217 IF ( PLTFIL(INDVAL,INDGRP,INDAVE).EQ.ANNPLT(I) .AND. & #218 & IPLUNT(INDVAL,INDGRP,INDAVE).EQ.IPPUNT(I) ) THEN #220 ELSEIF ( PLTFIL(INDVAL,INDGRP,INDAVE).EQ.ANNPLT(I) .AND. & #221 & IPLUNT(INDVAL,INDGRP,INDAVE).NE.IPPUNT(I) ) THEN #225 ELSEIF ( PLTFIL(INDVAL,INDGRP,INDAVE).NE.ANNPLT(I) .AND. & #226 & IPLUNT(INDVAL,INDGRP,INDAVE).EQ.IPPUNT(I) ) THEN OUPOST.f90#33 INTEGER :: I , J , IPRDT , IDAT , ISTR , ISTP #188 DO I = 1 , NUMGRP ! 0 #189 IF ( I.NE.INDGRP .OR. J.NE.INDAVE ) THEN ! 0 #190 IF ( PSTFIL(INDGRP,INDAVE).EQ.PSTFIL(I,J) .AND. & #191 & IPSUNT(INDGRP,INDAVE).EQ.IPSUNT(I,J) ) THEN #193 ELSEIF ( PSTFIL(INDGRP,INDAVE).EQ.PSTFIL(I,J) .AND. & #194 & IPSUNT(INDGRP,INDAVE).NE.IPSUNT(I,J) ) THEN #198 ELSEIF ( PSTFIL(INDGRP,INDAVE).NE.PSTFIL(I,J) .AND. & #199 & IPSUNT(INDGRP,INDAVE).EQ.IPSUNT(I,J) ) THEN #209 DO I = 1 , NUMGRP ! 0 #210 IF ( PSTFIL(INDGRP,INDAVE).EQ.ANNPST(I) .AND. & #211 & IPSUNT(INDGRP,INDAVE).EQ.IAPUNT(I) ) THEN #213 ELSEIF ( PSTFIL(INDGRP,INDAVE).EQ.ANNPST(I) .AND. & #214 & IPSUNT(INDGRP,INDAVE).NE.IAPUNT(I) ) THEN #218 ELSEIF ( PSTFIL(INDGRP,INDAVE).NE.ANNPST(I) .AND. & #219 & IPSUNT(INDGRP,INDAVE).EQ.IAPUNT(I) ) THEN #301 WRITE (IPSUNT(INDGRP,INDAVE),9007) (MODOPS(I),I=1,18) OURANK.f90#25 INTEGER :: I , J , K , IPRDT #151 DO I = 1 , NUMAVE #152 IF ( I.NE.INDAVE ) THEN ! 0 #153 IF ( RNKFIL(INDAVE).EQ.RNKFIL(I) .AND. IRKUNT(INDAVE) & #154 & .EQ.IRKUNT(I) ) THEN #156 ELSEIF ( RNKFIL(INDAVE).EQ.RNKFIL(I) .AND. IRKUNT(INDAVE) & #157 & .NE.IRKUNT(I) ) THEN #161 ELSEIF ( RNKFIL(INDAVE).NE.RNKFIL(I) .AND. IRKUNT(INDAVE) & #162 & .EQ.IRKUNT(I) ) THEN OUSEAS.f90#25 INTEGER :: I , J , K #106 DO I = 1 , NUMGRP #107 IF ( I.NE.INDGRP ) THEN ! 0 #108 IF ( SEAHRS(INDGRP).EQ.SEAHRS(I) .AND. ISHUNT(INDGRP) & #109 & .EQ.ISHUNT(I) ) THEN #111 ELSEIF ( SEAHRS(INDGRP).EQ.SEAHRS(I) .AND. ISHUNT(INDGRP) & #112 & .NE.ISHUNT(I) ) THEN #116 ELSEIF ( SEAHRS(INDGRP).NE.SEAHRS(I) .AND. ISHUNT(INDGRP) & #117 & .EQ.ISHUNT(I) ) THEN OUTOXX.f90#25 INTEGER :: I , J , IPRDT , NIDUM , NRDUM , NUMPER #148 DO I = 1 , NUMAVE #149 IF ( I.NE.INDAVE ) THEN ! 0 #150 IF ( TOXFIL(INDAVE).EQ.TOXFIL(I) .AND. ITXUNT(INDAVE) & #151 & .EQ.ITXUNT(I) ) THEN #153 ELSEIF ( TOXFIL(INDAVE).EQ.TOXFIL(I) .AND. ITXUNT(INDAVE) & #154 & .NE.ITXUNT(I) ) THEN #158 ELSEIF ( TOXFIL(INDAVE).NE.TOXFIL(I) .AND. ITXUNT(INDAVE) & #159 & .EQ.ITXUNT(I) ) THEN #178 & ITAB , NXTOX , NYTOX , (IDUM,I=1,NIDUM) #179 WRITE (ITXUNT(INDAVE)) TOXTHR(INDAVE) , (RDUM,I=1,NRDUM) OUTQA.f90#30 INTEGER :: I , IVAL , IDCST1 #44 DO I = 2 , 8 #45 IF ( IOSTAT(I).GT.0 ) OUTOPT = .TRUE. ! 21 PCALC.f90#87 INTEGER :: I , KITER , NDXHE , NDXZMID , NDXZPL , NDXHE3 PDEP.f90#35 DO I = 1 , NPD #36 DQCOR(I) = 1.0 ! 0 #37 WQCOR(I) = 1.0 #40 VSETL = VGRAV(I) ! 0 #41 CALL DEPLETE(VDEP(I),XARG,ROMBERG,DQCOR(I)) #45 IF ( WDPLETE .AND. PSCVRT(I).GT.0.0 ) WQCOR(I) & #46 & = EXP(-PSCVRT(I)*XARG/US) PERPLT.f90#31 INTEGER :: I , J , K #112 DO I = 1 , NUMGRP #113 IF ( I.NE.INDGRP ) THEN ! 0 #114 IF ( ANNPLT(INDGRP).EQ.ANNPLT(I) .AND. IPPUNT(INDGRP) & #115 & .EQ.IPPUNT(I) ) THEN #117 ELSEIF ( ANNPLT(INDGRP).EQ.ANNPLT(I) .AND. IPPUNT(INDGRP) & #118 & .NE.IPPUNT(I) ) THEN #122 ELSEIF ( ANNPLT(INDGRP).NE.ANNPLT(I) .AND. IPPUNT(INDGRP) & #123 & .EQ.IPPUNT(I) ) THEN #134 DO I = 1 , NHIVAL ! 0 #135 IF ( ANNPLT(INDGRP).EQ.PLTFIL(I,J,K) .AND. IPPUNT(INDGRP)& #136 & .EQ.IPLUNT(I,J,K) ) THEN #138 ELSEIF ( ANNPLT(INDGRP).EQ.PLTFIL(I,J,K) .AND. & #139 & IPPUNT(INDGRP).NE.IPLUNT(I,J,K) ) THEN #143 ELSEIF ( ANNPLT(INDGRP).NE.PLTFIL(I,J,K) .AND. & #144 & IPPUNT(INDGRP).EQ.IPLUNT(I,J,K) ) THEN PERPST.f90#28 INTEGER :: I , J #120 DO I = 1 , NUMGRP #121 IF ( I.NE.INDGRP ) THEN ! 0 #122 IF ( ANNPST(INDGRP).EQ.ANNPST(I) .AND. IAPUNT(INDGRP) & #123 & .EQ.IAPUNT(I) ) THEN #125 ELSEIF ( ANNPST(INDGRP).EQ.ANNPST(I) .AND. IAPUNT(INDGRP) & #126 & .NE.IAPUNT(I) ) THEN #130 ELSEIF ( ANNPST(INDGRP).NE.ANNPST(I) .AND. IAPUNT(INDGRP) & #131 & .EQ.IAPUNT(I) ) THEN #141 DO I = 1 , NUMGRP ! 0 #142 IF ( ANNPST(INDGRP).EQ.PSTFIL(I,J) .AND. IAPUNT(INDGRP) & #143 & .EQ.IPSUNT(I,J) ) THEN #145 ELSEIF ( ANNPST(INDGRP).EQ.PSTFIL(I,J) .AND. IAPUNT(INDGRP) & #146 & .NE.IPSUNT(I,J) ) THEN #150 ELSEIF ( ANNPST(INDGRP).NE.PSTFIL(I,J) .AND. IAPUNT(INDGRP) & #151 & .EQ.IPSUNT(I,J) ) THEN PFLINI.f90#26 INTEGER :: I , J #38 DO I = 1 , MXPLVL #39 IFLAG(I) = 0 ! 216000 #40 PFLHT(I) = -99.0 #41 PFLWS(I) = -99.0 #42 PFLWD(I) = -99.0 #43 PFLTA(I) = -99.0 #44 PFLSA(I) = -99.0 #45 PFLSW(I) = -99.0 #46 PFLSV(I) = -99.0 #47 PFLTG(I) = -99.0 #48 PFLTGZ(I) = -99.0 PFRACT.f90#52 DO I = 1 , 100 #55 TWOIZI = 2.*I*HSBL PITEFF.f90#28 INTEGER :: I , II , IUPWND #144 DO I = 1 , 5 ! 0 #145 XVERT(I) = XTEMP(I) ! 0 #146 YVERT(I) = YTEMP(I) PLOTFL.f90#31 INTEGER :: I , IVAL #35 DATA (NCHR2(I),I=1,10)/'1ST' , '2ND' , '3RD' , '4TH' , '5TH' , & #61 WRITE (IPLUNT(IVAL,IGRP,IAVE),9007) (MODOPS(I),I=1,18) PLTANN.f90#60 WRITE (IPPUNT(IGRP),9007) (MODOPS(I),I=1,18) PLUMEF.f90#45 INTEGER :: I , J PLUME_VOL.f90#30 INTEGER :: I , ISDX #248 DO I = 1 , NTAB ! 0 #250 XTAB = XTABLE(I) ! 0 #251 IF ( I.GT.1 ) THEN #252 DELTAX = XTABLE(I) - XTABLE(I-1) ! 0 #294 CALL RELDISP(XTAB,SR(I),SR3(I)) ! 0 #296 IF ( I.GT.1 ) THEN #297 SIGR = MAX(SRMIN,0.5*(SR(I)+SR(I-1))) ! 0 #299 SIGR = MAX(SRMIN,0.5*SR(I)) ! 0 #325 IF ( I.GT.1 ) THEN ! 0 #326 SIGR3 = MAX(SRMIN,0.5*(SR3(I)+SR3(I-1))) ! 0 #328 SIGR3 = MAX(SRMIN,0.5*SR3(I)) ! 0 #375 CALL RELDISP(XARG,SR(I),SR3(I)) ! 0 #377 DELTAX = XARG - XTABLE(I-1) #379 IF ( I.GT.1 ) THEN #380 SIGR = MAX(SRMIN,0.5*(SR(I)+SR(I-1))) ! 0 #382 SIGR = MAX(SRMIN,0.5*SR(I)) ! 0 #407 IF ( I.GT.1 ) THEN ! 0 #408 SIGR3 = MAX(SRMIN,0.5*(SR3(I)+SR3(I-1))) ! 0 #410 SIGR3 = MAX(SRMIN,0.5*SR3(I)) ! 0 PNPOLY.f90#64 INTEGER I , J , N , INOUT #73 DO I = 1 , N #74 XI = X(I) - PX ! 0 #75 YI = Y(I) - PY #82 J = 1 + MOD(I,N) ! 0 POLDST.f90#25 INTEGER :: I , J #31 DO I = 1 , IFC #32 IF ( FIELD(I).EQ.'DIST' ) ISC = I + 1 ! 21 #45 DO I = ISC , IFC #46 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 12 POLORG.f90#25 INTEGER :: I , ISDX #34 DO I = 1 , IFC #35 IF ( FIELD(I).EQ.'ORIG' ) ISC = I + 1 ! 12 PREGENCAR.f90#27 INTEGER :: I , K #36 DO I = 1 , IFC #37 IF ( FIELD(I).EQ.'XYINC' ) ISC = I + 1 ! 0 PREGENPOL.f90#26 INTEGER :: I , K #35 DO I = 1 , IFC #36 IF ( FIELD(I).EQ.'GDIR' ) ISC = I + 1 ! 18 PREINCLUD.f90#29 INTEGER :: I , ILREAL #78 READ (INCUNT,RDFRM,END=999) RUNST1 , (RUNST(I),I=1,ISTRG) PREPOLDST.f90#31 DO I = 1 , IFC #32 IF ( FIELD(I).EQ.'DIST' ) ISC = I + 1 ! 21 #44 DO I = ISC , IFC #45 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 12 PRERADRNG.f90#31 DO I = 1 , IFC #32 IF ( FIELD(I).EQ.'DDIR' ) ISC = I + 1 ! 0 #44 DO I = ISC , IFC #45 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 PRESET.f90#29 INTEGER :: I , J , ISPRD , IEPRD #45 DATA (NCHR1(I),I=1,10)/'FIRST' , 'SECOND' , 'THIRD' , 'FOURTH' , & #48 DATA (NCHR2(I),I=1,10)/'1ST' , '2ND' , '3RD' , '4TH' , '5TH' , & #82 READ (INUNIT,RDFRM,END=999) RUNST1 , (RUNST(I),I=1,ISTRG) #135 DO I = 3 , IFC ! 3 #136 IF ( FIELD(I).EQ.'CONC' .OR. FIELD(I).EQ.'DEPOS' .OR. & #137 & FIELD(I).EQ.'DDEP' .OR. FIELD(I).EQ.'WDEP' ) & #140 IF ( FIELD(I).EQ.'PVMRM' ) THEN #142 ELSEIF ( FIELD(I).EQ.'OLM' ) THEN #149 DO I = 3 , IFC ! 3 #150 IF ( FIELD(I).NE.'PERIOD' .AND. FIELD(I).NE.'ANNUAL' ) & #176 DO I = 4 , IFC ! 3 #178 CALL FSPLIT(PATH,KEYWRD,FIELD(I),ILEN_FLD,'-',RMARK,LPRD,& PREXYPNTS.f90#25 INTEGER :: I , JSET #32 DO I = 1 , IFC ! 0 #33 IF ( FIELD(I).EQ.'XPNTS' ) ISC = I + 1 ! 0 #44 DO I = ISC , IFC #45 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 #55 DO I = 1 , IFC ! 0 #56 IF ( FIELD(I).EQ.'YPNTS' ) ISC = I + 1 ! 0 #67 DO I = ISC , IFC #68 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 PRM_PDEP.f90#36 DO I = 1 , NPD #37 DQCOR(I) = 1.0 ! 0 #38 WQCOR(I) = 1.0 #41 VSETL = VGRAV(I) ! 0 #42 CALL PRM_DEPLETE(VDEP(I),XARG,ROMBERG,DQCOR(I)) #44 IF ( WDPLETE .AND. PSCVRT(I).GT.0.0 ) THEN ! 0 #47 WQCOR(I) = EXP(-PSCVRT(I)*XARG/US) ! 0 PRTANN.f90#30 INTEGER :: I , J , K , II , NX , NY , INDZ , INDC , ISRF , INDEXW #66 DO I = 1 , INNET #68 NPPX = 1 + INT((NUMXPT(I)-1)/NCPP) ! 2 #69 NPPY = 1 + INT((NUMYPT(I)-1)/NRPP) #82 WRITE (IOUNIT,9037) NTID(I) , NTTYP(I) ! 2 #89 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 2 #92 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #93 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #96 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #99 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #102 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #103 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #106 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #112 DO K = 1 + NRPP*(NY-1) , NUMYPT(I) ! 2 #113 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 72 #114 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #115 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #116 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #117 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 72 #118 YCOVAL = YCOORD(K,I) #123 & J=1+NCPP*(NX-1),NUMXPT(I)) #132 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #133 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #134 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #135 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #136 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #137 YCOVAL = YCOORD(K,I) #142 & J=1+NCPP*(NX-1),NUMXPT(I)) PRTDAY.f90#29 INTEGER :: I , J , K , II , NX , NY , INDZ , INDC , ISRF , INDEXW #60 DO I = 1 , INNET #62 NPPX = 1 + INT((NUMXPT(I)-1)/NCPP) ! 0 #63 NPPY = 1 + INT((NUMYPT(I)-1)/NRPP) #71 WRITE (IOUNIT,9037) NTID(I) , NTTYP(I) #78 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #81 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #82 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #85 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #88 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #91 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #92 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #95 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #101 DO K = 1 + NRPP*(NY-1) , NUMYPT(I) ! 0 #102 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #103 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #104 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #105 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #106 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #107 YCOVAL = YCOORD(K,I) #112 & J=1+NCPP*(NX-1),NUMXPT(I)) #121 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #122 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #123 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #124 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #125 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #126 YCOVAL = YCOORD(K,I) #131 & J=1+NCPP*(NX-1),NUMXPT(I)) PRTDET.f90#27 INTEGER :: I , J , N , NSPP , NPAGE , IHR , INGRP #64 & (WORKID(I),I=1+NSPP*(N-1),INGRP) #68 DO I = ISTAHR , IENDHR #69 WRITE (IOUNIT,9062) I , GRPVAL(I) , & #70 & (WHRVAL(I,J),J=1+NSPP*(N-1),INGRP) #72 WRITE (IOUNIT,9064) GRPAVE , (WAVEV(I),I=1+NSPP*(N-1),INGRP)! 0 #76 & (WORKID(I),I=1+NSPP*(N-1),NSPP*N) #80 DO I = ISTAHR , IENDHR #81 WRITE (IOUNIT,9062) I , GRPVAL(I) , & #82 & (WHRVAL(I,J),J=1+NSPP*(N-1),NSPP*N) #85 & (WAVEV(I),I=1+NSPP*(N-1),NSPP*N) PRTMET.f90#29 INTEGER :: I , J , K #38 WRITE (IOUNIT,9037) (IPROC(I),I=1,366) #56 WRITE (IOUNIT,9001) (UCAT(I),I=1,5) PRTOPT.f90#43 INTEGER :: I , ILMAX #192 WRITE (IOUNIT,9042) NUMAVE , (CHRAVE(I),I=1,NUMAVE) ! 3 #201 WRITE (IOUNIT,9042) NUMAVE , (CHRAVE(I),I=1,NUMAVE) ! 0 #209 WRITE (IOUNIT,9042) NUMAVE , (CHRAVE(I),I=1,NUMAVE) ! 0 PRTPM10.f90#27 INTEGER :: I , J , K , II , INDZ , INDC , NX , NY , ISRF , INDEXW #65 DO I = 1 , INNET #67 NPPX = 1 + INT((NUMXPT(I)-1)/NCPP) ! 0 #68 NPPY = 1 + INT((NUMYPT(I)-1)/NRPP) #75 WRITE (IOUNIT,9037) NTID(I) , NTTYP(I) #82 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #85 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #86 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #89 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #92 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #95 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #96 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #99 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #105 DO K = 1 + NRPP*(NY-1) , NUMYPT(I) ! 0 #106 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #107 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #108 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #109 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #110 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #111 YCOVAL = YCOORD(K,I) #116 & ,NUMXPT(I)) #125 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #126 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #127 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #128 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #129 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #130 YCOVAL = YCOORD(K,I) #135 & ,NUMXPT(I)) PRTPM10SUM.f90#25 INTEGER :: I , IVAL , INDMX #30 DATA (RANK(I),I=1,10)/' 1ST' , ' 2ND' , ' 3RD' , ' 4TH' , ' 5TH' ,& PRTREC.f90#34 INTEGER :: I , J , K , INDZ , NX , NY , INDC , ISRF #41 DO I = 1 , INNET #45 WRITE (IOUNIT,9037) NTID(I) , NTTYP(I) #46 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN #51 WRITE (IOUNIT,9036) XORIG(I) , YORIG(I) ! 3 #59 WRITE (IOUNIT,9040) (XCOORD(J,I),J=1,NUMXPT(I)) ! 3 #60 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN #69 WRITE (IOUNIT,9040) (YCOORD(J,I),J=1,NUMYPT(I)) ! 3 #78 NPPX = 1 + INT((NUMXPT(I)-1)/NCPP) #79 NPPY = 1 + INT((NUMYPT(I)-1)/NRPP) #83 WRITE (IOUNIT,9037) NTID(I) , NTTYP(I) #88 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #91 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #92 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #95 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #98 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #101 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #102 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #105 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #110 DO K = 1 + NRPP*(NY-1) , NUMYPT(I) ! 0 #111 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #112 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #113 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #114 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #115 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #116 YCOVAL = YCOORD(K,I) #121 & NUMXPT(I)) #130 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #131 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #132 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #133 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #134 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #135 YCOVAL = YCOORD(K,I) #140 & NUMXPT(I)) #157 NPPX = 1 + INT((NUMXPT(I)-1)/NCPP) #158 NPPY = 1 + INT((NUMYPT(I)-1)/NRPP) #162 WRITE (IOUNIT,9037) NTID(I) , NTTYP(I) #166 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #169 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #170 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #173 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #176 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #179 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #180 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #183 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #188 DO K = 1 + NRPP*(NY-1) , NUMYPT(I) ! 0 #189 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #190 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #191 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #192 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #193 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #194 YCOVAL = YCOORD(K,I) #199 & NUMXPT(I)) #208 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #209 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #210 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #211 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #212 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #213 YCOVAL = YCOORD(K,I) #218 & NUMXPT(I)) #237 NPPX = 1 + INT((NUMXPT(I)-1)/NCPP) #238 NPPY = 1 + INT((NUMYPT(I)-1)/NRPP) #242 WRITE (IOUNIT,9037) NTID(I) , NTTYP(I) #247 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #250 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #251 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #254 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #257 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #260 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #261 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #264 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #269 DO K = 1 + NRPP*(NY-1) , NUMYPT(I) ! 0 #270 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #271 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #272 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #273 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #274 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #275 YCOVAL = YCOORD(K,I) #280 & NUMXPT(I)) #289 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #290 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #291 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #292 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #293 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #294 YCOVAL = YCOORD(K,I) #299 & NUMXPT(I)) #318 DO I = 1 , NUMREC #319 IF ( RECTYP(I).EQ.'DC' ) THEN ! 0 #329 WRITE (BUF132(1:65),9045) AXR(I) , AYR(I) , AZELEV(I) & #330 & , AZHILL(I) , AZFLAG(I) #332 WRITE (BUF132(66:130),9045) AXR(I) , AYR(I) , & #333 & AZELEV(I) , AZHILL(I) , AZFLAG(I) #348 DO I = 1 , NUMREC #349 IF ( RECTYP(I).EQ.'DP' ) THEN ! 0 #351 XRMS = AXR(I) - AXS(IREF(I)) #352 YRMS = AYR(I) - AYS(IREF(I)) #366 WRITE (BUF132(1:65),9047) SRCID(IREF(I)) , RANGE , & #367 & RADIAL , AZELEV(I) , AZHILL(I) , AZFLAG(I) #369 WRITE (BUF132(66:130),9047) SRCID(IREF(I)) , RANGE , & #370 & RADIAL , AZELEV(I) , AZHILL(I) , AZFLAG(I) PRTSOC.f90#31 INTEGER :: I , N , NROW , NPAGE , INGRP , IOGRP #75 WRITE (IOUNIT,9068) EVGRP(IEVENT) , (HEADID(I),I=1,IOGRP) #86 WRITE (IOUNIT,9066) (WORKID(I),WAVEV(I),I=1+120*(N-1),INGRP)! 0 #88 WRITE (IOUNIT,9066) (WORKID(I),WAVEV(I),I=1+120*(N-1),120*N)! 0 PRTSRC.f90#49 INTEGER :: I , J , K , NL , I1 , I2 , I3 , IFR , IDW , ITO , & #72 DO I = 1 , NUMSRC #73 IF ( SRCTYP(I).EQ.'POINT' ) THEN ! 27 #75 IF ( URBSRC(I).EQ.'Y' ) THEN #87 IF ( ADSBH(J,I).NE.0.0 .AND. ADSBW(J,I).NE.0.0 .AND. & #88 & ADSBL(J,I).NE.0.0 ) BLDING = 'YES' #104 IF ( L_METHOD2(I) ) THEN ! 9 #107 WRITE (CNPD,'(I4,1X)') INPD(I) ! 9 #109 WRITE (IOUNIT,9047) SRCID(I) , CNPD , AQS(I) , AXS(I) , & #110 & AYS(I) , AZS(I) , AHS(I) , ATS(I) , & #111 & AVS(I) , ADS(I) , BLDING , URB , & #112 & QFLAG(I) #120 DO I = 1 , NUMSRC #121 IF ( SRCTYP(I).EQ.'VOLUME' ) THEN ! 27 #123 IF ( URBSRC(I).EQ.'Y' ) THEN #140 WRITE (IOUNIT,9075) SRCID(I) , INPD(I) , AQS(I) , AXS(I) , & #141 & AYS(I) , AZS(I) , AHS(I) , ASYINI(I) , & #142 & ASZINI(I) , URB , QFLAG(I) #150 DO I = 1 , NUMSRC #151 IF ( SRCTYP(I).EQ.'AREA' ) THEN ! 27 #153 IF ( URBSRC(I).EQ.'Y' ) THEN #172 WRITE (IOUNIT,9077) SRCID(I) , INPD(I) , AQS(I) , AXS(I) , & #173 & AYS(I) , AZS(I) , AHS(I) , AXINIT(I) , & #174 & AYINIT(I) , AANGLE(I) , ASZINI(I) , & #175 & URB , QFLAG(I) #184 DO I = 1 , NUMSRC #185 IF ( SRCTYP(I).EQ.'AREACIRC' ) THEN ! 27 #187 IF ( URBSRC(I).EQ.'Y' ) THEN #204 WRITE (IOUNIT,9079) SRCID(I) , INPD(I) , AQS(I) , AXS(I) , & #205 & AYS(I) , AZS(I) , AHS(I) , RADIUS(I) , & #206 & NVERTS(I) , ASZINI(I) , URB , QFLAG(I) #214 DO I = 1 , NUMSRC #215 IF ( SRCTYP(I).EQ.'AREAPOLY' ) THEN ! 27 #217 IF ( URBSRC(I).EQ.'Y' ) THEN #235 WRITE (IOUNIT,9081) SRCID(I) , INPD(I) , AQS(I) , AXS(I) , & #236 & AYS(I) , AZS(I) , AHS(I) , NVERTS(I) , & #237 & ASZINI(I) , URB , QFLAG(I) #245 DO I = 1 , NUMSRC #246 IF ( SRCTYP(I).EQ.'OPENPIT' ) THEN ! 27 #248 IF ( URBSRC(I).EQ.'Y' ) THEN #268 WRITE (IOUNIT,9083) SRCID(I) , INPD(I) , AQS(I) , AXS(I) , & #269 & AYS(I) , AZS(I) , AHS(I) , AXINIT(I) , & #270 & AYINIT(I) , AANGLE(I) , AVOLUM(I) , & #271 & URB , QFLAG(I) #298 WRITE (IOUNIT,9068) GRPID(J) , (WORKID(I),I=1,INGRP) ! 3 #300 WRITE (IOUNIT,9068) GRPID(J) , (WORKID(I),I=1,12*K) ! 0 #302 WRITE (IOUNIT,9067) (WORKID(I),I=1+12*(K-1),INGRP) ! 0 #304 WRITE (IOUNIT,9067) (WORKID(I),I=1+12*(K-1),12*K) ! 0 #333 WRITE (IOUNIT,9068) OLMID(J) , (WORKID(I),I=1,INGRP) ! 0 #335 WRITE (IOUNIT,9068) OLMID(J) , (WORKID(I),I=1,12*K) ! 0 #337 WRITE (IOUNIT,9067) (WORKID(I),I=1+12*(K-1),INGRP) ! 0 #339 WRITE (IOUNIT,9067) (WORKID(I),I=1+12*(K-1),12*K) ! 0 #347 DO I = 1 , NUMSRC , 4 #356 IF ( I+3.LE.NUMSRC ) THEN ! 0 #357 WRITE (IOUNIT,9061) SRCID(I) , ANO2_RATIO(I) , SRCID(I+1)& #358 & , ANO2_RATIO(I+1) , SRCID(I+2) , & #359 & ANO2_RATIO(I+2) , SRCID(I+3) , & #360 & ANO2_RATIO(I+3) #361 ELSEIF ( I+2.LE.NUMSRC ) THEN #362 WRITE (IOUNIT,9061) SRCID(I) , ANO2_RATIO(I) , SRCID(I+1)& #363 & , ANO2_RATIO(I+1) , SRCID(I+2) , & #364 & ANO2_RATIO(I+2) #365 ELSEIF ( I+1.LE.NUMSRC ) THEN #366 WRITE (IOUNIT,9061) SRCID(I) , ANO2_RATIO(I) , SRCID(I+1)& #367 & , ANO2_RATIO(I+1) #369 WRITE (IOUNIT,9061) SRCID(I) , ANO2_RATIO(I) ! 0 #376 DO I = 1 , NUMSRC #377 NPD = INPD(I) ! 27 #378 IF ( NPD.NE.0 .AND. .NOT.L_METHOD2(I) ) THEN #384 WRITE (IOUNIT,9050) SRCID(I) , SRCTYP(I) ! 0 #385 WRITE (IOUNIT,9051) (APHI(J,I),J=1,NPD) #387 WRITE (IOUNIT,9052) (APDIAM(J,I),J=1,NPD) #388 WRITE (IOUNIT,9053) (APDENS(J,I),J=1,NPD) #389 ELSEIF ( NPD.NE.0 .AND. L_METHOD2(I) ) THEN #396 WRITE (IOUNIT,9050) SRCID(I) , SRCTYP(I) ! 0 #397 WRITE (IOUNIT,99051) (FINEMASS(I),J=1,NPD) #401 WRITE (IOUNIT,9052) (APDIAM(J,I),J=1,NPD) #402 WRITE (IOUNIT,9053) (APDENS(J,I),J=1,NPD) #411 WRITE (IOUNIT,9050) SRCID(I) , SRCTYP(I) ! 0 #413 WRITE (IOUNIT,99090) PDIFF(I) ! 0 #415 WRITE (IOUNIT,99091) PDIFFW(I) #417 WRITE (IOUNIT,99093) RCLI(I) #419 WRITE (IOUNIT,9094) HENRY(I) #428 DO I = 1 , NUMSRC #437 IF ( ADSBH(J,I).NE.0.0 .AND. ADSBW(J,I).NE.0.0 .AND. & #438 & ADSBL(J,I).NE.0.0 ) BLDING = 'YES' #457 WRITE (IOUNIT,9062) SRCID(I) , & #458 & (J,ABS(ADSBH(J,I)),ADSBW(J,I),ADSBL(J,I)& #459 & ,ADSXADJ(J,I),ADSYADJ(J,I),J=1,NSEC) #471 DO I = 1 , NUMSRC #472 IF ( QFLAG(I).EQ.'SEASON' ) THEN ! 27 #484 WRITE (IOUNIT,9005) SRCID(I) , SRCTYP(I) ! 0 #487 WRITE (IOUNIT,9006) (QFACT(I1,I),I1=1,4) #493 DO I = 1 , NUMSRC #494 IF ( QFLAG(I).EQ.'MONTH' ) THEN ! 27 #510 WRITE (IOUNIT,9009) SRCID(I) , SRCTYP(I) ! 0 #511 WRITE (IOUNIT,9010) (QFACT(I1,I),I1=1,12) #517 DO I = 1 , NUMSRC #518 IF ( QFLAG(I).EQ.'HROFDY' ) THEN ! 27 #529 WRITE (IOUNIT,9009) SRCID(I) , SRCTYP(I) ! 0 #530 WRITE (IOUNIT,9014) (I1,QFACT(I1,I),I1=1,24) #535 DO I = 1 , NUMSRC #536 IF ( QFLAG(I).EQ.'WSPEED' ) THEN ! 27 #546 WRITE (IOUNIT,9009) SRCID(I) , SRCTYP(I) ! 0 #549 WRITE (IOUNIT,9024) (QFACT(I2,I),I2=1,6) #555 DO I = 1 , NUMSRC #556 IF ( QFLAG(I).EQ.'SEASHR' ) THEN ! 27 #564 WRITE (IOUNIT,9009) SRCID(I) , SRCTYP(I) #569 WRITE (IOUNIT,9014) (I2,QFACT(I2+IFR,I),I2=1,24) #575 DO I = 1 , NUMSRC #576 IF ( QFLAG(I).EQ.'SHRDOW' ) THEN ! 27 #583 WRITE (IOUNIT,99009) SRCID(I) , SRCTYP(I) #591 WRITE (IOUNIT,99014) (I3,QFACT(I3+IFR+IDW,I),I3=1,24) #598 DO I = 1 , NUMSRC #599 IF ( QFLAG(I).EQ.'SHRDOW7' ) THEN ! 27 #606 WRITE (IOUNIT,99009) SRCID(I) , SRCTYP(I) #614 WRITE (IOUNIT,99014) (I3,QFACT(I3+IFR+IDW,I),I3=1,24) PRTSUM.f90#33 INTEGER :: IWHP(NVAL) , I , IVAL , INDMX , IHST , INDLOC #39 DATA (RANK(I),I=1,10)/' 1ST' , ' 2ND' , ' 3RD' , ' 4TH' , ' 5TH' ,& #140 DO I = 1 , IHST ! 4 #141 INDLOC = HMLOC(I,IGRP,IAVE,ITYP) ! 8 #142 IF ( I.EQ.1 ) THEN #164 WRITE (IOUNIT,1004) GRPID(IGRP) , RANK(IWHP(I)) , & #165 & HMAX(I,IGRP,IAVE,ITYP) , & #166 & HMCLM(I,IGRP,IAVE,ITYP) , & #167 & HMDATE(I,IGRP,IAVE,ITYP) , XR2 , YR2 , ZE2 ,& #178 WRITE (IOUNIT,1002) GRPID(IGRP) , RANK(IWHP(I)) , & #179 & HMAX(I,IGRP,IAVE,ITYP) , & #180 & HMCLM(I,IGRP,IAVE,ITYP) , & #181 & HMDATE(I,IGRP,IAVE,ITYP) , XR2 , YR2 , ZE2 ,& #193 WRITE (IOUNIT,1005) RANK(IWHP(I)) , & #194 & HMAX(I,IGRP,IAVE,ITYP) , & #195 & HMCLM(I,IGRP,IAVE,ITYP) , & #196 & HMDATE(I,IGRP,IAVE,ITYP) , XR2 , YR2 , ZE2 ,& #207 WRITE (IOUNIT,1003) RANK(IWHP(I)) , & #208 & HMAX(I,IGRP,IAVE,ITYP) , & #209 & HMCLM(I,IGRP,IAVE,ITYP) , & #210 & HMDATE(I,IGRP,IAVE,ITYP) , XR2 , YR2 , ZE2 ,& PSIDE.f90#73 INTEGER :: I , KS , IMINUS , IPLUS , NOUT , ICON #85 DO I = 1 , 2 #86 KS = IVERT + I - 1 !2102900 #87 U(I) = UVERT(KS) #88 V1(I) = VVERT(KS) #89 VN(I) = VNVERT(KS) #90 W(I) = WVERT(KS) #97 DO I = 1 , 2 #98 IF ( VN(I).LT.-3.9 ) IMINUS = I + IMINUS !2102900 #99 IF ( VN(I).GT.3.9 ) IPLUS = I + IPLUS PSIDE2.f90#38 INTEGER :: I , J , ISEG , NPTS , NOUT , ICON #49 DO I = 1 , NSEGS #50 ULIST(J) = UASEGS(I) ! 58748 #52 ULIST(J) = UBSEGS(I) #59 DO I = 1 , NSEGS #60 USIGN(I) = 1 ! 58748 #61 IF ( UASEGS(I).GT.UBSEGS(I) ) THEN #62 USIGN(I) = -1 ! 17946 #63 TEMP = UASEGS(I) #64 UASEGS(I) = UBSEGS(I) #65 UBSEGS(I) = TEMP #67 IF ( UASEGS(I).EQ.UBSEGS(I) ) USIGN(I) = 0 ! 58748 #71 DO I = 2 , NPTS #72 U1 = ULIST(I-1) ! 96734 #73 U2 = ULIST(I) PSIDE2_TOX.f90#47 INTEGER :: I , J , ISEG , NPTS , NOUT , ICON #58 DO I = 1 , NSEGS #59 ULIST(J) = UASEGS(I) ! 0 #61 ULIST(J) = UBSEGS(I) #68 DO I = 1 , NSEGS #69 USIGN(I) = 1 ! 0 #70 IF ( UASEGS(I).GT.UBSEGS(I) ) THEN #71 USIGN(I) = -1 ! 0 #72 TEMP = UASEGS(I) #73 UASEGS(I) = UBSEGS(I) #74 UBSEGS(I) = TEMP #76 IF ( UASEGS(I).EQ.UBSEGS(I) ) USIGN(I) = 0 ! 0 #80 DO I = 2 , NPTS #81 U1 = ULIST(I-1) ! 0 #82 U2 = ULIST(I) PSIDE_TOX.f90#81 INTEGER :: I , KS , IMINUS , IPLUS , NOUT , ICON #93 DO I = 1 , 2 #94 KS = IVERT + I - 1 ! 0 #95 U(I) = UVERT(KS) #96 V1(I) = VVERT(KS) #97 VN(I) = VNVERT(KS) #98 W(I) = WVERT(KS) #105 DO I = 1 , 2 #106 IF ( VN(I).LT.-3.9 ) IMINUS = I + IMINUS ! 0 #107 IF ( VN(I).GT.3.9 ) IPLUS = I + IPLUS PSTANN.f90#75 WRITE (IAPUNT(IGRP),9007) (MODOPS(I),I=1,18) QATR2.f90#3 SUBROUTINE QATR2(XL,XU,EPS,NDIM,FCT,Y,IER,I,AUX) #64 INTEGER I , II , JI , J , JJ #82 DO I = 2 , NDIM #98 AUX(I) = HALF*AUX(I-1) + P*SM ! 0 #103 JI = I - 1 #105 II = I - J ! 0 #115 IF ( I.GE.3 ) THEN QATR3.f90#3 SUBROUTINE QATR3(XL,XU,EPS,EPST,NDIM,IMIN,Y,YD,IER,I,AUX) #42 INTEGER :: I , II , JI , J , JJ #44 DO I = 1 , NDIM !1037136 #45 AUX(I) = 0.0 ! 12445K #46 AUXD(I) = 0.0 #68 DO I = 2 , NDIM #87 AUX(I) = 0.5*AUX(I-1) + P*SM !4135462 #88 AUXD(I) = 0.5*AUXD(I-1) + P*SMD #93 JI = I - 1 #95 II = I - J ! 10341K #107 IF ( I.GE.IMIN ) THEN RADRNG.f90#25 INTEGER :: I , J #31 DO I = 1 , IFC #32 IF ( FIELD(I).EQ.'DDIR' ) ISC = I + 1 ! 0 #45 DO I = ISC , IFC #46 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 RANKFL.f90#25 INTEGER :: I , J , IRANK , IDATSV(NMAX) #39 WRITE (IRKUNT(IAVE),9007) (MODOPS(I),I=1,18) #56 DO I = 1 , NMAX ! 0 #57 IDATSV(I) = 0 ! 0 #61 DO I = 1 , IMXVAL(IAVE) #64 IF ( MXDATE(I,IGRP,IAVE,ITYP).EQ.IDATSV(J) ) & #69 IREC = MXLOCA(I,IGRP,IAVE,ITYP) #70 IDATSV(IRANK) = MXDATE(I,IGRP,IAVE,ITYP) #72 & RMXVAL(I,IGRP,IAVE,ITYP) , & #73 & MXDATE(I,IGRP,IAVE,ITYP) , AXR(IREC) , & RATE.f90#56 DO I = 2 , NENTP1 #57 IF ( X.LE.XCAT(I) ) THEN ! 0 #58 ALPHA = ALPHAP(I-1) ! 0 #59 BETA = BETAP(I-1) RECARD.f90#28 INTEGER :: I , ILSAVE #51 DO I = 1 , NREC #52 ZETMP1(I) = 0.0 ! 432 #53 ZETMP2(I) = 0.0 #54 ZFTMP1(I) = 0.0 #55 ZFTMP2(I) = 0.0 RECART.f90#87 DO I = 1 , INNET - 1 #89 IF ( FIELD(3).EQ.NTID(I) ) & #90 & CALL ERRHDL(PATH,MODNAM,'W','252',NTID(I)) REPOLR.f90#25 INTEGER :: I , IORSET , IXRSET , IDRSET , IGRSET #85 DO I = 1 , INNET - 1 #87 IF ( FIELD(3).EQ.NTID(I) ) & #88 & CALL ERRHDL(PATH,MODNAM,'W','252',NTID(I)) RESINI.f90#34 INTEGER :: I , J , K , L , M #49 DO I = 1 , NHIVAL #50 HIVALU(J,I,K,L,M) = 0.0 ! 1152 #51 NHIDAT(J,I,K,L,M) = 0 #52 HCLMSG(J,I,K,L,M) = ' ' #53 HMAX(I,K,L,M) = 0.0 #54 HMDATE(I,K,L,M) = 0 #55 HMLOC(I,K,L,M) = 0 #56 HMCLM(I,K,L,M) = ' ' #85 DO I = 1 , NUMREC ! 2 #86 ANNVAL(I,J,K) = 0.0 ! 288 #87 SUMANN(I,J,K) = 0.0 #89 DO I = 1 , NHIANN ! 2 #90 AMXVAL(I,J,K) = 0.0 ! 20 #91 IMXLOC(I,J,K) = 0 #98 DO I = 1 , NUMREC ! 2 #99 SUMH4H(I,J) = 0.0 ! 288 #101 DO I = 1 , NMXPM ! 2 #102 MXPMVAL(I,J) = 0.0 ! 20 #103 MXPMLOC(I,J) = 0 #112 DO I = 1 , NUMREC ! 192 #114 SHVALS(I,J,K,L,M) = 0.0 ! 27648 #123 DO I = 1 , 4 ! 48 #125 NSEAHR(I,J) = 0 ! 192 #126 NSEACM(I,J) = 0 RSDUMP.f90#43 INTEGER :: I , J , K , L , M #60 WRITE (IDPUNT) (((((HIVALU(I,J,K,L,M),I=1,NUMREC),J=1,NHIVAL& #62 WRITE (IDPUNT) (((((NHIDAT(I,J,K,L,M),I=1,NUMREC),J=1,NHIVAL& #64 WRITE (IDPUNT) (((((HCLMSG(I,J,K,L,M),I=1,NUMREC),J=1,NHIVAL& #69 WRITE (IDPUNT) ((SUMH4H(I,J),I=1,NUMREC),J=1,NUMGRP) #75 WRITE (IDPUNT) ((((RMXVAL(I,J,K,L),I=1,NMXVAL),J=1,NUMGRP),K& #77 WRITE (IDPUNT) ((((MXDATE(I,J,K,L),I=1,NMXVAL),J=1,NUMGRP),K& #79 WRITE (IDPUNT) ((((MXLOCA(I,J,K,L),I=1,NMXVAL),J=1,NUMGRP),K& #81 WRITE (IDPUNT) ((((MCLMSG(I,J,K,L),I=1,NMXVAL),J=1,NUMGRP),K& #86 WRITE (IDPUNT) (((((SHVALS(I,J,K,L,M),I=1,NUMREC),J=1,NUMGRP& #88 WRITE (IDPUNT) ((NSEAHR(I,J),I=1,4),J=1,24) #89 WRITE (IDPUNT) ((NSEACM(I,J),I=1,4),J=1,24) #94 WRITE (IDPUNT) (((ANNVAL(I,J,K),I=1,NUMREC),J=1,NUMGRP),K=1,& #97 WRITE (IDPUNT) (((AMXVAL(I,J,K),I=1,NHIANN),J=1,NUMGRP), & #99 WRITE (IDPUNT) (((IMXLOC(I,J,K),I=1,NHIANN),J=1,NUMGRP), & #104 WRITE (IDPUNT) (((ANNVAL(I,J,K),I=1,NUMREC),J=1,NUMGRP),K=1,& #106 WRITE (IDPUNT) (((SUMANN(I,J,K),I=1,NUMREC),J=1,NUMGRP),K=1,& #120 WRITE (IDPUN2) (((((HIVALU(I,J,K,L,M),I=1,NUMREC),J=1,NHIVAL& #122 WRITE (IDPUN2) (((((NHIDAT(I,J,K,L,M),I=1,NUMREC),J=1,NHIVAL& #124 WRITE (IDPUN2) (((((HCLMSG(I,J,K,L,M),I=1,NUMREC),J=1,NHIVAL& #129 WRITE (IDPUN2) ((SUMH4H(I,J),I=1,NUMREC),J=1,NUMGRP) #135 WRITE (IDPUN2) ((((RMXVAL(I,J,K,L),I=1,NMXVAL),J=1,NUMGRP),K& #137 WRITE (IDPUN2) ((((MXDATE(I,J,K,L),I=1,NMXVAL),J=1,NUMGRP),K& #139 WRITE (IDPUN2) ((((MXLOCA(I,J,K,L),I=1,NMXVAL),J=1,NUMGRP),K& #141 WRITE (IDPUN2) ((((MCLMSG(I,J,K,L),I=1,NMXVAL),J=1,NUMGRP),K& #146 WRITE (IDPUN2) (((((SHVALS(I,J,K,L,M),I=1,NUMREC),J=1,NUMGRP& #148 WRITE (IDPUN2) ((NSEAHR(I,J),I=1,4),J=1,24) #149 WRITE (IDPUN2) ((NSEACM(I,J),I=1,4),J=1,24) #154 WRITE (IDPUN2) (((ANNVAL(I,J,K),I=1,NUMREC),J=1,NUMGRP),K=1,& #157 WRITE (IDPUN2) (((AMXVAL(I,J,K),I=1,NHIANN),J=1,NUMGRP), & #159 WRITE (IDPUN2) (((IMXLOC(I,J,K),I=1,NHIANN),J=1,NUMGRP), & #164 WRITE (IDPUN2) (((ANNVAL(I,J,K),I=1,NUMREC),J=1,NUMGRP),K=1,& #166 WRITE (IDPUN2) (((SUMANN(I,J,K),I=1,NUMREC),J=1,NUMGRP),K=1,& RSINIT.f90#33 INTEGER :: I , J , K , L , M #43 READ (IRSUNT,ERR=99,END=999) (((((HIVALU(I,J,K,L,M),I=1,NUMREC)& #46 READ (IRSUNT,ERR=99,END=999) (((((NHIDAT(I,J,K,L,M),I=1,NUMREC)& #49 READ (IRSUNT,ERR=99,END=999) (((((HCLMSG(I,J,K,L,M),I=1,NUMREC)& #56 & ((SUMH4H(I,J),I=1,NUMREC),J=1,NUMGRP) #62 READ (IRSUNT,ERR=99,END=999) ((((RMXVAL(I,J,K,L),I=1,NMXVAL),J=& #64 READ (IRSUNT,ERR=99,END=999) ((((MXDATE(I,J,K,L),I=1,NMXVAL),J=& #66 READ (IRSUNT,ERR=99,END=999) ((((MXLOCA(I,J,K,L),I=1,NMXVAL),J=& #68 READ (IRSUNT,ERR=99,END=999) ((((MCLMSG(I,J,K,L),I=1,NMXVAL),J=& #74 READ (IRSUNT,ERR=99,END=999) (((((SHVALS(I,J,K,L,M),I=1,NUMREC)& #77 READ (IRSUNT,ERR=99,END=999) ((NSEAHR(I,J),I=1,4),J=1,24) ! 0 #78 READ (IRSUNT,ERR=99,END=999) ((NSEACM(I,J),I=1,4),J=1,24) ! 0 #83 READ (IRSUNT,ERR=99,END=999) (((ANNVAL(I,J,K),I=1,NUMREC),J=1, & #87 READ (IRSUNT,ERR=99,END=999) (((ANNVAL(I,J,K),I=1,NUMREC),J=1, & #89 READ (IRSUNT,ERR=99,END=999) (((SUMANN(I,J,K),I=1,NUMREC),J=1, & #97 DO I = 1 , NUMREC ! 0 #98 ANNVAL(I,J,K) = 0.0 ! 0 #106 READ (IRSUNT,ERR=99,END=999) (((AMXVAL(I,J,K),I=1,NHIANN),J=1, & #108 READ (IRSUNT,ERR=99,END=999) (((IMXLOC(I,J,K),I=1,NHIANN),J=1, & SBYVAL.f90#26 INTEGER :: I , INX , JC , IMIN #39 DO I = JC , INX #40 IF ( ARRIN1(I).LT.RMIN ) THEN ! 0 #41 IMIN = I ! 0 #42 RMIN = ARRIN1(I) SCAVRAT.f90#47 INTEGER :: I , N , ILQ , IMISS #66 DO I = 1 , NPD ! 0 #67 PSCVRT(I) = 0.0 ! 0 #68 WASHOUT(I) = 0.0 #69 ECOLL(I) = 0.0 #90 DO I = 1 , NPD #98 TERM1 = (4./(REYNOLD*SCHMIDT(I))) & #99 & *(1.+0.4*SQRT(REYNOLD)*SCHMIDT(I)**0.333333+ & #100 & 0.16*SQRT(REYNOLD*SCHMIDT(I))) #104 KAPPA = (PDIAM(I)*1.E-6)/(RDROP*0.02) #110 STOKE = TSTOP(I)*(VFALL-VGRAV(I))/(RDROP*0.01) #118 TERM3 = TERM3*SQRT(1./PDENS(I)) #120 ECOLL(I) = MIN(1.0,TERM1+TERM2+TERM3) #124 WASHOUT(I) = 1.5*(ZSUBP*ECOLL(I))/(2.*RDROP*0.01) #129 PSCVRT(I) = 1.5*ECOLL(I)*PRATE/(2.*RDROP*3.6E4) ! 0 #131 PSCVRT(I) = 0.0 ! 0 SETCAR.f90#27 INTEGER :: I , J , JSET #38 DO I = 1 , ICOUNT ! 0 #59 AXR(ISET) = XCOORD(I,INNET) ! 0 #85 DO I = 1 , IZE #93 AZELEV(ISET) = ZETMP2(I) ! 0 #100 DO I = 1 , IZH #108 AZHILL(ISET) = ZHTMP2(I) ! 0 #115 DO I = 1 , IZF #123 AZFLAG(ISET) = ZFTMP2(I) ! 0 #127 DO I = IDC1 + 1 , IRXR #128 NETID(I) = NETIDT ! 0 #129 RECTYP(I) = 'GC' SETIDG.f90#34 INTEGER :: I , II , ISTR , IDNUM #41 I = 8 #49 DO WHILE ( .NOT.HIT .AND. I.GE.1 ) #50 CHKI = INID(I:I) ! 0 #52 II = I ! 0 #55 I = I - 1 ! 0 #59 I = 1 ! 0 #61 CHKI = INID(I:I) #64 IDCHR1 = INID(ISTR:I) ! 0 #65 I = I + 1 #66 IF ( I.GT.II ) THEN #69 CHKI = INID(I:I) ! 0 #74 ISTR = I ! 0 #76 NUMID = INID(ISTR:I) ! 0 #77 I = I + 1 #78 IF ( I.GT.II ) THEN #81 CHKI = INID(I:I) ! 0 #86 ISTR = I ! 0 #87 DO WHILE ( I.LE.II ) #88 IDCHR2 = INID(ISTR:I) ! 0 #89 I = I + 1 #90 IF ( I.GT.II ) THEN #93 CHKI = INID(I:I) ! 0 SETPOL.f90#27 INTEGER :: I , J , JSET #39 DO I = 1 , ICOUNT ! 108 #61 AXR(ISET) = XINT + XCOORD(I,INNET)*SIN(YTEMP) #62 AYR(JSET) = YINT + XCOORD(I,INNET)*COS(YTEMP) #89 DO I = 1 , IZE #97 AZELEV(ISET) = ZETMP2(I) ! 0 #104 DO I = 1 , IZH #112 AZHILL(ISET) = ZHTMP2(I) ! 0 #119 DO I = 1 , IZF #127 AZFLAG(ISET) = ZFTMP2(I) ! 0 #131 DO I = IDC1 + 1 , IRXR #132 NETID(I) = NETIDT ! 432 #133 RECTYP(I) = 'GP' SETUP.f90#44 INTEGER :: I , IFSTAT #80 READ (INUNIT,RDFRM,END=999) RUNST1 , (RUNST(I),I=1,ISTRG) SET_METDATA.f90#38 INTEGER :: I , IA , IY , IM , ID , NL , GINDEX , NUMSW #228 DO I = 1 , NPLVLS #229 IF ( PFLHT(I).GE.ZI .AND. PFLSW(I).GE.0.0 ) THEN ! 3044 #231 SUMSW = SUMSW + PFLSW(I) SHOUT.f90#52 WRITE (ISHUNT(IGRP),9007) (MODOPS(I),I=1,18) SINDEX.f90#22 INTEGER :: I , IDIM , INDEXS #30 I = 1 #33 DO WHILE ( .NOT.FIND .AND. I.LE.IDIM ) #34 IF ( ELEM.EQ.ARRIN(I) ) THEN ! 831 #38 I = I + 1 ! 831 SOCARD.f90#35 INTEGER :: I , J , ILSAVE #53 DO I = 1 , NSRC ! 3 #55 IWRK2(I,J) = 0 ! 351 SOGRP.f90#26 INTEGER :: I , K , IH , IL #47 DO I = 1 , NUMGRP #48 IF ( TEMPID.EQ.GRPID(I) ) CONT = .TRUE. ! 0 #67 DO I = 1 , NUMSRC ! 3 #68 IGROUP(I,IGRP) = 1 ! 27 #72 DO I = 4 , IFC ! 0 #73 CALL FSPLIT(PATH,KEYWRD,FIELD(I),ILEN_FLD,'-',RMARK,LOWID, & SOPARM.f90#31 INTEGER :: I , ISDX #63 DO I = 4 , IFC ! 27 #64 CALL STONUM(FIELD(I),ILEN_FLD,TEMP(I-3),IMIT) ! 117 SPRTHT.f90#30 INTEGER :: IHNUM , I , J , K , II , INDZ , INDC , NX , NY , ISRF ,& #36 DATA (CHRVAL(I),I=1,10)/' 1ST' , ' 2ND' , ' 3RD' , ' 4TH' , & #63 DO I = 1 , INNET #65 NPPX = 1 + INT((NUMXPT(I)-1)/NCPP) ! 8 #66 NPPY = 1 + INT((NUMYPT(I)-1)/NRPP) #74 WRITE (IOUNIT,9037) NTID(I) , NTTYP(I) #81 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 8 #84 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #85 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #88 & (XCOORD(J,I),J=1+NCPP*(NX-1),NUMXPT(I)) #91 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #94 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #95 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #98 & (XCOORD(J,I),J=1+NCPP*(NX-1),NCPP*NX) #104 DO K = 1 + NRPP*(NY-1) , NUMYPT(I) ! 8 #105 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 288 #106 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #107 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #108 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #109 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 288 #110 YCOVAL = YCOORD(K,I) #117 & ,J=1+NCPP*(NX-1),NUMXPT(I)) #128 IF ( NTTYP(I).EQ.'GRIDCART' ) THEN ! 0 #129 INDZ = NETEND(I) - K*NUMXPT(I) + 1 ! 0 #130 YCOVAL = YCOORD(NUMYPT(I)-K+1,I) #131 ELSEIF ( NTTYP(I).EQ.'GRIDPOLR' ) THEN #132 INDZ = NETSTA(I) + (K-1)*NUMXPT(I) ! 0 #133 YCOVAL = YCOORD(K,I) #140 & ,J=1+NCPP*(NX-1),NUMXPT(I)) SRCQA.f90#41 INTEGER :: I , J , K , N , ITOTSRC , ITOTGRP #48 DO I = 1 , NUMSRC #52 IF ( IWRK2(I,1).GT.0 .OR. IWRK2(I,2).GT.0 .OR. IWRK2(I,3) & #53 & .GT.0 .OR. IWRK2(I,11).GT.0 .OR. IWRK2(I,12).GT.0 .OR. & #54 & IWRK2(I,13).GT.0 ) THEN #57 IF ( IWRK2(I,1).LT.NSEC ) & #58 & CALL ERRHDL(PATH,MODNAM,'E','236',SRCID(I)) #60 IF ( IWRK2(I,2).LT.NSEC ) & #61 & CALL ERRHDL(PATH,MODNAM,'E','237',SRCID(I)) #65 IF ( IWRK2(I,11).LT.NSEC ) & #66 & CALL ERRHDL(PATH,MODNAM,'E','241',SRCID(I)) #68 IF ( IWRK2(I,12).LT.NSEC ) & #69 & CALL ERRHDL(PATH,MODNAM,'E','246',SRCID(I)) #71 IF ( IWRK2(I,13).LT.NSEC ) & #72 & CALL ERRHDL(PATH,MODNAM,'E','247',SRCID(I)) #76 IF ( QFLAG(I).NE.' ' ) THEN ! 27 #77 IF ( QFLAG(I).EQ.'SEASON' .AND. IWRK2(I,4).LT.4 ) THEN ! 0 #79 CALL ERRHDL(PATH,MODNAM,'E','239',SRCID(I)) ! 0 #80 ELSEIF ( QFLAG(I).EQ.'MONTH' .AND. IWRK2(I,4).LT.12 ) THEN #82 CALL ERRHDL(PATH,MODNAM,'E','239',SRCID(I)) ! 0 #83 ELSEIF ( QFLAG(I).EQ.'HROFDY' .AND. IWRK2(I,4).LT.24 ) THEN #85 CALL ERRHDL(PATH,MODNAM,'E','239',SRCID(I)) ! 0 #86 ELSEIF ( QFLAG(I).EQ.'WSPEED' .AND. IWRK2(I,4).LT.6 ) THEN #88 CALL ERRHDL(PATH,MODNAM,'E','239',SRCID(I)) ! 0 #89 ELSEIF ( QFLAG(I).EQ.'SEASHR' .AND. IWRK2(I,4).LT.96 ) THEN #91 CALL ERRHDL(PATH,MODNAM,'E','239',SRCID(I)) ! 0 #92 ELSEIF ( QFLAG(I).EQ.'SHRDOW' .AND. IWRK2(I,4).LT.288 ) THEN #94 CALL ERRHDL(PATH,MODNAM,'E','239',SRCID(I)) ! 0 #95 ELSEIF ( QFLAG(I).EQ.'SHRDOW7' .AND. IWRK2(I,4).LT.672 ) & #98 CALL ERRHDL(PATH,MODNAM,'E','239',SRCID(I)) ! 0 #103 IF ( IWRK2(I,5).NE.0 .OR. IWRK2(I,6).NE.0 .OR. IWRK2(I,7) & #106 INPD(I) = IWRK2(I,5) ! 0 #109 IF ( IWRK2(I,5).NE.IWRK2(I,6) .OR. IWRK2(I,5).NE.IWRK2(I,7) & #110 & ) CALL ERRHDL(PATH,MODNAM,'E','240',SRCID(I)) #113 N = INPD(I) #116 ATOT = ATOT + APHI(J,I) ! 0 #120 & CALL ERRHDL(PATH,MODNAM,'W','330',SRCID(I)) #123 CALL ERRHDL(PATH,MODNAM,'E','244',SRCID(I)) ! 0 #127 ELSEIF ( SRCTYP(I).EQ.'OPENPIT' ) THEN #129 CALL ERRHDL(PATH,MODNAM,'E','323',SRCID(I)) ! 0 #133 IF ( INPD(I).EQ.0 ) THEN ! 27 #136 IF ( (DEPOS .OR. DDEP) .AND. SOGAS(I).EQ.'N' .AND. & #138 & CALL ERRHDL(PATH,MODNAM,'E','242',SRCID(I)) #142 IF ( SRCTYP(I).EQ.'AREAPOLY' ) THEN ! 27 #143 IF ( IWRK2(I,10).LT.NVERTS(I) ) THEN ! 0 #145 CALL ERRHDL(PATH,MODNAM,'E','265',SRCID(I)) ! 0 #148 AXVERT(NVERTS(I)+1,I) = AXVERT(1,I) ! 0 #149 AYVERT(NVERTS(I)+1,I) = AYVERT(1,I) #154 DO J = 1 , NVERTS(I) #155 SUM = SUM + (AXVERT(J,I)*AYVERT(J+1,I)-AXVERT(J+1,I) & #156 & *AYVERT(J,I)) #163 AXINIT(I) = SQRT(ABS(AREA)) #164 AYINIT(I) = SQRT(ABS(AREA)) #169 DO J = 1 , NVERTS(I) #170 SUMX = SUMX + (AXVERT(J,I)+AXVERT(J+1,I)) & #171 & *(AXVERT(J,I)*AYVERT(J+1,I)-AXVERT(J+1,I) & #172 & *AYVERT(J,I)) #173 SUMY = SUMY + (AYVERT(J,I)+AYVERT(J+1,I)) & #174 & *(AXVERT(J,I)*AYVERT(J+1,I)-AXVERT(J+1,I) & #175 & *AYVERT(J,I)) #178 AXCNTR(I) = SUMX/(6.*AREA) ! 0 #179 AYCNTR(I) = SUMY/(6.*AREA) #185 IF ( URBSRC(I).EQ.'Y' ) NUMURB = NUMURB + 1 ! 27 #192 CALL LOCATE(GRIDHT,1,29,AHS(I),NDXSTK(I)) #200 DO I = 1 , NUMSRC #201 IF ( IGROUP(I,J).EQ.1 ) ITOTSRC = ITOTSRC + 1 ! 27 #212 DO I = 1 , NUMSRC #215 IF ( IGRP_OLM(I,J).EQ.1 ) ITOTGRP = ITOTGRP + 1 ! 0 #218 IF ( ITOTGRP.GT.1 ) CALL ERRHDL(PATH,MODNAM,'E','282',SRCID(I))! 27 #225 DO I = 1 , NUMSRC ! 0 #227 IF ( AQS(I).LT.0.0 ) & #228 & CALL ERRHDL(PATH,MODNAM,'E','338',SRCID(I)) STODBL.f90#36 I = 1 #49 DO WHILE ( .NOT.MEND .AND. I.LE.LEN ) #50 CHK = STRVAR(I:I) ! 0 #97 I = I + 1 ! 0 STONUM.f90#24 INTEGER :: I , IMUTI , LENGTH #31 I = 1 #44 DO WHILE ( .NOT.MEND .AND. I.LE.LENGTH ) #45 CHK = STRVAR(I:I) ! 13521 #92 I = I + 1 ! 13521 TERHGT.f90#32 INTEGER :: I , J , IZE1 #40 DO I = 1 , IFC #41 IF ( FIELD(I).EQ.'ELEV' ) ISC = I + 1 ! 0 #67 DO I = ISC + 1 , IFC #68 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 #87 DO I = IZE1 , IZE ! 0 #88 IF ( I.LE.NREC ) THEN ! 0 #89 ZETMP2(I) = 0.3048*ZETMP2(I) ! 0 #90 ZHTMP2(I) = 0.3048*ZHTMP2(I) TGQA.f90#50 DO I = 2 , NUMSRC #51 IF ( AXS(I).GT.XSMAX ) THEN ! 0 #52 XSMAX = AXS(I) ! 0 #53 ELSEIF ( AXS(I).LT.XSMIN ) THEN #54 XSMIN = AXS(I) ! 0 #56 IF ( AYS(I).GT.YSMAX ) THEN ! 0 #57 YSMAX = AYS(I) ! 0 #58 ELSEIF ( AYS(I).LT.YSMIN ) THEN #59 YSMIN = AYS(I) ! 0 #68 DO I = 2 , NUMREC #69 IF ( AXR(I).GT.XRMAX ) THEN ! 0 #70 XRMAX = AXR(I) ! 0 #71 ELSEIF ( AXR(I).LT.XRMIN ) THEN #72 XRMIN = AXR(I) ! 0 #74 IF ( AYR(I).GT.YRMAX ) THEN ! 0 #75 YRMAX = AYR(I) ! 0 #76 ELSEIF ( AYR(I).LT.YRMIN ) THEN #77 YRMIN = AYR(I) ! 0 #109 DO I = 1 , NUMSRC ! 0 #112 ZINT = ZINTERP(AXS(I),AYS(I)) ! 0 #114 DIFF = ABS(AZS(I)-ZINT) #117 & CALL ERRHDL(PATH,MODNAM,'W','393',SRCID(I)) #122 DO I = 1 , NUMREC ! 0 #125 ZINT = ZINTERP(AXR(I),AYR(I)) ! 0 #127 DIFF = ABS(AZELEV(I)-ZINT) #130 WRITE (DUMMY,'("RE#",i5.5)') I ! 0 TOXXFL.f90#27 INTEGER :: I , IG , ICODE #58 & (IDCONC(IAVE,I),I=1,NPAIR) #60 & (TXCONC(IAVE,I),I=1,NPAIR) URBANS.f90#28 INTEGER :: I , IL , IH , K #42 DO I = 3 , IFC ! 0 #43 CALL FSPLIT(PATH,KEYWRD,FIELD(I),ILEN_FLD,'-',RMARK,LOWID, & VARINI.f90#47 INTEGER :: I , J , K , IERRST , IWARN , NDAY #158 DO I = 1 , NSRC #159 AXS(I) = 0.0 ! 27 #160 AYS(I) = 0.0 #161 AZS(I) = 0.0 #162 AQS(I) = 0.0 #163 AHS(I) = 0.0 #164 ADS(I) = 0.0 #165 AVS(I) = 0.0 #166 ATS(I) = 0.0 #167 ASYINI(I) = 0.0 #168 ASZINI(I) = 0.0 #169 AXINIT(I) = 0.0 #170 AYINIT(I) = 0.0 #171 AANGLE(I) = 0.0 #172 RADIUS(I) = 0.0 #173 AXCNTR(I) = 0.0 #174 AYCNTR(I) = 0.0 #175 AALPHA(I) = 0.0 #176 APDEFF(I) = 0.0 #177 AVOLUM(I) = 0.0 #178 SOPCRD(I) = 'N' #179 SOGAS(I) = 'N' #180 URBSRC(I) = 'N' #181 SRCID(I) = ' ' #182 SRCTYP(I) = ' ' #183 QFLAG(I) = ' ' #184 INPD(I) = 0 #185 NVERTS(I) = 0 #186 IELUNT(I) = 0 #187 EVAL(I) = .FALSE. #188 EVLFIL(I) = ' ' #190 PDIFF(I) = 0.0 #191 ALPHAS(I) = 0.0 #192 REACT(I) = 0.0 #193 HENRY(I) = 0.0 #194 L_METHOD2(I) = .FALSE. #196 ADSBH(J,I) = 0.0 ! 972 #197 ADSBW(J,I) = 0.0 #200 ADSBL(J,I) = 0.0 #201 ADSXADJ(J,I) = 0.0 #202 ADSYADJ(J,I) = 0.0 #207 QFACT(J,I) = 0.0 ! 27 #210 IGROUP(I,J) = 0 ! 27 #213 APDIAM(J,I) = 0.0 ! 27 #214 APHI(J,I) = 0.0 #215 APDENS(J,I) = 0.0 #218 AXVERT(J,I) = 0.0 ! 648 #219 AYVERT(J,I) = 0.0 #223 IGRP_OLM(I,J) = 0 ! 0 #225 L_OLMGRP(I) = .FALSE. ! 0 #228 ANO2_RATIO(I) = 0.10 ! 0 #231 CHI(K,I,J) = 0.0 ! 0 #237 DO I = 1 , NPDMAX ! 3 #238 PDIAM(I) = 0.0 ! 3 #239 PHI(I) = 0.0 #240 PDENS(I) = 0.0 #241 VGRAV(I) = 0.0 #242 TSTOP(I) = 0.0 #243 VDEP(I) = 0.0 #244 PSCVRT(I) = 0.0 #245 SCF(I) = 0.0 #246 WASHOUT(I) = 0.0 #252 DO I = 1 , NPDMAX #253 WQCOR(I) = 1.0 ! 3 #254 DQCOR(I) = 1.0 #327 DO I = 1 , NEVE ! 0 #328 AXR(I) = 0.0 ! 0 #329 AYR(I) = 0.0 #330 AZELEV(I) = 0.0 #331 AZHILL(I) = 0.0 #332 AZFLAG(I) = 0.0 #333 EVAPER(I) = 0 #334 EVDATE(I) = 0 #335 EVJDAY(I) = 0 #336 IDXEV(I) = 0 #337 EVNAME(I) = ' ' #338 EVGRP(I) = ' ' #364 DO I = 1 , MXPLVL #365 AIFLAG(J,I) = 0 ! 0 #366 APFLHT(J,I) = 0.0 #367 APFLWD(J,I) = 0.0 #368 APFLWS(J,I) = 0.0 #369 APFLTA(J,I) = 0.0 #370 APFLSA(J,I) = 0.0 #371 APFLSV(J,I) = 0.0 #372 APFLSW(J,I) = 0.0 #374 DO I = 1 , NSRC ! 0 #375 EV_HRQS(I,J) = 0.0 ! 0 #376 EV_HRTS(I,J) = 0.0 #377 EV_HRVS(I,J) = 0.0 #381 DO I = 1 , NREC ! 3 #382 AXR(I) = 0.0 ! 432 #383 AYR(I) = 0.0 #384 AZELEV(I) = 0.0 #385 AZHILL(I) = 0.0 #386 AZFLAG(I) = 0.0 #387 IREF(I) = 0 #388 NETID(I) = ' ' #389 RECTYP(I) = ' ' #390 NDXARC(I) = 0 #392 DO I = 1 , NNET ! 3 #393 XORIG(I) = 0.0 ! 3 #394 YORIG(I) = 0.0 #395 NUMXPT(I) = 0 #396 NUMYPT(I) = 0 #397 NETSTA(I) = 0 #398 NETEND(I) = 0 #399 NTID(I) = ' ' #400 NTTYP(I) = ' ' #402 XCOORD(J,I) = 0.0 ! 12 #405 YCOORD(J,I) = 0.0 ! 108 #412 DO I = 1 , NAVE ! 3 #413 KAVE(I) = 0 ! 6 #439 DO I = 1 , NPAIR #440 IDCONC(J,I) = 0 ! 600 #441 TXCONC(J,I) = 0.0 #443 DO I = 1 , NVAL ! 6 #444 NHIAVE(I,J) = 0 ! 12 #446 DO I = 1 , NGRP ! 6 #447 MAXFLE(I,J) = 0 ! 6 #448 IPSTFL(I,J) = 0 #449 IMXUNT(I,J) = 0 #450 IPSUNT(I,J) = 0 #451 IPSFRM(I,J) = 0 #452 THRESH(I,J) = 0.0 #453 THRFIL(I,J) = ' ' #454 PSTFIL(I,J) = ' ' #456 IPLTFL(K,I,J) = 0 ! 12 #457 IPLUNT(K,I,J) = 0 #458 PLTFIL(K,I,J) = ' ' #462 DO I = 1 , NGRP ! 3 #463 GRPID(I) = ' ' ! 3 #464 IANPST(I) = 0 #465 IANFRM(I) = 0 #466 IANPLT(I) = 0 #467 IAPUNT(I) = 0 #468 IPPUNT(I) = 0 #469 ISHUNT(I) = 0 #470 ISEAHR(I) = 0 #471 ANNPST(I) = ' ' #472 ANNPLT(I) = ' ' #475 DO I = 1 , NOLM ! 0 #476 OLMID(I) = ' ' ! 0 VDP.f90#117 INTEGER :: I , N #127 DATA (((RESTAB(I,J,K),I=1,9),J=1,6),K=1,5)/1.E07 , 60. , 120. , & #383 DO I = 1 , NPD #385 ST = TSTOP(I)*T1 ! 0 #391 D_SUB_B = 8.09E-14*(TA*SCF(I)/PDIAM(I)) #392 SCHMIDT(I) = XNU/D_SUB_B #403 RD(I) = 1.0/(GUST_ADJ*USTAR*(SCHMIDT(I)**(-2./3.)+XINERT)& #407 VDEP(I) = 1.0/(RA+RD(I)+RA*RD(I)*VGRAV(I)) + VGRAV(I) & #452 DO I = 1 , NPD ! 0 #453 WRITE (101,1002) KURDAT , ISRC , I , RA , RD(I) , VGRAV(I) ,& #454 & VDEP(I) VDP1.f90#55 INTEGER :: I , J , N , IO6 #63 DO I = 1 , NUMSRC #64 WRITE (IO6,*) 'SOURCE = ' , I ! 0 #65 WRITE (IO6,*) 'INPD = ' , INPD(I) #66 IF ( INPD(I).GT.0 ) THEN #68 & (APDIAM(N,I),N=1,INPD(I)) #70 & (APHI(N,I),N=1,INPD(I)) #72 & (APDENS(N,I),N=1,INPD(I)) #123 DO I = 1 , INPD(J) ! 0 #125 DIAMCM = 1.E-4*APDIAM(I,J) ! 0 #126 SCF(I) = 1. + 2.0*XMFP*(A1+A2*EXP(-A3*DIAMCM/XMFP)) & #133 AVGRAV(I,J) = (APDENS(I,J)-RHOAIR)*G*DIAMCM**2*SCF(I) & #139 ATSTOP(I,J) = AVGRAV(I,J)/G #146 DO I = 1 , NUMSRC #147 WRITE (IO6,*) 'SOURCE = ' , I ! 0 #149 & (AVGRAV(N,I),N=1,INPD(I)) #151 & (ATSTOP(N,I),N=1,INPD(I)) VRTCBL.f90#63 DO I = 0 , 1000 #66 TWOIZI = 2.*I*ZI*FACT #68 IF ( FACT.LT.0. .AND. I.EQ.0 ) GOTO 50 #100 DO I = 0 , 1000 #103 TWOIZI = 2.*I*ZITMP*FACT #105 IF ( FACT.LT.0. .AND. I.EQ.0 ) GOTO 150 VRTSBL.f90#53 DO I = 1 , 100 #56 TWOIZI = 2.*I*ZIARG #82 DO I = 1 , 100 #84 TWOIZI = 2.*I*ZITMP WAKE_DFSN.f90#536 DO I = 1 , NWAK #537 WRITE (IO6,32) I , XWAK(I) , SZWAK(I) , SYWAK(I) , DRWAK(I) ! 0 #544 DO I = 1 , NCAV #545 WRITE (IO6,33) I , XCAV(I) , SZCAV(I) , SYCAV(I) ! 0 WAKE_DRDX.f90#41 DO I = NWKM1 , 1 , -1 #42 IF ( X.GE.XWAK(I) ) THEN ! 28645K #43 IP1 = I + 1 !1358844 #44 DRDX = DRWAK(IP1) - (DRWAK(IP1)-DRWAK(I))*(XWAK(IP1)-X) & #45 & /(XWAK(IP1)-XWAK(I)) WAKE_XSIG.f90#89 DO I = NWKM1 , 1 , -1 #90 IF ( X.GE.XWAK(I) ) THEN ! 10189K #91 IP1 = I + 1 ! 412294 #92 FAC = (XWAK(IP1)-X)/(XWAK(IP1)-XWAK(I)) #93 SZ = SZWAK(IP1) - (SZWAK(IP1)-SZWAK(I))*FAC #94 SY = SYWAK(IP1) - (SYWAK(IP1)-SYWAK(I))*FAC #122 DO I = NCVM1 , 1 , -1 #123 IF ( X.GE.XCAV(I) ) THEN !9691098 #124 IP1 = I + 1 ! 375616 #125 FAC = (XCAV(IP1)-X)/(XCAV(IP1)-XCAV(I)) #126 SZC = SZCAV(IP1) - (SZCAV(IP1)-SZCAV(I))*FAC #127 SYC = SYCAV(IP1) - (SYCAV(IP1)-SYCAV(I))*FAC WSCATS.f90#30 INTEGER :: I , IWS #37 DO I = 3 , IFC ! 0 #38 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 #46 IWS = I - 2 ! 0 XYPNTS.f90#28 INTEGER :: I , J , JSET #35 DO I = 1 , IFC ! 0 #36 IF ( FIELD(I).EQ.'XPNTS' ) ISC = I + 1 ! 0 #48 DO I = ISC , IFC #49 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 #74 DO I = 1 , IFC ! 0 #75 IF ( FIELD(I).EQ.'YPNTS' ) ISC = I + 1 ! 0 #88 DO I = ISC , IFC #89 CALL STONUM(FIELD(I),ILEN_FLD,FNUM,IMIT) ! 0 ZIAVER.f90#57 DO I = 2 , NDXBLW #58 SUM = SUM + (HTS(I)-HTS(I-1))*0.5*(PARRAY(I)+PARRAY(I-1)) !2001510
HyperKWIC - Version 1.00DD executed at 20:00 on 1 Mar 2018 | Personal or Academic or Evaluation User | Free for Non-Commercial, Non-Government Use