49 char dirname[
LEN], filename[2 *
LEN];
51 int ntask = -1, rank = 0, size = 1;
55 MPI_Init(&argc, &argv);
56 MPI_Comm_rank(MPI_COMM_WORLD, &rank);
57 MPI_Comm_size(MPI_COMM_WORLD, &size);
62 ERRMSG(
"Give parameters: <dirlist> <ctl> <atm_in>");
65 if (!(dirlist = fopen(argv[1],
"r")))
66 ERRMSG(
"Cannot open directory list!");
69 while (fscanf(dirlist,
"%4999s", dirname) != EOF) {
72 if ((++ntask) % size != rank)
76 LOG(1,
"Parallelization: ntask= %d | rank= %d | size= %d",
87 mptrac_alloc(&ctl, &cache, &clim, &met0, &met1, &atm, &dd);
90 sprintf(filename,
"%s/%s", dirname, argv[2]);
97 sprintf(filename,
"%s/%s", dirname, argv[3]);
99 ERRMSG(
"Cannot open file!");
121 if (ctl->
dt_mod > fabs(met0->
lon[1] - met0->
lon[0]) * 111132. / 150.)
122 WARN(
"Violation of CFL criterion! Check DT_MOD!");
146 LOG(1,
"SIZE_NP = %d", atm->
np);
147 LOG(1,
"SIZE_MPI_TASKS = %d", size);
148 LOG(1,
"SIZE_OMP_THREADS = %d", omp_get_max_threads());
151 LOG(1,
"MEMORY_ATM = %g MByte",
sizeof(
atm_t) / 1024. / 1024.);
152 LOG(1,
"MEMORY_CACHE = %g MByte",
sizeof(
cache_t) / 1024. / 1024.);
153 LOG(1,
"MEMORY_CLIM = %g MByte",
sizeof(
clim_t) / 1024. / 1024.);
154 LOG(1,
"MEMORY_METEO = %g MByte",
sizeof(
met_t) / 1024. / 1024.);
157 mptrac_free(ctl, cache, clim, met0, met1, atm, dd);
void mptrac_alloc(ctl_t **ctl, cache_t **cache, clim_t **clim, met_t **met0, met_t **met1, atm_t **atm, dd_t **dd)
Allocates and initializes memory resources for MPTRAC.
void mptrac_run_timestep(ctl_t *ctl, cache_t *cache, clim_t *clim, met_t **met0, met_t **met1, atm_t *atm, double t, dd_t *dd)
Executes a single timestep of the MPTRAC model simulation.
void mptrac_init(ctl_t *ctl, cache_t *cache, clim_t *clim, atm_t *atm, const int ntask)
Initializes the MPTRAC model and its associated components.
void mptrac_get_met(ctl_t *ctl, clim_t *clim, const double t, met_t **met0, met_t **met1, dd_t *dd)
Retrieves meteorological data for the specified time.
void mptrac_read_clim(const ctl_t *ctl, clim_t *clim)
Reads various climatological data and populates the given climatology structure.
int mptrac_read_atm(const char *filename, const ctl_t *ctl, atm_t *atm)
Reads air parcel data from a specified file into the given atmospheric structure.
void mptrac_write_output(const char *dirname, const ctl_t *ctl, met_t *met0, met_t *met1, atm_t *atm, const double t)
Writes various types of output data to files in a specified directory.
void mptrac_read_ctl(const char *filename, int argc, char *argv[], ctl_t *ctl)
Reads control parameters from a configuration file and populates the given structure.
void mptrac_free(ctl_t *ctl, cache_t *cache, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm, dd_t *dd)
Frees memory resources allocated for MPTRAC.
MPTRAC library declarations.
int dd_init(ctl_t *ctl, dd_t *dd, atm_t *atm)
Initializes domain decomposition for parallel processing.
#define LEN
Maximum length of ASCII data lines.
#define ERRMSG(...)
Print an error message with contextual information and terminate the program.
#define WARN(...)
Print a warning message with contextual information.
#define START_TIMERS
Starts a timer for tracking.
#define PRINT_TIMERS
Print the current state of all timers.
#define LOG(level,...)
Print a log message with a specified logging level.
#define STOP_TIMERS
Stop the current timer.
int np
Number of air parcels.
int direction
Direction flag (1=forward calculation, -1=backward calculation).
int dd
Domain decomposition (0=no, 1=yes, with 2x2 if not specified).
double t_stop
Stop time of simulation [s].
double dt_mod
Time step of simulation [s].
double t_start
Start time of simulation [s].
Domain decomposition data structure.
int init
Shows if domain decomposition was initialized.
double lon[EX]
Longitudes [deg].
int main(int argc, char *argv[])