ACPARM.f90#2 SUBROUTINE ACPARM(ISDX,TEMP) #25 INTEGER :: ISDX #31 AQS(ISDX) = TEMP(1) #32 AHS(ISDX) = TEMP(2) #33 RADIUS(ISDX) = TEMP(3) #35 NVERTS(ISDX) = NINT(TEMP(4)) ! 0 #37 NVERTS(ISDX) = NVMAX - 4 ! 0 #40 ASZINI(ISDX) = TEMP(5) ! 0 #42 ASZINI(ISDX) = 0. ! 0 #48 IF ( AQS(ISDX).EQ.0.0 ) CALL ERRHDL(PATH,MODNAM,'W','320',' QS ') ! 0 #50 IF ( AHS(ISDX).LT.0.0 ) THEN #53 ELSEIF ( AHS(ISDX).GT.100.0 ) THEN #58 IF ( RADIUS(ISDX).LE.0.0 ) THEN ! 0 #61 ELSEIF ( RADIUS(ISDX).GT.10000.0 ) THEN #66 IF ( ASZINI(ISDX).LT.0.0 ) THEN ! 0 #69 ELSEIF ( ASZINI(ISDX).LT.1.0E-5 ) THEN #71 ASZINI(ISDX) = 1.0E-5 ! 0 #72 ELSEIF ( ASZINI(ISDX).GT.200.0 ) THEN #77 IF ( NVERTS(ISDX).LT.3 ) THEN ! 0 #81 ELSEIF ( NVERTS(ISDX).GT.NVMAX-4 ) THEN #88 CALL GENCIR(ISDX) ! 0 #91 AXCNTR(ISDX) = AXS(ISDX) #92 AYCNTR(ISDX) = AYS(ISDX) APARM.f90#2 SUBROUTINE APARM(ISDX,TEMP) #34 INTEGER :: I , J , ISDX #40 AQS(ISDX) = TEMP(1) #41 AHS(ISDX) = TEMP(2) #43 AXINIT(ISDX) = TEMP(3) ! 0 #44 AYINIT(ISDX) = AXINIT(ISDX) #45 AANGLE(ISDX) = 0. #46 ASZINI(ISDX) = 0. #48 AXINIT(ISDX) = TEMP(3) ! 9 #49 AYINIT(ISDX) = TEMP(4) #50 AANGLE(ISDX) = 0. #51 ASZINI(ISDX) = 0. #53 AXINIT(ISDX) = TEMP(3) ! 0 #54 AYINIT(ISDX) = TEMP(4) #55 AANGLE(ISDX) = TEMP(5) #56 ASZINI(ISDX) = 0. #58 AXINIT(ISDX) = TEMP(3) ! 0 #59 AYINIT(ISDX) = TEMP(4) #60 AANGLE(ISDX) = TEMP(5) #61 ASZINI(ISDX) = TEMP(6) #67 IF ( AQS(ISDX).EQ.0.0 ) CALL ERRHDL(PATH,MODNAM,'W','320',' QS ') ! 9 #69 IF ( AHS(ISDX).LT.0.0 ) THEN #72 ELSEIF ( AHS(ISDX).GT.100.0 ) THEN #78 IF ( AXINIT(ISDX).LT.0.0 ) & #81 IF ( AYINIT(ISDX).LT.0.0 ) & #83 IF ( AXINIT(ISDX).LT.1.0E-5 ) THEN #87 AXINIT(ISDX) = 1.0E-5 #89 IF ( AYINIT(ISDX).LT.1.0E-5 ) THEN ! 9 #93 AYINIT(ISDX) = 1.0E-5 #96 IF ( AXINIT(ISDX).GT.2000.0 ) & #99 IF ( AYINIT(ISDX).GT.2000.0 ) & #102 IF ( ABS(AANGLE(ISDX)).GT.180. ) & #106 IF ( ASZINI(ISDX).LT.0.0 ) THEN #109 ELSEIF ( ASZINI(ISDX).LT.1.0E-5 ) THEN #111 ASZINI(ISDX) = 1.0E-5 ! 9 #112 ELSEIF ( ASZINI(ISDX).GT.200.0 ) THEN #120 IF ( AYINIT(ISDX)/AXINIT(ISDX).GT.10.00001 .OR. AXINIT(ISDX) & #121 & /AYINIT(ISDX).GT.10.00001 ) & #122 & CALL ERRHDL(PATH,MODNAM,'W','391',SRCID(ISDX)) #125 NVERTS(ISDX) = 4 #131 AXVERT(1,ISDX) = AXS(ISDX) #132 AYVERT(1,ISDX) = AYS(ISDX) #134 AXVERT(2,ISDX) = AXVERT(1,ISDX) & #135 & + (AYINIT(ISDX)*SIN(AANGLE(ISDX)*DTORAD)) #136 AYVERT(2,ISDX) = AYVERT(1,ISDX) & #137 & + (AYINIT(ISDX)*COS(AANGLE(ISDX)*DTORAD)) #139 AXVERT(3,ISDX) = AXVERT(2,ISDX) & #140 & + (AXINIT(ISDX)*COS(AANGLE(ISDX)*DTORAD)) #141 AYVERT(3,ISDX) = AYVERT(2,ISDX) & #142 & - (AXINIT(ISDX)*SIN(AANGLE(ISDX)*DTORAD)) #144 AXVERT(4,ISDX) = AXVERT(3,ISDX) & #145 & - (AYINIT(ISDX)*SIN(AANGLE(ISDX)*DTORAD)) #146 AYVERT(4,ISDX) = AYVERT(3,ISDX) & #147 & - (AYINIT(ISDX)*COS(AANGLE(ISDX)*DTORAD)) #149 AXVERT(5,ISDX) = AXS(ISDX) #150 AYVERT(5,ISDX) = AYS(ISDX) #153 AXCNTR(ISDX) = 0.0 #154 AYCNTR(ISDX) = 0.0 #155 DO I = 1 , NVERTS(ISDX) #156 AXCNTR(ISDX) = AXCNTR(ISDX) + AXVERT(I,ISDX) ! 36 #157 AYCNTR(ISDX) = AYCNTR(ISDX) + AYVERT(I,ISDX) #159 AXCNTR(ISDX) = AXCNTR(ISDX)/NVERTS(ISDX) ! 9 #160 AYCNTR(ISDX) = AYCNTR(ISDX)/NVERTS(ISDX) APPARM.f90#2 SUBROUTINE APPARM(ISDX,TEMP) #27 INTEGER :: ISDX #33 AQS(ISDX) = TEMP(1) #34 AHS(ISDX) = TEMP(2) #35 NVERTS(ISDX) = NINT(TEMP(3)) #37 ASZINI(ISDX) = TEMP(4) ! 0 #39 ASZINI(ISDX) = 0. ! 0 #45 IF ( AQS(ISDX).EQ.0.0 ) CALL ERRHDL(PATH,MODNAM,'W','320',' QS ') ! 0 #47 IF ( AHS(ISDX).LT.0.0 ) THEN #50 ELSEIF ( AHS(ISDX).GT.100.0 ) THEN #55 IF ( ASZINI(ISDX).LT.0.0 ) THEN ! 0 #58 ELSEIF ( ASZINI(ISDX).LT.1.0E-5 ) THEN #60 ASZINI(ISDX) = 1.0E-5 ! 0 #61 ELSEIF ( ASZINI(ISDX).GT.200.0 ) THEN #66 IF ( NVERTS(ISDX).LT.3 ) THEN ! 0 #69 ELSEIF ( NVERTS(ISDX).GT.NVMAX-2 ) THEN #72 ELSEIF ( NVERTS(ISDX).GT.NVMAX-4 ) THEN ARVERT.f90#25 INTEGER :: K , ISDX #38 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) #40 ISET = IWRK2(ISDX,10) ! 0 #59 IF ( FNUMX.NE.AXS(ISDX) ) & #60 & CALL ERRHDL(PATH,MODNAM,'E','262',SRCID(ISDX)) #61 IF ( FNUMY.NE.AYS(ISDX) ) & #62 & CALL ERRHDL(PATH,MODNAM,'E','262',SRCID(ISDX)) #65 IF ( ISET.LE.NVERTS(ISDX) ) THEN ! 0 #67 AXVERT(ISET,ISDX) = FNUMX ! 0 #68 AYVERT(ISET,ISDX) = FNUMY #71 CALL ERRHDL(PATH,MODNAM,'E','264',SRCID(ISDX)) ! 0 #74 IWRK2(ISDX,10) = ISET ! 0 DISPOL.f90#30 INTEGER :: I1 , I2 , I3 , I4 , I5 , ISDX #166 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 0 #171 AXR(I1+1) = AXS(ISDX) + RANGE*SIN(DIRECT*DTORAD) ! 0 #172 AYR(I2+1) = AYS(ISDX) + RANGE*COS(DIRECT*DTORAD) #190 IREF(IRXR) = ISDX DSBLDG.f90#26 INTEGER :: I , IH , IL , ISDX #54 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 225 #56 IF ( SRCTYP(ISDX).EQ.'POINT' ) THEN ! 225 #58 CALL DSFILL(ISDX) ! 225 #61 CALL ERRHDL(PATH,MODNAM,'W','233',SRCID(ISDX)) ! 0 #80 ISDX = I ! 0 #82 CALL DSFILL(ISDX) DSFILL.f90#2 SUBROUTINE DSFILL(ISDX) #31 INTEGER :: J , K , ISDX #37 ISET = IWRK2(ISDX,1) ! 9 #50 ADSBH(ISET,ISDX) = FNUM ! 324 #60 IWRK2(ISDX,1) = ISET ! 9 #62 ISET = IWRK2(ISDX,2) ! 54 #75 ADSBW(ISET,ISDX) = FNUM ! 324 #85 IWRK2(ISDX,2) = ISET ! 54 #90 ISET = IWRK2(ISDX,11) ! 54 #103 ADSBL(ISET,ISDX) = FNUM ! 324 #113 IWRK2(ISDX,11) = ISET ! 54 #117 ISET = IWRK2(ISDX,12) ! 54 #130 ADSXADJ(ISET,ISDX) = FNUM ! 324 #137 IWRK2(ISDX,12) = ISET ! 54 #141 ISET = IWRK2(ISDX,13) ! 54 #154 ADSYADJ(ISET,ISDX) = FNUM ! 324 #161 IWRK2(ISDX,13) = ISET ! 54 EFFILL.f90#2 SUBROUTINE EFFILL(ISDX,IQMAX) #25 INTEGER :: J , K , ISDX , IQMAX #30 ISET = IWRK2(ISDX,4) #44 QFACT(ISET,ISDX) = FNUM ! 0 #55 IWRK2(ISDX,4) = ISET ! 0 EMVARY.f90#32 INTEGER :: I , IH , IL , ISDX , IQMAX #64 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 0 #66 QFLAG(ISDX) = FIELD(4) ! 0 #67 IF ( QFLAG(ISDX).EQ.'SEASON' ) THEN #69 ELSEIF ( QFLAG(ISDX).EQ.'MONTH' ) THEN #71 ELSEIF ( QFLAG(ISDX).EQ.'HROFDY' ) THEN #73 ELSEIF ( QFLAG(ISDX).EQ.'WSPEED' ) THEN #75 ELSEIF ( QFLAG(ISDX).EQ.'SEASHR' ) THEN #77 ELSEIF ( QFLAG(ISDX).EQ.'SHRDOW' ) THEN #79 ELSEIF ( QFLAG(ISDX).EQ.'SHRDOW7' ) THEN #86 CALL EFFILL(ISDX,IQMAX) ! 0 #109 ISDX = I ! 0 #110 QFLAG(ISDX) = FIELD(4) #111 IF ( QFLAG(ISDX).EQ.'SEASON' ) THEN #113 ELSEIF ( QFLAG(ISDX).EQ.'MONTH' ) THEN #115 ELSEIF ( QFLAG(ISDX).EQ.'HROFDY' ) THEN #117 ELSEIF ( QFLAG(ISDX).EQ.'WSPEED' ) THEN #119 ELSEIF ( QFLAG(ISDX).EQ.'SEASHR' ) THEN #121 ELSEIF ( QFLAG(ISDX).EQ.'SHRDOW' ) THEN #123 ELSEIF ( QFLAG(ISDX).EQ.'SHRDOW7' ) THEN #130 CALL EFFILL(ISDX,IQMAX) ! 0 EVLOC.f90#25 INTEGER :: ISDX #51 CALL SINDEX(EVNAME,NEVE,USEVN,ISDX,FIND) #84 AZELEV(ISDX) = FNUM ! 0 #91 AZHILL(ISDX) = FNUM ! 0 #94 AZELEV(ISDX) = 0. ! 0 #95 AZHILL(ISDX) = 0. #104 AZFLAG(ISDX) = FNUM ! 0 #107 AZFLAG(ISDX) = 0. ! 0 #111 AXR(ISDX) = SETAXR ! 0 #112 AYR(ISDX) = SETAYR #114 AXR(ISDX) = SETAXR*SIN(SETAYR*DTORAD) ! 0 #115 AYR(ISDX) = SETAXR*COS(SETAYR*DTORAD) EVPER.f90#31 INTEGER :: ISDX , IMN , IDY , IEVYR2 , IEVYR4 #69 CALL SINDEX(EVNAME,NEVE,USEVN,ISDX,FIND) #74 CALL ERRHDL(PATH,MODNAM,'E','313',EVNAME(ISDX)) ! 0 #105 CALL SINDEX(GRPID,NGRP,EVGRP(IEVENT),ISDX,FIND) #110 IDXEV(IEVENT) = ISDX ! 0 GASDEP.f90#28 INTEGER :: I , IH , IL , ISDX #63 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 0 #65 SOGAS(ISDX) = 'Y' ! 0 #78 PDIFF(ISDX) = FNUM ! 0 #89 PDIFFW(ISDX) = FNUM ! 0 #100 RCLI(ISDX) = FNUM ! 0 #111 HENRY(ISDX) = FNUM ! 0 #130 ISDX = I ! 0 #131 SOGAS(ISDX) = 'Y' #144 PDIFF(ISDX) = FNUM ! 0 #155 PDIFFW(ISDX) = FNUM ! 0 #166 RCLI(ISDX) = FNUM ! 0 #177 HENRY(ISDX) = FNUM ! 0 GENCIR.f90#2 SUBROUTINE GENCIR(ISDX) #27 INTEGER :: I , ISDX , NSIDES , NEWRAD #33 NSIDES = NVERTS(ISDX) #38 AREA = PI*RADIUS(ISDX)*RADIUS(ISDX) #47 AXVERT(I,ISDX) = (NEWRAD*SIN(ANG/RTODEG)) + AXS(ISDX) #48 AYVERT(I,ISDX) = (NEWRAD*COS(ANG/RTODEG)) + AYS(ISDX) #52 AXVERT(NSIDES+1,ISDX) = AXVERT(1,ISDX) ! 0 #53 AYVERT(NSIDES+1,ISDX) = AYVERT(1,ISDX) #57 AXINIT(ISDX) = SQRT(AREA) #58 AYINIT(ISDX) = SQRT(AREA) INPPDM.f90#28 INTEGER :: I , J , K , IH , IL , ISDX , IING #44 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 0 #46 ISET = IWRK2(ISDX,5) ! 0 #59 APDIAM(ISET,ISDX) = FNUM ! 0 #67 IWRK2(ISDX,5) = ISET ! 0 INPPDN.f90#28 INTEGER :: I , J , K , IH , IL , ISDX , IING #44 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 0 #46 ISET = IWRK2(ISDX,7) ! 0 #57 & CALL ERRHDL(PATH,MODNAM,'E','334',SRCID(ISDX)) #62 APDENS(ISET,ISDX) = FNUM ! 0 #70 IWRK2(ISDX,7) = ISET ! 0 #100 & CALL ERRHDL(PATH,MODNAM,'E','334',SRCID(ISDX)) INPPHI.f90#26 INTEGER :: I , J , K , IH , IL , ISDX , IING #42 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 0 #44 ISET = IWRK2(ISDX,6) ! 0 #55 & CALL ERRHDL(PATH,MODNAM,'E','332',SRCID(ISDX)) #60 APHI(ISET,ISDX) = FNUM ! 0 #68 IWRK2(ISDX,6) = ISET ! 0 #98 & CALL ERRHDL(PATH,MODNAM,'E','332',SRCID(ISDX)) METH_2.f90#25 INTEGER :: I , IH , IL , ISDX #56 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 0 #58 L_METHOD2(ISDX) = .TRUE. ! 0 #69 FINEMASS(ISDX) = FNUM ! 0 #79 APDIAM(1,ISDX) = FNUM ! 0 #82 APHI(1,ISDX) = 1.0 #83 APDENS(1,ISDX) = 1.0 #86 INPD(ISDX) = 1 #105 ISDX = I ! 0 #106 L_METHOD2(ISDX) = .TRUE. #117 FINEMASS(ISDX) = FNUM ! 0 #127 APDIAM(1,ISDX) = FNUM ! 0 #130 APHI(1,ISDX) = 1.0 #131 APDENS(1,ISDX) = 1.0 #134 INPD(ISDX) = 1 NO2RAT.f90#26 INTEGER :: I , IH , IL , ISDX , IING , IPINDEX #57 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) ! 0 #68 & CALL ERRHDL(PATH,MODNAM,'E','336',SRCID(ISDX)) #70 ANO2_RATIO(ISDX) = FNUM OPARM.f90#2 SUBROUTINE OPARM(ISDX,TEMP) #30 INTEGER :: I , ISDX #36 AQS(ISDX) = TEMP(1) #37 AHS(ISDX) = TEMP(2) #38 AXINIT(ISDX) = TEMP(3) #39 AYINIT(ISDX) = TEMP(4) #40 AVOLUM(ISDX) = TEMP(5) #41 AANGLE(ISDX) = 0. #42 IF ( IFC.EQ.9 ) AANGLE(ISDX) = TEMP(6) #47 IF ( AQS(ISDX).EQ.0.0 ) CALL ERRHDL(PATH,MODNAM,'W','320',' QS ') #49 IF ( AHS(ISDX).LT.0.0 ) THEN #52 ELSEIF ( AHS(ISDX).GT.200.0 ) THEN #58 IF ( AXINIT(ISDX).LT.0.0 ) & #61 IF ( AYINIT(ISDX).LT.0.0 ) & #63 IF ( AXINIT(ISDX).LT.1.0E-5 ) THEN #67 AXINIT(ISDX) = 1.0E-5 #69 IF ( AYINIT(ISDX).LT.1.0E-5 ) THEN ! 0 #73 AYINIT(ISDX) = 1.0E-5 #76 IF ( AXINIT(ISDX).GT.2000.0 ) & #79 IF ( AYINIT(ISDX).GT.2000.0 ) & #82 IF ( ABS(AANGLE(ISDX)).GT.180. ) & #86 IF ( AVOLUM(ISDX).LE.0.0 ) & #91 IF ( AYINIT(ISDX)/AXINIT(ISDX).GT.10.0 .OR. AXINIT(ISDX) & #92 & /AYINIT(ISDX).GT.10.0 ) & #93 & CALL ERRHDL(PATH,MODNAM,'W','392',SRCID(ISDX)) #96 EFFDEP = AVOLUM(ISDX)/(AXINIT(ISDX)*AYINIT(ISDX)) #98 IF ( AHS(ISDX).GT.EFFDEP ) & #99 & CALL ERRHDL(PATH,MODNAM,'E','322',SRCID(ISDX)) #108 AXVERT(1,ISDX) = AXS(ISDX) #109 AYVERT(1,ISDX) = AYS(ISDX) #111 AXVERT(2,ISDX) = AXVERT(1,ISDX) & #112 & + (AYINIT(ISDX)*SIN(AANGLE(ISDX)*DTORAD)) #113 AYVERT(2,ISDX) = AYVERT(1,ISDX) & #114 & + (AYINIT(ISDX)*COS(AANGLE(ISDX)*DTORAD)) #116 AXVERT(3,ISDX) = AXVERT(2,ISDX) & #117 & + (AXINIT(ISDX)*COS(AANGLE(ISDX)*DTORAD)) #118 AYVERT(3,ISDX) = AYVERT(2,ISDX) & #119 & - (AXINIT(ISDX)*SIN(AANGLE(ISDX)*DTORAD)) #121 AXVERT(4,ISDX) = AXVERT(3,ISDX) & #122 & - (AYINIT(ISDX)*SIN(AANGLE(ISDX)*DTORAD)) #123 AYVERT(4,ISDX) = AYVERT(3,ISDX) & #124 & - (AYINIT(ISDX)*COS(AANGLE(ISDX)*DTORAD)) #126 AXVERT(5,ISDX) = AXS(ISDX) #127 AYVERT(5,ISDX) = AYS(ISDX) #130 IF ( AYINIT(ISDX).GE.AXINIT(ISDX) ) THEN #131 AALPHA(ISDX) = AANGLE(ISDX) ! 0 #132 ELSEIF ( AXINIT(ISDX).GT.AYINIT(ISDX) ) THEN #133 AALPHA(ISDX) = AANGLE(ISDX) + 90.0 ! 0 #137 APDEFF(ISDX) = AVOLUM(ISDX)/(AXINIT(ISDX)*AYINIT(ISDX)) ! 0 #140 ASZINI(ISDX) = APDEFF(ISDX)/4.3 #143 AXCNTR(ISDX) = 0.0 #144 AYCNTR(ISDX) = 0.0 #146 AXCNTR(ISDX) = AXCNTR(ISDX) + AXVERT(I,ISDX) ! 0 #147 AYCNTR(ISDX) = AYCNTR(ISDX) + AYVERT(I,ISDX) #149 AXCNTR(ISDX) = AXCNTR(ISDX)/NVERT ! 0 #150 AYCNTR(ISDX) = AYCNTR(ISDX)/NVERT PLUME_VOL.f90#2 SUBROUTINE PLUME_VOL(XARG,ISDX,BVARG,BVARG3,BHARG,VOLOUT) #30 INTEGER :: I , ISDX #55 ISRC = ISDX POLORG.f90#25 INTEGER :: I , ISDX #54 CALL SINDEX(SRCID,NSRC,SOID,ISDX,FIND) #60 XINT = AXS(ISDX) ! 3 #61 YINT = AYS(ISDX) PPARM.f90#2 SUBROUTINE PPARM(ISDX,TEMP) #25 INTEGER :: ISDX #31 AQS(ISDX) = TEMP(1) #32 AHS(ISDX) = TEMP(2) #33 ATS(ISDX) = TEMP(3) #34 AVS(ISDX) = TEMP(4) #35 ADS(ISDX) = TEMP(5) #40 IF ( AQS(ISDX).EQ.0.0 ) CALL ERRHDL(PATH,MODNAM,'W','320',' QS ') #42 IF ( AHS(ISDX).LT.0.0 ) THEN #45 ELSEIF ( AHS(ISDX).GT.600.0 ) THEN #50 IF ( ATS(ISDX).EQ.0.0 ) THEN ! 9 #52 ATS(ISDX) = -1.0E-5 ! 0 #53 ELSEIF ( ATS(ISDX).GT.2000.0 ) THEN #58 IF ( AVS(ISDX).LT.0.0 ) THEN ! 9 #62 AVS(ISDX) = 1.0E-5 #63 ELSEIF ( AVS(ISDX).LT.1.0E-5 ) THEN #65 AVS(ISDX) = 1.0E-5 ! 0 #66 ELSEIF ( AVS(ISDX).GT.50.0 ) THEN #71 IF ( ADS(ISDX).LT.0.0 ) THEN ! 9 #74 ELSEIF ( ADS(ISDX).LT.1.0E-5 ) THEN #76 ADS(ISDX) = 1.0E-5 ! 0 #77 ELSEIF ( ADS(ISDX).GT.20.0 ) THEN SOPARM.f90#31 INTEGER :: I , ISDX #51 CALL SINDEX(SRCID,NSRC,FIELD(3),ISDX,FIND) ! 27 #55 IF ( SOPCRD(ISDX).EQ.'Y' ) THEN ! 27 #57 CALL ERRHDL(PATH,MODNAM,'E','315',SRCID(ISDX)) ! 0 #60 SOPCRD(ISDX) = 'Y' ! 27 #71 IF ( SRCTYP(ISDX).EQ.'POINT' ) THEN ! 27 #85 CALL PPARM(ISDX,TEMP) ! 9 #86 ELSEIF ( SRCTYP(ISDX).EQ.'VOLUME' ) THEN #100 CALL VPARM(ISDX,TEMP) ! 9 #101 ELSEIF ( SRCTYP(ISDX).EQ.'AREA' ) THEN #115 CALL APARM(ISDX,TEMP) ! 9 #116 ELSEIF ( SRCTYP(ISDX).EQ.'AREAPOLY' ) THEN #130 CALL APPARM(ISDX,TEMP) ! 0 #131 ELSEIF ( SRCTYP(ISDX).EQ.'AREACIRC' ) THEN #145 CALL ACPARM(ISDX,TEMP) ! 0 #147 ELSEIF ( SRCTYP(ISDX).EQ.'OPENPIT' ) THEN #161 CALL OPARM(ISDX,TEMP) ! 0 VPARM.f90#2 SUBROUTINE VPARM(ISDX,TEMP) #25 INTEGER :: ISDX #31 AQS(ISDX) = TEMP(1) #32 AHS(ISDX) = TEMP(2) #33 ASYINI(ISDX) = TEMP(3) #34 ASZINI(ISDX) = TEMP(4) #39 IF ( AQS(ISDX).EQ.0.0 ) CALL ERRHDL(PATH,MODNAM,'W','320',' QS ') #41 IF ( AHS(ISDX).LT.0.0 ) THEN #44 ELSEIF ( AHS(ISDX).GT.100.0 ) THEN #49 IF ( ASYINI(ISDX).LT.0.0 ) THEN ! 9 #52 ELSEIF ( ASYINI(ISDX).LT.1.0E-5 ) THEN #56 ASYINI(ISDX) = 1.0E-5 #57 ELSEIF ( ASYINI(ISDX).GT.200.0 ) THEN #62 IF ( ASZINI(ISDX).LT.0.0 ) THEN ! 9 #65 ELSEIF ( ASZINI(ISDX).LT.1.0E-5 ) THEN #69 ASZINI(ISDX) = 1.0E-5 #70 ELSEIF ( ASZINI(ISDX).GT.200.0 ) THEN
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