gcmprocpy 1.0.0__tar.gz

gcmprocpy-1.0.0/PKG-INFO

@@ -0,0 +1,23 @@
+Metadata-Version: 2.1
+Name: gcmprocpy
+Version: 1.0.0
+Summary: A Python3 post processing tool for TIE-GCM and WACCM-X
+Home-page: https://github.com/NCAR/gcmprocpy
+Author: Nikhil Rao
+Author-email: nikhilr@ucar.edu
+Requires-Python: >=3.8
+Description-Content-Type: text/markdown
+Requires-Dist: cartopy
+Requires-Dist: matplotlib
+Requires-Dist: numpy
+Requires-Dist: xarray
+Requires-Dist: ipython
+Requires-Dist: geomag
+
+# GCMPROCPY
+
+GCMprocpy is a tool used for TIE-GCM and WACCM-X post processing and plot generation.
+
+## Documentation on readthedocs
+## https://gcmprocpy.readthedocs.io
+

gcmprocpy-1.0.0/README.md

@@ -0,0 +1,7 @@
+# GCMPROCPY
+
+GCMprocpy is a tool used for TIE-GCM and WACCM-X post processing and plot generation.
+
+## Documentation on readthedocs
+## https://gcmprocpy.readthedocs.io
+

gcmprocpy-1.0.0/setup.cfg

@@ -0,0 +1,4 @@
+[egg_info]
+tag_build = 
+tag_date = 0
+

gcmprocpy-1.0.0/setup.py

@@ -0,0 +1,37 @@
+from setuptools import setup, find_packages
+
+with open("README.md", "r", encoding="utf-8") as fh:
+    long_description = fh.read()
+
+setup(
+    name='gcmprocpy',
+    version='1.0.0',
+    author = "Nikhil Rao",
+    author_email = "nikhilr@ucar.edu",
+    description='A Python3 post processing tool for TIE-GCM and WACCM-X',
+    long_description=long_description,
+    long_description_content_type="text/markdown",
+    url='https://github.com/NCAR/gcmprocpy', 
+    python_requires='>=3.8',
+    install_requires=[
+        'cartopy',
+        'matplotlib',
+        'numpy',
+        'xarray',
+        'ipython',
+        'geomag',
+    ],
+    package_dir={'': 'src'},  
+    packages=find_packages(where='src'), 
+    entry_points={
+        'console_scripts': [
+            'lat_lon= gcmprocpy.cmd.cmd_lat_lon:cmd_plt_lat_lon',
+            'lev_var= gcmprocpy.cmd.cmd_lev_var:cmd_plt_lev_var',
+            'lev_lat= gcmprocpy.cmd.cmd_lev_lat:cmd_plt_lev_lat',
+            'lev_lon= gcmprocpy.cmd.cmd_lev_lon:cmd_plt_lev_lon',
+            'lev_time= gcmprocpy.cmd.cmd_lev_time:cmd_plt_lev_time',
+            'lat_time= gcmprocpy.cmd.cmd_lat_time:cmd_plt_lat_time',
+            
+        ]
+    }
+)

gcmprocpy-1.0.0/src/gcmprocpy/__init__.py

@@ -0,0 +1,5 @@
+from .plot_gen import plt_lat_lon, plt_lev_var, plt_lev_lon, plt_lev_lat, plt_lev_time, plt_lat_time
+from .data_parse import arr_lat_lon,arr_lev_var,arr_lev_lon, arr_lev_lat,arr_lev_time,arr_lat_time, arr_var, get_mtime, get_time, time_list, var_list, level_list, lon_list, lat_list
+from .io import load_datasets
+from .mov_gen import mov_lat_lon
+from .data_emissions import mkeno53, mkeco215, mkeoh83, arr_mkeno53, arr_mkeco215, arr_mkeoh83

gcmprocpy-1.0.0/src/gcmprocpy/cmd/cmd_lat_lon.py

@@ -0,0 +1,50 @@
+#!/usr/bin/env python3
+from ..plot_gen import plt_lat_lon
+from ..io import load_datasets, save_output
+import argparse
+import os
+
+def cmd_parser():
+    parser = argparse.ArgumentParser(description="Parser for loading, plotting, and saving")
+
+    # Loading datasets
+    parser.add_argument('-dir','--directory', type=str, help='Path to the directory containing the datasets')
+    parser.add_argument('-dsf','--dataset_filter', type=str, help='Filter for the dataset file names', default=None)
+    
+    # Saving output
+    parser.add_argument('-o_dir','--output_directory', type=str, help='Directory where the plot will be saved.', default=os.getcwd())
+    parser.add_argument('-o_file','--filename', type=str, required=True, help='Filename for the saved plot.')
+    parser.add_argument('-o_format','--output_format', type=str, help='Format of the output plot, e.g., "png", "pdf".', default='jpg')
+
+    # Plotting parameters
+    parser.add_argument('-var','--variable_name', type=str, required=True, help='The name of the variable with latitude, longitude, and lev/ilev dimensions.')
+    parser.add_argument('-t','--time', type=str, help='The selected time, e.g., "2022-01-01T12:00:00".')
+    parser.add_argument('-mt','--mtime', nargs='+', type=int, help='The selected time as a list, e.g., [1, 12, 0] for 1st day, 12 hours, 0 mins.')
+    parser.add_argument('-lvl','--level', type=float, help='The selected lev/ilev value.')
+    parser.add_argument('-unit','--variable_unit', type=str, help='The desired unit of the variable.')
+    parser.add_argument('-ci','--contour_intervals', type=int, help='The number of contour intervals. Defaults to 20.')
+    parser.add_argument('-cv','--contour_value', type=int, help='The value between each contour interval.')
+    parser.add_argument('-si','--symmetric_interval', action='store_true', help='If True, the contour intervals will be symmetric around zero. Defaults to False.')
+    parser.add_argument('-cmc','--cmap_color', type=str, help='The color map of the contour. Defaults to "viridis" for Density, "inferno" for Temp, "bwr" for Wind, "viridis" for undefined.')
+    parser.add_argument('-lc','--line_color', type=str, default='white', help='The color for all lines in the plot. Defaults to "white".')
+    parser.add_argument('-cst','--coastlines', action='store_true', help='Shows coastlines on the plot. Defaults to False.')
+    parser.add_argument('-nsh','--nightshade', action='store_true', help='Shows nightshade on the plot. Defaults to False.')
+    parser.add_argument('-gm','--gm_equator', action='store_true', help='Shows geomagnetic equator on the plot. Defaults to False.')
+    parser.add_argument('-lat_min','--latitude_minimum', type=float, help='Minimum latitude to slice plots. Defaults to -87.5.')
+    parser.add_argument('-lat_max','--latitude_maximum', type=float, help='Maximum latitude to slice plots. Defaults to 87.5.')
+    parser.add_argument('-lon_min','--longitude_minimum', type=float, help='Minimum longitude to slice plots. Defaults to -180.')
+    parser.add_argument('-lon_max','--longitude_maximum', type=float, help='Maximum longitude to slice plots. Defaults to 175.')
+    return (parser)
+
+
+
+
+def cmd_plt_lat_lon():
+    parser = cmd_parser()
+    args = parser.parse_args()
+    datasets = load_datasets(args.directory,args.dataset_filter)
+    plot =  plt_lat_lon(datasets, variable_name=args.variable_name, time=args.time, mtime=args.mtime, level=args.level, variable_unit=args.variable_unit, contour_intervals=args.contour_intervals, contour_value=args.contour_value, symmetric_interval=args.symmetric_interval, cmap_color=args.cmap_color, line_color=args.line_color, coastlines=args.coastlines, nightshade=args.nightshade, gm_equator=args.gm_equator, latitude_minimum=args.latitude_minimum, latitude_maximum=args.latitude_maximum, longitude_minimum=args.longitude_minimum, longitude_maximum=args.longitude_maximum)
+    save_output(args.output_directory,args.filename,args.output_format,plot)
+
+
+

gcmprocpy-1.0.0/src/gcmprocpy/cmd/cmd_lat_time.py

@@ -0,0 +1,45 @@
+#!/usr/bin/env python3
+from ..plot_gen import plt_lat_time
+from ..io import load_datasets, save_output
+import argparse
+import os
+
+def cmd_parser():
+    parser = argparse.ArgumentParser(description="Parser for loading, plotting, and saving")
+
+    # Loading datasets
+    parser.add_argument('-dir','--directory', type=str, help='Path to the directory containing the datasets')
+    parser.add_argument('-dsf','--dataset_filter', type=str, help='Filter for the dataset file names', default=None)
+    
+    # Saving output
+    parser.add_argument('-o_dir','--output_directory', type=str, help='Directory where the plot will be saved.', default=os.getcwd())
+    parser.add_argument('-o_file','--filename', type=str, required=True, help='Filename for the saved plot.')
+    parser.add_argument('-o_format','--output_format', type=str, required=True, help='Format of the output plot, e.g., "png", "pdf".', default='jpg')
+
+    # Plotting parameters
+    parser.add_argument('-var','--variable_name', type=str, help="The name of the variable with latitude, longitude, time, and ilev dimensions.")
+    parser.add_argument('-lvl','--level', type=float, help="The specific level value for the plot.", default=None)
+    parser.add_argument('-lon','--longitude', type=float, help="The specific longitude value for the plot.", default=None)
+    parser.add_argument('-unit','--variable_unit', type=str, help="The desired unit of the variable.", default=None)
+    parser.add_argument('-ci','--contour_intervals', type=int, help="The number of contour intervals. Defaults to 10.", default=10)
+    parser.add_argument('-cv','--contour_value', type=int, help="The value between each contour interval.", default=None)
+    parser.add_argument('-si','--symmetric_interval', action='store_true', help="If True, the contour intervals will be symmetric around zero. Defaults to False.")
+    parser.add_argument('-cmc','--cmap_color', type=str, help="The color map of the contour. Defaults to 'viridis' for Density, 'inferno' for Temp, 'bwr' for Wind, 'viridis' for undefined.", default=None)
+    parser.add_argument('-lc','--line_color', type=str, help="The color for all lines in the plot. Defaults to 'white'.", default='white')
+    parser.add_argument('-lat_min','--latitude_minimum', type=float, help="Minimum latitude value for the plot. Defaults to -87.5.", default=-87.5)
+    parser.add_argument('-lat_max','--latitude_maximum', type=float, help="Maximum latitude value for the plot. Defaults to 87.5.", default=87.5)
+    parser.add_argument('--mtime_minimum', type=float, help="Minimum time value for the plot. Defaults to None.", default=None)
+    parser.add_argument('--mtime_maximum', type=float, help="Maximum time value for the plot. Defaults to None.", default=None)
+    
+    return (parser)
+
+
+
+
+def cmd_plt_lat_time():
+    parser = cmd_parser()
+    args = parser.parse_args()
+    datasets = load_datasets(args.directory,args.dataset_filter)
+    plot = plt_lat_time(datasets, variable_name=args.variable_name, level=args.level, longitude=args.longitude, variable_unit=args.variable_unit, contour_intervals=args.contour_intervals, contour_value=args.contour_value, symmetric_interval=args.symmetric_interval, cmap_color=args.cmap_color, line_color=args.line_color, latitude_minimum=args.latitude_minimum, latitude_maximum=args.latitude_maximum, mtime_minimum=args.mtime_minimum, mtime_maximum=args.mtime_maximum)
+    save_output(args.output_directory,args.filename,args.output_format,plot)
+

gcmprocpy-1.0.0/src/gcmprocpy/cmd/cmd_lev_lat.py

@@ -0,0 +1,45 @@
+#!/usr/bin/env python3
+from ..plot_gen import plt_lev_lat
+from ..io import load_datasets, save_output
+import argparse
+import os
+
+def cmd_parser():
+    parser = argparse.ArgumentParser(description="Parser for loading, plotting, and saving")
+
+    # Loading datasets
+    parser.add_argument('-dir','--directory', type=str, help='Path to the directory containing the datasets')
+    parser.add_argument('-dsf','--dataset_filter', type=str, help='Filter for the dataset file names', default=None)
+    
+    # Saving output
+    parser.add_argument('-o_dir','--output_directory', type=str, help='Directory where the plot will be saved.', default=os.getcwd())
+    parser.add_argument('-o_file','--filename', type=str, required=True, help='Filename for the saved plot.')
+    parser.add_argument('-o_format','--output_format', type=str, required=True, help='Format of the output plot, e.g., "png", "pdf".', default='jpg')
+
+    # Plotting parameters
+    parser.add_argument('-var','--variable_name', type=str, help="The name of the variable with latitude, longitude, and lev/ilev dimensions.")
+    parser.add_argument('-t','--time', type=str, help="The selected time, e.g., '2022-01-01T12:00:00'.", default=None)
+    parser.add_argument('-mt','--mtime', type=int, nargs=3, help="The selected time as a list, e.g., [1, 12, 0] for 1st day, 12 hours, 0 mins.", default=None)
+    parser.add_argument('-lon','--longitude', type=float, help="The specific longitude value for the plot.", default=None)
+    parser.add_argument('-unit','--variable_unit', type=str, help="The desired unit of the variable.", default=None)
+    parser.add_argument('-ci','--contour_intervals', type=int, help="The number of contour intervals. Defaults to 20.", default=20)
+    parser.add_argument('-cv','--contour_value', type=int, help="The value between each contour interval.", default=None)
+    parser.add_argument('-si','--symmetric_interval', action='store_true', help="If True, the contour intervals will be symmetric around zero. Defaults to False.")
+    parser.add_argument('-cmc','--cmap_color', type=str, help="The color map of the contour. Defaults to 'viridis' for Density, 'inferno' for Temp, 'bwr' for Wind, 'viridis' for undefined.", default=None)
+    parser.add_argument('-lc','--line_color', type=str, help="The color for all lines in the plot. Defaults to 'white'.", default='white')
+    parser.add_argument('-lvl_min','--level_minimum', type=float, help="Minimum level value for the plot. Defaults to None.", default=None)
+    parser.add_argument('-lvl_max','--level_maximum', type=float, help="Maximum level value for the plot. Defaults to None.", default=None)
+    parser.add_argument('-lat_min','--latitude_minimum', type=float, help="Minimum latitude value for the plot. Defaults to -87.5.", default=-87.5)
+    parser.add_argument('-lat_max','--latitude_maximum', type=float, help="Maximum latitude value for the plot. Defaults to 87.5.", default=87.5)
+    
+    return (parser)
+
+
+
+
+def cmd_plt_lev_lat():
+    parser = cmd_parser()
+    args = parser.parse_args()
+    datasets = load_datasets(args.directory,args.dataset_filter)
+    plot = plt_lev_lat(datasets, variable_name=args.variable_name, time=args.time, mtime=args.mtime, longitude=args.longitude, variable_unit=args.variable_unit, contour_intervals=args.contour_intervals, contour_value=args.contour_value, symmetric_interval=args.symmetric_interval, cmap_color=args.cmap_color, line_color=args.line_color, level_minimum=args.level_minimum, level_maximum=args.level_maximum, latitude_minimum=args.latitude_minimum, latitude_maximum=args.latitude_maximum)
+    save_output(args.output_directory,args.filename,args.output_format,plot)

gcmprocpy-1.0.0/src/gcmprocpy/cmd/cmd_lev_lon.py

@@ -0,0 +1,45 @@
+#!/usr/bin/env python3
+from ..plot_gen import plt_lev_lon
+from ..io import load_datasets, save_output
+import argparse
+import os
+
+def cmd_parser():
+    parser = argparse.ArgumentParser(description="Parser for loading, plotting, and saving")
+
+    # Loading datasets
+    parser.add_argument('-dir','--directory', type=str, help='Path to the directory containing the datasets')
+    parser.add_argument('-dsf','--dataset_filter', type=str, help='Filter for the dataset file names', default=None)
+    
+    # Saving output
+    parser.add_argument('-o_dir','--output_directory', type=str, help='Directory where the plot will be saved.', default=os.getcwd())
+    parser.add_argument('-o_file','--filename', type=str, required=True, help='Filename for the saved plot.')
+    parser.add_argument('-o_format','--output_format', type=str, required=True, help='Format of the output plot, e.g., "png", "pdf".', default='jpg')
+
+    # Plotting parameters
+    parser.add_argument('-var','--variable_name', type=str, help="The name of the variable with latitude, longitude, and lev/ilev dimensions.")
+    parser.add_argument('-lat','--latitude', type=float, help="The specific latitude value for the plot.")
+    parser.add_argument('-t','--time', type=str, help="The selected time, e.g., '2022-01-01T12:00:00'.", default=None)
+    parser.add_argument('-mt','--mtime', type=int, nargs=3, help="The selected time as a list, e.g., [1, 12, 0] for 1st day, 12 hours, 0 mins.", default=None)
+    parser.add_argument('-unit','--variable_unit', type=str, help="The desired unit of the variable.", default=None)
+    parser.add_argument('-ci','--contour_intervals', type=int, help="The number of contour intervals. Defaults to 20.", default=20)
+    parser.add_argument('-cv','--contour_value', type=int, help="The value between each contour interval.", default=None)
+    parser.add_argument('-si','--symmetric_interval', action='store_true', help="If True, the contour intervals will be symmetric around zero. Defaults to False.")
+    parser.add_argument('-cmc','--cmap_color', type=str, help="The color map of the contour. Defaults to 'viridis' for Density, 'inferno' for Temp, 'bwr' for Wind, 'viridis' for undefined.", default=None)
+    parser.add_argument('-lc','--line_color', type=str, help="The color for all lines in the plot. Defaults to 'white'.", default='white')
+    parser.add_argument('-lvl_min','--level_minimum', type=float, help="Minimum level value for the plot. Defaults to None.", default=None)
+    parser.add_argument('-lvl_max','--level_maximum', type=float, help="Maximum level value for the plot. Defaults to None.", default=None)
+    parser.add_argument('-lon_min','--longitude_minimum', type=float, help="Minimum longitude value for the plot. Defaults to -180.", default=-180)
+    parser.add_argument('-lon_max','--longitude_maximum', type=float, help="Maximum longitude value for the plot. Defaults to 175.", default=175)
+    return (parser)
+
+
+
+
+def cmd_plt_lev_lon():
+    parser = cmd_parser()
+    args = parser.parse_args()
+    datasets = load_datasets(args.directory,args.dataset_filter)
+    plot = plt_lev_lon(datasets, variable_name=args.variable_name, latitude=args.latitude, time=args.time, mtime=args.mtime, variable_unit=args.variable_unit, contour_intervals=args.contour_intervals, contour_value=args.contour_value, symmetric_interval=args.symmetric_interval, cmap_color=args.cmap_color, line_color=args.line_color, level_minimum=args.level_minimum, level_maximum=args.level_maximum, longitude_minimum=args.longitude_minimum)
+    save_output(args.output_directory,args.filename,args.output_format,plot)
+

gcmprocpy-1.0.0/src/gcmprocpy/cmd/cmd_lev_time.py

@@ -0,0 +1,45 @@
+#!/usr/bin/env python3
+from ..plot_gen import plt_lev_time
+from ..io import load_datasets, save_output
+import argparse
+import os
+
+def cmd_parser():
+    parser = argparse.ArgumentParser(description="Parser for loading, plotting, and saving")
+
+    # Loading datasets
+    parser.add_argument('-dir','--directory', type=str, help='Path to the directory containing the datasets')
+    parser.add_argument('-dsf','--dataset_filter', type=str, help='Filter for the dataset file names', default=None)
+    
+    # Saving output
+    parser.add_argument('-o_dir','--output_directory', type=str, help='Directory where the plot will be saved.', default=os.getcwd())
+    parser.add_argument('-o_file','--filename', type=str, required=True, help='Filename for the saved plot.')
+    parser.add_argument('-o_format','--output_format', type=str, required=True, help='Format of the output plot, e.g., "png", "pdf".', default='jpg')
+
+    # Plotting parameters
+    parser.add_argument('-var','--variable_name', type=str, help="The name of the variable with latitude, longitude, time, and ilev dimensions.")
+    parser.add_argument('-lat','--latitude', type=float, help="The specific latitude value for the plot.")
+    parser.add_argument('-lon','--longitude', type=float, help="The specific longitude value for the plot.", default=None)
+    parser.add_argument('-unit','--variable_unit', type=str, help="The desired unit of the variable.", default=None)
+    parser.add_argument('-ci','--contour_intervals', type=int, help="The number of contour intervals. Defaults to 10.", default=10)
+    parser.add_argument('-cv','--contour_value', type=int, help="The value between each contour interval.", default=None)
+    parser.add_argument('-si','--symmetric_interval', action='store_true', help="If True, the contour intervals will be symmetric around zero. Defaults to False.")
+    parser.add_argument('-cmc','--cmap_color', type=str, help="The color map of the contour. Defaults to 'viridis' for Density, 'inferno' for Temp, 'bwr' for Wind, 'viridis' for undefined.", default=None)
+    parser.add_argument('-lc','--line_color', type=str, help="The color for all lines in the plot. Defaults to 'white'.", default='white')
+    parser.add_argument('-lvl_min','--level_minimum', type=float, help="Minimum level value for the plot. Defaults to None.", default=None)
+    parser.add_argument('-lvl_max','--level_maximum', type=float, help="Maximum level value for the plot. Defaults to None.", default=None)
+    parser.add_argument('--mtime_minimum', type=float, help="Minimum time value for the plot. Defaults to None.", default=None)
+    parser.add_argument('--mtime_maximum', type=float, help="Maximum time value for the plot. Defaults to None.", default=None)
+    return (parser)
+
+
+
+
+def cmd_plt_lev_time():
+    parser = cmd_parser()
+    args = parser.parse_args()
+    datasets = load_datasets(args.directory,args.dataset_filter)
+    plot = plt_lev_time(datasets, variable_name=args.variable_name, latitude=args.latitude, longitude=args.longitude, variable_unit=args.variable_unit, contour_intervals=args.contour_intervals, contour_value=args.contour_value, symmetric_interval=args.symmetric_interval, cmap_color=args.cmap_color, line_color=args.line_color, level_minimum=args.level_minimum, level_maximum=args.level_maximum, mtime_minimum=args.mtime_minimum, mtime_maximum=args.mtime_maximum)
+    save_output(args.output_directory,args.filename,args.output_format,plot)
+
+

gcmprocpy-1.0.0/src/gcmprocpy/cmd/cmd_lev_var.py

@@ -0,0 +1,39 @@
+#!/usr/bin/env python3
+from ..plot_gen import plt_lev_var
+from ..io import load_datasets, save_output
+import argparse
+import os
+
+def cmd_parser():
+    parser = argparse.ArgumentParser(description="Parser for loading, plotting, and saving")
+
+    # Loading datasets
+    parser.add_argument('-dir','--directory', type=str, help='Path to the directory containing the datasets')
+    parser.add_argument('-dsf','--dataset_filter', type=str, help='Filter for the dataset file names', default=None)
+    
+    # Saving output
+    parser.add_argument('-o_dir','--output_directory', type=str, help='Directory where the plot will be saved.', default=os.getcwd())
+    parser.add_argument('-o_file','--filename', type=str, required=True, help='Filename for the saved plot.')
+    parser.add_argument('-o_format','--output_format', type=str, required=True, help='Format of the output plot, e.g., "png", "pdf".', default='jpg')
+
+    # Plotting parameters
+    parser.add_argument('-var','--variable_name', type=str, help='The name of the variable with latitude, longitude, and lev/ilev dimensions')
+    parser.add_argument('-lat','--latitude', type=float, help='The specific latitude value for the plot')
+    parser.add_argument('-t','--time', type=str, help='The selected time, e.g., "2022-01-01T12:00:00"', default=None)
+    parser.add_argument('-mt','--mtime', nargs=3, type=int, help='The selected time as a list, e.g., [1, 12, 0] for 1st day, 12 hours, 0 mins', default=None)
+    parser.add_argument('-lon','--longitude', type=float, help='The specific longitude value for the plot', default=None)
+    parser.add_argument('-unit','--variable_unit', type=str, help='The desired unit of the variable', default=None)
+    parser.add_argument('-lvl_min','--level_minimum', type=float, help='Minimum level value for the plot', default=None)
+    parser.add_argument('-lvl_max','--level_maximum', type=float, help='Maximum level value for the plot', default=None)
+    return (parser)
+
+
+
+
+def cmd_plt_lev_var():
+    parser = cmd_parser()
+    args = parser.parse_args()
+    datasets = load_datasets(args.directory,args.dataset_filter)
+    plot = plt_lev_var(datasets,args.variable_name,args.latitude,args.time,args.mtime,args.longitude,args.variable_unit,args.level_minimum,args.level_maximum)
+    save_output(args.output_directory,args.filename,args.output_format,plot)
+

gcmprocpy-1.0.0/src/gcmprocpy/convert_units.py

@@ -0,0 +1,56 @@
+conversion_units = {
+    'cm-3': {'m-3': 1000000.0},
+    'km': {'m': 1000.0},
+    'cm/s': {'m/s': 0.01, 'km/s': 1e-5},
+    'erg/g/s': {'J/kg/s': 1e-7},
+    'ergs/cm2/s': {'J/m2/s': 0.001},
+    'millibars': {'pascals': 100},
+    'microbars': {'pascals': 0.1},
+    'g/cm3': {'kg/m3': 1000.0},
+    'degrees': {'radians': 0.0174533},
+    'degrees_east': {'radians': 0.0174533},
+    'degrees_north': {'radians': 0.0174533},
+    'km/s': {'m/s': 1000.0, 'cm/s': 100000.0},
+    'GW': {'MW': 1000.0, 'kW': 1000000.0, 'W': 1000000000.0},
+    'keV': {'eV': 1000.0, 'J': 1.60218e-16},
+    'nT': {'µT': 0.001, 'T': 1e-9},
+    'cm': {'m': 0.01, 'km': 1e-5},
+    'K': {
+        'C': {'factor': 1, 'offset': -273.15},  # Celsius = Kelvin - 273.15
+        'F': {'factor': 9/5, 'offset': -459.67}  # Fahrenheit = Kelvin * 9/5 - 459.67
+    }
+}
+
+def convert_units(data, from_unit, to_unit):
+    """
+    Convert data from one unit to another based on predefined conversion factors.
+
+    Args:
+        data (float): Numeric data to be converted.
+        from_unit (str): The current unit of the data.
+        to_unit (str): The desired unit to convert to.
+
+    Returns:
+        float: Converted data.
+    """
+    
+    if from_unit == to_unit:
+        return data, from_unit
+    
+    # Check if conversion is possible
+    if from_unit in conversion_units and to_unit in conversion_units[from_unit]:
+        conversion = conversion_units[from_unit][to_unit]
+        print(f"Converting from {from_unit} to {to_unit}")
+
+        # Check if conversion is a dictionary or direct float
+        if isinstance(conversion, dict):
+            factor = conversion.get('factor', 1)
+            offset = conversion.get('offset', 0)
+
+            return data * factor + offset, to_unit
+        else:
+            return data * conversion, to_unit
+    else:
+        #raise ValueError(f"No conversion factor found for {from_unit} to {to_unit}")
+        print(f"No conversion factor found for {from_unit} to {to_unit}")
+        return data, from_unit