(*<pre>*)
(*
          Section 1: Header

This is a Mathematica notebook file. This file is created for Hua
on Monday, 14 March 1994.

The main model is caps (version 1).

The sub-models are:
     fuelCost (version 1)
     fuelConsM (version 1)
     fuelConsA (version 1).

The scenario is caps(1).

The report can be obtained by running Mathmatica with this file.

*)

(*
          Section 2   Table of Contents

     Section 1: Header
     Section 2: Table of Contents
     Section 3: Model Descriptions
     Section 4: Variable Descriptions
     Section 5: Variable Declarations
     Section 6: Model Declarations
     Section 7: Model Evaluations
     Section 8: Post-Evaluation Analyses
     Section 9: End of Input
*)

(*          Section 3: Model Descriptions          *)

(*
 Model: caps (version 1).
caps: A model for determining fuel costs t enable Comparison 
of Alternative Propulsion Systems.

Model version 1 is created by Hua Hua
 on February 11, 1994

Model source: 
Clark Pritchett

Model is maintained by:
Hua Hua OPIM Department University of Pennsylvania PA 19104-6366 
(215)-898-5715
*)

caps::usage = " caps: A model for determining fuel costs t enable 
Comparison of Alternative Propulsion Systems. "

(* The sub-models are following: *)

(*
 Model: fuelCost (version 1).
fuelCoast: A model determins fuel costs as a function of time.

Model version 1 is created by Hua Hua
 on February 11, 1994

Model source: 
Clark Pritchett

Model is maintained by:
Hua Hua OPIM Department University of Pennsylvania PA 19104-6366 
(215)-898-5715
*)

fuelCost::usage = " fuelCoast: A model determins fuel 
costs as a function of time. "


(*
 Model: fuelConsM (version 1).
fuelConsM: A model determins fuel consumption for Mains.

Model version 1 is created by Hua Hua
 on February 11, 1994

Model source: 
Clark Pritchett

Model is maintained by:
Hua Hua OPIM Department University of Pennsylvania PA 19104-6366 
(215)-898-5715
*)

fuelConsM::usage = " fuelConsM: A model determins fuel 
consumption for Mains. "


(*
 Model: fuelConsA (version 1).
fuelConsA: A model determins fuel consumption for Auxiliaries.

Model version 1 is created by Hua Hua
 on February 11, 1994

Model source: 
Clark Pritchett

Model is maintained by:
Hua Hua OPIM Department University of Pennsylvania PA 19104-6366 
(215)-898-5715
*)

fuelConsA::usage = " fuelConsA: A model determins fuel 
consumption for Auxiliaries. "



(*           Section 4: Variable Descriptions           *)

fracL::usage = " fracL is defined as a local variable for the global 
variable fracTimeload . fracTimeload[i] a variable: 
The fraction of time at (electrical) load i, 
where i is an index; i = 1, 2, ... "

fracS::usage = " fracS is defined as a local variable for the global 
variable fracTimespeed . fracTimespeed[i] a variable: 
The fraction of time a vessel spends at speed i, 
where i is an index; i = 1, 2, ... "

fcA::usage = " fcA is defined as a local variable for the global 
variable fuelConsumedAuxiliaries . total fuel consumed by Auxiliaries 
"

fcM::usage = " fcM is defined as a local variable for the global 
variable fuelConsumedMains . total fuel consumed by Mains "

c0::usage = " c0 is defined as a local variable for the global 
variable fuelCostNow . c0, avariable: The present cost of fuel 
in dollars per gallon. "

g1::usage = " g1 is defined as a local variable for the global 
variable gamma1 . gamma1, a variable: 
pounds per gallon of diesel fuel. "

g2::usage = " g2 is defined as a local variable for the global 
variable gamma2 . gamma2, a variable: kilowatts per horse power. "

hpr::usage = " hpr is defined as a local variable for the global 
variable hpRequiredspeed . hpRequiredspeed[i] a variable: 
The horsepower required for speed i, where i is 
an index; i = 1, 2, ... "

kwr::usage = " kwr is defined as a local variable for the global 
variable kwRequiredspeed . kwRequiredspeed[i] a variable: 
The electrical power required at load i, 
where i is an index; i = 1, 2, ... "

sfcL::usage = " sfcL is defined as a local variable for the global 
variable specFuelConsDieselLoad . specFuelConsDieselLoad[i] 
a variable: The specific load of a diesel engine (lb / hp * hr) 
at vessel speed i, where i is an index; i = 1, 2, ... "

sfcS::usage = " sfcS is defined as a local variable for the global 
variable specFuelConsDieselspeed . specFuelConsDieselspeed[i] 
a variable: The specific fuel consumption of a diesel 
engine (lb / hp * hr) at vessel speed i, where i is an index; 
i = 1, 2, ... "

ufc::usage = " ufc is defined as a local variable for the global 
variable unitFuelCost . cost of fuel in the year t "


(*         Section 5: Variable Declarations         *)

(*
This is a scenario for model caps
Scenario source: 
Hemant K. Bhargava
*)


c0=0.93
fracL[1]=0.6
fracL[2]=0.3
fracL[3]=0.1
fracS[1]=0.35
fracS[2]=0.4
fracS[3]=0.2
fracS[4]=0.05
g1=0.75
g2=0.7457
hpr[1]=100
hpr[2]=1200
hpr[3]=2100
hpr[4]=6500
kwr[1]=20
kwr[2]=60
kwr[3]=125
sfcL[1]=0.42
sfcL[2]=0.41
sfcL[3]=0.4
sfcS[1]=0.38
sfcS[2]=0.37
sfcS[3]=0.37
sfcS[4]=0.34
t=0

(*         Section 6: Model Declarations         *)

fc[t_] := ufc[t] * (fcM + fcA)

ufc[t_] := c0 * (1 + 17 / 60 * t - 1 / 120 * t ^ 2)

fcM := Sum[fracS[i] * hpr[i] * sfcS[i], { i, 1, 4 } ] / g1

fcA := Sum[fracL[i] * kwr[i] * sfcL[i], { i, 1, 3 } ] / g1 * g2

(*         Section 7: Model Evaluations         *)

Print["Model Output:"]


Print["According to the model(s) and the scenario,
   the variable fc is equal to  ", fc]

(*         Section 8: Post-Evaluation Analyses         *)



(*         Section 9: End of Input         *)

(*</pre>*)