PROGRAM FIRST C Version 1.1.1 C C THIS IS A PROGRAM TO ESTIMATE BOTH ACUTE AND CHRONIC TIER ONE, C UPPER LEVEL (HIGHER EXPOSURE) DRINKING WATER CONCENTRATIONS FOR C FOOD QUALITY PROTECTION ACT (FQPA) ASSESSMENTS. IT ESTIMATES C PESTICIDE CONCENTRATIONS IN A VULNERABLE INDEX RESERVOIR LOCATED C A HIGH USE AREA FOR THE PESTICIDE BEING ASSESSED. THE PROGRAM C CONSIDERS REDUCTIONS IN DISSOLVED CONCENTRATION DUE TO THE C PERCENTAGE OF THE WATERSHED WHICH IS CROPPED (PERCENT CROPPED AREA), C REDUCTION IN DISSOLVED PESTICIDE CONCENTRATION DUE TO ADSORPTION C OF PESTICIDE TO SOIL OR SEDIMENT, INCORPORATION, DEGRADATION IN C SOIL BEFORE WASHOFF TO A WATER BODY, DIRECT DEPOSITION OF SPRAY C DRIFT INTO THE WATER BODY, AND DEGRADATION OF THE PESTICIDE WITHIN C THE WATER BODY -IT IS DESIGNED TO MIMIC A PRZM-EXAMS SIMULATION FOR C A HIGH EXPOSURE SITE C Release notes version 1.1.1: This release was created primarily for the purpose C of standardizing the compiler versions in EFED. This version has been compiled with C Intel Fortran 10.0 along with Visual Studio 2005. No calulations-type changes were made C from version 1.1.0. In this version only grammar was corrected in the output text. C DIMENSION X0(153),X1(153),X2(153),X3(153),X4(153),X5(153) DIMENSION XV(153,6) C EQUIVALENCE (X0(1),XV(1,1)),(X1(1),XV(1,2)),(X2(1),XV(1,3)) EQUIVALENCE (X3(1),XV(1,4)),(X4(1),XV(1,5)),(X5(1),XV(1,6)) C REAL CHRONIC(401),SDCONC(401),ROCONC(401),DEGFRF(10),ADSFRR(401), 2 ADSFRS(401),PSTMSF(600),PSTMSP(600) REAL INCORP,APPRAT,KOC,SOL,KADS1,KADSUS,KADSUR,SDINIT,SDFIN, 2 ROINIT,ROFIN,METHAP,KMETP,KDEGF,APPEFF,PCTSRO,ROAREA,WBAREA, 3 CONC0,CONC4,CONC21,CONC60,CONC90,CON365,DEGHAP,KDEGP,METHAF, 4 METRAT,HYDHAP,DRIFT,KHYDP,FOTHAP,KFOTP,SUM4,SUM21,SUM60, 5 SUM90,SUM365,KD,YLOC,YLOCEN,PONDEP,KMETF,KDFRAC C REAL CHECK1,CHECK2,CHECK3,CHECK4,CHECK5,CHECK6,CHECK7 REAL SCON1,SCON2,SCON3,SCON4,SCON5,SCON6,SCON7,SCON8,SCON9,SCON10 REAL ROCO1,ROCO2,ROCO3,ROCO4,ROCO5,ROCO6,ROCO7,ROCO8,ROCO9,ROCO10 REAL DEGF1,DEGF2,DEGF3,DEGF4,DEGF5,DEGF6,DEGF7,DEGF8,DEGF9,DEGF10 REAL LOOK1,LOOK2 REAL PCA,FLADJP,FLADJA,KDADJP,KDADJA,HLADJP,HLADJA,PKFAC,AVFAC C INTEGER CODE,STORM,APFLAG,APPNUM,APSPAC,NBCPC CHARACTER*1 METHOD,AGAIN,WETTED,ADSORP,AIRFLG,GRNFLG,ORCFLG CHARACTER*4 SOLUNI CHARACTER*6 SPTYPE CHARACTER*10 COLON,COLOFF,CLEAR,CROP CHARACTER*16 CHMNAM CHARACTER*20 OUTFIL CHARACTER*22 UNITS C C DESCRIPTION OF VARIABLES C C ADSFRR FRACTION OF RUNOFF ADSORTION TO SEDIMENT COMPLETED C ADSFRS FRACTION OF SPRAY DRIFT ADSORTION TO SEDIMENT COMPLETED C ADSORP FLAG TO CHOSE BETWEEN KOC AND Kd PLUS ORGANIC CARBON AS C INPUT VALUES C APFLAG FLAG TO INDICATE SURFACE APPLICATION C APPEFF APPLICATION EFFICIENCY C APPNUM MAXIMUM NUMBER OF APPLICATION PERMITTED ON LABEL C APPRAT APPLICATION RATE C APSPAC INTERVAL IN DAYS BETWEEN PESTICIDE APPLICATIONS C CONC PEAK CONCENTRATION IN THE RESERVOIR C CON365 AVERAGE CONCENTRATION IN THE RESERVOIR DURING THE FIRST 365 DAYS C DEGFRF FRACTION OF PESTICIDE REMAINING IN FIELD AT TIME OF RAIN C DEGHAP CALCULATED OVERALL HALFLIFE IN THE RESERVOIR C FLADJA ADJUSTMENT FACTOR TO REDUCE THE ANNUAL AVERAGE PESTICIDE C CONCENTRATION BASED ON FLOW THROUGH THE RESERVOIR C FLADJP ADJUSTMENT FACTOR TO REDUCE THE PEAK ANNUAL PESTICIDE C CONCENTRATION BASED ON FLOW THROUGH THE RESERVOIR C FOC FRACTION ORGANIC CARBON C FOTHAP PHOTOLYSIS HALFLIFE IN THE RESERVOIR = NOMINAL HALF-LIFE / 124 C HYDHAP HYDROLYSIS HALFLIFE IN THE RESERVOIR C INCORP DEPTH OF INCORPORATION C INITCONC PARTIAL CALCULATION FOR UNITS CONVERSION AND RUNOFF DEPTH C KADS1 BINDING RATE FOR SPRAY DRIFT ON DAY 1 C KADSUR ULTIMATE BINDING RATE FOR RUNOFF ON DAYS TWO C THROUGH SIXTY C KADSUS ULTIMATE BINDING RATE FOR SPRAY DRIFT ON DAYS TWO C THROUGH SIXTY C KD SOIL ADSORPTION COEFFICIENT C KDADJA ADJUSTMENT FACTOR TO REDUCE THE ANNUAL AVERAGE PESTICIDE C CONCENTRATION BASED ON FLOW THROUGH THE RESERVOIR AS A C FUNCTION OF KD C KDADJP ADJUSTMENT FACTOR TO REDUCE THE PEAK ANNUAL PESTICIDE C CONCENTRATION BASED ON FLOW THROUGH THE RESERVOIR AS A C FUNCTION OF KD C KMETF METABOLIC DEGRADATION RATE IN THE FIELD C KDEGF OVERALL DEGRADATION RATE IN THE FIELD C KDEGP OVERALL DEGRADATION RATE IN THE RESERVOIR C KFOTP PHOTOLYSIS DEGRADATION RATE IN THE RESERVOIR C KHYDP HYDROLYSIS DEGRADATION RATE IN THE RESERVOIR C KMETP AEROBIC METABOLIC DEGRADATION RATE IN THE RESERVOIR C KOC ORGANIC CARBON PARTITION COEFFICIENT C KDFRAC DISSOLVED FRACTION OF THE PESTICIDE AFTER ADSORPTION C METHAF AEROBIC METABOLIC SOIL HALFLIFE C METHAP AEROBIC METABOLIC HALFLIFE IN THE RESERVOIR C METRAT METRIC APPLICATION RATE IN KILOGRAMS PER HECTARE C ORCFLG FLAG TO INDICATE TYPE OF ORCHARD AIRBLAST APPLICATION C PCA PERCENT CROPPED AREA (USED AS A DECIMAL) C PCTSRO PERCENT SURFACE RUNOFF FROM THE FIELD C PSTMSF ARRAY WITH VALUES FOR MASS OF PESTICIDE IN THE FIELD ON C THE DAY OF APPLICATION JUST AFTER APPLICATION C PSTMSP ARRAY WITH VALUES FOR MASS OF PESTICIDE IN THE RESERVOIR ON C THE DAY OF APPLICATION JUST AFTER APPLICATION C ROAREA AREA OF FIELD FROM WHICH RUNOFF OCCURS C ROFIN FINAL CONCENTRATION IN THE RESERVOIR DUE TO RUNOFF C ROINIT INITIAL CONCENTRATION IN THE RESERVOIR DUE TO RUNOFF C SDINIT INITIAL CONCENTRATION IN THE RESERVOIR DUE TO SPRAY DRIFT C SDFIN FINAL CONCENTRATION IN THE RESERVOIR DUE TO SPRAY DRIFT C SOL SOLUBILITY C SUM4 SUM OF THE FIRST 4 DAYS RESERVOIR CONCENTRATIONS C SUM21 SUM OF THE FIRST 21 DAYS RESERVOIR CONCENTRATIONS C SUM60 SUM OF THE FIRST 60 DAYS RESERVOIR CONCENTRATIONS C SUM90 SUM OF THE FIRST 90 DAYS RESERVOIR CONCENTRATIONS C SUM365 SUM OF THE FIRST 365 DAYS RESERVOIR CONCENTRATIONS C TDEGF DEGRADATION TIME IN THE FIELD C WBAREA SURFACE AREA OF THE WATERBODY (INDEX RESERVOIR) C WETTED FLAG TO INDICATE THE PESTICIDE IS WETTED-IN AND RUNOFF C OCCURS ONE THE DAY OF APPLICATION C YLOC DISTANCE FROM EDGE OF DOWNWIND SWATH TO NEAR EDGE OF C RESERVOIR IN METRIC UNITS (METERS) C YLOCEN DISTANCE FROM EDGE OF DOWNWIND SWATH TO NEAR EDGE OF C RESERVOIR IN ENGLISH UNITS (FEET) C COLON=CHAR(27)//'[34;47m' WRITE(*,2) COLON 2 FORMAT(A10) C CLEAR=CHAR(27)//'[2J' WRITE(*,2) CLEAR C WRITE(*,5) 5 FORMAT(///,3X,' FIRST ',// 2 3X,' (F)QPA (I)NDEX (R)ESERVOIR (S)CREENING (T)OOL ',/// 3 3X,' ENVIRONMENTAL FATE AND EFFECTS DIVISION ',/ 4 3X,' OFFICE OF PESTICIDE PROGRAMS ',/ 5 3X,' U.S. ENVIRONMENTAL PROTECTION AGENCY ',// 6 3X,' TIER ONE SCREENING MODEL ',/ 7 3X,' FOR DRINKING WATER PESTICIDE EXPOSURE ',// 8 3X,' VERSION 1.1.1 ',/ 9 3X,' March 25, 2008 ') C WRITE(*,10) 10 FORMAT(///,3X,'THIS IS A PROGRAM TO ESTIMATE TIER ONE ACUTE AND C 2HRONIC',/ 3 3X,'CONCENTRATION VALUES FOR PESTICIDES IN DRINKING WATER BASED', 4/ 5 3X,'UPON AN INDEX RESERVOIR LOCATED IN SHIPMAN, ILLINOIS ',// 6 3X,'THE PROGRAM IS USED TO ESTIMATE CONSERVATIVE EXPOSURE ',/ 7 3X,'VALUES FOR PESTICIDES IN ANY AREA OF THE UNITED STATES ',// 8 3X,'PLEASE ENTER A RUN NUMBER TO CONTINUE ---> ',$) READ(*,*) CODE C C OPEN FILES FOR PROGRAM OUTPUT C WRITE(*,11) 11 FORMAT(///,3X,'PLEASE SELECT AN OUTPUT FILE NAME ---> ',$) READ(*,12) OUTFIL 12 FORMAT(A20) C OPEN(UNIT=6,FILE=OUTFIL,STATUS='UNKNOWN') C C ENTER THE NAME OF THE CHEMICAL AND CROP FOR THE OUTPUT FILE C 99 WRITE(*,13) 13 FORMAT(///,3X,'PLEASE ENTER THE CHEMICAL NAME ---> ',$) READ(*,14) CHMNAM 14 FORMAT(A16) C WRITE(*,15) 15 FORMAT(///,3X,'PLEASE ENTER THE CROP NAME ---> ',$) READ(*,16) CROP 16 FORMAT(A10) C C ENTER THE USAGE INFORMATION C WRITE(*,20) 20 FORMAT(//////////////////// 2 3X,'THE PROGRAM ASSUMES THAT RAINFALL AND RESULTING RUNOFF ARE',/ 3 3X,'SUFFICIENT TO REMOVE UP TO 8% OF THE PESTICIDE ',/ 4 3X,'FROM THE PORTION OF THE 427 ACRES (172.8 HECTARES) OF ',/ 5 3X,'FIELDS IN THE RESERVOIR WATERSHED WHERE THE CROP IS GROWN',// 6/// 7 3X,'THE PORTION OF THE CHEMICAL THAT IS REMOVED FROM THE FIELDS' 8,/ 9 3X,'IN THIS WAY, FLOWS INTO THE RESERVOIR AND IS DISSOLVED IN ',/ A 3X,'THE RESERVOIR WATER ',// B 3X,'THE CHEMICAL CONCENTATION IN THE RESERVOIR REPRESENTS THE ',/ C 3X,'PART THAT IS DISSOLVED AND NOT BOUND TO FIELD SOIL OR TO ',/ D 3X,'RESERVOIR BOTTOM SEDIMENTS ',///// E 3X,'THE FOLLOWING INFORMATION SHOULD BE TAKEN FROM THE MOST ',/ F 3X,'CURRENT, ACCEPTED LABEL FOR THE USE SITE IN QUESTION ',// G 3X,'PLEASE ENTER APPLICATION RATE (IN POUNDS a.i. PER ACRE) ---> C H 3X,'NOTE: TO ENTER THE RATE IN KG/HA, PLEASE ENTER ZERO (0) ---> I',$) C READ(*,21) APPRAT 21 FORMAT(F10.0) C IF(APPRAT.LE.0.0)THEN WRITE(*,45) 45 FORMAT(//////////////////////// 2 3X,'PLEASE ENTER THE APPLICATION RATE (IN KG/HA) ---> ',$) READ(*,21) METRAT APPRAT = METRAT * 0.89218 ENDIF C WRITE(*,22) 22 FORMAT(/// 2 3X,'ENTER MAXIMUM NUMBER OF APPLICATIONS PERMITTED PER YEAR---> ' 3,$) C READ(*,32) APPNUM C 32 FORMAT(I6) C APPTOT = APPRAT * APPNUM C IF(APPNUM.GT.1)THEN WRITE(*,36) 36 FORMAT(/// 2 3X,'PLEASE ENTER INTERVAL BETWEEN APPLICATIONS (DAYS)---> ',$) C READ(*,32) APSPAC ELSE APSPAC=1 ENDIF C C CALCULATE THE TOTAL TIME OF DEGRADATION IN THE FIELD FROM THE FIRST C APPLICATION UNTIL THE STORM C TDEGF = APSPAC * (APPNUM-1) + 2 C C ENTER THE PERCENT CROPPED AREA (FRACTION OF THE WATERSHED AREA WHICH IS C PLANTED IN THE DESIGNATED CROP - AS A DECIMAL) C WRITE(*,29) 29 FORMAT(//////// 2 3X,'THE AMOUNT OF PESTICIDE IN THE WATERSHED AVAILABLE FOR ',/ 3 3X,'WASHOFF BY RAINFALL INTO THE RESERVOIR IS DEPENDENT ON ',/ 4 3X,'ON EXTENT OF THE WATERSHED ON WHICH THE CROP IS GROWN - ',/// 5 3X,'THE PROGRAM REPRESENTS THIS AREA BY A PERCENT CROPPED AREA',/ 6 3X,'(PCA) FACTOR FOR THE CROP AS FOUND IN THE FOLLOWING LIST:',// 7 3X,'CORN 0.46 CORN-SOY 0.83 CORN-SOY-WHEAT 0.83',/ 8 3X,'SOYBEANS 0.41 CORN-WHEAT 0.56 CORN-SOY-COTTON 0.83',/ 9 3X,'WHEAT 0.56 CORN-COTTON 0.46 SOY-WHEAT-COTTON 0.58',/ A 3X,'COTTON 0.20 SOY- WHEAT 0.56 SOYBEANS-COTTON 0.49',/ B 3X,'WHEAT-COT 0.20 ',// C 3X,'ALL OTHERS 0.87 ',///// C D 3X,'CITRUS 0.13 GRAPES 0.13 APPLES 0.13',/ C B 3X,'ALFALFA 0.25 SUGARCANE 0.25 STRAWBERRIES 0.25 ',/ C C 3X,'SORGHUM 0.25 SUGAR BEETS 0.25 LETTUCE 0.25 ',/ C D 3X,'PASTURE 0.25 SUNFLOWERS 0.25 SWEET CORN 0.25 ',/ C E 3X,'DRY BEANS 0.25 PEANUTS 0.25 TOMATOES 0.25 ',/ C F 3X,'OATS 0.25 PECANS ',/ C G 3X,'BARLEY 0.25 ',/ C H 3X,'HAY 0.25 ',//// I 3X,'PLEASE ENTER THE APPROPRIATE PCA FACTOR (DECIMAL) ---> ',$) C READ(*,21) PCA C IF(PCA.GT.1.0)THEN 229 WRITE(*,38) C 38 FORMAT(/// 2,' WARNING: THE PCA VALUE MUST NOT EXCEED 1.00 ') C WRITE(*,39) 39 FORMAT(///// 2 3X,'THE AMOUNT OF PESTICIDE IN THE WATERSHED AVAILABLE FOR ',/ 3 3X,'WASHOFF BY RAINFALL INTO THE RESERVOIR IS DEPENDENT ON ',/ 4 3X,'ON EXTENT OF THE WATERSHED ON WHICH THE CROP IS GROWN - ',/// 5 3X,'THE PROGRAM REPRESENTS THIS AREA BY A PERCENT CROPPED AREA',/ 6 3X,'(PCA) FACTOR FOR THE CROP AS FOUND IN THE FOLLOWING LIST:',// 7 3X,'CORN 0.46 CORN-SOY 0.83 CORN-SOY-WHEAT 0.83',/ 8 3X,'SOYBEANS 0.41 CORN-WHEAT 0.56 CORN-SOY-COTTON 0.83',/ 9 3X,'WHEAT 0.56 CORN-COTTON 0.46 SOY-WHEAT-COTTON 0.58',/ A 3X,'COTTON 0.20 SOY- WHEAT 0.56 SOYBEANS-COTTON 0.49',/ B 3X,'WHEAT-COT 0.20 ',// C 3X,'ALL OTHERS 0.87 ',///// C D 3X,'ALFALFA 0.25 SUGARCANE 0.25 STRAWBERRIES 0.25 ',/ C E 3X,'SORGHUM 0.25 SUGAR BEETS 0.25 LETTUCE 0.25 ',/ C F 3X,'PASTURE 0.25 SUNFLOWERS 0.25 SWEET CORN 0.25 ',/ C G 3X,'DRY BEANS 0.25 PEANUTS 0.25 TOMATOES 0.25 ',/ C H 3X,'OATS 0.25 PECANS ',/ C I 3X,'BARLEY 0.25 ',/ C J 3X,'HAY 0.25 ',//// K 3X,'PLEASE ENTER THE APPROPRIATE PCA FACTOR (DECIMAL) ---> ',$) C READ(*,21) PCA C IF(PCA.GT.1.0)THEN GOTO 229 ELSE WRITE(*,23) ENDIF C ELSE WRITE(*,23) ENDIF C C ENTER THE BINDING FACTOR (Kd OR Koc) C 23 FORMAT(//////////////////// 2 3X,'THE DISSOLVED PESTICIDE CONCENTRATION IN THE RESERVOIR IS ',/ 3 3X,'CALCULATED BY SUBTRACTION OF THE PORTION OF THE CHEMICAL ',/ 4 3X,'THAT IS BOUND TO FIELD SOIL, TO FIELD ORGANIC MATTER OR ',/ 5 3X,'TO RESERVOIR BOTTOM OR SUSPENDED SEDIMENT ',/// 6 3X,'THIS BOUND FRACTION IS ESTIMATED BY USE OF THE SOIL/WATER ',/ 7 3X,'EQUILIBRIUM PARTITION COEFFICIENT (Kd) OR THE ORGANIC ',/ 8 3X,'CARBON NORMALIZED SOIL/WATER EQUILIBRIUM PARTITION ',/ 9 3X,'COEFFICIENT (Koc) ',// A 3X,'SEE THE FIRST PROGRAM USERS MANUAL FOR THE APPROPRIATE Kd ',/ B 3X,'OR Koc VALUE TO USE ',/// C 3X,'TO USE A Kd VALUE, PLEASE ENTER IT HERE. TO USE A Koc ',/ D 3X,'VALUE, PLEASE ENTER ZERO (0) ---> ',$) C READ(*,21) KD KOC = KD * 86.207 ADSORP = 'B' C IF(KD.LE.0.0)THEN WRITE(*,37) C 37 FORMAT(//////////////////////// 2 3X,'PLEASE ENTER THE APPROPRIATE Koc VALUE ---> ',$) READ(*,21) KOC ADSORP = 'A' KD = 0.0116 * KOC ENDIF C C SET THE DAY ON WHICH THE RAINFALL/RUNOFF EVENT OCCURS C STORM = 2 C C ENTER THE DEGREDATION HALF-LIFE IN THE AGRICULTURAL FIELD (USUALLY THE C AEROBIC SOIL METABOLISM HALF-LIFE) C WRITE(*,24) 24 FORMAT(/////////////////////// 2 3X,'THE DISSOLVED PESTICIDE CONCENTRATION IS ALSO REDUCED BY ',/ 3 3X,'DEGRADATION IN THE FIELD PRIOR TO A RAINFALL/RUNOFF EVENT',// 4 3X,'THIS PROGRAM ASSUMES DEGRADATION BY AEROBIC METABOLISM ',/ 5 3X,'BETWEEN APPLICATIONS AS WELL AS FOR TWO DAYS AFTER THE LAST', 6/ 7 3X,'OF THE APPLICATIONS ',// 8 3X,'(IF STABLE TO AEROBIC SOIL METABOLISM OR IF DATA IS ',/ 9 3X,'UNAVAILABLE, PLEASE ENTER ZERO (0)) ',/// A 3X,'PLEASE ENTER SOIL AEROBIC METABOLIC HALFLIFE (IN DAYS) ---> ' B,$) C READ(*,21)METHAF C C ALLOW THE USER TO SELECT FOR A CHEMICAL THAT IS WETTED IN AT THE TIME C OF APPLICATION C WRITE(*,25) 25 FORMAT(//////////////////// 2 3X,'SOME PESTICIDE LABELS REQUIRE THAT THE PESTICIDE BE ',/ 3 3X,'ACTIVATED BY "WETTING-IN" EITHER THROUGH RAINFALL OR ',/ 4 3X,'THROUGH IRRIGATION ',/// 5 3X,'IN THIS CASE, RUNOFF TO THE RESERVOIR IS ASSUMED TO ',/ 6 3X,'OCCUR IMMEDIATELY RATHER THAN AFTER TWO DAYS ',// 7 3x,'IS THIS PESTICIDE TO BE WETTED-IN ? (Y or N) ---> ',$) C READ(*,26) WETTED 26 FORMAT(A1) C IF(WETTED.EQ.'Y'.OR.WETTED.EQ.'y') STORM = 0 C C ENTER THE SPRAY DRIFT AND BUFFER INFORMATION C WRITE(*,27) 27 FORMAT(//////////////////////// 2 3X,'THE DISSOLVED PESTICIDE CONCENTRATION MAY BE INCREASED BY ',/ 3 3X,'BY DEPOSITION OF SPRAY DRIFT EITHER DIRECTLY INTO THE RESERVO 4IR',/ 5 3X,'OR INTO THE STREAMS THAT FLOW INTO THE RESERVOIR ',///// 63X,'PLEASE ENTER A, B, C, or D TO SELECT METHOD OF APPLICATION:',/ 7 3X,'A: AERIAL SPRAY (DRIFT=16.0%; APPLICATION EFFICIENCY=95%) ',/ 8 3X,'B: GROUND SPRAY (DRIFT=6.4%; APPLICATION EFFICIENCY=99%) ',/ 9 3X,'C: ORCHARD-VINYARD AIRBLAST SPRAY (DRIFT=6.3%;APP EFF=99%)',/ A 3X,'D: GRANULAR (DRIFT=0.0%; APPL. EFFICIENCY=100%) ---> ',$) C READ(*,30)METHOD 30 FORMAT(A1) C IF(METHOD.EQ.'A'.OR.METHOD.EQ.'a')THEN C DRIFT = 0.16 APPEFF = 0.95 SPTYPE = 'AERIAL' C ELSEIF(METHOD.EQ.'B'.OR.METHOD.EQ.'b')THEN C DRIFT = 0.064 APPEFF = 0.99 SPTYPE = 'GROUND' C IF 6 ELSEIF(METHOD.EQ.'C'.OR.METHOD.EQ.'c')THEN C DRIFT = 0.063 APPEFF = 0.99 SPTYPE = 'ABLAST' C ELSEIF(METHOD.EQ.'D'.OR.METHOD.EQ.'d')THEN C YLOCEN = 0.0 YLOC = 0.0 DRIFT = 0.0 APPEFF = 1.0 SPTYPE = 'GRANUL' ENDIF C C FOR AERIAL OR GROUND APPLICATION CALL FOR BUFFER WIDTH C IF(METHOD.EQ.'A'.OR.METHOD.EQ.'a'.OR.METHOD.EQ.'B'.OR.METHOD.EQ.'b 2'.OR.METHOD.EQ.'C'.OR.METHOD.EQ.'c')THEN C C YLOC = YLOCEN * 0.3048 YLOCEN = 0.0 YLOC = 0.0 ENDIF C C ENTER DEPTH OF INCORPORATION FOR GRANULAR PESTICIDES C (This is a correction from Version 1.0 which mistakenly permitted C incorporation for granular application only) C IF(METHOD.EQ.'B'.OR.METHOD.EQ.'b'.OR.METHOD.EQ.'D'.OR.METHOD.EQ.'d 2')THEN C WRITE(*,28) 28 FORMAT(//////////////////// 2 3X,'THE DISSOLVED PESTICIDE CONCENTRATION MAY BE REDUCED BY',/ 3 3X,'INCORPORATION OF THE PESTICIDE AT THE TIME OF APPLICATION.',/ 4 3x,'IN THIS CASE, THE FOLLOWING DEPTHS ARE SUGGESTED: ',// 5 3X,' METHOD OF APPLICATION INCORPORATION DEPTH (IN) ',/ 6 3X,' _________________________ ________________________ ',/ 7 3X,' ',/ 8 3X,' BROADCAST 0.0 ',/ 9 3X,' DISKED IN AFTER BROADCAST 4.0 ',/ A 3X,'CHISEL PLOWED AFTER BROADCAST 6.0 ',/ B 3X,' SURFACE BANDED 0.0 ',/ C 3X,' BANDED - INCORPORATED 1.2 ',/ D 3X,' T - BANDED 1.5 ',/ E 3X,' IN FURROW 2.0 ',/ F 3X,' AERIAL or AIRBLAST SPRAY 0.0 ',/ G 3X,' GROUND SPRAY DEPENDS ON METHOD ',// H 3X,'PLEASE ENTER DEPTH OF INCORPORATION (IN INCHES) ---> ',$) C READ(*,21) INCORP APFLAG = 0 ENDIF C C ENTER THE CHEMICAL SOLUBILITY AS AN UPPER LIMIT ON THE DISSOLVED C CONCENTRATION C WRITE(*,31) 31 FORMAT(/////////////// 2 3X,'THE SOLUBILITY OF A PESTICIDE IN WATER IS A MEASURE OF THE',/ 3 3X,'MAXIMUM AMOUNT OF THE CHEMICAL THAT CAN BE DISSOLVED. ',///// 4 3X,'THUS, THE AQUEOUS PESTICIDE CONCENTRATION IN A WATER BODY',/ 5 3X,'CANNOT EXCEED THE SOLUBILITY OF THE CHEMICAL.',///// 6 3X,'PLEASE ENTER THE SOLUBILITY (IN PPM) ---> ',$) C READ(*,21) SOL C C CALCULATE THE DISSOLVED CONCENTRATION IN THE RESERVOIR AS A FUNCTION OF C THE CHEMICAL Kd BASED ON THE RELATIVE PERCENT DISSOLVED FROM A PRZM-EXAMS C SIMULATION - THE INTERPOLATION IS BASED ON THE FOLLOWING FORMULA: C C IF(KD.LE.A) C KDFRAC = 1.0 C ELSEIF(KD.LE.B) C KDFRAC = C + (D-C) * (B-Kd) / (B-A) C IF(KD.LE.5.00e-3)THEN KDFRAC = 1.0 ELSEIF(KD.LE.1.00e-2.AND.KD.GT.5.00e-3)THEN KDFRAC = 0.9991715 + (1.0 - 0.9991715) * (1.00e-2 - KD) / 2 (1.00e-2 - 5.00e-3) ELSEIF(KD.LE.5.00e-2.AND.KD.GT.1.00e-2)THEN KDFRAC = 0.9933720 + (0.9991715 - 0.9933720) * (5.00e-2 - KD) / 2 (5.00e-2 - 1.00e-2) ELSEIF(KD.LE.1.00e-1.AND.KD.GT.5.00e-2)THEN KDFRAC = 0.9859155 + (0.9933720 - 0.9859155) * (1.00e-1 - KD) / 2 (1.00e-1 - 5.00e-2) ELSEIF(KD.LE.3.00e-1.AND.KD.GT.1.00e-1)THEN KDFRAC = 0.9569180 + (0.9859155 - 0.9569180) * (3.00e-1 - KD) / 2 (3.00e-1 - 1.00e-1) ELSEIF(KD.LE.5.00e-1.AND.KD.GT.3.00e-1)THEN KDFRAC = 0.9295775 + (0.9569180 - 0.9295775) * (5.00e-1 - KD) / 2 (5.00e-1 - 3.00e-1) ELSEIF(KD.LE.7.50e-1.AND.KD.GT.5.00e-1)THEN KDFRAC = 0.8980944 + (0.9295775 - 0.8980944) * (7.50e-1 - KD) / 2 (7.50e-1 - 5.00e-1) ELSEIF(KD.LE.1.00e00.AND.KD.GT.7.50e-1)THEN KDFRAC = 0.8682684 + (0.8980944 - 0.8682684) * (1.00e00 - KD) / 2 (1.00e00 - 7.50e-1) ELSEIF(KD.LE.1.25e00.AND.KD.GT.1.00e00)THEN KDFRAC = 0.8409279 + (0.8682684 - 0.8409279) * (1.25e00 - KD) / 2 (1.25e00 - 1.00e00) ELSEIF(KD.LE.1.50e00.AND.KD.GT.1.25e00)THEN KDFRAC = 0.8147307 + (0.8409279 - 0.8147307) * (1.50e00 - KD) / 2 (1.50e00 - 1.25e00) ELSEIF(KD.LE.1.75e00.AND.KD.GT.1.50e00)THEN KDFRAC = 0.7904060 + (0.8147307 - 0.7904060) * (1.75e00 - KD) / 2 (1.75e00 - 1.50e00) ELSEIF(KD.LE.2.00e00.AND.KD.GT.1.75e00)THEN KDFRAC = 0.7675973 + (0.7904060 - 0.7675973) * (2.00e00 - KD) / 2 (2.00e00 - 1.75e00) ELSEIF(KD.LE.2.25e00.AND.KD.GT.2.00e00)THEN KDFRAC = 0.7461475 + (0.7675973 - 0.7461475) * (2.25e00 - KD) / 2 (2.25e00 - 2.00e00) ELSEIF(KD.LE.2.50e00.AND.KD.GT.2.25e00)THEN KDFRAC = 0.7260066 + (0.7461475 - 0.7260066) * (2.50e00 - KD) / 2 (2.50e00 - 2.25e00) ELSEIF(KD.LE.2.75e00.AND.KD.GT.2.50e00)THEN KDFRAC = 0.7070340 + (0.7260066 - 0.7070340) * (2.75e00 - KD) / 2 (2.75e00 - 2.50e00) ELSEIF(KD.LE.3.00e00.AND.KD.GT.2.75e00)THEN KDFRAC = 0.6891881 + (0.7070340 - 0.6891881) * (3.00e00 - KD) / 2 (3.00e00 - 2.75e00) ELSEIF(KD.LE.3.50e00.AND.KD.GT.3.00e00)THEN KDFRAC = 0.6562883 + (0.6891881 - 0.6562883) * (3.50e00 - KD) / 2 (3.50e00 - 3.00e00) ELSEIF(KD.LE.4.00e00.AND.KD.GT.3.50e00)THEN KDFRAC = 0.6269097 + (0.6562883 - 0.6269097) * (4.00e00 - KD) / 2 (4.00e00 - 3.50e00) ELSEIF(KD.LE.4.50e00.AND.KD.GT.4.00e00)THEN KDFRAC = 0.6004060 + (0.6269097 - 0.6004060) * (4.50e00 - KD) / 2 (4.50e00 - 4.00e00) ELSEIF(KD.LE.5.00e00.AND.KD.GT.4.50e00)THEN KDFRAC = 0.5765700 + (0.6004060 - 0.5765700) * (5.00e00 - KD) / 2 (5.00e00 - 4.50e00) ELSEIF(KD.LE.5.50e00.AND.KD.GT.5.00e00)THEN KDFRAC = 0.5548384 + (0.5765700 - 0.5548384) * (5.50e00 - KD) / 2 (5.50e00 - 5.00e00) ELSEIF(KD.LE.6.00e00.AND.KD.GT.5.50e00)THEN KDFRAC = 0.5352196 + (0.5548384 - 0.5352196) * (6.00e00 - KD) / 2 (6.00e00 - 5.50e00) ELSEIF(KD.LE.7.00e00.AND.KD.GT.6.00e00)THEN KDFRAC = 0.5007954 + (0.5352196 - 0.5007954) * (7.00e00 - KD) / 2 (7.00e00 - 6.00e00) ELSEIF(KD.LE.8.00e00.AND.KD.GT.7.00e00)THEN KDFRAC = 0.4717896 + (0.5007954 - 0.4717896) * (8.00e00 - KD) / 2 (8.00e00 - 7.00e00) ELSEIF(KD.LE.9.00e00.AND.KD.GT.8.00e00)THEN KDFRAC = 0.4471002 + (0.4717896 - 0.4471002) * (9.00e00 - KD) / 2 (9.00e00 - 8.00e00) ELSEIF(KD.LE.1.00e01.AND.KD.GT.9.00e00)THEN KDFRAC = 0.4257415 + (0.4471002 - 0.4257415) * (1.00e01 - KD) / 2 (1.00e01 - 9.00e00) ELSEIF(KD.LE.1.25e01.AND.KD.GT.1.00e01)THEN KDFRAC = 0.3837614 + (0.4257415 - 0.3837614) * (1.25e01 - KD) / 2 (1.25e01 - 1.00e01) ELSEIF(KD.LE.1.50e01.AND.KD.GT.1.25e01)THEN KDFRAC = 0.3528086 + (0.3837614 - 0.3528086) * (1.50e01 - KD) / 2 (1.50e01 - 1.25e01) ELSEIF(KD.LE.1.75e01.AND.KD.GT.1.50e01)THEN KDFRAC = 0.3292129 + (0.3528086 - 0.3292129) * (1.75e01 - KD) / 2 (1.75e01 - 1.50e01) ELSEIF(KD.LE.2.00e01.AND.KD.GT.1.75e01)THEN KDFRAC = 0.3107208 + (0.3292129 - 0.3107208) * (2.00e01 - KD) / 2 (2.00e01 - 1.75e01) ELSEIF(KD.LE.2.50e01.AND.KD.GT.2.00e01)THEN KDFRAC = 0.2834880 + (0.3107208 - 0.2834880) * (2.50e01 - KD) / 2 (2.50e01 - 2.00e01) ELSEIF(KD.LE.3.00e01.AND.KD.GT.2.50e01)THEN KDFRAC = 0.2646396 + (0.2834880 - 0.2646396) * (3.00e01 - KD) / 2 (3.00e01 - 2.50e01) ELSEIF(KD.LE.4.00e01.AND.KD.GT.3.00e01)THEN KDFRAC = 0.2400580 + (0.2646396 - 0.2400580) * (4.00e01 - KD) / 2 (4.00e01 - 3.00e01) ELSEIF(KD.LE.5.00e01.AND.KD.GT.4.00e01)THEN KDFRAC = 0.2249793 + (0.2400580 - 0.2249793) * (5.00e01 - KD) / 2 (5.00e01 - 4.00e01) ELSEIF(KD.LE.1.00e02.AND.KD.GT.5.00e01)THEN KDFRAC = 0.1939188 + (0.2249793 - 0.1939188) * (1.00e02 - KD) / 2 (1.00e02 - 5.00e01) ELSEIF(KD.LE.5.00e02.AND.KD.GT.1.00e02)THEN KDFRAC = 0.1788732 + (0.1939188 - 0.1788732) * (5.00e02 - KD) / 2 (5.00e02 - 1.00e02) ELSEIF(KD.LE.1.00e03.AND.KD.GT.5.00e02)THEN KDFRAC = 0.1615742 + (0.1788732 - 0.1615742) * (1.00e03 - KD) / 2 (1.00e03 - 5.00e02) ELSEIF(KD.LE.5.00e03.AND.KD.GT.1.00e03)THEN KDFRAC = 0.1425352 + (0.1615742 - 0.1425352) * (5.00e03 - KD) / 2 (5.00e03 - 1.00e03) ELSEIF(KD.LE.1.00e04.AND.KD.GT.5.00e03)THEN KDFRAC = 0.1258409 + (0.1425352 - 0.1258409) * (1.00e04 - KD) / 2 (1.00e04 - 5.00e03) ELSEIF(KD.LE.2.00e04.AND.KD.GT.1.00e04)THEN KDFRAC = 0.1021458 + (0.1258409 - 0.1021458) * (2.00e04 - KD) / 2 (2.00e04 - 1.00e04) ELSEIF(KD.LE.3.00e04.AND.KD.GT.2.00e04)THEN KDFRAC = 0.0859983 + (0.1021458 - 0.0859983) * (3.00e04 - KD) / 2 (3.00e04 - 2.00e04) ELSEIF(KD.LE.5.00e04.AND.KD.GT.3.00e04)THEN KDFRAC = 0.0653521 + (0.0859983 - 0.0653521) * (5.00e04 - KD) / 2 (5.00e04 - 3.00e04) ELSEIF(KD.LE.1.00e05.AND.KD.GT.5.00e04)THEN KDFRAC = 0.0408318 + (0.0653521 - 0.0408318) * (1.00e05 - KD) / 2 (1.00e05 - 5.00e04) ELSEIF(KD.LE.5.00e05.AND.KD.GT.1.00e05)THEN KDFRAC = 0.0102055 + (0.0408318 - 0.0102055) * (5.00e05 - KD) / 2 (5.00e05 - 1.00e05) ELSEIF(KD.LE.1.00e06.AND.KD.GT.5.00e05)THEN KDFRAC = 0.0052672 + (0.0102055 - 0.0052672) * (1.00e06 - KD) / 2 (1.00e06 - 5.00e05) ELSEIF(KD.GT.1.00e06)THEN KDFRAC = 0.001 ENDIF C C VALUES FOR Kd OVER 1,000,000 WILL NOT RUN FOR LONG TERM PRZM-EXAMS C SIMULATIONS AND ARE SET TO 0.001 C C ELSEIF(KD.LE.5.00e06.AND.KD.GT.1.00e06)THEN C KDFRAC = 0.0010813 + (0.0052672 - 0.0010813) * (5.00e06 - KD) / C 2 (5.00e06 - 1.00e06) C ELSEIF(KD.LE.1.00e07.AND.KD.GT.5.00e06)THEN C KDFRAC = 0.0005425 + (0.0010813 - 0.0005425) * (1.00e07 - KD) / C 2 (1.00e07 - 5.00e06) C ELSEIF(KD.LE.5.00e07.AND.KD.GT.1.00e07)THEN C KDFRAC = 0.0001088 + (0.0005425 - 0.0001088) * (5.00e07 - KD) / C 2 (5.00e07 - 1.00e07) C ELSEIF(KD.LE.1.00e08.AND.KD.GT.5.00e07)THEN C KDFRAC = 0.0000544 + (0.0001088 - 0.0000544) * (1.00e08 - KD) / C 2 (1.00e08 - 5.00e07) C ELSEIF(KD.LE.5.00e08.AND.KD.GT.1.00e08)THEN C KDFRAC = 0.0000109 + (0.0000544 - 0.0000109) * (5.00e08 - KD) / C 2 (5.00e08 - 1.00e08) C ELSEIF(KD.LE.1.00e09.AND.KD.GT.5.00e08)THEN C KDFRAC = 0.0000054 + (0.0000109 - 0.0000054) * (1.00e09 - KD) / C 2 (1.00e09 - 5.00e08) C ELSEIF(KD.GT.1.00e09)THEN C KDFRAC = 0.0000054 C ENDIF C C CALCULATE THE AMOUNT OF REDUCTION DUE TO DEGREDATION IN THE FIELD C IF(METHAF.LE.0.0) THEN KMETF = 0.0 ELSE KMETF = LOG(2.0) / METHAF ENDIF C C LIMIT THE AMOUNT OF REDUCTION DUE TO INCORPORATION C IF(INCORP.LE.0.0001) APFLAG = 1 IF(INCORP.LE.1.0) INCORP = 1.0 IF(INCORP.GE.6.0) INCORP = 6.0 C C SET THE PERCENT RUNOFF FROM THE FIELD AND THE AREA OF THE FIELD C PCTSRO = 0.08 ROAREA = 172.8 WBAREA = 5.28 C C CALCULATE THE FRACTION OF PESTICIDE REMAINING FOR TWO ANNUAL TURNOVERS C AS A FUNCTION OF Kd (THIS IS AN EMPIRICAL FUNCTION BASED ON PRZM-EXAMS C SIMULATION FOR THIS AMOUNT OF FLOW THROUGH THE RESERVOIR) C IF(KD.LE.1.0)THEN KDADJA = 0.705 KDADJP = 0.974 ELSEIF(KD.LE.10)THEN KDADJA = 0.705 - KD * ((0.705-0.664)/10.0) KDADJP = 0.974 - KD * ((0.974-0.930)/10.0) ELSEIF(KD.LE.100)THEN KDADJA = 0.664 - KD * ((0.664-0.732)/100.0) KDADJP = 0.930 - KD * ((0.930-0.960)/100.0) ELSEIF(KD.LE.1000)THEN KDADJA = 0.732 - KD * ((0.732-0.918)/1000.0) KDADJP = 0.960 - KD * ((0.960-0.992)/1000.0) ELSEIF(KD.LE.10000)THEN KDADJA = 0.918 - KD * ((0.918-0.987)/10000.0) KDADJP = 0.992 - KD * ((0.992-0.999)/10000.0) ELSE KDADJA = 1.0 KDADJP = 1.0 ENDIF C C SET THE FLOW ADJUSTMENT FACTOR EQUAL TO THE Kd-FLOW ADJUSTMENT FACTOR C FLADJA = KDADJA FLADJP = KDADJP C CALCULATE THE FRACTION OF THE APPLICATION RATE REMAINING IN THE FIELD C ON EACH OF THE SEVEN DAYS AFTER APPLICATION - SHIFT ALLOWS CALCULATION C OF A ZERO DAY VALUE C SHIFT = 0 C DO 100 I = 1,8 SHIFT = I-1 DEGFRF(I) = EXP(-KMETF*SHIFT) DEGF1 = DEGFRF(1) DEGF2 = DEGFRF(2) DEGF3 = DEGFRF(3) DEGF4 = DEGFRF(4) DEGF5 = DEGFRF(5) DEGF6 = DEGFRF(6) DEGF7 = DEGFRF(7) DEGF8 = DEGFRF(8) DEGF9 = DEGFRF(9) DEGF10 = DEGFRF(10) 100 CONTINUE C C INCREASE THE VALUE OF 'STORM' TO FIT THIS SHIFTED VALUE C STORM = STORM + 1 C C CALCULATE CHRONIC EEC,s C C ASK FOR RESERVOIR DEGRADATION HALFLIVES C WRITE(*,33) 33 FORMAT(//////////////////////// 2 3X,'THE CHRONIC DRINKING WATER PESTICIDE CONCENTRATION VALUE ',/ 3 3X,'IS ESTIMATED USING A DEGRADATION RATE THAT IS CALCULATED BY' 4,/ 5 3X,'SUMMING THE INDIVIDUAL AQUATIC DEGRADATION RATES (THE AEROBIC 6',/ 7 3X,'AQUATIC METABOLIC RATE IS ASSUMED TO INCLUDE HYDROLYSIS)',/// 8 3X,'ENTER ANY OR ALL OF THE FOLLOWING THAT ARE AVAILABLE: ',// 9 3X,'(PLEASE ENTER ZERO (0) FOR ANY THAT ARE STABLE OR ',/ A 3X,'FOR WHICH VALUES ARE UNAVAILABLE) ',// B 3X,'ENTER AEROBIC AQUATIC METABOLIC HALFLIFE IN DAYS (IF ',/ C 3X,'UNAVAILABLE, THE RECOMMENDED EFED DEFAULT IS 2 TIMES',/ R 3X,'THE AEROBIC SOIL INPUT VALUE) ---> D ',$) READ(*,21) METHAP C IF(METHAP.LE.0.0) THEN KMETP = 0.0 METHAP = 0.00 C WRITE(*,34) 34 FORMAT(///,3x,'ENTER pH 7 HYDROLYSIS HALFLIFE (DAYS) ---> ',$) READ(*,21) HYDHAP C IF(HYDHAP.LE.0.0) THEN KHYDP = 0.0 ELSE KHYDP = LOG(2.0) / HYDHAP ENDIF ELSE KMETP = LOG(2.0) / METHAP ENDIF C WRITE(*,35) 35 FORMAT(///,3x,'ENTER PHOTOLYSIS HALFLIFE (DAYS) ---> ',$) READ(*,21) FOTHAP C IF(FOTHAP.LE.0.0) THEN KFOTP = 0.0 ELSE KFOTP = (LOG(2.0) / FOTHAP) / 124 ENDIF C C CALCULATE THE OVERALL DEGRADATION RATE IN THE RESERVOIR C KDEGP = KHYDP + KFOTP + KMETP IF(KDEGP.LE.0.0) THEN DEGHAP = 1000000 ELSE DEGHAP = LOG(2.0) / KDEGP ENDIF C C ZERO OUT PESTICIDE MASS LEFT FROM PREVIOUS RUNS C DO 80 I = 1 , 600 PSTMSF(I) = 0.0 80 CONTINUE C DO 81 I = 1 , 600 PSTMSP(I) = 0.0 81 CONTINUE C C SET THE INITIAL PESTICIDE IN THE FIELD AT THE TIME OF APPLICATION TO C THE APPLICATION RATE C I = 1 PSTMSF(1) = APPRAT C C SET THE INITIAL PESTICIDE IN THE RESERVOIR AT THE TIME OF APPLICATION TO C THE APPLICATION RATE C I = 1 PSTMSP(1) = APPRAT C C SET OVERALL FIELD DEGRADATION RATE TO SOIL METABOLISM RATE C IF(METHAF.LE.0.0) THEN KDEGF = 0.0 ELSE KDEGF = LOG(2.0) / METHAF ENDIF C C CALCULATE THE AMOUNT OF PESTICIDE IN THE FIELD ON EACH APPLICATION DATE C DO 85 I = 2 , APPNUM PSTMSF(I) = PSTMSF(I-1) * EXP(-KDEGF*APSPAC) + APPRAT 85 CONTINUE C C CALCULATE THE AMOUNT OF PESTICIDE IN THE RESERVOIR ON EACH APPLICATION DATE C DO 86 I = 2 , APPNUM PSTMSP(I) = PSTMSP(I-1) * EXP(-KDEGP*APSPAC) + APPRAT 86 CONTINUE C C CALCULATE THE DISSOLVED CONCENTRATION DUE TO SPRAY DRIFT C C CALCULATE RATE OF ADSORPTION TO SEDIMENT ON DAY 1 AND THE INITIAL C DISSOLVED CONCENTRATION RESULTING FROM SPRAY DRIFT ALONE C KADS1 = (9.2529+1.751*KOC) / (1.341E6+KOC) C SDINIT=(1.12085*PSTMSP(APPNUM)*DRIFT*WBAREA*EXP(-KADS1))/144.146 C C CALCULATE THE ULTIMATE RATE OF ADSORPTION OF PESTICIDE IN SPRAY C DRIFT TO SEDIMENT C KADSUS = (9366.5+12.4*KOC) / (655000+KOC) C C CALCULATE THE DAILY FRACTION OF PESTICIDE IN SPRAY DRIFT WHICH IS C ADSORBED C DO 105 I = 1,375 ADSFRS(I) = EXP(-KADSUS*I) 105 CONTINUE C C CALCULATE THE FINAL DISSOLVED CONCENTRATION DUE TO SPRAY DRIFT LOAD * C THE VALUE 0.13875 IS THE RATIO BETWEEN THE POND AND INDEX RESERVOIR C (20.0/144.146) C SDFIN = 1.12085 * PSTMSP(APPNUM) * DRIFT * WBAREA * 0.13875 * 2 ((37.0388+9E-6*KOC) / (750+KOC)) C C CALCULATE THE DAILY DISSOLVED CONCENTRATION IN THE RESERVOIR DUE TO SPRAY C DRIFT C DO 110 I = 1,375 SDCONC(I) = (SDFIN + ADSFRS(I) * (SDINIT-SDFIN)) * EXP(-KDEGP*I) 110 CONTINUE C C CALCULATE THE INITIAL DISSOLVED CONCENTRATION IN THE RESERVOIR DUE TO THE C RUNOFF EVENT AND THE RESULTING DAILY DISSOLVED CONCENTRATION VALUES C ROINIT = (1.12085 * PSTMSF(APPNUM) * APPEFF * ROAREA * PCTSRO * 2 KDFRAC * DEGFRF(STORM) / INCORP) / 144.146 C C CALCULATE THE ULTIMATE RATE OF ADSORPTION OF PESTICIDE IN RUNOFF TO C SEDIMENT C KADSUR = (5742.9+7.6*KOC) / (405000+KOC) C C CALCULATE THE DAILY FRACTION OF PESTICIDE IN RUNOFF WHICH IS ADSORBED C DO 205 I = 1,375 ADSFRR(I) = EXP(-KADSUR*I) 205 CONTINUE C C CALCULATE THE FINAL DISSOLVED CONCENTRATION DUE TO RUNOFF LOAD - 6.262 C IS RUNOFF LOAD IN PRZM-EXAMS SIMULATION C THE VALUE 7.2073 IS THE RATIO BETWEEN THE RESERVOIR AND THE FARM POND C (144.146/20.0)) C ROFIN = 1.12085 * PSTMSF(APPNUM) * APPEFF * PCTSRO * ROAREA * 2 DEGFRF(STORM) * ((157.845+4.3E-6*KOC**1.215) / 3 (510+KOC**1.215)) / INCORP / (6.262 * 7.2073) C C CALCULATE THE DAILY DISSOLVED CONCENTRATION IN THE RESERVOIR DUE TO RUNOFF C DO 210 I=1,375 ROCONC(I) = (ROFIN+ADSFRR(I) * (ROINIT-ROFIN)) * EXP(-KDEGP*I) 210 CONTINUE C C CALCULATE TOTAL DAILY DISSOLVED CONCENTRATION VALUES BY SUMMING THE C DISSOLVED CONCENTRATION DUE TO RUNOFF AND THE CONCENTRATION DUE TO SPRAY C DRIFT (IF APPLICABLE) C IF(METHOD.EQ.'A'.OR.METHOD.EQ.'a'.OR.METHOD.EQ.'B'.OR. 2 METHOD.EQ.'b'.OR.METHOD.EQ.'C'.OR.METHOD.EQ.'c')THEN CONC0 = ROINIT + SDCONC(STORM) DO 215 I = 1,375 CHRONIC(I) = ROCONC(I) + SDCONC(I+STORM-1) 215 CONTINUE C ELSEIF(METHOD.EQ.'D'.OR.METHOD.EQ.'d')THEN CONC0 = ROINIT DO 220 I = 1,375 CHRONIC(I) = ROCONC(I) 220 CONTINUE ENDIF C C ZERO OUT ANY PESTICIDE REMAINING FROM PREVIOUS RUNS C SUM4 = 0.0 SUM21 = 0.0 SUM60 = 0.0 SUM90 = 0.0 SUM365 = 0.0 C C CALCULATE THE 96 HOUR AVERAGE CONC C DO 225 I = 1,3 SUM4 = SUM4 + CHRONIC(I) 225 CONTINUE CONC4 = (CONC0 + SUM4) / 4 C C CALCULATE THE 21 DAY AVERAGE CONC C DO 230 I = 1,20 SUM21 = SUM21 + CHRONIC(I) 230 CONTINUE CONC21 = (CONC0 + SUM21) / 21 C C CALCULATE THE 60 DAY AVERAGE CONC C DO 240 I=1,59 SUM60 = SUM60 + CHRONIC(I) 240 CONTINUE CONC60 = (CONC0 + SUM60) / 60 C C CALCULATE THE 90 DAY AVERAGE CONC C DO 250 I=1,89 SUM90 = SUM90 + CHRONIC(I) 250 CONTINUE CONC90 = (CONC0 + SUM90) / 90 C C CALCULATE THE 365 DAY AVERAGE CONC C DO 260 I=1,364 SUM365 = SUM365 + CHRONIC(I) 260 CONTINUE CON365 = (CONC0 + SUM365) / 365 C C ADJUST CONCENTRATION VALUES FOR PCA AND FLOW THROUGH THE RESERVOIR C CONC0 = PCA * CONC0 CON365 = PCA * CON365 C CONC0 = FLADJP * CONC0 CON365 = FLADJA * CON365 C C RECALIBRATE THE ANNUAL AVERAGE CONCENTRATION VALUES TO EXCEED THE C LOUISIANA SWEET POTATO PRZM-EXAMS SCENARIO C IF(KD.GE.10.0) THEN AVFAC = 2.0 * LOG10(KD) CON365 = CON365 * AVFAC ENDIF C C BE SURE THE FINAL DISSOLVED CONCENTRATION DOES NOT EXCEED THE C SOLUBILITY OF THE CHEMICAL C IF(CONC0.GE.SOL) CONC0 = SOL IF(CONC4.GE.SOL) CONC4 = SOL IF(CONC21.GE.SOL) CONC21 = SOL IF(CONC60.GE.SOL) CONC60 = SOL IF(CONC90.GE.SOL) CONC90 = SOL IF(CON365.GE.SOL) CON365 = SOL C C CALCULATE THE UNITS OF THE ANSWER C IF(CONC0.GE.1.0)THEN UNITS = 'MILLIGRAMS/LITER (PPM)' ENDIF C IF(CONC0.LT.1.0.AND.CONC0.GE.0.001)THEN CONC0 = CONC0 * 1000 CONC4 = CONC4 * 1000 CONC21 = CONC21 * 1000 CONC60 = CONC60 * 1000 CONC90 = CONC90 * 1000 CON365 = CON365 * 1000 UNITS = 'MICROGRAMS/LITER (PPB)' ENDIF C IF(CONC0.LT.0.001)THEN CONC0 = CONC0 * 1000000 CONC4 = CONC4 * 1000000 CONC21 = CONC21 * 1000000 CONC60 = CONC60 * 1000000 CONC90 = CONC90 * 1000000 CON365 = CON365 * 1000000 UNITS = 'NANOGRAMS/LITER (PPTr)' ENDIF C IF(SOL.GE.1.0)THEN SOLUNI = 'PPM' ENDIF C IF(SOL.LT.1.0.AND.SOL.GE.0.001) THEN SOL = SOL * 1000 SOLUNI = 'PPB' ENDIF C IF(SOL.LT.0.001) THEN SOL = SOL * 1000000 SOLUNI = 'PPTr' ENDIF C C WRITE OUTPUT TO THE SCREEN AND TO THE OUTPUT FILE C IF(ADSORP.EQ.'A'.OR.ADSORP.EQ.'a')THEN C WRITE(*,50)CODE,CHMNAM,CROP,SOLUNI WRITE(6,350)CODE,CHMNAM,CROP,SOLUNI 50 FORMAT(//,3X,'RUN No.',I4,' FOR ',A16,' ON ',A12,' * INPUT VALU 2ES * ',/ 3 3X,'------------------------------------------------------------- 4-------',/ 5 3X,' RATE (#/AC) No.APPS & SOIL SOLUBIL APPL TYPE %CROPPED 6 INCORP',/ 7 3X,' ONE(MULT) INTERVAL Koc (',A4,') (%DRIFT) AREA 8 (IN)',/ 9 3X,'------------------------------------------------------------- A-------') C 350 FORMAT(///,3X,'RUN No.',I4,' FOR ',A16,' ON ',A12,' * INPUT VAL 2UES * ',/ 3 3X,'------------------------------------------------------------- 4-------',/ 5 3X,' RATE (#/AC) No.APPS & SOIL SOLUBIL APPL TYPE %CROPPED 6 INCORP',/ 7 3X,' ONE(MULT) INTERVAL Koc (',A4,') (%DRIFT) AREA 8 (IN)',/ 9 3X,'------------------------------------------------------------- A-------') C IF(APFLAG.EQ.1) INCORP = 0.0 C WRITE(*,51)APPRAT,'(',PSTMSF(APPNUM),')',APPNUM,APSPAC,KOC,SOL, 2 SPTYPE,'(',DRIFT*100,')',PCA*100.0,INCORP WRITE(6,51)APPRAT,'(',PSTMSF(APPNUM),')',APPNUM,APSPAC,KOC,SOL, 2 SPTYPE,'(',DRIFT*100,')',PCA*100.0,INCORP 51 FORMAT(F7.3,A1,F7.3,A1,1X,I3,1X,I3,F10.1,F7.1,3X,A6,A1,F4.1,A1, 2 1X,F5.1,2X,F4.1) C ELSEIF(ADSORP.EQ.'B'.OR.ADSORP.EQ.'b')THEN C WRITE(*,57)CODE,CHMNAM,CROP,SOLUNI WRITE(6,257)CODE,CHMNAM,CROP,SOLUNI C 57 FORMAT(//,3X,'RUN No.',I4,' FOR ',A16,' ON ',A12,' * INPUT VALU 2ES * ',/ 3 3X,'------------------------------------------------------------- 4-------',/ 5 3X,'RATE (#/AC) No.APPS & SOIL SOLUBIL APPL TYPE %CROPPED 6 INCORP',/ 7 3X,' ONE(MULT) INTERVAL Kd (',A4,') (%DRIFT) AREA 8 (IN)',/ 9 3X,'------------------------------------------------------------- A-------') C 257 FORMAT(///,3X,'RUN No.',I4,' FOR ',A16,' ON ',A12,' * INPUT VAL 2UES * ',/ 3 3X,'------------------------------------------------------------- 4-------',/ 5 3X,'RATE (#/AC) No.APPS & SOIL SOLUBIL APPL TYPE %CROPPED 6 INCORP',/ 7 3X,' ONE(MULT) INTERVAL Kd (',A4,') (%DRIFT) AREA 8 (IN)',/ 9 3X,'------------------------------------------------------------- A-------') C IF(APFLAG.EQ.1) INCORP = 0.0 C WRITE(*,58)APPRAT,'(',PSTMSF(APPNUM),')',APPNUM,APSPAC,KD,SOL, 2 SPTYPE,'(',DRIFT*100,')',PCA*100.0,INCORP WRITE(6,58)APPRAT,'(',PSTMSF(APPNUM),')',APPNUM,APSPAC,KD,SOL, 2 SPTYPE,'(',DRIFT*100,')',PCA*100.0,INCORP 58 FORMAT(F7.3,A1,F7.3,A1,1X,I3,1X,I3,F10.1,F7.1,3X,A6,A1,F4.1,A1, 2 1X,F5.1,4X,F4.1) C ENDIF C IF(DEGHAP.GT.999999) DEGHAP = 0.0 IF(HYDHAP.LE.0.0)THEN WRITE(6,60) WRITE(*,60) WRITE(6,61)METHAF,STORM-1,'N/A ',FOTHAP,'-',FOTHAP*124,METHAP, 2 DEGHAP WRITE(*,61)METHAF,STORM-1,'N/A ',FOTHAP,'-',FOTHAP*124,METHAP, 2 DEGHAP ELSE WRITE(6,60) WRITE(*,60) WRITE(6,62)METHAF,STORM-1,HYDHAP,FOTHAP,'-',FOTHAP*124,METHAP, 2 DEGHAP WRITE(*,62)METHAF,STORM-1,HYDHAP,FOTHAP,'-',FOTHAP*124,METHAP, 2 DEGHAP ENDIF C 60 FORMAT(//,1X,' FIELD AND RESERVOIR HALFLIFE VALUES (DAYS) ',/ 1 3X,'----------------------------------------------------------- 2---------',/ 3 3X,'METABOLIC DAYS UNTIL HYDROLYSIS PHOTOLYSIS METABOLIC 4 COMBINED',/ 5 3X,' (FIELD) RAIN/RUNOFF (RESERVOIR) (RES.-EFF) (RESER.) 6 (RESER.) ',/ 7 3X,'----------------------------------------------------------- 8---------') C 61 FORMAT(3X,F7.2,7X,I2,7X,A7,3X,F6.2,A1,F8.2,3X,F6.2,3X,F7.2) C 62 FORMAT(3X,F7.2,7X,I2,8X,F7.2,2X,F6.2,A1,F8.2,2X,F6.2,4X,F7.2) C WRITE(6,70)UNITS WRITE(*,70)UNITS C WRITE(6,72)CONC0,CON365 WRITE(*,72)CONC0,CON365 C 70 FORMAT(//,3X,'UNTREATED WATER CONC (',A22,') Ver 1.1.1 MAR 26, 20 208',/ 3 3X,'------------------------------------------------------------- 4-------',/ 5 3X,' PEAK DAY (ACUTE) ANNUAL AVERAGE (CHRONIC) ',/ 6 3X,' CONCENTRATION CONCENTRATION ',/ 7 3X,'------------------------------------------------------------- 8-------') C C 71 FORMAT(5X,F7.2,5X,F7.2,7X,F7.2,7X,F7.2,7X,F7.2) C 72 FORMAT(12X,F7.3,20X,F7.3) C WRITE(*,53) 53 FORMAT(/,3X,'DO YOU WANT TO DO ANOTHER RUN (Y OR N) ---> ',$) READ(*,30) AGAIN C IF(AGAIN.EQ.'Y'.OR.AGAIN.EQ.'y')THEN WRITE(*,54) 54 FORMAT(///3X,'PLEASE ENTER A NEW RUN NUMBER ---> ',$) READ(*,*) CODE C C SET ALL INPUT VALUES BACK TO ZERO FOR SUBSEQUENT RUNS C APPRAT = 0 APPNUM = 0 APPTOT = 0 KD = 0 KOC = 0 KDFRAC = 0 SOL = 0 METHAF = 0 METHAP = 0 HYDHAP = 0 FOTHAP = 0 DEGHAP = 0 INCORP = 0 PCTSRO = 0 ROAREA = 0 ROINIT = 0 ROFIN = 0 SDINIT = 0 SDFIN = 0 SHIFT = 0 STORM = 0 KADS1 = 0 KADSUR = 0 KADSUS = 0 KDEGP = 0 KMETF = 0 KHYDP = 0 KMETP = 0 KFOTP = 0 C CONC0 = 0 CONC4 = 0 CONC21 = 0 CONC60 = 0 CONC90 = 0 CON365 = 0 SUM4 = 0 SUM21 = 0 SUM60 = 0 SUM90 = 0 SUM365 = 0 C DO 1000 K=1,375 ROCONC(K) = 0 SDCONC(K) = 0 CHRONIC(K) = 0 1000 CONTINUE C DO 1001 K = 1,100 ADSFRR(K) = 0 ADSFRS(K) = 0 1001 CONTINUE C DO 1002 K = 1,8 DEGFRF(K) = 0 1002 CONTINUE C I = 0 K = 0 C GOTO 99 ENDIF C COLOFF = CHAR(27)//'[0m' WRITE(*,2) COLOFF WRITE(*,2) CLEAR C STOP END