READ_CTD.R

read_CTD <- function ( filenc_core ){

##################################################
#### DATA MODE 
##################################################

DATA_MODE=unlist(strsplit(ncvar_get(filenc_core,"DATA_MODE"),split=""))

DATA_MODE_CTD=DATA_MODE[1]

##################################################
#### Get the Best CTD
##################################################

if ( DATA_MODE_CTD == "R" ) {

	PSAL=ncvar_get(filenc_core,"PSAL")
	TEMP=ncvar_get(filenc_core,"TEMP")
	PRES=ncvar_get(filenc_core,"PRES")

} else {
	
	PSAL=ncvar_get(filenc_core,"PSAL_ADJUSTED")
	TEMP=ncvar_get(filenc_core,"TEMP_ADJUSTED")
	PRES=ncvar_get(filenc_core,"PRES_ADJUSTED")

}

#### CTD is always the 1st profile of the nc
if (filenc_core$dim$N_PROF$len == 1 ) {

        PSAL_CTD=PSAL
        TEMP_CTD=TEMP
        PRES_CTD=PRES

} else {

        PSAL_CTD=PSAL[,1]
        TEMP_CTD=TEMP[,1]
        PRES_CTD=PRES[,1]

}

result=(list(PRES=PRES_CTD,PSAL=PSAL_CTD,TEMP=TEMP_CTD))
#result=(list(PRES=PRES,PSAL=PSAL,TEMP=TEMP))
return(result)

}