R/load_dssat.r

In palderman/dssatR: R Interface to DSSAT-CSM models

load.dssat <- function(model='cropgro',sqno=NULL,type='OUT'){
    files = c('PlantGro','PlantN','PlantC','PlantP','ETPhot','Weather','SoilWat','SoilNi','SoilPi','Dormancy','SOMLITC','SOMLITN')
    dirfiles = list.files()
    if(exists('output')) rm('output',inherits=T)
    for (i in 1:length(files)) {
        temp = NULL
        if (type=='OUT'&&any(dirfiles==paste(files[i],'.',type,sep=''))) {
            read.dssat(fileout=paste(files[i],'.OUT',sep=''),
                       model,sqno=sqno)
            temp = eval(as.name(files[i]))
        }else if(type=='R'&&any(grepl(tolower(files[i]),dirfiles))){
#            read.dssat.binary(files[i])
            if(!is.null(sqno)) colnames(temp)[1] = 'SQNO'
        }
        if(!is.null(temp)&&!exists('output')){
            assign('output',temp)
        }else if(!is.null(temp)){
            if(is.null(sqno)){
                temp = temp[,!colnames(temp)%in%colnames(output)|
                    colnames(temp)%in%c('TRNO','DATE')]
                output = merge(output,temp,by=c('TRNO','DATE'),sort=T,all=T)
            }else{
                temp = temp[,!colnames(temp)%in%colnames(output)|
                    colnames(temp)%in%'DATE']
                output = merge(output,temp,by='DATE',sort=T,all=T)
            }
        }
    }
    if ('TRNO'%in%colnames(output)){
        output=output[order(output$TRNO,output$DATE),]
    }else{
        output=output[order(output$DATE),]
    }
    assign('output',output,envir=globalenv())
}