forked from Unidata/netcdf-c
/
tst_atts.c
2495 lines (2192 loc) · 80.5 KB
/
tst_atts.c
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
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
686
687
688
689
690
691
692
693
694
695
696
697
698
699
700
701
702
703
704
705
706
707
708
709
710
711
712
713
714
715
716
717
718
719
720
721
722
723
724
725
726
727
728
729
730
731
732
733
734
735
736
737
738
739
740
741
742
743
744
745
746
747
748
749
750
751
752
753
754
755
756
757
758
759
760
761
762
763
764
765
766
767
768
769
770
771
772
773
774
775
776
777
778
779
780
781
782
783
784
785
786
787
788
789
790
791
792
793
794
795
796
797
798
799
800
801
802
803
804
805
806
807
808
809
810
811
812
813
814
815
816
817
818
819
820
821
822
823
824
825
826
827
828
829
830
831
832
833
834
835
836
837
838
839
840
841
842
843
844
845
846
847
848
849
850
851
852
853
854
855
856
857
858
859
860
861
862
863
864
865
866
867
868
869
870
871
872
873
874
875
876
877
878
879
880
881
882
883
884
885
886
887
888
889
890
891
892
893
894
895
896
897
898
899
900
901
902
903
904
905
906
907
908
909
910
911
912
913
914
915
916
917
918
919
920
921
922
923
924
925
926
927
928
929
930
931
932
933
934
935
936
937
938
939
940
941
942
943
944
945
946
947
948
949
950
951
952
953
954
955
956
957
958
959
960
961
962
963
964
965
966
967
968
969
970
971
972
973
974
975
976
977
978
979
980
981
982
983
984
985
986
987
988
989
990
991
992
993
994
995
996
997
998
999
1000
/* This is part of the netCDF package.
Copyright 2006 University Corporation for Atmospheric Research/Unidata.
See COPYRIGHT file for conditions of use.
This is a very simple example which writes a netCDF file with
Unicode names encoded with UTF-8. It is the NETCDF3 equivalent
of tst_unicode.c
$Id: tst_atts.c 2792 2014-10-27 06:02:59Z wkliao $
*/
#include <nc_tests.h>
#include <stdlib.h>
#include <string.h>
/* The data file we will create. */
#define FILE_NAME "tst_atts.nc"
void
check_err(const int stat, const int line, const char *file) {
if (stat != NC_NOERR) {
(void)fprintf(stderr,"line %d of %s: %s\n", line, file, nc_strerror(stat));
fflush(stderr);
exit(1);
}
}
int
create_file()
{
int stat; /* return status */
int ncid; /* netCDF id */
/* dimension ids */
int Dr_dim;
int D1_dim;
int D2_dim;
int D3_dim;
int D4_dim;
/* dimension lengths */
size_t Dr_len = NC_UNLIMITED;
size_t D1_len = 1;
size_t D2_len = 2;
size_t D3_len = 3;
size_t D4_len = 4;
/* variable ids */
int c_id;
int b_id;
int s_id;
int i_id;
int f_id;
int d_id;
int cr_id;
int br_id;
int sr_id;
int ir_id;
int fr_id;
int dr_id;
int c1_id;
int b1_id;
int s1_id;
int i1_id;
int f1_id;
int d1_id;
int c2_id;
int b2_id;
int s2_id;
int i2_id;
int f2_id;
int d2_id;
int c3_id;
int b3_id;
int s3_id;
int i3_id;
int f3_id;
int d3_id;
int c4_id;
int b4_id;
int s4_id;
int i4_id;
int f4_id;
int d4_id;
int cr1_id;
int br2_id;
int sr3_id;
int ir4_id;
int f11_id;
int d12_id;
int c13_id;
int b14_id;
int s21_id;
int i22_id;
int f23_id;
int d24_id;
int c31_id;
int b32_id;
int s33_id;
int i34_id;
int f41_id;
int d42_id;
int c43_id;
int b44_id;
int sr11_id;
int ir12_id;
int fr13_id;
int dr14_id;
int cr21_id;
int br22_id;
int sr23_id;
int ir24_id;
int fr31_id;
int dr32_id;
int cr33_id;
int br34_id;
int sr41_id;
int ir42_id;
int fr43_id;
int dr44_id;
int c111_id;
int b112_id;
int s113_id;
int i114_id;
int f121_id;
int d122_id;
int c123_id;
int b124_id;
int s131_id;
int i132_id;
int f133_id;
int d134_id;
int c141_id;
int b142_id;
int s143_id;
int i144_id;
int f211_id;
int d212_id;
int c213_id;
int b214_id;
int s221_id;
int i222_id;
int f223_id;
int d224_id;
int c231_id;
int b232_id;
int s233_id;
int i234_id;
int f241_id;
int d242_id;
int c243_id;
int b244_id;
int s311_id;
int i312_id;
int f313_id;
int d314_id;
int c321_id;
int b322_id;
int s323_id;
int i324_id;
int f331_id;
int d332_id;
int c333_id;
int b334_id;
int s341_id;
int i342_id;
int f343_id;
int d344_id;
int c411_id;
int b412_id;
int s413_id;
int i414_id;
int f421_id;
int d422_id;
int c423_id;
int b424_id;
int s431_id;
int i432_id;
int f433_id;
int d434_id;
int c441_id;
int b442_id;
int s443_id;
int i444_id;
/* rank (number of dimensions) for each variable */
# define RANK_c 0
# define RANK_b 0
# define RANK_s 0
# define RANK_i 0
# define RANK_f 0
# define RANK_d 0
# define RANK_cr 1
# define RANK_br 1
# define RANK_sr 1
# define RANK_ir 1
# define RANK_fr 1
# define RANK_dr 1
# define RANK_c1 1
# define RANK_b1 1
# define RANK_s1 1
# define RANK_i1 1
# define RANK_f1 1
# define RANK_d1 1
# define RANK_c2 1
# define RANK_b2 1
# define RANK_s2 1
# define RANK_i2 1
# define RANK_f2 1
# define RANK_d2 1
# define RANK_c3 1
# define RANK_b3 1
# define RANK_s3 1
# define RANK_i3 1
# define RANK_f3 1
# define RANK_d3 1
# define RANK_c4 1
# define RANK_b4 1
# define RANK_s4 1
# define RANK_i4 1
# define RANK_f4 1
# define RANK_d4 1
# define RANK_cr1 2
# define RANK_br2 2
# define RANK_sr3 2
# define RANK_ir4 2
# define RANK_f11 2
# define RANK_d12 2
# define RANK_c13 2
# define RANK_b14 2
# define RANK_s21 2
# define RANK_i22 2
# define RANK_f23 2
# define RANK_d24 2
# define RANK_c31 2
# define RANK_b32 2
# define RANK_s33 2
# define RANK_i34 2
# define RANK_f41 2
# define RANK_d42 2
# define RANK_c43 2
# define RANK_b44 2
# define RANK_sr11 3
# define RANK_ir12 3
# define RANK_fr13 3
# define RANK_dr14 3
# define RANK_cr21 3
# define RANK_br22 3
# define RANK_sr23 3
# define RANK_ir24 3
# define RANK_fr31 3
# define RANK_dr32 3
# define RANK_cr33 3
# define RANK_br34 3
# define RANK_sr41 3
# define RANK_ir42 3
# define RANK_fr43 3
# define RANK_dr44 3
# define RANK_c111 3
# define RANK_b112 3
# define RANK_s113 3
# define RANK_i114 3
# define RANK_f121 3
# define RANK_d122 3
# define RANK_c123 3
# define RANK_b124 3
# define RANK_s131 3
# define RANK_i132 3
# define RANK_f133 3
# define RANK_d134 3
# define RANK_c141 3
# define RANK_b142 3
# define RANK_s143 3
# define RANK_i144 3
# define RANK_f211 3
# define RANK_d212 3
# define RANK_c213 3
# define RANK_b214 3
# define RANK_s221 3
# define RANK_i222 3
# define RANK_f223 3
# define RANK_d224 3
# define RANK_c231 3
# define RANK_b232 3
# define RANK_s233 3
# define RANK_i234 3
# define RANK_f241 3
# define RANK_d242 3
# define RANK_c243 3
# define RANK_b244 3
# define RANK_s311 3
# define RANK_i312 3
# define RANK_f313 3
# define RANK_d314 3
# define RANK_c321 3
# define RANK_b322 3
# define RANK_s323 3
# define RANK_i324 3
# define RANK_f331 3
# define RANK_d332 3
# define RANK_c333 3
# define RANK_b334 3
# define RANK_s341 3
# define RANK_i342 3
# define RANK_f343 3
# define RANK_d344 3
# define RANK_c411 3
# define RANK_b412 3
# define RANK_s413 3
# define RANK_i414 3
# define RANK_f421 3
# define RANK_d422 3
# define RANK_c423 3
# define RANK_b424 3
# define RANK_s431 3
# define RANK_i432 3
# define RANK_f433 3
# define RANK_d434 3
# define RANK_c441 3
# define RANK_b442 3
# define RANK_s443 3
# define RANK_i444 3
/* variable shapes */
int cr_dims[RANK_cr];
int br_dims[RANK_br];
int sr_dims[RANK_sr];
int ir_dims[RANK_ir];
int fr_dims[RANK_fr];
int dr_dims[RANK_dr];
int c1_dims[RANK_c1];
int b1_dims[RANK_b1];
int s1_dims[RANK_s1];
int i1_dims[RANK_i1];
int f1_dims[RANK_f1];
int d1_dims[RANK_d1];
int c2_dims[RANK_c2];
int b2_dims[RANK_b2];
int s2_dims[RANK_s2];
int i2_dims[RANK_i2];
int f2_dims[RANK_f2];
int d2_dims[RANK_d2];
int c3_dims[RANK_c3];
int b3_dims[RANK_b3];
int s3_dims[RANK_s3];
int i3_dims[RANK_i3];
int f3_dims[RANK_f3];
int d3_dims[RANK_d3];
int c4_dims[RANK_c4];
int b4_dims[RANK_b4];
int s4_dims[RANK_s4];
int i4_dims[RANK_i4];
int f4_dims[RANK_f4];
int d4_dims[RANK_d4];
int cr1_dims[RANK_cr1];
int br2_dims[RANK_br2];
int sr3_dims[RANK_sr3];
int ir4_dims[RANK_ir4];
int f11_dims[RANK_f11];
int d12_dims[RANK_d12];
int c13_dims[RANK_c13];
int b14_dims[RANK_b14];
int s21_dims[RANK_s21];
int i22_dims[RANK_i22];
int f23_dims[RANK_f23];
int d24_dims[RANK_d24];
int c31_dims[RANK_c31];
int b32_dims[RANK_b32];
int s33_dims[RANK_s33];
int i34_dims[RANK_i34];
int f41_dims[RANK_f41];
int d42_dims[RANK_d42];
int c43_dims[RANK_c43];
int b44_dims[RANK_b44];
int sr11_dims[RANK_sr11];
int ir12_dims[RANK_ir12];
int fr13_dims[RANK_fr13];
int dr14_dims[RANK_dr14];
int cr21_dims[RANK_cr21];
int br22_dims[RANK_br22];
int sr23_dims[RANK_sr23];
int ir24_dims[RANK_ir24];
int fr31_dims[RANK_fr31];
int dr32_dims[RANK_dr32];
int cr33_dims[RANK_cr33];
int br34_dims[RANK_br34];
int sr41_dims[RANK_sr41];
int ir42_dims[RANK_ir42];
int fr43_dims[RANK_fr43];
int dr44_dims[RANK_dr44];
int c111_dims[RANK_c111];
int b112_dims[RANK_b112];
int s113_dims[RANK_s113];
int i114_dims[RANK_i114];
int f121_dims[RANK_f121];
int d122_dims[RANK_d122];
int c123_dims[RANK_c123];
int b124_dims[RANK_b124];
int s131_dims[RANK_s131];
int i132_dims[RANK_i132];
int f133_dims[RANK_f133];
int d134_dims[RANK_d134];
int c141_dims[RANK_c141];
int b142_dims[RANK_b142];
int s143_dims[RANK_s143];
int i144_dims[RANK_i144];
int f211_dims[RANK_f211];
int d212_dims[RANK_d212];
int c213_dims[RANK_c213];
int b214_dims[RANK_b214];
int s221_dims[RANK_s221];
int i222_dims[RANK_i222];
int f223_dims[RANK_f223];
int d224_dims[RANK_d224];
int c231_dims[RANK_c231];
int b232_dims[RANK_b232];
int s233_dims[RANK_s233];
int i234_dims[RANK_i234];
int f241_dims[RANK_f241];
int d242_dims[RANK_d242];
int c243_dims[RANK_c243];
int b244_dims[RANK_b244];
int s311_dims[RANK_s311];
int i312_dims[RANK_i312];
int f313_dims[RANK_f313];
int d314_dims[RANK_d314];
int c321_dims[RANK_c321];
int b322_dims[RANK_b322];
int s323_dims[RANK_s323];
int i324_dims[RANK_i324];
int f331_dims[RANK_f331];
int d332_dims[RANK_d332];
int c333_dims[RANK_c333];
int b334_dims[RANK_b334];
int s341_dims[RANK_s341];
int i342_dims[RANK_i342];
int f343_dims[RANK_f343];
int d344_dims[RANK_d344];
int c411_dims[RANK_c411];
int b412_dims[RANK_b412];
int s413_dims[RANK_s413];
int i414_dims[RANK_i414];
int f421_dims[RANK_f421];
int d422_dims[RANK_d422];
int c423_dims[RANK_c423];
int b424_dims[RANK_b424];
int s431_dims[RANK_s431];
int i432_dims[RANK_i432];
int f433_dims[RANK_f433];
int d434_dims[RANK_d434];
int c441_dims[RANK_c441];
int b442_dims[RANK_b442];
int s443_dims[RANK_s443];
int i444_dims[RANK_i444];
/* enter define mode */
#ifdef TEST_PNETCDF
stat = nc_create_par(FILE_NAME, NC_CLOBBER|NC_PNETCDF, MPI_COMM_WORLD, MPI_INFO_NULL, &ncid);
#else
stat = nc_create(FILE_NAME, NC_CLOBBER, &ncid);
#endif
check_err(stat,__LINE__,__FILE__);
/* define dimensions */
stat = nc_def_dim(ncid, "Dr", Dr_len, &Dr_dim);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_dim(ncid, "D1", D1_len, &D1_dim);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_dim(ncid, "D2", D2_len, &D2_dim);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_dim(ncid, "D3", D3_len, &D3_dim);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_dim(ncid, "D4", D4_len, &D4_dim);
check_err(stat,__LINE__,__FILE__);
/* define variables */
stat = nc_def_var(ncid, "c", NC_CHAR, RANK_c, 0, &c_id);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_var(ncid, "b", NC_BYTE, RANK_b, 0, &b_id);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_var(ncid, "s", NC_SHORT, RANK_s, 0, &s_id);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_var(ncid, "i", NC_INT, RANK_i, 0, &i_id);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_var(ncid, "f", NC_FLOAT, RANK_f, 0, &f_id);
check_err(stat,__LINE__,__FILE__);
stat = nc_def_var(ncid, "d", NC_DOUBLE, RANK_d, 0, &d_id);
check_err(stat,__LINE__,__FILE__);
cr_dims[0] = Dr_dim;
stat = nc_def_var(ncid, "cr", NC_CHAR, RANK_cr, cr_dims, &cr_id);
check_err(stat,__LINE__,__FILE__);
br_dims[0] = Dr_dim;
stat = nc_def_var(ncid, "br", NC_BYTE, RANK_br, br_dims, &br_id);
check_err(stat,__LINE__,__FILE__);
sr_dims[0] = Dr_dim;
stat = nc_def_var(ncid, "sr", NC_SHORT, RANK_sr, sr_dims, &sr_id);
check_err(stat,__LINE__,__FILE__);
ir_dims[0] = Dr_dim;
stat = nc_def_var(ncid, "ir", NC_INT, RANK_ir, ir_dims, &ir_id);
check_err(stat,__LINE__,__FILE__);
fr_dims[0] = Dr_dim;
stat = nc_def_var(ncid, "fr", NC_FLOAT, RANK_fr, fr_dims, &fr_id);
check_err(stat,__LINE__,__FILE__);
dr_dims[0] = Dr_dim;
stat = nc_def_var(ncid, "dr", NC_DOUBLE, RANK_dr, dr_dims, &dr_id);
check_err(stat,__LINE__,__FILE__);
c1_dims[0] = D1_dim;
stat = nc_def_var(ncid, "c1", NC_CHAR, RANK_c1, c1_dims, &c1_id);
check_err(stat,__LINE__,__FILE__);
b1_dims[0] = D1_dim;
stat = nc_def_var(ncid, "b1", NC_BYTE, RANK_b1, b1_dims, &b1_id);
check_err(stat,__LINE__,__FILE__);
s1_dims[0] = D1_dim;
stat = nc_def_var(ncid, "s1", NC_SHORT, RANK_s1, s1_dims, &s1_id);
check_err(stat,__LINE__,__FILE__);
i1_dims[0] = D1_dim;
stat = nc_def_var(ncid, "i1", NC_INT, RANK_i1, i1_dims, &i1_id);
check_err(stat,__LINE__,__FILE__);
f1_dims[0] = D1_dim;
stat = nc_def_var(ncid, "f1", NC_FLOAT, RANK_f1, f1_dims, &f1_id);
check_err(stat,__LINE__,__FILE__);
d1_dims[0] = D1_dim;
stat = nc_def_var(ncid, "d1", NC_DOUBLE, RANK_d1, d1_dims, &d1_id);
check_err(stat,__LINE__,__FILE__);
c2_dims[0] = D2_dim;
stat = nc_def_var(ncid, "c2", NC_CHAR, RANK_c2, c2_dims, &c2_id);
check_err(stat,__LINE__,__FILE__);
b2_dims[0] = D2_dim;
stat = nc_def_var(ncid, "b2", NC_BYTE, RANK_b2, b2_dims, &b2_id);
check_err(stat,__LINE__,__FILE__);
s2_dims[0] = D2_dim;
stat = nc_def_var(ncid, "s2", NC_SHORT, RANK_s2, s2_dims, &s2_id);
check_err(stat,__LINE__,__FILE__);
i2_dims[0] = D2_dim;
stat = nc_def_var(ncid, "i2", NC_INT, RANK_i2, i2_dims, &i2_id);
check_err(stat,__LINE__,__FILE__);
f2_dims[0] = D2_dim;
stat = nc_def_var(ncid, "f2", NC_FLOAT, RANK_f2, f2_dims, &f2_id);
check_err(stat,__LINE__,__FILE__);
d2_dims[0] = D2_dim;
stat = nc_def_var(ncid, "d2", NC_DOUBLE, RANK_d2, d2_dims, &d2_id);
check_err(stat,__LINE__,__FILE__);
c3_dims[0] = D3_dim;
stat = nc_def_var(ncid, "c3", NC_CHAR, RANK_c3, c3_dims, &c3_id);
check_err(stat,__LINE__,__FILE__);
b3_dims[0] = D3_dim;
stat = nc_def_var(ncid, "b3", NC_BYTE, RANK_b3, b3_dims, &b3_id);
check_err(stat,__LINE__,__FILE__);
s3_dims[0] = D3_dim;
stat = nc_def_var(ncid, "s3", NC_SHORT, RANK_s3, s3_dims, &s3_id);
check_err(stat,__LINE__,__FILE__);
i3_dims[0] = D3_dim;
stat = nc_def_var(ncid, "i3", NC_INT, RANK_i3, i3_dims, &i3_id);
check_err(stat,__LINE__,__FILE__);
f3_dims[0] = D3_dim;
stat = nc_def_var(ncid, "f3", NC_FLOAT, RANK_f3, f3_dims, &f3_id);
check_err(stat,__LINE__,__FILE__);
d3_dims[0] = D3_dim;
stat = nc_def_var(ncid, "d3", NC_DOUBLE, RANK_d3, d3_dims, &d3_id);
check_err(stat,__LINE__,__FILE__);
c4_dims[0] = D4_dim;
stat = nc_def_var(ncid, "c4", NC_CHAR, RANK_c4, c4_dims, &c4_id);
check_err(stat,__LINE__,__FILE__);
b4_dims[0] = D4_dim;
stat = nc_def_var(ncid, "b4", NC_BYTE, RANK_b4, b4_dims, &b4_id);
check_err(stat,__LINE__,__FILE__);
s4_dims[0] = D4_dim;
stat = nc_def_var(ncid, "s4", NC_SHORT, RANK_s4, s4_dims, &s4_id);
check_err(stat,__LINE__,__FILE__);
i4_dims[0] = D4_dim;
stat = nc_def_var(ncid, "i4", NC_INT, RANK_i4, i4_dims, &i4_id);
check_err(stat,__LINE__,__FILE__);
f4_dims[0] = D4_dim;
stat = nc_def_var(ncid, "f4", NC_FLOAT, RANK_f4, f4_dims, &f4_id);
check_err(stat,__LINE__,__FILE__);
d4_dims[0] = D4_dim;
stat = nc_def_var(ncid, "d4", NC_DOUBLE, RANK_d4, d4_dims, &d4_id);
check_err(stat,__LINE__,__FILE__);
cr1_dims[0] = Dr_dim;
cr1_dims[1] = D1_dim;
stat = nc_def_var(ncid, "cr1", NC_CHAR, RANK_cr1, cr1_dims, &cr1_id);
check_err(stat,__LINE__,__FILE__);
br2_dims[0] = Dr_dim;
br2_dims[1] = D2_dim;
stat = nc_def_var(ncid, "br2", NC_BYTE, RANK_br2, br2_dims, &br2_id);
check_err(stat,__LINE__,__FILE__);
sr3_dims[0] = Dr_dim;
sr3_dims[1] = D3_dim;
stat = nc_def_var(ncid, "sr3", NC_SHORT, RANK_sr3, sr3_dims, &sr3_id);
check_err(stat,__LINE__,__FILE__);
ir4_dims[0] = Dr_dim;
ir4_dims[1] = D4_dim;
stat = nc_def_var(ncid, "ir4", NC_INT, RANK_ir4, ir4_dims, &ir4_id);
check_err(stat,__LINE__,__FILE__);
f11_dims[0] = D1_dim;
f11_dims[1] = D1_dim;
stat = nc_def_var(ncid, "f11", NC_FLOAT, RANK_f11, f11_dims, &f11_id);
check_err(stat,__LINE__,__FILE__);
d12_dims[0] = D1_dim;
d12_dims[1] = D2_dim;
stat = nc_def_var(ncid, "d12", NC_DOUBLE, RANK_d12, d12_dims, &d12_id);
check_err(stat,__LINE__,__FILE__);
c13_dims[0] = D1_dim;
c13_dims[1] = D3_dim;
stat = nc_def_var(ncid, "c13", NC_CHAR, RANK_c13, c13_dims, &c13_id);
check_err(stat,__LINE__,__FILE__);
b14_dims[0] = D1_dim;
b14_dims[1] = D4_dim;
stat = nc_def_var(ncid, "b14", NC_BYTE, RANK_b14, b14_dims, &b14_id);
check_err(stat,__LINE__,__FILE__);
s21_dims[0] = D2_dim;
s21_dims[1] = D1_dim;
stat = nc_def_var(ncid, "s21", NC_SHORT, RANK_s21, s21_dims, &s21_id);
check_err(stat,__LINE__,__FILE__);
i22_dims[0] = D2_dim;
i22_dims[1] = D2_dim;
stat = nc_def_var(ncid, "i22", NC_INT, RANK_i22, i22_dims, &i22_id);
check_err(stat,__LINE__,__FILE__);
f23_dims[0] = D2_dim;
f23_dims[1] = D3_dim;
stat = nc_def_var(ncid, "f23", NC_FLOAT, RANK_f23, f23_dims, &f23_id);
check_err(stat,__LINE__,__FILE__);
d24_dims[0] = D2_dim;
d24_dims[1] = D4_dim;
stat = nc_def_var(ncid, "d24", NC_DOUBLE, RANK_d24, d24_dims, &d24_id);
check_err(stat,__LINE__,__FILE__);
c31_dims[0] = D3_dim;
c31_dims[1] = D1_dim;
stat = nc_def_var(ncid, "c31", NC_CHAR, RANK_c31, c31_dims, &c31_id);
check_err(stat,__LINE__,__FILE__);
b32_dims[0] = D3_dim;
b32_dims[1] = D2_dim;
stat = nc_def_var(ncid, "b32", NC_BYTE, RANK_b32, b32_dims, &b32_id);
check_err(stat,__LINE__,__FILE__);
s33_dims[0] = D3_dim;
s33_dims[1] = D3_dim;
stat = nc_def_var(ncid, "s33", NC_SHORT, RANK_s33, s33_dims, &s33_id);
check_err(stat,__LINE__,__FILE__);
i34_dims[0] = D3_dim;
i34_dims[1] = D4_dim;
stat = nc_def_var(ncid, "i34", NC_INT, RANK_i34, i34_dims, &i34_id);
check_err(stat,__LINE__,__FILE__);
f41_dims[0] = D4_dim;
f41_dims[1] = D1_dim;
stat = nc_def_var(ncid, "f41", NC_FLOAT, RANK_f41, f41_dims, &f41_id);
check_err(stat,__LINE__,__FILE__);
d42_dims[0] = D4_dim;
d42_dims[1] = D2_dim;
stat = nc_def_var(ncid, "d42", NC_DOUBLE, RANK_d42, d42_dims, &d42_id);
check_err(stat,__LINE__,__FILE__);
c43_dims[0] = D4_dim;
c43_dims[1] = D3_dim;
stat = nc_def_var(ncid, "c43", NC_CHAR, RANK_c43, c43_dims, &c43_id);
check_err(stat,__LINE__,__FILE__);
b44_dims[0] = D4_dim;
b44_dims[1] = D4_dim;
stat = nc_def_var(ncid, "b44", NC_BYTE, RANK_b44, b44_dims, &b44_id);
check_err(stat,__LINE__,__FILE__);
sr11_dims[0] = Dr_dim;
sr11_dims[1] = D1_dim;
sr11_dims[2] = D1_dim;
stat = nc_def_var(ncid, "sr11", NC_SHORT, RANK_sr11, sr11_dims, &sr11_id);
check_err(stat,__LINE__,__FILE__);
ir12_dims[0] = Dr_dim;
ir12_dims[1] = D1_dim;
ir12_dims[2] = D2_dim;
stat = nc_def_var(ncid, "ir12", NC_INT, RANK_ir12, ir12_dims, &ir12_id);
check_err(stat,__LINE__,__FILE__);
fr13_dims[0] = Dr_dim;
fr13_dims[1] = D1_dim;
fr13_dims[2] = D3_dim;
stat = nc_def_var(ncid, "fr13", NC_FLOAT, RANK_fr13, fr13_dims, &fr13_id);
check_err(stat,__LINE__,__FILE__);
dr14_dims[0] = Dr_dim;
dr14_dims[1] = D1_dim;
dr14_dims[2] = D4_dim;
stat = nc_def_var(ncid, "dr14", NC_DOUBLE, RANK_dr14, dr14_dims, &dr14_id);
check_err(stat,__LINE__,__FILE__);
cr21_dims[0] = Dr_dim;
cr21_dims[1] = D2_dim;
cr21_dims[2] = D1_dim;
stat = nc_def_var(ncid, "cr21", NC_CHAR, RANK_cr21, cr21_dims, &cr21_id);
check_err(stat,__LINE__,__FILE__);
br22_dims[0] = Dr_dim;
br22_dims[1] = D2_dim;
br22_dims[2] = D2_dim;
stat = nc_def_var(ncid, "br22", NC_BYTE, RANK_br22, br22_dims, &br22_id);
check_err(stat,__LINE__,__FILE__);
sr23_dims[0] = Dr_dim;
sr23_dims[1] = D2_dim;
sr23_dims[2] = D3_dim;
stat = nc_def_var(ncid, "sr23", NC_SHORT, RANK_sr23, sr23_dims, &sr23_id);
check_err(stat,__LINE__,__FILE__);
ir24_dims[0] = Dr_dim;
ir24_dims[1] = D2_dim;
ir24_dims[2] = D4_dim;
stat = nc_def_var(ncid, "ir24", NC_INT, RANK_ir24, ir24_dims, &ir24_id);
check_err(stat,__LINE__,__FILE__);
fr31_dims[0] = Dr_dim;
fr31_dims[1] = D3_dim;
fr31_dims[2] = D1_dim;
stat = nc_def_var(ncid, "fr31", NC_FLOAT, RANK_fr31, fr31_dims, &fr31_id);
check_err(stat,__LINE__,__FILE__);
dr32_dims[0] = Dr_dim;
dr32_dims[1] = D3_dim;
dr32_dims[2] = D2_dim;
stat = nc_def_var(ncid, "dr32", NC_DOUBLE, RANK_dr32, dr32_dims, &dr32_id);
check_err(stat,__LINE__,__FILE__);
cr33_dims[0] = Dr_dim;
cr33_dims[1] = D3_dim;
cr33_dims[2] = D3_dim;
stat = nc_def_var(ncid, "cr33", NC_CHAR, RANK_cr33, cr33_dims, &cr33_id);
check_err(stat,__LINE__,__FILE__);
br34_dims[0] = Dr_dim;
br34_dims[1] = D3_dim;
br34_dims[2] = D4_dim;
stat = nc_def_var(ncid, "br34", NC_BYTE, RANK_br34, br34_dims, &br34_id);
check_err(stat,__LINE__,__FILE__);
sr41_dims[0] = Dr_dim;
sr41_dims[1] = D4_dim;
sr41_dims[2] = D1_dim;
stat = nc_def_var(ncid, "sr41", NC_SHORT, RANK_sr41, sr41_dims, &sr41_id);
check_err(stat,__LINE__,__FILE__);
ir42_dims[0] = Dr_dim;
ir42_dims[1] = D4_dim;
ir42_dims[2] = D2_dim;
stat = nc_def_var(ncid, "ir42", NC_INT, RANK_ir42, ir42_dims, &ir42_id);
check_err(stat,__LINE__,__FILE__);
fr43_dims[0] = Dr_dim;
fr43_dims[1] = D4_dim;
fr43_dims[2] = D3_dim;
stat = nc_def_var(ncid, "fr43", NC_FLOAT, RANK_fr43, fr43_dims, &fr43_id);
check_err(stat,__LINE__,__FILE__);
dr44_dims[0] = Dr_dim;
dr44_dims[1] = D4_dim;
dr44_dims[2] = D4_dim;
stat = nc_def_var(ncid, "dr44", NC_DOUBLE, RANK_dr44, dr44_dims, &dr44_id);
check_err(stat,__LINE__,__FILE__);
c111_dims[0] = D1_dim;
c111_dims[1] = D1_dim;
c111_dims[2] = D1_dim;
stat = nc_def_var(ncid, "c111", NC_CHAR, RANK_c111, c111_dims, &c111_id);
check_err(stat,__LINE__,__FILE__);
b112_dims[0] = D1_dim;
b112_dims[1] = D1_dim;
b112_dims[2] = D2_dim;
stat = nc_def_var(ncid, "b112", NC_BYTE, RANK_b112, b112_dims, &b112_id);
check_err(stat,__LINE__,__FILE__);
s113_dims[0] = D1_dim;
s113_dims[1] = D1_dim;
s113_dims[2] = D3_dim;
stat = nc_def_var(ncid, "s113", NC_SHORT, RANK_s113, s113_dims, &s113_id);
check_err(stat,__LINE__,__FILE__);
i114_dims[0] = D1_dim;
i114_dims[1] = D1_dim;
i114_dims[2] = D4_dim;
stat = nc_def_var(ncid, "i114", NC_INT, RANK_i114, i114_dims, &i114_id);
check_err(stat,__LINE__,__FILE__);
f121_dims[0] = D1_dim;
f121_dims[1] = D2_dim;
f121_dims[2] = D1_dim;
stat = nc_def_var(ncid, "f121", NC_FLOAT, RANK_f121, f121_dims, &f121_id);
check_err(stat,__LINE__,__FILE__);
d122_dims[0] = D1_dim;
d122_dims[1] = D2_dim;
d122_dims[2] = D2_dim;
stat = nc_def_var(ncid, "d122", NC_DOUBLE, RANK_d122, d122_dims, &d122_id);
check_err(stat,__LINE__,__FILE__);
c123_dims[0] = D1_dim;
c123_dims[1] = D2_dim;
c123_dims[2] = D3_dim;
stat = nc_def_var(ncid, "c123", NC_CHAR, RANK_c123, c123_dims, &c123_id);
check_err(stat,__LINE__,__FILE__);
b124_dims[0] = D1_dim;
b124_dims[1] = D2_dim;
b124_dims[2] = D4_dim;
stat = nc_def_var(ncid, "b124", NC_BYTE, RANK_b124, b124_dims, &b124_id);
check_err(stat,__LINE__,__FILE__);
s131_dims[0] = D1_dim;
s131_dims[1] = D3_dim;
s131_dims[2] = D1_dim;
stat = nc_def_var(ncid, "s131", NC_SHORT, RANK_s131, s131_dims, &s131_id);
check_err(stat,__LINE__,__FILE__);
i132_dims[0] = D1_dim;
i132_dims[1] = D3_dim;
i132_dims[2] = D2_dim;
stat = nc_def_var(ncid, "i132", NC_INT, RANK_i132, i132_dims, &i132_id);
check_err(stat,__LINE__,__FILE__);
f133_dims[0] = D1_dim;
f133_dims[1] = D3_dim;
f133_dims[2] = D3_dim;
stat = nc_def_var(ncid, "f133", NC_FLOAT, RANK_f133, f133_dims, &f133_id);
check_err(stat,__LINE__,__FILE__);
d134_dims[0] = D1_dim;
d134_dims[1] = D3_dim;
d134_dims[2] = D4_dim;
stat = nc_def_var(ncid, "d134", NC_DOUBLE, RANK_d134, d134_dims, &d134_id);
check_err(stat,__LINE__,__FILE__);
c141_dims[0] = D1_dim;
c141_dims[1] = D4_dim;
c141_dims[2] = D1_dim;
stat = nc_def_var(ncid, "c141", NC_CHAR, RANK_c141, c141_dims, &c141_id);
check_err(stat,__LINE__,__FILE__);
b142_dims[0] = D1_dim;
b142_dims[1] = D4_dim;
b142_dims[2] = D2_dim;
stat = nc_def_var(ncid, "b142", NC_BYTE, RANK_b142, b142_dims, &b142_id);
check_err(stat,__LINE__,__FILE__);
s143_dims[0] = D1_dim;
s143_dims[1] = D4_dim;
s143_dims[2] = D3_dim;
stat = nc_def_var(ncid, "s143", NC_SHORT, RANK_s143, s143_dims, &s143_id);
check_err(stat,__LINE__,__FILE__);
i144_dims[0] = D1_dim;
i144_dims[1] = D4_dim;
i144_dims[2] = D4_dim;
stat = nc_def_var(ncid, "i144", NC_INT, RANK_i144, i144_dims, &i144_id);
check_err(stat,__LINE__,__FILE__);
f211_dims[0] = D2_dim;
f211_dims[1] = D1_dim;
f211_dims[2] = D1_dim;
stat = nc_def_var(ncid, "f211", NC_FLOAT, RANK_f211, f211_dims, &f211_id);
check_err(stat,__LINE__,__FILE__);
d212_dims[0] = D2_dim;
d212_dims[1] = D1_dim;
d212_dims[2] = D2_dim;
stat = nc_def_var(ncid, "d212", NC_DOUBLE, RANK_d212, d212_dims, &d212_id);
check_err(stat,__LINE__,__FILE__);
c213_dims[0] = D2_dim;
c213_dims[1] = D1_dim;
c213_dims[2] = D3_dim;
stat = nc_def_var(ncid, "c213", NC_CHAR, RANK_c213, c213_dims, &c213_id);
check_err(stat,__LINE__,__FILE__);
b214_dims[0] = D2_dim;
b214_dims[1] = D1_dim;
b214_dims[2] = D4_dim;
stat = nc_def_var(ncid, "b214", NC_BYTE, RANK_b214, b214_dims, &b214_id);
check_err(stat,__LINE__,__FILE__);
s221_dims[0] = D2_dim;
s221_dims[1] = D2_dim;
s221_dims[2] = D1_dim;
stat = nc_def_var(ncid, "s221", NC_SHORT, RANK_s221, s221_dims, &s221_id);
check_err(stat,__LINE__,__FILE__);
i222_dims[0] = D2_dim;
i222_dims[1] = D2_dim;
i222_dims[2] = D2_dim;
stat = nc_def_var(ncid, "i222", NC_INT, RANK_i222, i222_dims, &i222_id);
check_err(stat,__LINE__,__FILE__);
f223_dims[0] = D2_dim;
f223_dims[1] = D2_dim;
f223_dims[2] = D3_dim;
stat = nc_def_var(ncid, "f223", NC_FLOAT, RANK_f223, f223_dims, &f223_id);
check_err(stat,__LINE__,__FILE__);
d224_dims[0] = D2_dim;
d224_dims[1] = D2_dim;
d224_dims[2] = D4_dim;
stat = nc_def_var(ncid, "d224", NC_DOUBLE, RANK_d224, d224_dims, &d224_id);
check_err(stat,__LINE__,__FILE__);
c231_dims[0] = D2_dim;
c231_dims[1] = D3_dim;
c231_dims[2] = D1_dim;
stat = nc_def_var(ncid, "c231", NC_CHAR, RANK_c231, c231_dims, &c231_id);
check_err(stat,__LINE__,__FILE__);
b232_dims[0] = D2_dim;
b232_dims[1] = D3_dim;
b232_dims[2] = D2_dim;
stat = nc_def_var(ncid, "b232", NC_BYTE, RANK_b232, b232_dims, &b232_id);
check_err(stat,__LINE__,__FILE__);
s233_dims[0] = D2_dim;
s233_dims[1] = D3_dim;
s233_dims[2] = D3_dim;
stat = nc_def_var(ncid, "s233", NC_SHORT, RANK_s233, s233_dims, &s233_id);
check_err(stat,__LINE__,__FILE__);
i234_dims[0] = D2_dim;
i234_dims[1] = D3_dim;
i234_dims[2] = D4_dim;
stat = nc_def_var(ncid, "i234", NC_INT, RANK_i234, i234_dims, &i234_id);
check_err(stat,__LINE__,__FILE__);
f241_dims[0] = D2_dim;
f241_dims[1] = D4_dim;
f241_dims[2] = D1_dim;
stat = nc_def_var(ncid, "f241", NC_FLOAT, RANK_f241, f241_dims, &f241_id);
check_err(stat,__LINE__,__FILE__);
d242_dims[0] = D2_dim;
d242_dims[1] = D4_dim;
d242_dims[2] = D2_dim;
stat = nc_def_var(ncid, "d242", NC_DOUBLE, RANK_d242, d242_dims, &d242_id);
check_err(stat,__LINE__,__FILE__);
c243_dims[0] = D2_dim;
c243_dims[1] = D4_dim;
c243_dims[2] = D3_dim;
stat = nc_def_var(ncid, "c243", NC_CHAR, RANK_c243, c243_dims, &c243_id);
check_err(stat,__LINE__,__FILE__);
b244_dims[0] = D2_dim;
b244_dims[1] = D4_dim;
b244_dims[2] = D4_dim;