sub_period_common.c
Go to the documentation of this file.00001
00002
00003
00004
00005
00006
00007
00012
00013
00014
00015
00016
00017
00018
00019
00020
00021
00022
00023
00024
00025
00026
00027
00028
00029
00030
00031
00032
00033
00034
00035
00036
00037
00038
00039
00040
00041
00042
00043
00044
00045
00046
00047
00048
00049
00050
00051
00052
00053
00054
00055 #include <utils.h>
00056
00058 int
00059 sub_period_common(double **buf_sub, int *ntime_sub, double *bufin, int *year, int *month, int *day,
00060 int *year_learn, int *month_learn, int *day_learn, int timedim, int ndima, int ndimb, int ntime, int ntime_learn) {
00078 int *buf_sub_i = NULL;
00079
00080 int dima;
00081 int dimb;
00082 int t;
00083 int tt;
00084
00085
00086 *ntime_sub = 0;
00087
00088
00089 for (t=0; t<ntime; t++) {
00090
00091 for (tt=0; tt<ntime_learn; tt++) {
00092 if (year[t] == year_learn[tt] &&
00093 month[t] == month_learn[tt] &&
00094 day[t] == day_learn[tt]) {
00095
00096 buf_sub_i = (int *) realloc(buf_sub_i, ((*ntime_sub)+1) * sizeof(int));
00097 if (buf_sub_i == NULL) alloc_error(__FILE__, __LINE__);
00098 buf_sub_i[(*ntime_sub)++] = t;
00099 }
00100 }
00101 }
00102
00103 if ( (*ntime_sub) == 0 ) {
00104 (void) fprintf(stderr, "%s: FATAL ERROR: No common subperiod! Maybe a problem in the time representation in the control run file.\nAborting.\n", __FILE__);
00105 (void) printf("MODEL TIMES ntime=%d\n", ntime);
00106
00107 for (t=0; t<ntime; t++)
00108 (void) printf("%d %d %d\n", year[t], month[t], day[t]);
00109 (void) printf("LEARNING TIMES ntime=%d\n", ntime_learn);
00110 for (t=0; t<ntime_learn; t++)
00111 (void) printf("%d %d %d\n", year_learn[t], month_learn[t], day_learn[t]);
00112
00113 return -1;
00114 }
00115
00116 (void) printf("%s: Sub-period: %d %d %d %d %d %d. Indexes: %d %d\n",__FILE__, year[buf_sub_i[0]], month[buf_sub_i[0]],
00117 day[buf_sub_i[0]], year[buf_sub_i[(*ntime_sub)-1]],month[buf_sub_i[(*ntime_sub)-1]],
00118 day[buf_sub_i[(*ntime_sub)-1]], buf_sub_i[0], buf_sub_i[(*ntime_sub)-1]);
00119
00120
00121 (*buf_sub) = (double *) malloc((*ntime_sub)*ndima*ndimb * sizeof(double));
00122 if ((*buf_sub) == NULL) alloc_error(__FILE__, __LINE__);
00123
00124 if (timedim == 3)
00125
00126 for (t=0; t<(*ntime_sub); t++)
00127 for (dimb=0; dimb<ndimb; dimb++)
00128 for (dima=0; dima<ndima; dima++)
00129 (*buf_sub)[dima+(dimb*ndima)+t*ndima*ndimb] = bufin[dima+dimb*ndima+buf_sub_i[t]*ndima*ndimb];
00130 else if (timedim == 1)
00131
00132 for (t=0; t<(*ntime_sub); t++)
00133 for (dimb=0; dimb<ndimb; dimb++)
00134 for (dima=0; dima<ndima; dima++)
00135 (*buf_sub)[t+dima*(*ntime_sub)+dimb*(ndima*(*ntime_sub))] = bufin[buf_sub_i[t]+dima*ntime+dimb*(ndima*ntime)];
00136 else
00137 (void) fprintf(stderr, "%s: Fatal error: timedim argument must be equal to 1 or 3.\n", __FILE__);
00138
00139
00140 (void) free(buf_sub_i);
00141
00142
00143 return 0;
00144 }