# Make a script to output data for CDIAC with 2000 data # Include: #wind speed #wind direction #pressure #soil temperature (5, 10, and 15 cm) (giving 5,10,30cm) #air temperature #precipitation #relative humidity #soil moisture (10 and 20 cm) (giving 10,30cm) #vapor pressure deficit # dayname_cbind( # all99[,"time"],round(all99[,"hour"],1), # all99[,"wsrmy"],all99[,"wdrmy"], # all99[,"pressure"], # all99[,"ts7"], all99[,"ts8"], all99[,"ts9"], # all99[,"ta1"], # all99[,"rain"], # all99[,"rh1"], # all99[,"sh2o1"],all99[,"sh2o2"], # all99[,"vpdef"],all99[,"par2"] # ) rnet <- all99[,13] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all99[j,14] } par <- all99[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all99[j,15]*1.0514 -1.0405 } dayname_cbind( all99[,1],round(all99[,2],1), all99[,3],all99[,8], all99[,13],all99[,15], all99[,17:51], all99[,57], all99[,61:68], all99[,72:74], all99[,84] ) # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC99.dat",sep=",",append=F, dimnames.write = F) # 2000 Data Here rnet <- all00[,13] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all00[j,14] } par <- all00[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all00[j,15]*1.0503+2.1769 } dayname_cbind( all00[,1],round(all00[,2],1), all00[,3],all00[,8], rnet, par, all00[,17:51], all00[,57], all00[,61:68], all00[,72:74], all00[,84] ) # To store the splus object file! Need to set up the names! newname_names(dayname) newname[1]_"time" newname[2]_"hour" newname[3]_"wsrmy" newname[4]_"wdrmy" newname[5]_"rnet" newname[6]_"par" newname[42]_"o3s" newname[54]_"vpdef names(dayname)_newname # cdiac00_dayname # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC00.dat",sep=",",append=F, dimnames.write = F) # 2001 Data Now rnet <- all01[,13] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all01[j,14] } par <- all01[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all01[j,15]*1.1002-2.9861 } dayname_cbind( all01[,1],round(all01[,2],1), all01[,3],all01[,8], rnet, par, all01[,17:51], all01[,57], all01[,61:68], all01[,72:74], all01[,84] ) # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC01.dat",sep=",",append=F, dimnames.write = F) # 2002 Data Here rnet <- all02[,13] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all02[j,14] } par <- all02[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all02[j,15]*1.0006+0.6580 } dayname_cbind( all02[,1],round(all02[,2],1), all02[,3],all02[,8], rnet, par, all02[,17:51], all02[,57], all02[,61:68], all02[,72:74], all02[,84] ) # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC02.dat",sep=",",append=F, dimnames.write = F) # 2003 data here rnet <- all03[,13] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all03[j,14] } par <- all03[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all03[j,15] # *1.0006+0.6580 CHECK THIS VALUE, or not use?? } dayname_cbind( all03[,1],round(all03[,2],1), all03[,3],all03[,8], rnet, par, all03[,17:51], all03[,57], all03[,61:68], all03[,72:74], all03[,84] ) # Clean out bad points for CDIAC archive dayname[10960,49]_as.numeric(NA) dayname[16283,52]_as.numeric(NA) dayname[11646:11653,52]_as.numeric(NA) dayname[10960,53]_as.numeric(NA) # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC03.dat",sep=",",append=F, dimnames.write = F) # 2004 data here rnet <- all04[,14] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all04[j,14] } par <- all04[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all04[j,15] # *1.0006+0.6580 } dayname_cbind( all04[,1],round(all04[,2],1), all04[,3],all04[,8], rnet, par, all04[,17:51], all04[,57], all04[,61:68], all04[,72:74], all04[,84] ) # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC04.dat",sep=",",append=F, dimnames.write = F) # 2005 data here rnet <- all05[,14] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all05[j,14] } par <- all05[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all05[j,15] # *1.0006+0.6580 } dayname_cbind( all05[,1],round(all05[,2],1), all05[,3],all05[,8], rnet, par, all05[,17:51], all05[,57], all05[,61:68], all05[,72:74], all05[,84] ) # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC05.dat",sep=",",append=F, dimnames.write = F) # 2006 data here rnet <- all06[,14] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all06[j,14] } par <- all06[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all06[j,15] # *1.0006+0.6580 } dayname_cbind( all06[,1],round(all06[,2],1), all06[,3],all06[,8], rnet, par, all06[,17:51], all06[,57], all06[,61:68], all06[,72:74], all06[,84] ) # Remove some 'bad' spots left in so you know that sonic, et al. not great dayname[4183:4190,17]_-99999 # sonic air temp dayname[2908,51]_-99999 #fmom dayname[15839,51]_-99999 #fmom dayname[15839,52]_-99999 #fsensht # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC06.dat",sep=",",append=F, dimnames.write = F) # 2007 data here: note stoppped in Jun, restarted Aug-Oct rnet <- all07[,14] i_length(rnet) for (j in 1:i) { if (is.na(rnet[j])) rnet[j]_all07[j,14] } par <- all07[,16] i_length(par) for (j in 1:i) { if (is.na(par[j])) par[j]_all07[j,15] # *1.0006+0.6580 } dayname_cbind( all07[,1],round(all07[,2],1), all07[,3],all07[,8], rnet, par, all07[,17:51], all07[,57], all07[,61:68], all07[,72:74], all07[,84] ) plot(dayname[,1],dayname[,3]) plot(dayname[,1],dayname[,4]) plot(dayname[,1],dayname[,5]) plot(dayname[,1],dayname[,6]) plot(dayname[,1],dayname[,7]) plot(dayname[,1],dayname[,8]) plot(dayname[,1],dayname[,9]) plot(dayname[,1],dayname[,10]) plot(dayname[,1],dayname[,11]) plot(dayname[,1],dayname[,12]) plot(dayname[,1],dayname[,13]) plot(dayname[,1],dayname[,14]) plot(dayname[,1],dayname[,15]) plot(dayname[,1],dayname[,16]) plot(dayname[,1],dayname[,17]) plot(dayname[,1],dayname[,18]) plot(dayname[,1],dayname[,19]) plot(dayname[,1],dayname[,20]) plot(dayname[,1],dayname[,21]) plot(dayname[,1],dayname[,22]) plot(dayname[,1],dayname[,23]) plot(dayname[,1],dayname[,24]) plot(dayname[,1],dayname[,25]) plot(dayname[,1],dayname[,26]) plot(dayname[,1],dayname[,27]) plot(dayname[,1],dayname[,28]) plot(dayname[,1],dayname[,29]) # plot(dayname[,1],dayname[,30]) # plot(dayname[,1],dayname[,31]) plot(dayname[,1],dayname[,32]) plot(dayname[,1],dayname[,33]) plot(dayname[,1],dayname[,34]) plot(dayname[,1],dayname[,35]) plot(dayname[,1],dayname[,36]) plot(dayname[,1],dayname[,37]) plot(dayname[,1],dayname[,38]) plot(dayname[,1],dayname[,39]) plot(dayname[,1],dayname[,40]) plot(dayname[,1],dayname[,41]) plot(dayname[,1],dayname[,42]) plot(dayname[,1],dayname[,43]) plot(dayname[,1],dayname[,44]) plot(dayname[,1],dayname[,45]) plot(dayname[,1],dayname[,46]) plot(dayname[,1],dayname[,47]) plot(dayname[,1],dayname[,48]) plot(dayname[,1],dayname[,49]) plot(dayname[,1],dayname[,50]) plot(dayname[,1],dayname[,51]) plot(dayname[,1],dayname[,52]) plot(dayname[,1],dayname[,53]) plot(dayname[,1],dayname[,54]) # Remove some 'bad' spots left in so you know that sonic, et al. not great dayname[1802:1803,52]_-99999 #fsensht dayname[1813,52]_-99999 #fsensht dayname[1973:1974,52]_-99999 #fsensht dayname[2723,52]_-99999 #fsensht dayname[5921,52]_-99999 #fsensht # Get rid of NAs values in the file dayname[ is.na(dayname)]_ -99999 dayname[dayname == NA]_ -99999 write.table(dayname, file="CDIAC07.dat",sep=",",append=F, dimnames.write = F)