file_1D_out.c

浏览该文件的文档.

 
#include <math.h>
#include <string.h>
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
 
#include "../include/var_struc.h"
#include "../include/file_io.h"
 
 
#define PRINT_NC(v, v_print)                        \
    do {                                \
    strcpy(file_data, add_out);                 \
    strcat(file_data, #v);                      \
    strcat(file_data, ".dat");                  \
    if((fp_write = fopen(file_data, "w")) == NULL)          \
        {                               \
        printf("Cannot open solution output file: %s!\n", #v);  \
        exit(1);                        \
        }                               \
    for(k = 0; k < N; ++k)                      \
        {                               \
        for(j = 0; j < m; ++j)                  \
            fprintf(fp_write, "%.10g\t", (v_print));        \
        fprintf(fp_write, "\n");                \
        }                               \
    fclose(fp_write);                       \
    } while (0)
 
void file_1D_write(const int m, const int N, const struct cell_var_stru CV, 
                    double * X[], const double * cpu_time, const char * problem, const double time_plot[])
{
  // Records the time when the program is running.
  /*
  struct tm * local_time;
  time_t t;
  t=time(NULL);
  local_time=localtime(&t);
  char str_time[100];
  sprintf(str_time, "_%02d%02d%02d%02d%02d%02d", local_time->tm_year-100, local_time->tm_mon+1, local_time->tm_mday, local_time->tm_hour, local_time->tm_min, local_time->tm_sec);
  */
    char add_out[FILENAME_MAX+40];
    // Get the address of the output data folder of the test example.
    example_io(problem, add_out, 0);
    
    char file_data[FILENAME_MAX+40];
    FILE * fp_write;
 
//===================Write Output Data File=========================
 
    int k, j;
    PRINT_NC(RHO, CV.RHO[k][j]);
    PRINT_NC(U,   CV.U[k][j]);
    PRINT_NC(P,   CV.P[k][j]);
    PRINT_NC(E,   CV.E[k][j]);
#ifdef RADIAL_BASICS
    PRINT_NC(R,   X[k][j]);
#else
    PRINT_NC(X, 0.5 * (X[k][j] + X[k][j+1]));
#endif
 
    strcpy(file_data, add_out);
    strcat(file_data, "time_plot.dat");
    printf("%s\n",file_data);
    if((fp_write = fopen(file_data, "w")) == NULL)
    {
        printf("Cannot open solution output file: time_plot!\n");
        exit(1);
    }
    for(k = 0; k < N; ++k)
    fprintf(fp_write, "%.10g\n", time_plot[k]);
    fclose(fp_write);
 
//======================Write Log File============================
    config_write(add_out, cpu_time, problem);
}