Python Interface to MAPL: Difference between revisions

Newer edit →

@@ Line 1: / Line 1: @@
-Here is a very high level example of a MAPL ocean/atmosphere couple model in Python:
+This document presents a very high level description of a python interface to MAPL.
- import MAPL
+= GEOS-5 Use Cases =
- atm = MAPL.GridComp(rc='atm.rc')
+The examples below are based on the actual GEOS-5 GCM Application, and a few components at the very top of the hierarchy. Here we assumed that MAPL Grid Components can be implemented entirely in C/Fortran or in Python.
- ocn = MAPL.GridComp(rc='ocn.rc')
- atm.initialize()
+== GEOS-5 Application with Fortran Root Component ===
- ocn.initialize()
- for t in clock.tick():
+In this example, the very top root component is implemented entirely in a shared library.
-    atm.run()
-    ocn.run()
-   atm.finalize()
+""""
-   ocn.finalize()
+The GEOSgcm application written in python, starting from a
+Root Component implemented in Fortran.
+"""
+import MAPL
+if __name == "__main__":
+#   Instantiate the root component
+#   ------------------------------
+    gcs = MAPL.GridComp(lib='/some/path/libGEOSana.so')
+#   Instantiate MAPL, using a "Fortran" as a root component
+#   -------------------------------------------------------
+    mapl = MAPL(root=gcs)
+#   Run the cap
+#   -----------
+    mapl.go()