51 char dirname[
LEN], filename[2 *
LEN];
53 int ntask = -1, rank = 0, size = 1;
57 MPI_Init(&argc, &argv);
58 MPI_Comm_rank(MPI_COMM_WORLD, &rank);
59 MPI_Comm_size(MPI_COMM_WORLD, &size);
70 ERRMSG(
"Give parameters: <dirlist> <ctl> <atm_in>");
73 if (!(dirlist = fopen(argv[1],
"r")))
74 ERRMSG(
"Cannot open directory list!");
77 while (fscanf(dirlist,
"%4999s", dirname) != EOF) {
80 if ((++ntask) % size != rank)
84 LOG(1,
"Parallelization: ntask= %d | rank= %d | size= %d",
98 sprintf(filename,
"%s/%s", dirname, argv[2]);
105 sprintf(filename,
"%s/%s", dirname, argv[3]);
107 ERRMSG(
"Cannot open file!");
129 if (ctl->
dt_mod > fabs(met0->
lon[1] - met0->
lon[0]) * 111132. / 150.)
130 WARN(
"Violation of CFL criterion! Check DT_MOD!");
135 static int dd_init_flg = 0;
136 if (t == ctl->
t_start || !dd_init_flg)
137 dd_init(ctl, &mpi_info, atm, &met0, t, &dd_init_flg);
159 LOG(1,
"SIZE_NP = %d", atm->
np);
160 LOG(1,
"SIZE_MPI_TASKS = %d", size);
161 LOG(1,
"SIZE_OMP_THREADS = %d", omp_get_max_threads());
164 LOG(1,
"MEMORY_ATM = %g MByte",
sizeof(
atm_t) / 1024. / 1024.);
165 LOG(1,
"MEMORY_CACHE = %g MByte",
sizeof(
cache_t) / 1024. / 1024.);
166 LOG(1,
"MEMORY_CLIM = %g MByte",
sizeof(
clim_t) / 1024. / 1024.);
167 LOG(1,
"MEMORY_METEO = %g MByte",
sizeof(
met_t) / 1024. / 1024.);
171 mptrac_free(ctl, cache, clim, met0, met1, atm, &mpi_info);
void mptrac_alloc(ctl_t **ctl, cache_t **cache, clim_t **clim, met_t **met0, met_t **met1, atm_t **atm)
Allocates and initializes memory resources for MPTRAC.
void mptrac_get_met(ctl_t *ctl, clim_t *clim, const double t, met_t **met0, met_t **met1)
Retrieves meteorological data for the specified time.
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_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_free(ctl_t *ctl, cache_t *cache, clim_t *clim, met_t *met0, met_t *met1, atm_t *atm)
Frees memory resources allocated 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)
Executes a single timestep of the MPTRAC model simulation.
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.
MPTRAC library declarations.
#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).
double t_stop
Stop time of simulation [s].
double dt_mod
Time step of simulation [s].
double t_start
Start time of simulation [s].
double lon[EX]
Longitudes [deg].
int main(int argc, char *argv[])