-
Notifications
You must be signed in to change notification settings - Fork 0
/
iau2wkt.c
358 lines (285 loc) · 33.7 KB
/
iau2wkt.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
#include<stdio.h>
#include<stdlib.h>
#include<string.h>
/*------------------------------------------------------------------------------
* Name: create_IAU2000_WKT_v2.py
* Author: Trent Hare
* Original: Jan 2006
* Last modified Feb 2016 to:
* ---- update to report IAU Mean from reports,
* ---- Asteroids and IAU reported Comets,
* ---- and two new projections (Mollweide and Robinson)
*
* Description: This Python script creates a IAU2000/2009 WKT projection strings for WMS services
*
* License: Public Domain
*
* INPUT: (naifcodes_radii_m_wAsteroids_IAU2000.csv or naifcodes_radii_m_wAsteroids_IAU2009.csv)
*
* Example file format:
* Naif_id,Body,IAU2000_Mean,IAU2000_Semimajor,IAU2000_Axisb,IAU2000_Semiminor
* 199,Mercury,2439700.00,2439700.00,2439700.00,2439700.00
* 299,Venus,6051800.00,6051800.00,6051800.00,6051800.00
* 399,Earth,6371000.00,6378140.00,6378140.00,6356750.00
* 301,Moon,1737400.00,1737400.00,1737400.00,1737400.00
* 499,Mars,3389500.00,3396190.00,3396190.00,3376200.00
* 401,Phobos,11100.00,13400.00,11200.00,9200.00
* 402,Deimos,6200.00,7500.00,6100.00,5200.00
* ...
* OUTPUT:
* Example: WMS#,GEOGCS["Mars 2000",DATUM["D_Mars_2000",SPHEROID["Mars_2000_IAU_IAG",3396190,169.8944472236118]],PRIMEM["Reference_Meridian",0],UNIT["Decimal_Degree",0.0174532925199433]]
*------------------------------------------------------------------------------
*/
void usage(){
printf("usage:\n\tiau2wkt naifcodes_radii_m_wAsteroids_IAU2000.csv outputFileName.wtk\n");
}
const char* getfield(char* line, int num)
{
const char* tok;
for (tok = strtok(line, ",");
tok && *tok;
tok = strtok(NULL, ",\n"))
{
if (!--num)
return tok;
}
return NULL;
}
int main(int argc, char * argv[])
{
if(argc < 2) usage();
FILE * stream;
stream = fopen(argv[1],"r");
/*grab year from file name*/
int theYear = 0;
char * pch;
pch = strtok(argv[1], "U.");
int count=0;
while (pch != NULL)
{
if(count==1) theYear=atoi(pch);
pch = strtok(NULL, "U.");
count++;
}
/*printf("theYear = %d\n",theYear);*/
/*Check if the Year has been found from filename*/
FILE *f1;
f1 = fopen(argv[2], "w") ;
if ( theYear ){
/*Open output file for writing*/
fprintf(f1, "# IAU%i WKT Codes\n", theYear);
} else {
printf("Can't parse the year from filename: %s \n",argv[2]);
usage();
exit(EXIT_FAILURE);
}
/*Read input file line by line*/
char line[1024];
int theNaifNum = 0;
char * theTarget = NULL;
char * theMeanChar = NULL;
float theMean = 0.0, flattening = 0.0;
float theA = 0.0, theB = 0.0, theC = 0.0;
/*Start parsing line by line*/
while (fgets(line, 1024, stream))
{
/*Need to replace original line in tmp each time used*/
char* tmp = strdup(line);
/*printf("\nline is %s",line);*/
theNaifNum = atoi(getfield(tmp,1));
/*printf("NAIF: %d\n",theNaifNum);*/
tmp = strdup(line);
if(getfield(tmp,3) != NULL){
tmp = strdup(line);
theMeanChar = (char *) getfield(tmp,3);
} else {
theMeanChar = NULL;
}
/*printf("MeanChar: %s\n",theMeanChar);*/
tmp = strdup(line);
if (theNaifNum != 0 && theMeanChar)/*#is it a number?*/
{
theTarget = (char *) getfield(tmp,2);
/*printf("Target: %s\t",theTarget);*/
tmp = strdup(line);
/*Check that third token in line is not empty*/
if (atof(getfield(tmp,3)) != 0)/*Then Radii values exist in input table*/
{
tmp = strdup(line);
theMean = atof(getfield(tmp,3));
tmp = strdup(line);
theA = atof(getfield(tmp,4));
tmp = strdup(line);
theB = atof(getfield(tmp,5));
tmp = strdup(line);
theC = atof(getfield(tmp,6));
/*printf("Mean: %f\t",theMean);*/
/*printf("A: %f\t",theA);*/
/*printf("B: %f\t",theB);*/
/*printf("C: %f\n",theC);*/
/*#Check to see if the Mean should be used, for triaxial bodies*/
if ((theA != theB) && (theA != theC))
{
theA = theMean;
theC = theMean;
}
flattening = ((theA - theC) / theA);
if (flattening < 0.0000000001
&& flattening > -0.0000000001
&& flattening != 0.0)
{
/*Inverse flattening if too small*/
flattening = 1.0 / flattening;
}
/*Even = Areocentric*/
/*Odd = Areographic*/
/*GEOIDS*/
fprintf(f1, "# IAU%d WKT Codes for %s\n", theYear, theTarget);
int gisCode = theNaifNum * 100;
fprintf(f1,"# GEOGCS %s Areocentric\n", theTarget);
fprintf(f1, "%d,GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 1;
fprintf(f1,"# GEOGCS %s Areographic\n", theTarget);
fprintf(f1, "%d,GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
/*#Static Projections*/
gisCode = theNaifNum * 100 + 10; /*# Equirectangular, ocentric, clon=0*/
fprintf(f1,"# Equirectangular %s Areocentric, clon=0\n", theTarget);
fprintf(f1,"%d,PROJCS[\"%s_Equidistant_Cylindrical\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Equirectangular\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 11; /*# Equirectangular, ographic, clon=0*/
fprintf(f1,"# Equirectangular %s Areographic, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Equidistant_Cylindrical\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Equirectangular\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 12; /*# Equirectangular, ocentric, clon=180*/
fprintf(f1,"# Equirectangular %s Areocentric, clon=180\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Equidistant_Cylindrical\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Equirectangular\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",180],PARAMETER[\"Standard_Parallel_1\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 13; /*# Equirectangular, ographic, clon=180*/
fprintf(f1,"# Equirectangular %s Areographic, clon=180\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Equidistant_Cylindrical\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Equirectangular\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",180],PARAMETER[\"Standard_Parallel_1\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 14; /*# Sinusoidal, ocentric, clon=0*/
fprintf(f1,"# Sinusoidal %s Areocentric, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Sinusoidal\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Sinusoidal\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 15; /*# Sinusoidal, ographic, clon=0*/
fprintf(f1,"# Sinusoidal %s Areographic, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Sinusoidal\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Sinusoidal\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 16; /*# Sinusoidal, ocentric, clon=180*/
fprintf(f1,"# Sinusoidal %s Areocentric, clon=180\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Sinusoidal\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Sinusoidal\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",180],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 17; /*# Sinusoidal, ographic, clon=180*/
fprintf(f1,"# Sinusoidal %s Areographic, clon=180\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Sinusoidal\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Sinusoidal\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",180],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 18; /*# North Polar, ocentric, clon=0*/
fprintf(f1,"# North Polar %s Areocentric, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_North_Pole_Stereographic\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Stereographic\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Scale_Factor\",1],PARAMETER[\"Latitude_Of_Origin\",90],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 19; /*# North Polar, ographic, clon=0*/
fprintf(f1,"# North Polar %s Areographic, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_North_Pole_Stereographic\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Stereographic\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Scale_Factor\",1],PARAMETER[\"Latitude_Of_Origin\",90],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 20; /*# South Polar, ocentric, clon=0*/
fprintf(f1,"# South Polar %s Areocentric, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_South_Pole_Stereographic\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Stereographic\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Scale_Factor\",1],PARAMETER[\"Latitude_Of_Origin\",-90],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 21; /*# South Polar, ographic, clon=0*/
fprintf(f1,"# South Polar %s Areographic, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_South_Pole_Stereographic\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Stereographic\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Scale_Factor\",1],PARAMETER[\"Latitude_Of_Origin\",-90],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 22; /*# Mollweide, ocentric, clon=0*/
fprintf(f1,"# Mollweide %s Areocentric, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Mollweide\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Mollweide\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 23; /*# Mollweide, ographic, clon=0*/
fprintf(f1,"# Mollweide %s Areographic, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Mollweide\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Mollweide\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 24; /*# Mollweide, ocentric, clon=180*/
fprintf(f1,"# Mollweide %s Areocentric, clon=180\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Mollweide\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Mollweide\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",180],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 25; /*# Mollweide, ocentric, clon=180*/
fprintf(f1,"# Mollweide %s Areographic, clon=180\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Mollweide\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Mollweide\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",180],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 26; /*# Robinson, ocentric, clon=0*/
fprintf(f1,"# Robinson %s Areocentric, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Robinson\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Robinson\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 27; /*# Robinson, ographic, clon=0*/
fprintf(f1,"# Robinson %s Areographic, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Robinson\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Robinson\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 28; /*# Robinson, ocentric, clon=180*/
fprintf(f1,"# Robinson %s Areocentric, clon=180\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Robinson\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Robinson\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",180],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 29; /*# Robinson, ographic, clon=180*/
fprintf(f1,"# Robinson %s Areographic, clon=180\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Robinson\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Robinson\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",180],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
/*#AUTO Projections*/
gisCode = theNaifNum * 100 + 60; /*# Sinusoidal, ocentric*/
fprintf(f1,"# Sinusoidal AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Sinusoidal_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Sinusoidal\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 61; /*# Sinusoidal, ographic*/
fprintf(f1,"# Sinusoidal AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Sinusoidal_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Sinusoidal\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 62; /*# Stereographic, ocentric, clon=0*/
fprintf(f1,"# Stereographic AUTO %s Areocentric, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Stereographic_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Stereographic\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Scale_Factor\",1],PARAMETER[\"Latitude_Of_Origin\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 63; /*# Stereographic, ographic, clon=0*/
fprintf(f1,"# Stereographic AUTO %s Areographic, clon=0\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Stereographic_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Stereographic\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Scale_Factor\",1],PARAMETER[\"Latitude_Of_Origin\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 64; /*# Transverse Mercator, ocentric*/
fprintf(f1,"# Transverse Mercator AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Transverse_Mercator_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Scale_Factor\",0.9996],PARAMETER[\"Latitude_Of_Origin\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 65; /*# Transverse Mercator, ographic*/
fprintf(f1,"# Transverse Mercator AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Transverse_Mercator_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Transverse_Mercator\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Scale_Factor\",0.9996],PARAMETER[\"Latitude_Of_Origin\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 66; /*# Orthographic, ocentric*/
fprintf(f1,"# Orthographic AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Orthographic_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Orthographic\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Longitude_Of_Center\",0.0],PARAMETER[\"Latitude_Of_Center\",90.0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 67; /*# Orthographic, ographic*/
fprintf(f1,"# Orthographic AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Orthographic_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Orthographic\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Longitude_Of_Center\",0.0],PARAMETER[\"Latitude_Of_Center\",90.0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 68; /*# Equidistant_Cylindrical, ocentric*/
fprintf(f1,"# Equidistant Cylindrical AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Equidistant_Cylindrical_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Equirectangular\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 69; /*# Equidistant_Cylindrical, ographic*/
fprintf(f1,"# Equidistant Cylindrical AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Equidistant_Cylindrical_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Equirectangular\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 70; /*# Lambert_Conformal_Conic, ocentric*/
fprintf(f1,"# Lambert Conformal Conic AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Lambert_Conformal_Conic_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Conformal_Conic_2SP\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",-20],PARAMETER[\"Standard_Parallel_2\",20],PARAMETER[\"Latitude_Of_Origin\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 71; /*# Lambert_Conformal_Conic, ographic*/
fprintf(f1,"# Lambert Conformal Conic AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Lambert_Conformal_Conic_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Conformal_Conic_2SP\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",-20],PARAMETER[\"Standard_Parallel_2\",20],PARAMETER[\"Latitude_Of_Origin\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 72; /*# Lambert_Azimuthal_Equal_Area, ocentric*/
fprintf(f1,"# Lambert Azimuthal Equal Area AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Lambert_Azimuthal_Equal_Area_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Latitude_Of_Origin\",90],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 73; /*# Lambert_Azimuthal_Equal_Area, ographic*/
fprintf(f1,"# Lambert Azimuthal Equal Area AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Lambert_Azimuthal_Equal_Area_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Latitude_Of_Origin\",90],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 74; /*# Mercator, ocentric*/
fprintf(f1,"# Mercator AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Mercator_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Mercator_1SP\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 75; /*# Mercator, ographic*/
fprintf(f1,"# Mercator AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Mercator_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Mercator_1SP\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],PARAMETER[\"Standard_Parallel_1\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 76; /*# Albers, ocentric*/
fprintf(f1,"# Albers AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Albers_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Albers_Conic_Equal_Area\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0.0],PARAMETER[\"Standard_Parallel_1\",60.0],PARAMETER[\"Standard_Parallel_2\",20.0],PARAMETER[\"Latitude_Of_Origin\",40.0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 77; /*# Albers, ographic*/
fprintf(f1,"# Albers AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Albers_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Albers_Conic_Equal_Area\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0.0],PARAMETER[\"Standard_Parallel_1\",60.0],PARAMETER[\"Standard_Parallel_2\",20.0],PARAMETER[\"Latitude_Of_Origin\",40.0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 78; /*# Oblique Cylindrical Equal Area, ocentric*/
fprintf(f1,"# Oblique Cylindrical Equal Area AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Oblique_Cylindrical_Equal_Area_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Oblique_Cylindrical_Equal_Area\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0.0],PARAMETER[\"Standard_Parallel_1\",0.0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 79; /*# Oblique Cylindrical Equal Area, ographic*/
fprintf(f1,"# Oblique Cylindrical Equal Area AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Oblique_Cylindrical_Equal_Area_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Oblique_Cylindrical_Equal_Area\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0.0],PARAMETER[\"Standard_Parallel_1\",0.0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 80; /*# Mollweide, ocentric*/
fprintf(f1,"# Mollweide AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Mollweide_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Mollweide\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 81; /*# Mollweide, ographic*/
fprintf(f1,"# Mollweide AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Mollweide_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Mollweide\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 82; /*# Robinson, ocentric*/
fprintf(f1,"# Robinson AUTO %s Areocentric\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Robinson_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Robinson\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
gisCode = theNaifNum * 100 + 83; /*# Robinson, ographic*/
fprintf(f1,"# Robinson AUTO %s Areographic\n", theTarget);
fprintf(f1, "%d,PROJCS[\"%s_Robinson_AUTO\",GEOGCS[\"%s %d\",DATUM[\"D_%s_%d\",SPHEROID[\"%s_%d_IAU_IAG\",%f,%f]],PRIMEM[\"Reference_Meridian\",0],UNIT[\"Decimal_Degree\",0.0174532925199433]],PROJECTION[\"Robinson\"],PARAMETER[\"False_Easting\",0],PARAMETER[\"False_Northing\",0],PARAMETER[\"Central_Meridian\",0],UNIT[\"Meter\",1],AUTHORITY[\"IAU%d\",\"%d\"]]\n",gisCode,theTarget,theTarget,theYear,theTarget,theYear,theTarget,theYear,theA,flattening,theYear,gisCode);
}
}
free(tmp);
}
fclose(stream);
fclose(f1);
return(EXIT_SUCCESS);
}