1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
|
#include <time.h>
#include <stdio.h>
#include<string.h>
#include "cmor.h"
#include <stdlib.h>
#include <math.h>
void read_time(it, time, time_bnds)
int it;
double time[];
double time_bnds[];
{
time[0] = (it - 0.5) * 30;
time_bnds[0] = (it - 1) * 30;
time_bnds[1] = it * 30;
}
#include "reader_2D_3D.h"
void read_coords(alats, alons, plevs, bnds_lat, bnds_lon, lon, lat, lev)
double *alats, *alons;
int *plevs;
double *bnds_lat, *bnds_lon;
int lon, lat, lev;
{
int i;
for (i = 0; i < lon; i++) {
alons[i] = i * 360. / lon;
bnds_lon[2 * i] = (i - 0.5) * 360. / lon;
bnds_lon[2 * i + 1] = (i + 0.5) * 360. / lon;
};
for (i = 0; i < lat; i++) {
alats[i] = (lat - i) * 10;
bnds_lat[2 * i] = (lat - i) * 10 + 5.;
bnds_lat[2 * i + 1] = (lat - i) * 10 - 5.;
};
for (i = 0; i < lev; i++) {
plevs[i] = (i + 1) * 100;
}
}
int main()
{
/* dimension parameters: */
/* --------------------------------- */
#define ntimes 2 /* number of time samples to process */
#define lon 3 /* number of longitude grid cells */
#define lat 4 /* number of latitude grid cells */
#define lev 4 /* number of standard pressure levels */
double iplevs[lev];
double lon_coords[lon];
double lat_coords[lat];
double lon_bounds[lon * 2];
double lat_bounds[lat * 2];
double data2d[lat * lon];
double data3d[lev * lat * lon];
int myvars[10];
int zfactor_id;
int tables[4];
int axes_ids[CMOR_MAX_DIMENSIONS];
int i, j, k;
int ierr = 0;
double Time[ntimes];
double bnds_time[ntimes * 2];
double tolerance = 1.e-4;
double lon0 = 280.;
double lat0 = 0.;
double delta_lon = 10.;
double delta_lat = 10.;
char id[CMOR_MAX_STRING];
double tmpf = 0.;
char returnvalue[CMOR_MAX_STRING];
printf("Testing variable with both singleton and generic dimensions.\n");
int exit_mode = CMOR_EXIT_ON_MAJOR;
int cmor_mode = CMOR_REPLACE;
ierr |= cmor_setup(NULL, &cmor_mode, NULL, &exit_mode, NULL, NULL);
ierr |= cmor_dataset_json("Test/CMOR_input_example.json");
ierr |= cmor_load_table("Tables/CMIP6_Emon.json", &tables[1]);
for (i = 0; i < ntimes; i++)
read_time(i, &Time[i], &bnds_time[2 * i]);
ierr |= cmor_axis(&axes_ids[0], "time", "days since 2018", ntimes, &Time[0], 'd',
&bnds_time[0], 2, NULL);
read_coords(&lat_coords[0], &lon_coords[0], &iplevs[0], &lat_bounds[0], &lon_bounds[0],
lon, lat, lev);
ierr |= cmor_axis(&axes_ids[1], "latitude", "degrees_north", lat, lat_coords, 'd',
lat_bounds, 2, "");
ierr |= cmor_axis(&axes_ids[2], "longitude", "degrees_east", lon, lon_coords, 'd',
lon_bounds, 2, "");
char landUse[4][CMOR_MAX_STRING];
strncpy(landUse[0], "primary_and_secondary_land", CMOR_MAX_STRING);
strncpy(landUse[1], "pastures", CMOR_MAX_STRING);
strncpy(landUse[2], "crops", CMOR_MAX_STRING);
strncpy(landUse[3], "urban", CMOR_MAX_STRING);
ierr |= cmor_axis(&axes_ids[3], "landUse", "", 4, landUse, 'c',
NULL, CMOR_MAX_STRING, "");
ierr |= cmor_variable(&myvars[0], "nwdFracLut", "%", 4, axes_ids, 'd', NULL,
&tolerance, NULL, "nwdFracLut", "no history", "no future");
for (i = 0; i < ntimes; i++) {
read_3d_input_files(i, "CLOUD", &data3d[0], lat, lon, lev);
ierr |= cmor_write(myvars[0], data3d, 'd', NULL, 1, NULL, NULL, NULL);
}
ierr |= cmor_close_variable(myvars[0], returnvalue, NULL);
printf("File: '%s' has been written\n", returnvalue);
if(cmor_nerrors != 0) {
printf("Error occured in test_cmor_nwdFracLut.\n");
exit(1);
}
ierr |= cmor_close();
return ierr;
}
|
