doc-calculator 0.2__py3-none-any.whl

doc_calculator/__init__.py

@@ -0,0 +1,2 @@
+from .core import DirectOperatingCost
+from .gemseo_discipline import GemseoDirectOperatingCost

doc_calculator/core/DOC_Calculator.py

@@ -0,0 +1,424 @@
+from .utils.params import Params
+from .utils.util_functions import _assign_input
+from typing import Dict, Tuple
+import math
+
+class DirectOperatingCost(object):
+
+    # Constants
+    MACH_NUMBER_FACTOR = 1.0  # 1.0 = Assumed subsonic cruise
+    FLIGHT_TIME_OFFSET = 0.25
+
+    def __init__(self, aircraft:dict, params:Params=Params()) -> None:
+        """
+        ### Description
+        This code enables to evaluate Direct and Total Operating Costs
+        for Short/Medium haul airliners and regional aircraft (jet and
+        propeller driven ).
+
+        ### Input Dict Keys
+        List of variables composing the aircraft dict (case insensitive):
+        - adp    (USD M)   Aircraft Delivery Price
+
+        - mtow  (Tonns)    Max Take-off Weight
+        - pld   (Tonns)    Payload
+        - mew   (Tonns)    Manufacturer Empty Weight
+        - bengw (Tonns)    Bare engine weight
+        - enpri (USD M)    Engine Price
+        - en               Thermal Engines number
+        - crewc            Attendants number
+        - crewtech         Number of Tech. crew members (pilots)
+        - bt      (HR)     Sector Block Time
+        - bf      (KG)     Sector Block Fuel
+        - sector  (NM)     Sector assumed for DOCs evaluation
+
+        - ieng             Engine maintenance cost flag (1=Calculate, 2=Assigned in eoc parameter)
+        - eoc    (USD/BHR) Engine overhaul cost (Necessary for ieng=2)
+        - shp    (HP)      Thermal Engine Shaft Horse Power (Necessary for ieng=1)
+
+        - afspare (fraction of adp)   Airframe spares in [0.0  1.0]
+        - enspare (fraction of enpri) Spare engines and spares in [0.0  1.0]
+        - dyrs    (YRS)               Depreciation period
+        - rval                        Aircraft Residual value in [0.0  1.0]
+        - rinsh   (fraction of adp)   Insurance rate in [0.0  1.0]
+        - crtechr  (USD/BHR)          Tech. crew members (pilots) hourly tot.cost
+        - crcabhr  (USD/BHR/Att )     Cab. crew member hourly cost
+        - labor_rate       (USD/MHR)  Maintenance Labour rate
+        - fuelpri  (USD/US GAL)       Fuel Price
+        - ioc_fact                    Coeff. for IOC evaluation in [0.0  Inf)
+        - util     (BHR/Annum)        Aircraft Annual Utilisation
+        - lifespan                    Lifespan of the aircraft (years)
+        - l_app (EPNdB)               Certified noise level at the approach measure point
+        - l_lat (EPNdB)               Certified noise level at the lateral measure point
+        - l_flyov (EPNdB)             Certified noise level at the fly-over measure point
+
+        - cnox     (USD)           Unit rate for NOX (generally referred to nitrogen oxides) (Emission tariff that depends on the airport)
+        - nox_value (Kg)           Emission value of NOX. Equivalent of nitrogen oxide exhausted by an aircraft,
+                                   in kilogram, in the "Landing and Take-Off Cycle, LTO
+        - cco     (USD)            Unit rate for CO (generally referred to nitrogen oxides) (Emission tariff that depends on the airport)
+        - co_value (Kg)            Emission value of CO. Equivalent of carbon monoxide exhausted by an aircraft, 
+                                   in kilogram, in the "Landing and Take-Off Cycle, LTO
+        - co2_value      (kg)      Mass of Emitted CO2, in kilograms
+        - prico2   (USD/kg)        co2_value price for unit of emitted co2_value mass
+
+        - n_bat                    Number of batteries pack     
+        - n_fc                     Number of fuel cells
+        - n_repbat                 Number of battery(ies) replacement during the aircraft lifespan
+        - batprice  (USD)          Battery(ies) price [USD]
+        - rvbat     (USD)          Battery(ies) residual value (at the end of its own life) [USD]
+        - lrbat     (USD/MMH)      Maintenance labor rate for battery(ies)
+        - tlbat     (MMH)          Maintenance man-hour for battery(ies) [hours]
+        - f_bat                    Maintenance frequency for battery(ies) - as the reciprocal of the number of flights before the next check
+        - n_repfc                  Number of fuel cell(s) replacement during the aircraft lifespan
+        - fcprice   (USD)          Fuel Cell price [USD]
+        - rvfc      (USD)          Fuel Cell(s) residual value (at the end of its own life) [USD]
+        - lrfc      (USD/MMH)      Maintenance labor rate for fuel cell(s)
+        - tlfc      (MMH)          Maintenance man-hour for fuel cell(s) [hours]
+        - f_fc                     Maintenance frequency for fuel cell(s) - as the reciprocal of the number of flights before the next check
+        - n_reppe                  Number of power electronics replacement during the aircraft lifespan
+        - peprice   (USD)          Power electronics price [USD]
+        - rvpe      (USD)          Power electronics residual value (at the end of its own life) [USD]
+        - lrpe      (USD/MMH)      Maintenance labor rate for power electronics
+        - tlpe      (MMH)          Maintenance man-hour for power electronics [hours]
+        - f_pe                     Maintenance frequency for power electronics - as the reciprocal of the number of flights before the next check
+
+        - n_em                     Number of electric machines
+        - emprice   (USD)          Price of a single electric machine
+        - lrem      (USD/MMH)      Maintenance labor rate for electric machine(s)
+        - speml     (USD)          Spare Parts Cost Line Maintenance electric machine
+        - spemb     (USD)          Spare Parts Cost Base Maintenance electric machine (If not known -> Assumption: spemb = 9.5*speml)
+        - tleml     (MMH)          Maintenance man-hour for line maintenance electric machine(s) [hours]
+        - tlemb     (MMH)          Maintenance man-hour for base maintenance electric machine(s) [hours]
+        - f_eml      (times/YR)    Maintenance frequency for electric machine(s) Line Maint.
+        - f_emb      (times/BH)    Maintenance frequency for electric machine(s) Base Maint.
+
+        - enerpri      (USD/kWh)   Electricity price 
+        - ener_req  (kWh)          Electricity Requirement (from battery)
+        - h2_pri         (USD/kg)  H2 price
+        - h2_req (kg)              H2 requirements
+        
+        """
+        self.aircraft = _assign_input(input=aircraft)
+        self._params = params
+
+        return None
+    
+    def calculate_doc(self) -> Dict[str, float]:
+
+        bt = self.aircraft["bt"]
+        financial = self._calculate_financial_cost()
+        operating = self._calculate_cash_operating_cost()
+
+        doc_total = sum(financial.values()) + sum(operating.values())
+
+        return {
+            **financial,
+            **operating,
+            "DOC [USD/BHR]": doc_total,
+            "DOC [USD/flight]": bt * doc_total,
+        }
+        
+    def calculate_ioc(self) -> Dict[str, float]:
+    
+        ioc_factor = self.aircraft["ioc_fact"]
+        bt = self.aircraft["bt"]
+        operating_cost = self._calculate_cash_operating_cost()
+        ioc_bhr = ioc_factor * sum(operating_cost.values())
+
+        return {
+            "IOC [USD/BHR]": ioc_bhr,
+            "IOC [USD/flight]": bt * ioc_bhr,
+        }
+    
+    def _calculate_cash_operating_cost(self) -> Dict[str, float]:
+
+        fuel                            = self._calculate_fuel_cost()
+        electric_energy                 = self._calculate_electric_energy_price()
+        h2                              = self._calculate_h2_price()
+        cockpit_crew                    = self._calculate_cockpit_crew_cost()
+        cabin_crew                      = self._calculate_cabin_crew_cost()
+        landing_fees                    = self._calculate_landing_fees()
+        nav_charges                     = self._calculate_navigation_charges()
+        ground_charges                  = self._calculate_ground_handling_charges()
+        noise_charges                   = self._calculate_noise_charges()
+        airframe_maintenance            = self._calculate_airframe_maintenance_cost()
+        thermal_engine_maintenance      = self._calculate_thermal_engine_maintenance_cost()
+        nox_emission_charges            = self._calculate_nox_emission_charges()   
+        co_emission_charges             = self._calculate_co_emission_charges()   
+        co2_emission_charges            = self._calculate_co2_emission_charges()
+
+        electric_machine_maint_line, electric_machine_maint_base    = self._calculate_electric_machine_maintenance_cost()
+        battery_maint_line, battery_maint_base                      = self._calculate_battery_maintenance_cost()
+        fuelcell_maint_line, fuelcell_maint_base                    = self._calculate_fuel_cell_maintenance_cost()
+        power_elec_maint_line, power_elec_maint_base                = self._calculate_power_electronic_maintenance_cost() 
+
+        return {
+            "FUEL [USD/BHR]": fuel,
+            "ELECTRYCITY [USD/BHR]": electric_energy,
+            "H2 [USD/BHR]": h2,
+            "COCKPIT CREW [USD/BHR]": cockpit_crew,
+            "CABIN CREW [USD/BHR]": cabin_crew,
+            "LANDING FEES [USD/BHR]": landing_fees,
+            "NAVIGATION CHARGES [USD/BHR]": nav_charges,
+            "GROUND HANDLING [USD/BHR]": ground_charges,
+            "NOISE CHARGES [USD/BHR]": noise_charges,
+            "NOX EMISSION CHARGES [USD/BHR]": nox_emission_charges,
+            "CO EMISSION CHARGES [USD/BHR]": co_emission_charges,
+            "CO2 EMISSION CHARGES [USD/BHR]": co2_emission_charges,
+            "AIRFRANE MAINTENANCE [USD/BHR]": airframe_maintenance,
+            "THERM. ENG. MAINTENANCE [USD/BH]": thermal_engine_maintenance,
+            "ELECTRIC MACHINE LINE MAINT. [USD/BH]": electric_machine_maint_line,
+            "ELECTRIC MACHINE BASE MAINT. [USD/BH]": electric_machine_maint_base,
+            "BATTERY LINE MAINT. [USD/BH]": battery_maint_line,
+            "BATTERY BASE MAINT. [USD/BH]": battery_maint_base,
+            "FUEL CELL LINE MAINT. [USD/BH]": fuelcell_maint_line,
+            "FUEL CELL BASE MAINT. [USD/BH]": fuelcell_maint_base,
+            "POWER ELECTR. LINE MAINT. [USD/BH]": power_elec_maint_line, 
+            "POWER ELECTR. BASE MAINT. [USD/BH]": power_elec_maint_base     
+            }
+    
+    def _calculate_financial_cost(self) -> Dict[str, float]:
+        return {
+            "INSURANCE [USD/BHR]": self._calculate_insurance_cost(),
+            "DEPRECIATION [USD/BHR]": self._calculate_depreciation(),
+            "INTEREST [USD/BHR]": self._calculate_interest()
+            }
+    
+    def _calculate_co2_emission_charges(self) -> float:
+        co2_value = self.aircraft["co2_value"]    
+        prico2 = self.aircraft["prico2"]
+        bt     = self.aircraft["bt"] 
+
+        return (1.0-self._params.AEC)*co2_value*prico2/bt
+    
+    def _calculate_battery_maintenance_cost(self) -> Tuple[float, float]:
+        n_bat     = self.aircraft["n_bat"]
+        n_repbat  = self.aircraft["n_repbat"]
+        batprice = self.aircraft["batprice"]
+        rvbat    = self.aircraft["rvbat"]
+        lrbat    = self.aircraft["lrbat"]
+        tlbat    = self.aircraft["tlbat"]
+        f_bat     = self.aircraft["f_bat"]
+        lifespan = self.aircraft["lifespan"]
+        util     = self.aircraft["util"]
+
+        battery_maint_line = n_bat*lrbat*tlbat*f_bat
+        battery_maint_base = n_bat*(n_repbat*(batprice-rvbat))/(lifespan*util) # replace
+        return battery_maint_line, battery_maint_base
+    
+    def _calculate_fuel_cell_maintenance_cost(self) -> Tuple[float, float]:
+        n_fc      = self.aircraft["n_fc"]
+        n_repfc   = self.aircraft["n_repfc"]
+        fcprice  = self.aircraft["fcprice"]
+        rvfc     = self.aircraft["rvfc"]
+        lrfc     = self.aircraft["lrfc"]
+        tlfc     = self.aircraft["tlfc"]
+        f_fc      = self.aircraft["f_fc"]
+        lifespan = self.aircraft["lifespan"]
+        util     = self.aircraft["util"]
+
+        fuel_cell_maint_line = n_fc*lrfc*tlfc*f_fc
+        fuel_cell_maint_base = n_fc*(n_repfc*(fcprice-rvfc))/(lifespan*util) # replace
+        return fuel_cell_maint_line, fuel_cell_maint_base
+    
+    def _calculate_power_electronic_maintenance_cost(self) -> Tuple[float, float]:
+        n_reppe   = self.aircraft["n_reppe"]
+        peprice  = self.aircraft["peprice"]
+        rvpe     = self.aircraft["rvpe"]
+        lrpe     = self.aircraft["lrpe"]
+        tlpe     = self.aircraft["tlpe"]
+        f_pe      = self.aircraft["f_pe"]
+        lifespan = self.aircraft["lifespan"]
+        util     = self.aircraft["util"]
+
+        power_electronic_maint_line = lrpe*tlpe*f_pe
+        power_electronic_maint_base = (n_reppe*(peprice-rvpe))/(lifespan*util) # replace
+        return power_electronic_maint_line, power_electronic_maint_base
+
+    def _calculate_electric_machine_maintenance_cost(self) -> Tuple[float, float]:
+        n_em      = self.aircraft["n_em"]     
+        speml    = self.aircraft["speml"]       # spare parts cost line maintenance
+        spemb    = self.aircraft["spemb"]       # spare parts cost base maintenance
+        lrem     = self.aircraft["lrem"]
+        tleml    = self.aircraft["tleml"]       # maintenance man hour line maint.
+        tlemb    = self.aircraft["tlemb"]       # maintenance man hour base maint.
+        f_eml     = self.aircraft["f_eml"]        # maintenance frequency line maint
+        f_emb     = self.aircraft["f_emb"]        # maintenance frequency base maint
+        lifespan = self.aircraft["lifespan"]
+        util     = self.aircraft["util"]
+        
+        electric_machine_maint_line = n_em*(speml + lrem*tleml)*lifespan*f_eml/util
+        electric_machine_maint_base = n_em*(spemb + lrem*tlemb)*f_emb*0.80
+        return electric_machine_maint_line, electric_machine_maint_base
+    
+    def _calculate_airframe_maintenance_cost(self) -> float:
+        n_bat     = self.aircraft["n_bat"]         
+        n_em      = self.aircraft["n_em"]
+        n_fc      = self.aircraft["n_fc"]
+        batprice  = self.aircraft["batprice"]
+        fcprice   = self.aircraft["fcprice"]
+        emprice   = self.aircraft["emprice"]
+        bt        = self.aircraft["bt"]
+        adp       = self.aircraft["adp"]
+        en        = self.aircraft["en"]
+        enpri     = self.aircraft["enpri"]
+        mew       = self.aircraft["mew"]
+        bengw     = self.aircraft["bengw"]
+        labor_rate = self.aircraft["labor_rate"]
+
+        FT  = bt - self.FLIGHT_TIME_OFFSET
+        AFW = mew-(bengw*en)
+
+        C_A_FH = 3.08*(adp - en*enpri - n_bat*batprice/1.0e6 - n_em*emprice/1.0e6 - n_fc*fcprice/1.0e6)
+        C_A_FC = 6.24*(adp - en*enpri - n_bat*batprice/1.0e6 - n_em*emprice/1.0e6 - n_fc*fcprice/1.0e6)
+        K_A_FC = 0.05*AFW*2.2 + 6 - 630.0/(AFW*2.2 + 120.0)
+        K_A_FH = 0.59*K_A_FC
+        
+        airframe_labor_cost    = (K_A_FH*FT + K_A_FC)*labor_rate*math.sqrt(self.MACH_NUMBER_FACTOR)/bt
+        airframe_material_cost = (C_A_FH*FT + C_A_FC)/bt
+
+        return airframe_material_cost + airframe_labor_cost
+
+    def _calculate_thermal_engine_maintenance_cost(self) -> float:
+        ieng = self.aircraft["ieng"]
+        en   = self.aircraft["en"]
+        
+        if ieng == 1:
+            bt    = self.aircraft["bt"]
+            labor_rate    = self.aircraft["labor_rate"]
+            shp   = self.aircraft["shp"]
+            enpri = self.aircraft["enpri"]
+
+            FT       = bt - self.FLIGHT_TIME_OFFSET
+            K_ICE_FC = (0.3 + 0.03*shp/1000.0)*en
+            K_ICE_FH = (0.65 + 0.03*shp/1000.0)*en
+            C_ICE_FC = 2.0*en*enpri*10.0
+            C_ICE_FH = 2.5*en*enpri*10.0
+
+            thermal_engine_labor_cost    = (K_ICE_FH*FT + K_ICE_FC)*labor_rate/bt
+            thermal_engine_material_cost = (C_ICE_FH*FT + C_ICE_FC)/bt
+
+            thermal_engine_maintenance_cost = thermal_engine_material_cost + thermal_engine_labor_cost
+
+        elif ieng == 2:
+            eoc = self.aircraft["eoc"]
+
+            thermal_engine_maintenance_cost = eoc*en
+            
+        else:
+            raise ValueError(f"ieng Value {ieng} not valid")
+        
+        return thermal_engine_maintenance_cost
+    
+    def _calculate_nox_emission_charges(self) -> float:
+        cnox        = self.aircraft["cnox"]
+        nox_value = self.aircraft["nox_value"]
+        bt           = self.aircraft["bt"]
+
+        return (cnox*nox_value)/bt
+
+    def _calculate_co_emission_charges(self) -> float:
+        cco         = self.aircraft["cco"]
+        co_value  = self.aircraft["co_value"]
+        bt           = self.aircraft["bt"]
+
+        return (cco*co_value)/bt
+    
+    def _calculate_noise_charges(self) -> float:
+        ta      = self._params.TA
+        td      = self._params.TD
+        cnoise  = self._params.CNOISE
+
+        l_app   = self.aircraft["l_app"]
+        l_flyov = self.aircraft["l_flyov"]
+        l_lat   = self.aircraft["l_lat"]
+        bt      = self.aircraft["bt"]
+
+        DELTAA = (l_app-ta)/10.0
+        DELTAD = (((l_flyov+l_lat)/2.0)-td)/10.0
+        
+        return (cnoise*(10.0**(DELTAA)+10.0**(DELTAD)))/bt
+
+    def _calculate_ground_handling_charges(self) -> float:
+        pld   = self.aircraft["pld"]
+        bt    = self.aircraft["bt"]
+
+        return (self._params.HTONN*pld)/bt
+    
+    def _calculate_navigation_charges(self) -> float:
+        mtow   = self.aircraft["mtow"]
+        bt     = self.aircraft["bt"]
+        sector = self.aircraft["sector"]
+
+        return (self._params.ENR*sector*1.853/100.0)*math.sqrt(mtow/50.0)/bt
+    
+    def _calculate_landing_fees(self) -> float:
+        mtow = self.aircraft["mtow"]
+        bt   = self.aircraft["bt"]
+
+        return (self._params.LANDINGUR*mtow)/bt
+
+    def _calculate_cabin_crew_cost(self) -> float:
+        crcabhr = self.aircraft["crcabhr"]
+        crewc   = self.aircraft["crewc"]
+
+        return crcabhr*crewc
+    
+    def _calculate_cockpit_crew_cost(self) -> float:
+        crtechr  = self.aircraft["crtechr"]
+        crewtech = self.aircraft["crewtech"]
+
+        return crtechr*crewtech
+    
+    def _calculate_h2_price(self) -> None:
+        h2_pri = self.aircraft["h2_pri"]
+        h2_req = self.aircraft["h2_req"]
+        bt     = self.aircraft["bt"]
+
+        return h2_pri*h2_req/bt
+    
+    def _calculate_electric_energy_price(self) -> None:
+        enerpri  = self.aircraft["enerpri"]
+        ener_req = self.aircraft["ener_req"]
+        bt       = self.aircraft["bt"]
+
+        return ener_req*enerpri/bt
+    
+    def _calculate_fuel_cost(self) -> float:
+        fuelpri = self.aircraft["fuelpri"]
+        bf      = self.aircraft["bf"]
+        bt      = self.aircraft["bt"]
+
+        return (0.328*fuelpri*bf)/bt
+    
+    def _calculate_insurance_cost(self) -> float:
+        rinsh = self.aircraft["rinsh"]
+        adp   = self.aircraft["adp"]*1.0e6
+        util  = self.aircraft["util"]
+
+        return (rinsh*adp)/util
+    
+    def _calculate_investment(self) -> float:
+        adp     = self.aircraft["adp"]*1.0e6
+        afspare = self.aircraft["afspare"]
+        enpri   = self.aircraft["enpri"]
+        en      = self.aircraft["en"]
+        enspare = self.aircraft["enspare"]
+
+        return adp+(afspare*(adp-enpri*en))+(enspare*enpri*en)
+    
+    def _calculate_interest(self) -> float:
+        util = self.aircraft["util"]
+
+        INVEST = self._calculate_investment()
+        
+        return (self._params.INTEREST_RATE*INVEST)/util
+        
+    def _calculate_depreciation(self) -> float:
+        rval   = self.aircraft["rval"]
+        dyrs   = self.aircraft["dyrs"]
+        util   = self.aircraft["util"]
+
+        INVEST = self._calculate_investment()
+
+        return ((1-rval)*INVEST)/(dyrs*util)

doc_calculator/core/__init__.py

@@ -0,0 +1 @@
+from .DOC_Calculator import DirectOperatingCost

doc_calculator/core/utils/params.py

@@ -0,0 +1,97 @@
+from dataclasses import dataclass
+
+@dataclass
+class Params():
+    """
+    ### Description
+    The dataclass stores the main unit rates and constants necessary to conduct
+    a cost analysis
+    
+    - AEC: Portion of free allocated certificate for CO2 emissions.
+    - ENR: Enroute navigation charges unit rate.
+    - INTEREST_RATE: annual interest rate
+    - LANDINGUR: Landing charge unit rate.
+    - HTONN: Coefficient of cost of handling per tonn of payload
+    - CNOISE: Unit noise rate (Noise tariff that depends on the airport) (USD)
+    - TA : Arrival airport threshold noise (EPNdB)
+    - TD: Departure airport threshold noise (EPNdB)
+    """
+
+    AEC: float = 0.15
+    ENR: float = 68.5
+    INTEREST_RATE = 0.053
+    LANDINGUR: float = 10.0
+    HTONN: float = 45.0
+    CNOISE: float = 4.15
+    TA: float = 89.0
+    TD: float = 92.0
+
+default_dict = {"ADP": 0.0,
+        "MTOW": 0.0,
+        "PLD": 0.,
+        "MEW": 0.0,
+        "BENGW": 0.0,
+        "ENPRI": 0.0,
+        "EN": 0.0,
+        "CREWTECH": 0.0,
+        "CREWC": 0.0,
+        "BT": 1.0,
+        "BF": 0.0,
+        "SECTOR": 0.0,
+        "IENG": 1.0,
+        "SHP": 0.0,
+        "AFSPARE": 0.0,
+        "ENSPARE": 0.0,
+        "DYRS": 1.0,
+        "RVAL": 0.0,
+        "RINSH": 0.0,
+        "CRTECHR": 0.0,
+        "CRCABHR": 0.0,
+        "LABOR_RATE": 0.0,
+        "FUELPRI": 0.0,
+        "IOC_FACT": 0.0,
+        "UTIL": 1.0,
+        "LIFESPAN": 1.0,
+        "L_APP": 0.0,
+        "L_LAT": 0.0,
+        "L_FLYOV": 0.0,
+        "CNOX": 0.0,
+        "NOX_VALUE":0.0,
+        "CCO": 0.0,
+        "CO_VALUE": 0.0,
+        "PRICO2": 0.0,
+        "CO2_VALUE": 0.0, 
+        "ENERPRI": 0.0,
+        "ENER_REQ": 0.0,
+        "H2_PRI": 0.0,
+        "H2_REQ": 0.0,
+        "N_BAT": 0.0,
+        "N_FC": 0.0,
+        "N_REPBAT": 0.0,
+        "BATPRICE": 0.0,
+        "RVBAT": 0.0,
+        "LRBAT": 0.0,
+        "TLBAT": 0.0,
+        "F_BAT": 0.0,
+        "N_REPFC": 0.0,
+        "FCPRICE": 0.0,
+        "RVFC": 0.0,
+        "LRFC": 0.0,
+        "TLFC": 0.0,
+        "F_FC": 0.0,
+        "N_REPPE": 0.0,
+        "PEPRICE": 0.0,
+        "RVPE": 0.0,
+        "LRPE": 0.0,
+        "TLPE": 0.0,
+        "F_PE": 0.0,
+        "N_EM": 0.0,
+        "EMPRICE": 0.0,
+        "LREM": 0.0,
+        "SPEML": 0.0,
+        "SPEMB": 0.0,
+        "TLEML": 0.0,
+        "TLEMB": 0.0,
+        "F_EML": 0.0,
+        "F_EMB": 0.0
+    }

doc_calculator/core/utils/util_functions.py

@@ -0,0 +1,11 @@
+from .params import default_dict as dt
+
+def _assign_input(input:dict) -> dict:
+    default_dict = {k.lower(): v for k, v in dt.items()} 
+    input        = {k.lower(): v for k, v in input.items()} 
+
+    for key, value in default_dict.items():
+        if key in input:
+            default_dict[key] = input[key]
+
+    return default_dict

doc_calculator/gemseo_discipline/DOC_Calculator.py

@@ -0,0 +1,54 @@
+from gemseo.core.discipline.discipline import Discipline
+from ..core import DirectOperatingCost
+from ..core.utils.params import Params
+from .utils.utils_functions import create_default_gemseo_grammar
+import numpy as np
+
+
+class GemseoDirectOperatingCost(Discipline):
+
+    def __init__(self, name="DOC_Calculator", **kwargs):
+        super().__init__(name)
+
+        # define input grammar
+        self.input_grammar.update_from_data(create_default_gemseo_grammar())
+
+        # define output grammar
+        self.output_grammar.update_from_names(["DOC", "IOC", "TOC"])
+
+        # define default data
+        self.default_input_data = create_default_gemseo_grammar()
+
+        # read kwargs params
+        if kwargs:
+            self._params = kwargs["params"]
+        else:
+            self._params = Params()
+
+    def _run(self, input_data):
+
+        # create DOC class aircraft dict
+        aircraft = {}
+
+        for key, value in input_data.items():
+            aircraft[key] = value[0]
+        
+        # instance of the DOC class
+        doc_calc_object = DirectOperatingCost(aircraft, params=self._params)
+
+        # DOC [USD/flight]
+        direct_operating_cost = doc_calc_object.calculate_doc()
+
+        # IOC [USD/flight]
+        indirect_operating_cost = doc_calc_object.calculate_ioc()
+
+        direct_operating_cost_per_flight   = direct_operating_cost["DOC [USD/flight]"]
+        indirect_operating_cost_per_flight = indirect_operating_cost["IOC [USD/flight]"]
+        total_operating_cost_per_flight    = direct_operating_cost_per_flight + indirect_operating_cost_per_flight
+
+        # write output
+        return {
+                "DOC": np.array([direct_operating_cost_per_flight]), 
+                "IOC": np.array([indirect_operating_cost_per_flight]), 
+                "TOC": np.array([total_operating_cost_per_flight])
+                }

doc_calculator/gemseo_discipline/__init__.py

@@ -0,0 +1 @@
+from .DOC_Calculator import GemseoDirectOperatingCost

doc_calculator/gemseo_discipline/utils/utils_functions.py

@@ -0,0 +1,10 @@
+from doc_calculator.core.utils.params import default_dict
+import numpy as np
+
+def create_default_gemseo_grammar() -> dict:
+    default_grammar_dict = {}
+
+    for key, value in default_dict.items():
+        default_grammar_dict[key] = np.array([value])
+
+    return default_grammar_dict

doc_calculator-0.2.dist-info/LICENSE

@@ -0,0 +1,21 @@
+MIT License
+
+Copyright (c) 2025 Michele Tuccillo
+
+Permission is hereby granted, free of charge, to any person obtaining a copy
+of this software and associated documentation files (the "Software"), to deal
+in the Software without restriction, including without limitation the rights
+to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+copies of the Software, and to permit persons to whom the Software is
+furnished to do so, subject to the following conditions:
+
+The above copyright notice and this permission notice shall be included in all
+copies or substantial portions of the Software.
+
+THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+SOFTWARE.

doc_calculator-0.2.dist-info/METADATA

@@ -0,0 +1,6 @@
+Metadata-Version: 2.2
+Name: doc_calculator
+Version: 0.2
+License-File: LICENSE
+Requires-Dist: gemseo
+Dynamic: requires-dist

doc_calculator-0.2.dist-info/RECORD

@@ -0,0 +1,15 @@
+doc_calculator/__init__.py,sha256=l24hn9jSD9FBynMc19l8v8OmxNLKQdKiawL8nbqe8hA,95
+doc_calculator/core/DOC_Calculator.py,sha256=7SPq6GayIm5EP6xwiGnaJ2ReqfdF290yLqCZRg7xb7k,20068
+doc_calculator/core/__init__.py,sha256=ZceJPuJOipP9v84b1sloJGHXppCEcPzElwHw7liG7yA,47
+doc_calculator/core/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+doc_calculator/core/utils/params.py,sha256=FS99-wK_O4-8vXKozPfsk67F6tLVRg_ekAxKehvU5dU,2472
+doc_calculator/core/utils/util_functions.py,sha256=jjPIs-shUAye2hgalqB1rancJ4FPtyRzl0XvXRB_7YI,345
+doc_calculator/gemseo_discipline/DOC_Calculator.py,sha256=hgnZ7eOqNNh34NeehY6c5VEmjzdoVFlryxTi_dUl_XY,1909
+doc_calculator/gemseo_discipline/__init__.py,sha256=lwTMrOr2oUeNTElP14atth9NdxyG6Qr1tUV2231YqBM,53
+doc_calculator/gemseo_discipline/utils/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+doc_calculator/gemseo_discipline/utils/utils_functions.py,sha256=kSaSP2FFllAjw7FI87KEb_IUfJ1uKVVVeREPVJufUKQ,295
+doc_calculator-0.2.dist-info/LICENSE,sha256=0hLLMdm78mvfi_SZni23FLSVq0vYbvmRuLPUwy6Gc7I,1094
+doc_calculator-0.2.dist-info/METADATA,sha256=x8Q2LEF1h0Ozi6DYcKh87GhVv3apkziOjOjYihDt6Jg,129
+doc_calculator-0.2.dist-info/WHEEL,sha256=In9FTNxeP60KnTkGw7wk6mJPYd_dQSjEZmXdBdMCI-8,91
+doc_calculator-0.2.dist-info/top_level.txt,sha256=dTASglOF0CjM5b8VIFNUVVa3wXtrMNXgEwQ4l8wtrss,15
+doc_calculator-0.2.dist-info/RECORD,,

doc_calculator-0.2.dist-info/WHEEL

@@ -0,0 +1,5 @@
+Wheel-Version: 1.0
+Generator: setuptools (75.8.0)
+Root-Is-Purelib: true
+Tag: py3-none-any
+

doc_calculator-0.2.dist-info/top_level.txt

@@ -0,0 +1 @@
+doc_calculator