eelab adds

xemc3/core/load.py

@@ -1519,6 +1519,106 @@ def write_info_file(fn: str, ds: xr.Dataset) -> None:
             "TOTAL_RAD": dict(long_name="Total radiation", units="W"),
         },
     ),
+    "ENERGY_TRANSPORT": dict(
+        type="mapping",
+        vars={
+            "qconve": dict(
+                units="Wm^{-2}", long_name="Parallel electron convective heat flux"
+            ),
+            "qconvi": dict(
+                units="Wm^{-2}", long_name="Parallel ion convective heat flux"
+            ),
+            "qconde": dict(
+                units="Wm^{-2}", long_name="Parallel electron conductive heat flux"
+            ),
+            "qcondi": dict(
+                units="Wm^{-2}", long_name="Parallel ion conductive heat flux"
+            ),
+        },
+    ),
+    "ENERGY_BALANCE": dict(
+        type="mapping",
+        vars={
+            "div_qconve": dict(
+                units="Wm^{-3}",
+                long_name="Divergence of parallel electron convective heat flux",
+                scale=1e6
+            ),
+            "div_qconvi": dict(
+                units="Wm^{-3}",
+                long_name="Divergence of parallel ion convective heat flux",
+                scale=1e6
+            ),
+            "div_qconde": dict(
+                units="Wm^{-3}",
+                long_name="Divergence of parallel ion conductive heat flux",
+                scale=1e6
+            ),
+            "div_qcondi": dict(
+                units="Wm^{-3}",
+                long_name="Divergence of parallel ion conductive heat flux",
+                scale=1e6
+            ),
+            "div_qcross": dict(
+                units="Wm^{-3}", long_name="Divergence of cross-field heat flux"
+                scale=1e6
+            ),
+        },
+    ),
+    "PARTICLE_TRANSPORT": dict(
+        type="mapping",
+        vars={
+            "Gpara": dict(
+                units="m^{-2}s^{-1}", long_name="Particle flux along field lines",
+                scale=1e4
+            )
+        },
+    ),
+    "PARTICLE_BALANCE": dict(
+        type="mapping",
+        vars={
+            "div_Gpara": dict(
+                units="m^{-3}s^{-1}",
+                long_name="Divergence of parallel particle flux",
+                scale=1e6
+            ),
+            "div_Gcross": dict(
+                units="m^{-3}s^{-1}",
+                long_name="Divergence of cross-field particle flux",
+                scale=1e6
+            ),
+        },
+    ),
+    "SOURCE_P": dict(type="mapping", 
+                     vars = dict(units="m^{-3}s^{-1}A^{-1}",
+                                 long_name="Particle sources (normalized)",
+                                 scale=1e6
+                                )
+                    ),
+    "SOURCE_M":dict(type="mapping", 
+                    vars=dict(units="kg m^{-2}s^{-1}A^{-1}",
+                              long_name="Momentum sources (normalized)",
+                              scale=1e7
+                             )
+                   ),
+    "SOURCE_E_E":dict(type="mapping", 
+                      dict(units="W m^{-3}s^{-1}A^{-1}",
+                           long_name="Electron energy sources (normalized)",
+                           scale=1e6
+                          )
+                     ),
+    "SOURCE_E_I":dict(type="mapping", 
+                      dict(units="W m^{-3}s^{-1}A^{-1}",
+                           long_name="Particle sources (normalized)",
+                           scale=1e6
+                          )
+                     ),
+    "TOTAL_FLUX": dict(type="scalar",
+                       vars = {"pflux_total": dict(units="A", 
+                                                   long_name="Total particle flux")
+                               "sp_main" : dict(long_name="Ionization fraction in core")
+                              }
+                      )
 }
 
 if False: