#ifndef GEMINI3D_H
#define GEMINI3D_H

// needed when compiling C files
#if __STDC_VERSION__ < 202311L
#include <stdbool.h>
#endif

#ifdef __cplusplus
extern "C" {
#endif

enum { LMAX = 1000 };

struct params {
  // order and lengths must match in Fortran and C
  // see gemini_main.f90 "cparams"
  int fortran_nml;
  int fortran_cli;
  int debug;
  int dryrun;
  char out_dir[LMAX];
  // .ini [base]
  int ymd[3];
  float UTsec0;
  float tdur;
  float dtout;
  float activ[3];
  float tcfl;
  float Teinf;
  // .ini
};

// cfgC for now is just passed around as void**, to modify it would have to be struct *


/* This is now housed in the main C program */
extern int gemini_main(struct params *, int*, int*);

extern void help_gemini_bin();

/* interfaces to libgemini */
extern void gemini_cfg_alloc_C(void**);
extern void gemini_cfg_dealloc_C(void**);
extern void cli_in_C(struct params*, int*, int*, void**);
extern void read_config_in_C(struct params*, void**);
extern void set_magnetic_pole_in_C(void**);
extern void grid_size_in_C(void**);
extern void get_fullgrid_size_C(int*, int*, int*);
extern void get_config_vars_C(void**, bool*, int*, double*, double*);
extern void get_subgrid_size_C(int*, int*, int*);
extern void get_species_size_C(int*);
extern void get_fullgrid_lims_C(double*, double*, double*, double*, double*, double*);
extern void gemini_work_alloc_C(void**,void**);
extern void gemini_work_dealloc_C(void**, void**);
//extern void memblock_from_C(double**, double**, double**);
extern void set_start_values_auxtimevars_C(double*, double*, double*);
extern void set_start_timefromcfg_C(void**, int*, double*, double*);
extern void set_start_values_auxvars_C(int*, void**, double**);
extern void get_cfg_timevars_C(void**,double*,int*,double*,int*,int*);
extern void msisinit_C(void**);
extern void init_neutralBG_input_C(void**, int*, void**, double*, double*, int*, double*, void**);
extern void set_update_cadence_C(int*);
extern void neutral_atmos_winds_C(void**, int*, void**, int*, double*, void**);
extern void check_finite_output_C(void**, double**, double**, double*);
extern void get_solar_indices_C(void**, double*, double*);
extern void v12rhov1_C(double**, double**);
extern void T2rhoe_C(double**, double**);
extern void interface_vels_allspec_C(int*, void**, double**, void**, int*);
extern void sweep3_allparams_C(double**, double**, void**, int*, void**, double*);
extern void sweep3_allspec_mass_C(double**, double**, void**, int*, void**, double*);
extern void sweep3_allspec_momentum_C(double**, double**, void**, int*, void**, double*);
extern void sweep3_allspec_energy_C(double**, double**, void**, int*, void**, double*);
extern void sweep1_allparams_C(double**, double**, void**, int*, void**, double*);
extern void sweep1_allspec_mass_C(double**, double**, void**, int*, void**, double*);
extern void sweep1_allspec_momentum_C(double**, double**, void**, int*, void**, double*);
extern void sweep1_allspec_energy_C(double**, double**, void**, int*, void**, double*);
extern void sweep2_allparams_C(double**, double**, void**, int*, void**, double*);
extern void sweep2_allspec_mass_C(double**, double**, void**, int*, void**, double*);
extern void sweep2_allspec_momentum_C(double**, double**, void**, int*, void**, double*);
extern void sweep2_allspec_energy_C(double**, double**, void**, int*, void**, double*);
extern void rhov12v1_C(double**, double**);
extern void clean_param_C(int*, int*, void**, double**);
extern void VNRicht_artvisc_C(double**, void**);
extern void compression_C(double**, double**, void**, int*, void**, double*);
extern void rhoe2T_C(double**, double**);
extern void energy_diffusion_C(void**, int*, void**, double**, double**, void**, double*);
extern void source_loss_allparams_C(void**, double**, double**, double**, void**, int*, void**,
                                      double*);
extern void source_loss_mass_C(void**, double**,double**,double**,void**,int*,void**,double*);
extern void source_loss_momentum_C(void**, double**,double**,double**,void**,int*,void**,double*);
extern void source_loss_energy_C(void**, double**,double**,double**,void**,int*,void**,double*);
extern void clear_ionization_arrays_C(void**);
extern void impact_ionization_C(void**,double**,void**,int*,void**,double*,double*,int*,double*);
extern void solar_ionization_C(void**,double**,void**,int*,void**,double*,int*,double*);
extern void check_dryrun_C(void**);
extern void maxcfl_C(double**, int*, void**, double*, double*);
extern void dateinc_C(double*, int*, double*);
extern void plasma_output_nompi_C(void**, int*, double*, double**, double**, int*, double*, double*, double*);
extern void read_fullsize_gridcenter_C(void**);
extern void gemini_grid_dealloc_C(int*,void**);
extern void gemini_grid_generate_C(int*, void**);
extern void gemini_grid_generate_altnull_C(int*, void**, double*);
extern void setv2v3_C(double*, double*);
extern void set_global_boundaries_allspec_C(int*, void**, double**, double**, void**, int*);
extern void set_electrodynamics_commtype_C(int*, void**);
extern void electrodynamics_test_C(void**,int*,void**,double**,double**,double**,void**);
extern void get_neutralperturb_interptype_C(void**,int*);
extern void init_precipinput_C(void**, int*, void**, double*, double*, int*, double*, void**);
extern void precip_perturb_C(void**,void**,int*,void**,double*,double*,int*,double*);
extern void init_efieldinput_nompi_C(void**, int*, void**, double*, double*, int*, double*, void**);
extern void efield_perturb_nompi_C(void**,void**,int*,void**,double*,double*,int*,double*);
extern void interp3_C(double**,double**,double**,int*,int*,int*,double**,double**,int*,int*,
		double*,double*,double*,double**,double**,int*);
extern void interp2_C(double**,double**,int*,int*,double**,double**,int*,int*,double*,double*,double**,
		double**);
extern void get_it_C(int*);
extern void itinc_C();
extern void init_solfluxinput_C(void**, int*, void**, double*, double*, int*, double*, void**);
extern void solflux_perturb_C(void**, void**, int*, void**, double*, double*, int*, double*);
extern void source_neut_C(void**,double**,void**,int*,void**);

/* interfaces for libgemini_mpi */
// some of these will very likely need to be rewritten when used with forestclaw
extern void mpisetup_C();
extern void mpiparms_C(int*, int*);
extern void init_procgrid_C(int*, int*, int*, int*);
extern void read_grid_C(void**, int*, void**);
extern void outdir_fullgridvaralloc_C(void**, void**, int*, int*, int*);
extern void calc_subgrid_size_in_C(int*, int*);
extern void get_initial_state_C(void**, double**, double**, void**, int*, void**, double*, int*, double*, double*, double*);
extern void init_inputdata_C(void**, int*, void**, double*, double*, int*, double*, void**);
extern void init_Efieldinput_C(void**, int*, void**, double*, double*, void**, int*, double*);
extern void pot2perpfield_C(int*, void**, double**);
extern void BGfield_Lagrangian_C(void**, int*, void**, double**, void**);
extern void init_neutralperturb_C(double*, void**, int*, void**, void**, int*, double*);
extern void get_initial_drifts_C(void**, int*, void**, double**, double**, double**, void**);
extern void neutral_atmos_wind_update_C(void**);
extern void neutral_perturb_C(void**, void**, int*, void**, double*, double*, int*, double*);
extern void efield_perturb_C(void**, void**, int*, void**, double*, double*, int*, double*);
extern void inputdata_perturb_C(void**, void**, int*, void**, double*, double*, int*, double*);
extern void electrodynamics_C(void**, double**, double**, double**, void**, int*, void**, double*, double*, int*, double*);
extern void get_gavg_Tinf_C(void**);
extern void clear_neutral_perturb_C(void**);
extern void clear_neutral_background_C(void**);
extern void halo_interface_vels_allspec_C(int*, void**, double**, int*);
extern void halo_allparams_C(int*, void**, double**, double**);
extern void halo_fluidvars_C(int*, void**, double**, double**);
extern void RK2_prep_mpi_allspec_C(int*, void**, double**);
extern void RK2_global_boundary_allspec_C(int*, void**, double**);
extern void dt_select_C(void**, int*, void**, double**, double**, double*, double*, double*, double*);
extern void check_fileoutput_C(void**, double**, double**, void**, double*, double*, double*, double*, int*, int*, double*);

#ifdef __cplusplus
}
#endif

#endif