;;creates a century site.100 file named outfname
;; input it is a large array of strings for weather data
;; c3 is an optional keyword determining the percent of c3 vs c4 grass

pro site100, input, outdir, outfname, c3=c3

	dat=input
	input=strarr(62)

	for i=0, 61 do begin
		input(i) = strcompress(dat(i), /remove_all)
		if dat(i) ge 1 then input(i) = input(i)+' '
	endfor
	cd, outdir, current=current
	
	openw, oun, /get, outfname

	if keyword_set(c3) then begin
		if c3 gt 100 then begin
			print, 'ERROR, c3=',c3,'... defaulting to CPR'
			printf, oun, 'CPR   Jornada'
		endif else if c3 lt 0 then begin
			print, 'ERROR, c3=',c3,'... defaulting to G5'
			printf, oun, 'G5    grass__mixed_75%warm'
		endif else if c3 lt 37 then begin
			print, 'G5 = ', outfname
			printf, oun, 'G5    grass__mixed_75%warm'
		endif else if c3 lt 63 then begin
			print, 'G3 = ', outfname
			printf, oun, 'G3    grass__mixed_50%_warm_50%_cool'
		endif else if c3 lt 85 then begin
			print, 'G4 = ', outfname			
			printf, oun, 'G4    grass__mixed_25%_warm'
		endif else if c3 le 100 then begin
			print, 'CPR = ', outfname
			printf, oun, 'CPR   Jornada'
		endif else begin
			print, 'UNCLASSIFIED C3 =',c3
			retall
		endelse
	
	endif else $
		printf, oun, 'CPR   Jornada'
	
	printf, oun, '*** Climate parameters'
	printf, oun, input(0),'           ''PRECIP(1)'' '
	printf, oun, input(1),'           ''PRECIP(2)'' '
	printf, oun, input(2),'           ''PRECIP(3)'' '
	printf, oun, input(3),'           ''PRECIP(4)'' '
	printf, oun, input(4),'           ''PRECIP(5)'' '
	printf, oun, input(5),'           ''PRECIP(6)'' '
	printf, oun, input(6),'           ''PRECIP(7)'' '
	printf, oun, input(7),'           ''PRECIP(8)'' '
	printf, oun, input(8),'           ''PRECIP(9)'' '
	printf, oun, input(9),'           ''PRECIP(10)'' '
	printf, oun, input(10),'           ''PRECIP(11)'' '
	printf, oun, input(11),'           ''PRECIP(12)'' '
	printf, oun, input(12),'           ''PRCSTD(1)'' '
	printf, oun, input(13),'           ''PRCSTD(2)'' '
	printf, oun, input(14),'           ''PRCSTD(3)'' '
	printf, oun, input(15),'           ''PRCSTD(4)'' '
	printf, oun, input(16),'           ''PRCSTD(5)'' '
	printf, oun, input(17),'           ''PRCSTD(6)'' '
	printf, oun, input(18),'           ''PRCSTD(7)'' '
	printf, oun, input(19),'           ''PRCSTD(8)'' '
	printf, oun, input(20),'           ''PRCSTD(9)'' '
	printf, oun, input(21),'           ''PRCSTD(10)'' '
	printf, oun, input(22),'           ''PRCSTD(11)'' '
printf, oun, input(23),'           ''PRCSTD(12)'' '
printf, oun, input(24),'           ''PRCSKW(1)'' '
printf, oun, input(25),'           ''PRCSKW(2)'' '
printf, oun, input(26),'           ''PRCSKW(3)'' '
printf, oun, input(27),'           ''PRCSKW(4)'' '
printf, oun, input(28),'           ''PRCSKW(5)'' '
printf, oun, input(29),'           ''PRCSKW(6)'' '
printf, oun, input(30),'           ''PRCSKW(7)'' '
printf, oun, input(31),'           ''PRCSKW(8)'' '
printf, oun, input(32),'           ''PRCSKW(9)'' '
printf, oun, input(33),'           ''PRCSKW(10)'' '
printf, oun, input(34),'           ''PRCSKW(11)'' '
printf, oun, input(35),'           ''PRCSKW(12)'' '
printf, oun, input(36),'           ''TMN2M(1)''  '
printf, oun, input(37),'           ''TMN2M(2)''  '
printf, oun, input(38),'           ''TMN2M(3)''  '
printf, oun, input(39),'           ''TMN2M(4)''  '
printf, oun, input(40),'           ''TMN2M(5)''  '
printf, oun, input(41),'           ''TMN2M(6)''  '
printf, oun, input(42),'           ''TMN2M(7)''  '
printf, oun, input(43),'           ''TMN2M(8)''  '
printf, oun, input(44),'           ''TMN2M(9)''  '
printf, oun, input(45),'           ''TMN2M(10)'' '
printf, oun, input(46),'           ''TMN2M(11)'' '
printf, oun, input(47),'           ''TMN2M(12)'' '
printf, oun, input(48),'           ''TMX2M(1)''  '
printf, oun, input(49),'           ''TMX2M(2)''  '
printf, oun, input(50),'           ''TMX2M(3)''  '
printf, oun, input(51),'           ''TMX2M(4)''  '
printf, oun, input(52),'           ''TMX2M(5)''  '
printf, oun, input(53),'           ''TMX2M(6)''  '
printf, oun, input(54),'           ''TMX2M(7)''  '
printf, oun, input(55),'           ''TMX2M(8)''  '
printf, oun, input(56),'           ''TMX2M(9)''  '
printf, oun, input(57),'           ''TMX2M(10)'' '
printf, oun, input(58),'           ''TMX2M(11)'' '
printf, oun, input(59),'           ''TMX2M(12)'' '
printf, oun, '*** Site and control parameters'
printf, oun, '0.00000           ''IVAUTO''    '
printf, oun, '1.00000           ''NELEM''     '
printf, oun, input(60),'          ''SITLAT''    '
printf, oun, input(61),'          ''SITLNG''    '
printf, oun, '0.80000           ''SAND''      '
printf, oun, '0.10000           ''SILT''      '
printf, oun, '0.10000           ''CLAY''      '
printf, oun, '1.20000           ''BULKD''     '
printf, oun, '6.00000           ''NLAYER''    '
printf, oun, '5.00000           ''NLAYPG''    '
printf, oun, '1.00000           ''DRAIN''     '
printf, oun, '0.00000           ''BASEF''     '
printf, oun, '0.00000           ''STORMF''    '
printf, oun, '1.00000           ''SWFLAG''    '
printf, oun, '0.20000           ''AWILT(1)''  '
printf, oun, '0.20000           ''AWILT(2)''  '
printf, oun, '0.20000           ''AWILT(3)''  '
printf, oun, '0.20000           ''AWILT(4)''  '
printf, oun, '0.20000           ''AWILT(5)''  '
printf, oun, '0.20000           ''AWILT(6)''  '
printf, oun, '0.20000           ''AWILT(7)''  '
printf, oun, '0.20000           ''AWILT(8)''  '
printf, oun, '0.20000           ''AWILT(9)''  '
printf, oun, '0.30000           ''AWILT(10)'' '
printf, oun, '0.30000           ''AFIEL(1)''  '
printf, oun, '0.30000           ''AFIEL(2)''  '
printf, oun, '0.30000           ''AFIEL(3)''  '
printf, oun, '0.30000           ''AFIEL(4)''  '
printf, oun, '0.30000           ''AFIEL(5)''  '
printf, oun, '0.30000           ''AFIEL(6)''  '
printf, oun, '0.30000           ''AFIEL(7)''  '
printf, oun, '0.30000           ''AFIEL(8)''  '
printf, oun, '0.30000           ''AFIEL(9)''  '
printf, oun, '0.00000           ''AFIEL(10)'' '
printf, oun, '7.50000           ''PH''        '
printf, oun, '1.00000           ''PSLSRB''    '
printf, oun, '2.00000           ''SORPMX''    '
printf, oun, '*** External nutrient input parameters'
printf, oun, '0.05000           ''EPNFA(1)''  '
printf, oun, '0.00700           ''EPNFA(2)''  '
printf, oun, '-0.90000          ''EPNFS(1)''  '
printf, oun, '0.01500           ''EPNFS(2)''  '
printf, oun, '0.00000           ''SATMOS(1)'' '
printf, oun, '0.00000           ''SATMOS(2)'' '
printf, oun, '0.00000           ''SIRRI''     '
printf, oun, '*** Organic matter initial values'
printf, oun, '20.00000          ''SOM1CI(1,1)'' '
printf, oun, '0.00000           ''SOM1CI(1,2)'' '
printf, oun, '60.00000          ''SOM1CI(2,1)'' '
printf, oun, '0.00000           ''SOM1CI(2,2)'' '
printf, oun, '1250.00000        ''SOM2CI(1)'' '
printf, oun, '0.00000           ''SOM2CI(2)'' '
printf, oun, '740.00000         ''SOM3CI(1)'' '
printf, oun, '0.00000           ''SOM3CI(2)'' '
printf, oun, '16.00000          ''RCES1(1,1)'' '
printf, oun, '50.00000          ''RCES1(1,2)'' '
printf, oun, '50.00000          ''RCES1(1,3)'' '
printf, oun, '12.00000          ''RCES1(2,1)'' '
printf, oun, '50.00000          ''RCES1(2,2)'' '
printf, oun, '50.00000          ''RCES1(2,3)'' '
printf, oun, '17.00000          ''RCES2(1)''  '
printf, oun, '117.00000         ''RCES2(2)''  '
printf, oun, '117.00000         ''RCES2(3)''  '
printf, oun, '8.00000           ''RCES3(1)''  '
printf, oun, '62.00000          ''RCES3(2)''  '
printf, oun, '62.00000          ''RCES3(3)''  '
printf, oun, '100.00000         ''CLITTR(1,1)'' '
printf, oun, '0.00000           ''CLITTR(1,2)'' '
printf, oun, '100.00000         ''CLITTR(2,1)'' '
printf, oun, '0.00000           ''CLITTR(2,2)'' '
printf, oun, '66.00000          ''RCELIT(1,1)'' '
printf, oun, '300.00000         ''RCELIT(1,2)'' '
printf, oun, '300.00000         ''RCELIT(1,3)'' '
printf, oun, '66.00000          ''RCELIT(2,1)'' '
printf, oun, '300.00000         ''RCELIT(2,2)'' '
printf, oun, '300.00000         ''RCELIT(2,3)'' '
printf, oun, '0.00000           ''AGLCIS(1)'' '
printf, oun, '0.00000           ''AGLCIS(2)'' '
printf, oun, '0.00000           ''AGLIVE(1)'' '
printf, oun, '0.00000           ''AGLIVE(2)'' '
printf, oun, '0.00000           ''AGLIVE(3)'' '
printf, oun, '200.00000         ''BGLCIS(1)'' '
printf, oun, '0.00000           ''BGLCIS(2)'' '
printf, oun, '3.00000           ''BGLIVE(1)'' '
printf, oun, '0.40000           ''BGLIVE(2)'' '
printf, oun, '0.40000           ''BGLIVE(3)'' '
printf, oun, '50.00000          ''STDCIS(1)'' '
printf, oun, '0.00000           ''STDCIS(2)'' '
printf, oun, '0.80000           ''STDEDE(1)'' '
printf, oun, '0.20000           ''STDEDE(2)'' '
printf, oun, '0.20000           ''STDEDE(3)'' '
printf, oun, '*** Forest organic matter initial parameters'
printf, oun, '0.00000           ''RLVCIS(1)'' '
printf, oun, '0.00000           ''RLVCIS(2)'' '
printf, oun, '0.00000           ''RLEAVE(1)'' '
printf, oun, '0.00000           ''RLEAVE(2)'' '
printf, oun, '0.00000           ''RLEAVE(3)'' '
printf, oun, '0.00000           ''FBRCIS(1)'' '
printf, oun, '0.00000           ''FBRCIS(2)'' '
printf, oun, '0.00000           ''FBRCHE(1)'' '
printf, oun, '0.00000           ''FBRCHE(2)'' '
printf, oun, '0.00000           ''FBRCHE(3)'' '
printf, oun, '0.00000           ''RLWCIS(1)'' '
printf, oun, '0.00000           ''RLWCIS(2)'' '
printf, oun, '0.00000           ''RLWODE(1)'' '
printf, oun, '0.00000           ''RLWODE(2)'' '
printf, oun, '0.00000           ''RLWODE(3)'' '
printf, oun, '0.00000           ''FRTCIS(1)'' '
printf, oun, '0.00000           ''FRTCIS(2)'' '
printf, oun, '0.00000           ''FROOTE(1)'' '
printf, oun, '0.00000           ''FROOTE(2)'' '
printf, oun, '0.00000           ''FROOTE(3)'' '
printf, oun, '0.00000           ''CRTCIS(1)'' '
printf, oun, '0.00000           ''CRTCIS(2)'' '
printf, oun, '0.00000           ''CROOTE(1)'' '
printf, oun, '0.00000           ''CROOTE(2)'' '
printf, oun, '0.00000           ''CROOTE(3)'' '
printf, oun, '0.00000           ''WD1CIS(1)'' '
printf, oun, '0.00000           ''WD1CIS(2)'' '
printf, oun, '0.00000           ''WD2CIS(1)'' '
printf, oun, '0.00000           ''WD2CIS(2)'' '
printf, oun, '0.00000           ''WD3CIS(1)'' '
printf, oun, '0.00000           ''WD3CIS(2)'' '
printf, oun, '0.30000           ''W1LIG''     '
printf, oun, '0.30000           ''W2LIG''     '
printf, oun, '0.30000           ''W3LIG''     '
printf, oun, '*** Mineral initial parameters'
printf, oun, '0.25000           ''MINERL(1,1)'' '
printf, oun, '0.00000           ''MINERL(2,1)'' '
printf, oun, '0.00000           ''MINERL(3,1)'' '
printf, oun, '0.00000           ''MINERL(4,1)'' '
printf, oun, '0.00000           ''MINERL(5,1)'' '
printf, oun, '0.00000           ''MINERL(6,1)'' '
printf, oun, '0.00000           ''MINERL(7,1)'' '
printf, oun, '0.00000           ''MINERL(8,1)'' '
printf, oun, '0.00000           ''MINERL(9,1)'' '
printf, oun, '0.00000           ''MINERL(10,1)'' '
printf, oun, '0.50000           ''MINERL(1,2)'' '
printf, oun, '0.00000           ''MINERL(2,2)'' '
printf, oun, '0.00000           ''MINERL(3,2)'' '
printf, oun, '0.00000           ''MINERL(4,2)'' '
printf, oun, '0.00000           ''MINERL(5,2)'' '
printf, oun, '0.00000           ''MINERL(6,2)'' '
printf, oun, '0.00000           ''MINERL(7,2)'' '
printf, oun, '0.00000           ''MINERL(8,2)'' '
printf, oun, '0.00000           ''MINERL(9,2)'' '
printf, oun, '0.00000           ''MINERL(10,2)'' '
printf, oun, '0.50000           ''MINERL(1,3)'' '
printf, oun, '0.00000           ''MINERL(2,3)'' '
printf, oun, '0.00000           ''MINERL(3,3)'' '
printf, oun, '0.00000           ''MINERL(4,3)'' '
printf, oun, '0.00000           ''MINERL(5,3)'' '
printf, oun, '0.00000           ''MINERL(6,3)'' '
printf, oun, '0.00000           ''MINERL(7,3)'' '
printf, oun, '0.00000           ''MINERL(8,3)'' '
printf, oun, '0.00000           ''MINERL(9,3)'' '
printf, oun, '0.00000           ''MINERL(10,3)'' '
printf, oun, '0.00000           ''PARENT(1)'' '
printf, oun, '50.00000          ''PARENT(2)'' '
printf, oun, '50.00000          ''PARENT(3)'' '
printf, oun, '0.00000           ''SECNDY(1)'' '
printf, oun, '15.00000          ''SECNDY(2)'' '
printf, oun, '2.00000           ''SECNDY(3)'' '
printf, oun, '0.00000           ''OCCLUD''    '
printf, oun, '*** Water initial parameters'
printf, oun, '0.00000           ''RWCF(1)''   '
printf, oun, '0.00000           ''RWCF(2)''   '
printf, oun, '0.00000           ''RWCF(3)''   '
printf, oun, '0.00000           ''RWCF(4)''   '
printf, oun, '0.00000           ''RWCF(5)''   '
printf, oun, '0.00000           ''RWCF(6)''   '
printf, oun, '0.00000           ''RWCF(7)''   '
printf, oun, '0.00000           ''RWCF(8)''   '
printf, oun, '0.00000           ''RWCF(9)''   '
printf, oun, '0.00000           ''RWCF(10)''  '
printf, oun, '0.00000           ''SNLQ''      '
printf, oun, '0.00000           ''SNOW''      '

close, oun
free_lun, oun
cd, current
end