Changes between Version 5 and Version 6 of Model Results Basque river basins

Timestamp:

Jun 13, 2012 4:15:42 PM (13 years ago)

Author:

aizkorri

Comment:

--

Model Results Basque river basins

@@ -58 +58 @@
 save(basquebasins,file=paste(datawd,"/dataEDAccm/basquebasins.RData",sep=""))
 
-#En 'EDA/EDACCM/prediction_1_dbeel_ok.R' run all 
+#En 'EDA/EDACCM/prediction_1_dbeel_ok.R' run all until
+# line 157 source("prediction_riverwidth_writedatabase_Basque.R")
 
 #En 'EDA/EDACCM/prediction_riverwidh_basque.R' run all 
 
-
-load(file=str_c(datawd,"/dataEDAccm/basquebasins",".Rdata")) 
-basquebasins@data
-sum(ccm$nb_eel_per_segment[ccm$wso_id==],na.rm=T)*0.05
+{{
+load(file=str_c(datawd,"/dataEDAccm/basquebasins",".Rdata"))
+bb<-basquebasins@data
+SEBasin<-array(0,c(length(levels(as.factor(bb$unidad_hid[bb$emp==1]))),5))
+colnames(SEBasin)<-c("Basin","Nb_YE","Nb_SE","Water_surface","Nb_YE/100m2")
+for(i in 1:length(levels(as.factor(bb$unidad_hid[bb$emp==1])))){
+        SEBasin[i,1]<-levels(as.factor(bb$unidad_hid[bb$emp==1]))[i]
+        SEBasin[i,2]<-round(sum(ccm$nb_eel_per_segment[ccm$wso1_id%in%bb$wso1_id[bb$unidad_hid==SEBasin[i,1]]],na.rm=T),0)
+        SEBasin[i,3]<-round(sum(ccm$nb_eel_per_segment[ccm$wso1_id%in%bb$wso1_id[bb$unidad_hid==SEBasin[i,1]]],na.rm=T)*0.05,0)
+        SEBasin[i,4]<-round(sum(ccm$riverareakm[ccm$wso1_id%in%bb$wso1_id[bb$unidad_hid==SEBasin[i,1]]]),3)
+        SEBasin[i,5]<-round(100*sum(ccm$nb_eel_per_segment[ccm$wso1_id%in%bb$wso1_id[bb$unidad_hid==SEBasin[i,1]]])/sum(ccm$riverarea[ccm$wso1_id%in%bb$wso1_id[bb$unidad_hid==SEBasin[i,1]]]),3)
+}
+SEBasin<-as.data.frame(SEBasin)
+SEBasin
+write.table(SEBasin,file=paste(datawd,"/dataEDAccm/SEBasin.csv",sep=""),sep=";",row.names=F,col.names=T)
+}}
 
 
+