12-14

NEFSC · Feb 14, 2025 · 3a24b5b · 3a24b5b
1 parent 672fbef
commit 3a24b5b
Show file tree

Hide file tree

Showing 14 changed files with 710 additions and 3 deletions.
diff --git a/R/get_thresholds_NRHA_atlantis.R b/R/get_thresholds_NRHA_atlantis.R
@@ -0,0 +1,41 @@
+# Use pulled survdat.rds file to get environmental thresholds for Atlantis groups
+library(dplyr)
+library(here)
+# Get output data
+NRHA_dt <-read.csv(here("data","NRHA_data.csv"))
+atlantis_groups_dt <- read.csv(here("inputs","atlantis_codes_svspp_survey_thresholds.csv"))
+# Modify data tables for ease of use
+colnames(atlantis_groups_dt)[2] <- "COMNAME"
+atlantis_groups_dt <- group_by(atlantis_groups_dt,Code)
+atlantis_NRHA_dt <- inner_join(NRHA_dt,atlantis_groups_dt,by="COMNAME")
+
+# The data includes a lot of 0's in the surface temperature data and some above 40 (104 F)
+# for the surface temperature data.  This constricts the ranges to exclude 0's and temperatures above 40.
+atlantis_NRHA_dt$BottTEMP[(atlantis_NRHA_dt$BottTEMP == 0) | (atlantis_NRHA_dt$BottTEMP > 40)] <- NA
+atlantis_NRHA_dt$SurfTEMP[(atlantis_NRHA_dt$SurfTEMP == 0) | (atlantis_NRHA_dt$SurfTEMP > 40)] <- NA
+
+atlantis_NRHA_summary_table <- group_by(atlantis_NRHA_dt,Code,SEASON)
+atlantis_NRHA_summary_table_SEASON <- summarise(atlantis_NRHA_summary_table, min_bottom_temp = min(BottTEMP, na.rm=T), max_bottom_temp = max(BottTEMP, na.rm=T),
+                                                min_surface_temp = min(SurfTEMP,na.rm=T), max_surface_temp = max(SurfTEMP,na.rm=T),
+                                                min_bottom_sal = min(BottSALIN,na.rm=T), max_bottom_sal = max(BottSALIN,na.rm=T),
+                                                min_surface_sal = min(SurfSALIN,na.rm=T), max_surface_sal = max(SurfSALIN,na.rm=T))
+
+atlantis_NRHA_summary_table_SPECIES <- group_by(atlantis_NRHA_dt,Code)
+atlantis_NRHA_summary_table_SPECIES <- summarise(atlantis_NRHA_summary_table_SPECIES, min_bottom_temp = min(BottTEMP, na.rm=T), max_bottom_temp = max(BottTEMP, na.rm=T),
+                                                 min_surface_temp = min(SurfTEMP,na.rm=T), max_surface_temp = max(SurfTEMP,na.rm=T),
+                                                 min_bottom_sal = min(BottSALIN,na.rm=T), max_bottom_sal = max(BottSALIN,na.rm=T),
+                                                 min_surface_sal = min(SurfSALIN,na.rm=T), max_surface_sal = max(SurfSALIN,na.rm=T))
+
+# Removes infinity values from the outputs and sets to NA
+atlantis_NRHA_summary_table_SEASON$min_bottom_sal[(atlantis_NRHA_summary_table_SEASON$min_bottom_sal == 'Inf') | (atlantis_NRHA_summary_table_SEASON$min_bottom_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SEASON$max_bottom_sal[(atlantis_NRHA_summary_table_SEASON$max_bottom_sal == 'Inf') | (atlantis_NRHA_summary_table_SEASON$max_bottom_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SEASON$min_surface_sal[(atlantis_NRHA_summary_table_SEASON$min_surface_sal == 'Inf') | (atlantis_NRHA_summary_table_SEASON$min_surface_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SEASON$max_surface_sal[(atlantis_NRHA_summary_table_SEASON$max_surface_sal == 'Inf') | (atlantis_NRHA_summary_table_SEASON$max_surface_sal == '-Inf')] <- 'NA'
+
+atlantis_NRHA_summary_table_SPECIES$min_bottom_sal[(atlantis_NRHA_summary_table_SPECIES$min_bottom_sal == 'Inf') | (atlantis_NRHA_summary_table_SPECIES$min_bottom_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SPECIES$max_bottom_sal[(atlantis_NRHA_summary_table_SPECIES$max_bottom_sal == 'Inf') | (atlantis_NRHA_summary_table_SPECIES$max_bottom_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SPECIES$min_surface_sal[(atlantis_NRHA_summary_table_SPECIES$min_surface_sal == 'Inf') | (atlantis_NRHA_summary_table_SPECIES$min_surface_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SPECIES$max_surface_sal[(atlantis_NRHA_summary_table_SPECIES$max_surface_sal == 'Inf') | (atlantis_NRHA_summary_table_SPECIES$max_surface_sal == '-Inf')] <- 'NA'
+
+write.csv(atlantis_NRHA_summary_table_SEASON,here("thresholds","atlantis_seasonal_thresholds_NRHA.csv"),row.names=FALSE)
+write.csv(atlantis_NRHA_summary_table_SPECIES,here("thresholds","atlantis_group_thresholds_NRHA.csv"),row.names=FALSE)
diff --git a/R/get_thresholds_survdat_atlantis.R b/R/get_thresholds_survdat_atlantis.R
@@ -4,12 +4,18 @@ library(here)
 
 survdat_dt <-readRDS(here("data","surveyPull.rds"))
 summary_table <- group_by(survdat_dt$survdat,SVSPP,SEASON)
+
+# The data includes a lot of 0's in the surface temperature data and some above 31 (the max reported temperature in the MAB)
+# for the surface temperature data.  This constricts the ranges to exclude 0's and temperatures above 31.
+
+summary_table$BOTTEMP[(summary_table$BOTTEMP == 0) | (summary_table$BOTTEMP > 31)] <- NA
+summary_table$SURFTEMP[(summary_table$SURFTEMP == 0) | (summary_table$SURFTEMP > 31)] <- NA
 summary_table_SEASON <- summarise(summary_table, min_bottom_temp = min(BOTTEMP, na.rm=T), max_bottom_temp = max(BOTTEMP, na.rm=T),
                                 min_surface_temp = min(SURFTEMP,na.rm=T), max_surface_temp = max(SURFTEMP,na.rm=T),
                                           min_bottom_sal = min(BOTSALIN,na.rm=T), max_bottom_sal = max(BOTSALIN,na.rm=T),
                                           min_surface_sal = min(SURFSALIN,na.rm=T), max_surface_sal = max(SURFSALIN,na.rm=T))
 
-summary_table_SPECIES <- group_by(survdat_dt$survdat,SVSPP)
+summary_table_SPECIES <- group_by(summary_table,SVSPP)
 summary_table_SPECIES <- summarise(summary_table_SPECIES, min_bottom_temp = min(BOTTEMP, na.rm=T), max_bottom_temp = max(BOTTEMP, na.rm=T),
                                   min_surface_temp = min(SURFTEMP,na.rm=T), max_surface_temp = max(SURFTEMP,na.rm=T),
                                   min_bottom_sal = min(BOTSALIN,na.rm=T), max_bottom_sal = max(BOTSALIN,na.rm=T),
@@ -50,5 +56,5 @@ summary_table_ATLANTIS_SPECIES$max_bottom_sal[(summary_table_ATLANTIS_SPECIES$ma
 summary_table_ATLANTIS_SPECIES$min_surface_sal[(summary_table_ATLANTIS_SPECIES$min_surface_sal == 'Inf') | (summary_table_ATLANTIS_SPECIES$min_surface_sal == '-Inf')] <- 'NA'
 summary_table_ATLANTIS_SPECIES$max_surface_sal[(summary_table_ATLANTIS_SPECIES$max_surface_sal == 'Inf') | (summary_table_ATLANTIS_SPECIES$max_surface_sal == '-Inf')] <- 'NA'
 
-write.csv(summary_table_ATLANTIS_SEASON,here("thresholds","atlantis_seasonal_thresholds.csv"),row.names=FALSE)
-write.csv(summary_table_ATLANTIS_SPECIES,here("thresholds","atlantis_group_thresholds.csv"),row.names=FALSE)
+write.csv(summary_table_ATLANTIS_SEASON,here("thresholds","atlantis_seasonal_thresholds_survdat.csv"),row.names=FALSE)
+write.csv(summary_table_ATLANTIS_SPECIES,here("thresholds","atlantis_group_thresholds_survdat.csv"),row.names=FALSE)
diff --git a/get_thresholds_NRHA_Ecomon.R b/get_thresholds_NRHA_Ecomon.R
@@ -0,0 +1,157 @@
+# Use pulled survdat.rds file to get environmental thresholds for Atlantis groups
+library(dplyr)
+library(here)
+library(geosphere)
+
+# Temporal and spatial bounds
+
+isWithinBounds <- function(CTD_loc,p1_loc,SPATIAL_BOUND) {
+
+  distance <- distm(CTD_loc,p1_loc,fun=distGeo) * 0.000621371
+  if (distance <= SPATIAL_BOUND) {
+    return(TRUE)
+  } else {
+    return(FALSE)
+  }
+}
+
+convertCTDdate <- function(date) {
+  split.points <- c(4,6,8)
+  converted_date <- substring(date,c(1,split.points + 1), c(split.points, nchar(date)))
+  converted_date <- as.Date(paste(converted_date[1],"-",converted_date[2],"-",converted_date[3],sep=""))
+  return(converted_date)
+}
+
+TEMPORAL_BOUND <- 7      # days
+SPATIAL_BOUND <- 1       # miles
+# Get output data
+NRHA_dt <-read.csv(here("data","NRHA_data.csv"))
+CTD_dt <-read.csv(here("data","CTD_data.csv"))
+
+# create empty dts to be filled for later analysts
+rows_CTD_dt <- nrow(CTD_dt)
+surfbot_dt <- CTD_dt[-c(1:rows_CTD_dt),]
+nearest_dt <- surfbot_dt
+
+# create list of uniqueIDs for CTD casts
+uniqueID_list <- unique(CTD_dt$UniqueID)[]
+
+
+# Make temporary data structure with just first (surface) and last (bottom) records, 
+# and add the first row to nearest_dt and the first and last to to surfbot_dt
+
+numCasts <- length(uniqueID_list)
+
+for(i in 1:numCasts) {
+  print(i)
+  temp_dt <- filter(CTD_dt, UniqueID == uniqueID_list[i])
+  nearest_dt <- rbind(nearest_dt,temp_dt[1,])
+  surfbot_dt <- rbind(surfbot_dt,temp_dt[1,])
+  surfbot_dt <- rbind(surfbot_dt,temp_dt[nrow(temp_dt),])
+}
+
+surfbot_dt_old <- surfbot_dt
+nearest_dt_old <- nearest_dt
+
+
+nearest_dt <- mutate(nearest_dt, DATE = convertCTDdate(Date))
+surfbot_dt$Date <- convertCTDdate(surfbot_dt$Date)
+
+surfbot_numRows <- nrow(surfbot_dt)
+
+newDates <- c(convertCTDdate(surfbot_dt$Date[1]))
+for (i in 2:surfbot_numRows) {
+  print(i)
+  newDate <- convertCTDdate(surfbot_dt$Date[i])
+  newDates <- append(newDates,newDate)
+}
+
+surfbot_dt$Date <- newDates
+atlantis_groups_dt <- read.csv(here("inputs","atlantis_codes_svspp_survey_thresholds.csv"))
+
+# Modify data tables for ease of use
+colnames(atlantis_groups_dt)[2] <- "COMNAME"
+atlantis_groups_dt <- group_by(atlantis_groups_dt,Code)
+
+atlantis_NRHA_dt <- inner_join(NRHA_dt,atlantis_groups_dt,by="COMNAME")
+
+atlantis_NRHA_tows_dt <- unique(select(atlantis_NRHA_dt,DATE,Start_Lat,Start_Lon))
+
+atlantis_NRHA_dt$SurfTEMP <- NA
+atlantis_NRHA_dt$BottTEMP <- NA
+atlantis_NRHA_dt$SurfSALIN <- NA
+atlantis_NRHA_dt$BottSALIN <- NA
+
+# Convert date columns to date object
+atlantis_NRHA_dt$DATE <- as.Date(atlantis_NRHA_dt$DATE)
+
+numRecords <- nrow(atlantis_NRHA_dt)
+
+for (n in 1:numRecords) {
+  loc_tow <- c(atlantis_NRHA_dt$Start_Lat[n], atlantis_NRHA_dt$Start_Lon[n])
+  date_tow <- as.Date(atlantis_NRHA_dt$DATE[n])
+
+  min_date <- date_tow - TEMPORAL_BOUND
+  max_date <- date_tow + TEMPORAL_BOUND
+
+  min_lat <- atlantis_NRHA_dt$Start_Lat[n] - 0.2
+  max_lat <- atlantis_NRHA_dt$Start_Lat[n] + 0.2
+  min_lon <- atlantis_NRHA_dt$Start_Lon[n] - 0.2
+  max_lon <- atlantis_NRHA_dt$Start_Lon[n] + 0.2
+
+  surfbot_dt_filtered <- filter(surfbot_dt, Date >= min_date & Date <= max_date)
+  surfbot_dt_filtered <- filter(surfbot_dt_filtered, Latitude >= min_lat & Latitude <= max_lat & Longitude >= min_lon & Longitude <= max_lon)
+  print(n)
+  numCasts <- nrow(surfbot_dt_filtered)
+
+  if (numCasts > 0) {
+    for (i in 1:numCasts) {
+      loc_ctd <- c(surfbot_dt_filtered$Latitude[i], surfbot_dt_filtered$Longitude[i])
+      if (isWithinBounds(loc_tow,loc_ctd,SPATIAL_BOUND)) {
+        print("Within Bounds")
+        atlantis_NRHA_dt$SurfTEMP[n] <- surfbot_dt_filtered$Temperature[i]
+        atlantis_NRHA_dt$SurfSALIN[n] <- surfbot_dt_filtered$Salinity[i]
+        i <- i + 1
+        atlantis_NRHA_dt$BottTEMP[n] <- surfbot_dt_filtered$Temperature[i]
+        atlantis_NRHA_dt$BottSALIN[n] <- surfbot_dt_filtered$Salinity[i]
+      } else {
+        i <- i + 1
+      }
+    }
+  }
+
+}
+
+# The data includes a lot of 0's in the surface temperature data and some above 40 (104 F)
+# for the surface temperature data.  This constricts the ranges to exclude 0's and temperatures above 40.
+
+atlantis_NRHA_dt$BottTEMP[(atlantis_NRHA_dt$BottTEMP == 0) | (atlantis_NRHA_dt$BottTEMP > 40)] <- NA
+atlantis_NRHA_dt$SurfTEMP[(atlantis_NRHA_dt$SurfTEMP == 0) | (atlantis_NRHA_dt$SurfTEMP > 40)] <- NA
+
+atlantis_NRHA_summary_table <- group_by(atlantis_NRHA_dt,Code,SEASON)
+atlantis_NRHA_summary_table_SEASON <- summarise(atlantis_NRHA_summary_table, min_bottom_temp = min(BottTEMP, na.rm=T), max_bottom_temp = max(BottTEMP, na.rm=T),
+                                  min_surface_temp = min(SurfTEMP,na.rm=T), max_surface_temp = max(SurfTEMP,na.rm=T),
+                                  min_bottom_sal = min(BottSALIN,na.rm=T), max_bottom_sal = max(BottSALIN,na.rm=T),
+                                  min_surface_sal = min(SurfSALIN,na.rm=T), max_surface_sal = max(SurfSALIN,na.rm=T))
+
+
+atlantis_NRHA_summary_table_SPECIES <- group_by(atlantis_NRHA_dt,Code)
+atlantis_NRHA_summary_table_SPECIES <- summarise(atlantis_NRHA_summary_table_SPECIES, min_bottom_temp = min(BottTEMP, na.rm=T), max_bottom_temp = max(BottTEMP, na.rm=T),
+                                   min_surface_temp = min(SurfTEMP,na.rm=T), max_surface_temp = max(SurfTEMP,na.rm=T),
+                                   min_bottom_sal = min(BottSALIN,na.rm=T), max_bottom_sal = max(BottSALIN,na.rm=T),
+                                   min_surface_sal = min(SurfSALIN,na.rm=T), max_surface_sal = max(SurfSALIN,na.rm=T))
+
+# Removes infinity values from the outputs and sets to NA
+atlantis_NRHA_summary_table_SEASON$min_bottom_sal[(atlantis_NRHA_summary_table_SEASON$min_bottom_sal == 'Inf') | (atlantis_NRHA_summary_table_SEASON$min_bottom_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SEASON$max_bottom_sal[(atlantis_NRHA_summary_table_SEASON$max_bottom_sal == 'Inf') | (atlantis_NRHA_summary_table_SEASON$max_bottom_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SEASON$min_surface_sal[(atlantis_NRHA_summary_table_SEASON$min_surface_sal == 'Inf') | (atlantis_NRHA_summary_table_SEASON$min_surface_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SEASON$max_surface_sal[(atlantis_NRHA_summary_table_SEASON$max_surface_sal == 'Inf') | (atlantis_NRHA_summary_table_SEASON$max_surface_sal == '-Inf')] <- 'NA'
+
+atlantis_NRHA_summary_table_SPECIES$min_bottom_sal[(atlantis_NRHA_summary_table_SPECIES$min_bottom_sal == 'Inf') | (atlantis_NRHA_summary_table_SPECIES$min_bottom_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SPECIES$max_bottom_sal[(atlantis_NRHA_summary_table_SPECIES$max_bottom_sal == 'Inf') | (atlantis_NRHA_summary_table_SPECIES$max_bottom_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SPECIES$min_surface_sal[(atlantis_NRHA_summary_table_SPECIES$min_surface_sal == 'Inf') | (atlantis_NRHA_summary_table_SPECIES$min_surface_sal == '-Inf')] <- 'NA'
+atlantis_NRHA_summary_table_SPECIES$max_surface_sal[(atlantis_NRHA_summary_table_SPECIES$max_surface_sal == 'Inf') | (atlantis_NRHA_summary_table_SPECIES$max_surface_sal == '-Inf')] <- 'NA'
+
+
+write.csv(atlantis_NRHA_summary_table_SEASON,here("thresholds","seasonal_thresholds_NRHA_ecomon_t7_d1.csv"),row.names=FALSE)
+write.csv(atlantis_NRHA_summary_table_SPECIES,here("thresholds","group_thresholds_NRHA_ecomon_t7_d1.csv"),row.names=FALSE)
diff --git a/inputs/atlantis_codes_svspp_survey_thresholds.csv b/inputs/atlantis_codes_svspp_survey_thresholds.csv
@@ -93,6 +93,7 @@ FDF,Cunner,176
 TAU,Tautog,177
 FDF,Oyster toadfish,185
 FDF,Wrymouth,191
+WOL,Atlantic wolffish,192
 OPT,Ocean pout,193
 FDF,Northern puffer,196
 GOO,Goosefish,197

diff --git a/thresholds/atlantis_group_thresholds_NRHA.csv b/thresholds/atlantis_group_thresholds_NRHA.csv
@@ -28,6 +28,7 @@
 "TYL",7.95,17.95,6.64,28.8,"33.118","36.359","32.057","36.43"
 "WHK",1.3,25.3,0.6,29.13,"20.25","38.78","0.05","36.389"
 "WIF",-1.4,28.15999985,-1.4,32.3,"5.7","38.78","0","36"
+"WOL",1.5,14.8,1.5,20.2,"31.294","35.051","28.3","33.877"
 "WPF",-1.2,38.1,-0.9,32,"0.5","41.5","0","36.94"
 "WSK",-1.4,23.6,-1.4,24.41538,"11.58998","36.431","11","36.43"
 "WTF",1.5,26.1,1.5,26.73,"20.25","35.894","11","36.478"

diff --git a/thresholds/atlantis_group_thresholds_survdat.csv b/thresholds/atlantis_group_thresholds_survdat.csv
@@ -46,6 +46,7 @@
 "TYL",3.25,19.66,2,28.8,"31.032","36.359","31.017","36.43"
 "WHK",1.6,25.3,0.6,29.13,"29.029","35.835","24.249","36.389"
 "WIF",-1.4,23.12,-1.4,27,"26.547","35.459","15.925","35.679"
+"WOL",1.3,14.8,0.5,21.4,"31.294","35.051","29.283","33.877"
 "WPF",-1.2,27.7,-0.9,28.7,"23.267","36.341","15.925","36.354"
 "WSK",-1.4,21.88,-1.4,25.2,"24.138","36.431","15.925","36.43"
 "WTF",1,18.53,1,26.73,"30.691","35.894","24.249","36.478"

diff --git a/thresholds/atlantis_seasonal_thresholds_NRHA.csv b/thresholds/atlantis_seasonal_thresholds_NRHA.csv
@@ -112,6 +112,10 @@
 "WIF","Spring",0.1,22.70000076,0.1,27.4,"14","35.459","4","35"
 "WIF","Summer",4.3,28.15999985,4.7,32.3,"5.7","36","2.3","36"
 "WIF","Winter",-1.4,18.9,-1.4,12.5,"7","36.19844","5","35.43703"
+"WOL","Fall",3.2,14.8,6.7,18.5,"31.935","35.051","30.345","33.395"
+"WOL","Spring",1.5,11.23,1.5,11.9,"31.294","35.048","28.3","33.877"
+"WOL","Summer",3.63,10,5.5,20.2,"32.191","34.692","30.105","32.33"
+"WOL","Winter",4.7,10.1,3.57,9,"32.112","33.603","31.815","32.08"
 "WPF","Fall",5,27.7,6.9,28.87,"0.5","41.5","0.05","36.7"
 "WPF","Spring",-0.9,38.1,-0.9,29.4,"2","36.78","0","36.6"
 "WPF","Summer",4.4,30.09,4.7,32,"3.6","36.6","4","36.94"

diff --git a/thresholds/atlantis_seasonal_thresholds_survdat.csv b/thresholds/atlantis_seasonal_thresholds_survdat.csv
@@ -214,6 +214,10 @@
 "WIF","SPRING",1.23,16,0.9,17.3,"26.547","35.459","15.925","33.888"
 "WIF","SUMMER",2.5,21,8.2,27,"30.904","32.774","30.057","32.053"
 "WIF","WINTER",-1.4,11.25,-1.4,8.51,"31.076","34.836","29.102","34.121"
+"WOL","FALL",2.9,14.8,6.7,18.5,"31.935","35.051","30.345","33.395"
+"WOL","SPRING",1.3,11.23,0.5,13.28,"31.294","35.048","29.283","33.877"
+"WOL","SUMMER",1.7,10.1,7.9,21.4,"32.191","34.692","30.105","32.33"
+"WOL","WINTER",1.8,7.91,0.8,5.1,"32.112","33.603","31.815","32.08"
 "WPF","FALL",5.34,27.7,6.3,28.7,"27.223","36.106","20.356","34.822"
 "WPF","SPRING",1.4,19.7,1.4,20.6,"23.267","36.341","15.925","36.354"
 "WPF","SUMMER",5.1,27.4,9.9,28.2,"30.904","31.594","30.87","31.561"

diff --git a/thresholds/group_thresholds_NRHA_ecomon_t14_d1.csv b/thresholds/group_thresholds_NRHA_ecomon_t14_d1.csv
@@ -0,0 +1,34 @@
+"Code","min_bottom_temp","max_bottom_temp","min_surface_temp","max_surface_temp","min_bottom_sal","max_bottom_sal","min_surface_sal","max_surface_sal"
+"BLF",5.78,28.08,5.94,27.05,"-1e+10","36.447","-1e+10","36.389"
+"BSB",3.15,28.8,3.25,27.53,"-1e+10","36.43","-1e+10","36.404"
+"BUT",2.53,29.23,3.8,27.5,"-1e+10","36.578","-1e+10","36.509"
+"COD",0.88,23.08,1.35,17.89,"-1e+10","34.085","-1e+10","35.5"
+"DOG",1.21,27.66,2.09,22.91,"-1e+10","36.578","-1e+10","36.255"
+"DRM",4.54,29.04,4.15,28.33,"-1e+10","36.578","-1e+10","36.455"
+"FDE",0.88,25.9,1.27,26.6,"-1e+10","36.031","-1e+10","36.168"
+"HAD",1.17,27.59,2.69,18.58,"-1e+10","36.089","-1e+10","35.785"
+"HAL",1.17,16.27,2.55,13.61,"-1e+10","33.64","31.663","35.001"
+"HER",0.88,25.16,1.23,21.233,"-1e+10","35.393","-1e+10","35.645"
+"LSK",0.88,28.36,0.32,23.45,"-1e+10","36.114","-1e+10","35.894"
+"MAK",1.21,23.42,2.48,22.53,"-1e+10","36.031","-1e+10","35.838"
+"OHK",3.87,29.21,4.63,15.59,"-1e+10","36.579","-1e+10","35.956"
+"OPT",1.3,27.59,1.27,19.85,"-1e+10","35.321","-1e+10","35.688"
+"PLA",1.21,22.04,1.35,17.174,"-1e+10","34.176","30.544","35.56"
+"POL",1.21,22.84,2.18,17.14,"-1e+10","35.26","-1e+10","35.512"
+"QHG",6.07,18.7199,3.57,13.05,"30.521","33.751","30.559","34.185"
+"RED",1.68,22.29,2.34,14.72,"-1e+10","34.783","31.43","35.58"
+"RHK",1.21,27.59,1.71,22.0707,"-1e+10","36.466","-1e+10","36.312"
+"SAL",Inf,-Inf,Inf,-Inf,"NA","NA","NA","NA"
+"SCU",4.25,29.23,4.12,28.33,"-1e+10","36.478","-1e+10","36.5"
+"SDF",6.05,22.92,5.33,14.08,"20.77","34.645","24.49","35.505"
+"SHK",1.21,29.21,1.74,24.46,"-1e+10","36.578","-1e+10","36.2"
+"SK",1.17,29.76,2.18,27.5,"-1e+10","36.579","-1e+10","36.5"
+"SUF",1.55,29.4,1.74,28.22,"-1e+10","36.558","-1e+10","36.554"
+"TAU",2.89,24.94985,2.84,23.03,"-1e+10","33.582","-1e+10","33.818"
+"TYL",6.42,28.51,7.95,17.95,"-1e+10","36.407","33.118","36.359"
+"WHK",1.64,26.68,2.92,18.01,"-1e+10","36.409","29.206","35.845"
+"WIF",0.88,24.47,0.32,23.12,"-1e+10","35.209","-1e+10","35.187"
+"WPF",0.88,28.58,0.32,26.7,"-1e+10","36.089","-1e+10","35.587"
+"WSK",0.88,23.77,0.32,21.88,"-1e+10","36.089","-1e+10","35.764"
+"WTF",1.55,27.17,2.09,18.33,"-1e+10","36.409","-1e+10","35.894"
+"YTF",1.17,25.21,0.32,20.02,"-1e+10","35.772","-1e+10","35.645"