Skip to content

Commit

Permalink
Changed PB and QB for OceanQuahog to match MAB model. Created resulti…
Browse files Browse the repository at this point in the history
…ng files.
  • Loading branch information
MaxGrezlik committed Aug 2, 2024
1 parent e60eabb commit 753dfe1
Show file tree
Hide file tree
Showing 8 changed files with 9 additions and 4 deletions.
13 changes: 9 additions & 4 deletions R/GB_alternate_balancing.R
Original file line number Diff line number Diff line change
Expand Up @@ -138,21 +138,26 @@ GB.params.adj$model[Group %in% c('AtlHerring', 'AtlMackerel'),
# Step 02 - Ocean Quahog -------------------------------------
GB.new <- rpath(GB.params.adj, eco.name = 'Georges Bank')
check.mort(GB.new, 'OceanQuahog')

#Changing PB and QB to match MAB
GB.params.adj$model[Group == 'OceanQuahog', PB := 0.05]
GB.params.adj$model[Group == 'OceanQuahog', QB := 0.3]

# Macrobenthos and Megabenthos. Some AmLobster and OtherDemersals
# Move predation on OceanQuahog to AtlScallop and SurfClam
GB.params.adj$diet[Group == 'OceanQuahog', Macrobenthos := 0.0002] # removed 0.005
GB.params.adj$diet[Group == 'AtlScallop', Macrobenthos := Macrobenthos + 0.0025]
GB.params.adj$diet[Group == 'SurfClam', Macrobenthos := Macrobenthos + 0.0025]

GB.params.adj$diet[Group == 'OceanQuahog', Megabenthos := Megabenthos - 0.028]
GB.params.adj$diet[Group == 'SurfClam', Megabenthos := Megabenthos + 0.028]
GB.params.adj$diet[Group == 'OceanQuahog', Megabenthos := Megabenthos - 0.03]
GB.params.adj$diet[Group == 'SurfClam', Megabenthos := Megabenthos + 0.03]

GB.params.adj$diet[Group == 'OceanQuahog', AmLobster := AmLobster - 0.02]
GB.params.adj$diet[Group == 'SurfClam', AmLobster := AmLobster + 0.01]
GB.params.adj$diet[Group == 'Detritus', AmLobster := AmLobster + 0.01]

GB.params.adj$diet[Group == 'OceanQuahog', OtherDemersals := OtherDemersals - 0.06]
GB.params.adj$diet[Group == 'OtherDemersals', OtherDemersals := OtherDemersals + 0.06]
# GB.params.adj$diet[Group == 'OceanQuahog', OtherDemersals := OtherDemersals - 0.06]
# GB.params.adj$diet[Group == 'OtherDemersals', OtherDemersals := OtherDemersals + 0.06]


#Step 03 - SouthernDemersals ----
Expand Down
Binary file modified data/alternate.GB.bal.rda
Binary file not shown.
Binary file modified data/alternate.GB.params.bal.rda
Binary file not shown.
Binary file modified outputs/GB_PB_slope.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified outputs/GB_QB_slope.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified outputs/GB_biomass_slope.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified outputs/GB_lifespan_slope.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
Binary file modified outputs/GB_prebal.png
Loading
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.

0 comments on commit 753dfe1

Please sign in to comment.