SAMPLE SIMULATION PROGRAM




 C THIS PROGRAM CALCULATES THE NET PRESENT VALUE OF INVESTING IN A SCLAR
 C PHCTOVCLTAIC POWERED IRRIGATION SYSTEM IN YEAR J. THIS VERSION OF THE
 C PFOGRAA CONSID=ES MAKING TH;- INVESTMENT IN EACH YEAR FROM 1990 TO
 C YEAR 2000. THERE ARE TWO SETS OF VARIABLES WHICH MUST BF SPECIFIED BY
 C THE USER: I RRI';ATION SUBMO EL VARILLES AND ECONOMIC VALUATIfN
 C SUSMCDEL VAfIA.LES
 C THE IRRIGATION SJROODEL VARIABLES ARE:
 C AI=INCH=S OF VATER APPLIED PER IRRIGATION (GROSS)
 C A=TOTAL ACRES TO 9E IRRIGATLD
 C H=HCURS PER CAY IN NHICH IRRIGATION OCCURS
 C D=0AYS TO COMPLETE ONE IRRIGATION (THIS SHOULD NOT EXCEED 75
 C PERCENT OF THE IRRIGATION FREQUENCY)
 C TOH=TCTAL DYNAMIC HEAD GF THE SYSTEM (FEET)
 C TIR=NUMqER OF IRRTATTIONS PER YEAR.
 C THE ECONOMIC SU3MODEL VARIABLES ARE:
 C PRB=9ASE OF PERIOD MARKET PRICE OF ELECTRICITY (DOLLARS PER KWH)
 C ASSUMING TIME OF DAY UTILITY RATES. IN THIS PRCGPAM. PRB IS A
 C WEIGHTED AVERAGE PRICE OF THE PEAK AND OFF PEAK RATES.
 C R=DISCOUNT FACTOR
 C BR=EJY BACK CONSTRAINT OF PRGP3OTIONALITY FOR ELECTRICITY RESCLD TO
 C THE UTILITY GRID
 C TU=EXPECTED LLNG TERM ESCALTION RATE (ANNJAL) OF KWH CHANGES.
 C Z=AERAY COST OUTLOOK. Z=I IMPLIES A PESSIMISTIC VIEW. Z=2
 C IMPLIES AN OPTIMISTIC VIEW(D.O.Ee PROJECTIONS). Z=3 OEONOTES
 C A MCCERATE COST DECLINE.
 1 INTEGER Z
 2 DIMENSION PK'P(21).PR(41)*P3(41),VRS(41A)TPE(41I)
 I FV(21)TAC(21)
 C ENTER TME USER SPECIFIED VARIABLES AT THIS POINT
 3 Z=2
 4 AI=2
 5 A=40
 6 H=13
 ? D=6
 E TDH=250
 TIR=6
 10 PrSB=.006
 11 R=.06
 12 BR=I.5
 13 TL=.06
 14 DC 5 1=1.21
 15 PV(I)=0..
 16 5 CONTINUE
 C THIS PCRTICN OF THE PROGRAM DEFINES THE PHYSICALL SYSTEM DESIGN CN THE
 C BASIS CF SIMPLE ENGINEERING EQUATIONS GIVEN IN HARRISON AND CHOATE
 C (196 ) AND UTILIZED THE ORLANOO/HERNOCN SOLMET DATA TAPE.
 C THE PUMPING RATE IN GALLONS PER MINUTE IS GIVEN HY
 17 GFM=(453*(AI)*(A))/((ll)*(D))
 C THE CCNTINUCUS 3RAKE HORSEPOWER REQUIRED IS CALCULATED BY
 1S BHFC=((GPM)(TDH))/2613.6
 C THIS CALCULATICE ASSUMES A PUMP EFFICIENCY CF .75 AND MOTOR
 C EFFICIENCY OF -55.
 C THE CCTINrUOUS Kd DEMAND OF THE MOTOR IS GIVEN BY
 A1 CKW=.7457*3H0C
 C THE FCLLC.ING EQUATIONS CALCULATE THE PV ARRAY SIZE AND KWP
 C RATING. THE CALCULATION IS: AZ(4**2)=(I/0.03143)-CKWe*
 C THIS IS EQUIVALENT TO AZ=31.7L(2e*CKW.
 C (NQTE:SYSTEM EFFICIENCY IS ASSUMED TO BE 0.069).
20 AZ=Z1.7662*CKW
21 PKW=AZ/12.5
 C THE FCLLC#ING EQUATIONS CALCULATE TOTAL ARRAY OUT'UT (TAO)
 C (KWH). ELECTRICITY UTILIZED FOR INRIGATION(EDIR) AND NET ARRAY
 C GENERATION FOR RESALE TO THE UTILITY GRID (AQN)
22 TAO=113.194?*AZ
 C NOTE: ANNUAL ARRAY GENERATION IS 116.1949KWH/(M*m*.).
23 EDI= (CKdI*r()f(D)*(TIRj
24 ACN=TA3-EDIR
25 P 'IT 1OS5
26 105 FCFMAT(5X. 3HGPMF.5X.HHtHPC.9X3MHCKWg9X2HAZ.IOX,3HPKW,9Xs3TAQ)
27 WRITE(6, 1 04) O;nM HPC.CKW AZ.PKW.TAQ
2d 104 FCFuAT(FL3.2,2X.F10.2.2Xe.F 1.2.ZA.FIO.2.2X.rlo.22X.FIO.2e2XI
21 P ['RI1TI05
39 106 FCl'.AT (4X, HE IRGX, 3HAQNI
31 FRITE(6.1 7) EDIR.AON
J2 IC7 FCRV4T(F103.2.2X.FL0.2.2X)
 C THE FCLLC'I:IG STATEMENTS GErNERATE THE ARRAYS UTILIZED IN THE
 C ECONOMIC VALUATION' PORTION iF THE POGCRAYM
33 GC TO(50.60.701.Z
34 50 DC 100 J=1.2t
35 IF(J.LE.7))=2
36 IF(J.GT*7ID=O
37 rpJ-1
38 PKI (J) a( D*10750EXP( -0. 047947*M