Fix bug in GDAC not being propogated

rscripts/split_by_type.R

@@ -113,7 +113,8 @@ for ( year in my_args[1]:my_args[2] ) {
                 mdata.gdac$pt <- 5 # all GCC data should be from ships
                 mdata.gdac$dups <- NA
                 mdata.gdac$dck <- 0 
-
+                mdata.gdac$irf <- 0
+                
                 ## Fill missing columns with NA
                 miss.cols <- colnames(mdata)[which(!colnames(mdata)%in%colnames(mdata.gdac))]
                 mdata.gdac[,miss.cols] <- NA
@@ -266,7 +267,7 @@ for ( year in my_args[1]:my_args[2] ) {
         mdata$data.type <- ifelse(mdata$dck == 151 & mdata$lon >= 260 & mdata$lon <= 270 & mdata$yr >= 1890 & mdata$yr <= 1919, "reject", mdata$data.type )
 
                                         #print(table(mdata$data.type))
-
+        
         if ( max(mdata$pt) == 99 ) {
             cat("PT problem")
             print(table(mdata$pt,mdata$dck))