from pycoupler.run import run_lpjml
                from pycoupler.coupler import LPJmLCoupler

                # Run lpjml simulation with socket coupling set
                run_lpjml(config_file=config_coupled_fn)
                # Establish coupler connection to LPJmL
                lpjml = LPJmLCoupler(config_file=config_coupled_fn)
                # Get initial data of previous years of simulation
                inputs = lpjml.read_input()
                outputs = lpjml.read_historic_output()

                for year in lpjml.get_sim_years():
                    # Placeholder to interact with inputs
                    lpjml.send_input(inputs)
                    lpjml.read_output(outputs)

                lpjml.close()