03_stitch_MOM_100sim.R

library(MSEtool)
library(dplyr)

source("99_MOM_fixes.R")
source("MOM_fn.R")

byc <- c("BSH", "SMA", "WHM", "BUM")
targ <- c("BET", "SWO")

multiHist_byc <- lapply(paste0('MOM/multiHist_', byc, '_100sim.rds'), readRDS)
multiHist_targ <- lapply(paste0('MOM/multiHist_', targ, '_100sim.rds'), readRDS)

MOM_byc <- lapply(paste0('MOM/MOM_', byc, '_100sim.rds'), readRDS)
MOM_targ <- lapply(paste0('MOM/MOM_', targ, '_100sim.rds'), readRDS)

LL_byc <- list(1:9, c(1:12)[-c(8, 10, 12)], 2, 2) %>% structure(names = byc)
LL_targ <- list(10:18, 1:10) %>% structure(names = targ)

M2_byc <- Map(function(x, y) aggregate_fleet(x = x, LL = y), x = MOM_byc, y = LL_byc)
M2_targ <- Map(function(x, y) aggregate_fleet(x = x, LL = y), x = MOM_targ, y = LL_targ)

MOM <- local({
  args <- c(M2_targ, M2_byc) %>% structure(names = c(targ, byc))
  do.call("MOM_stitch", args)
})


#lapply(args,function(x)x@cpars[[1]][[1]]$Mat_age[1,,1])
lapply(MOM@cpars,function(x)print(x[[1]]$Mat_age[1,,1]))
#lapply(args,function(x)x@cpars[[1]][[1]]$Len_age[1,,1])
lapply(MOM@cpars,function(x)print(x[[1]]$Len_age[1,,1]))

MOM2 = MOM_simplify(MOM)
MOM3 = fix_selectivity_1(MOM2) #  MOM3@cpars[[1]][[1]]$V[1,1:4,1:10] # first 10 years of sim 1 & CAL_bins and CAL_binsmid fix

class(MOM3@cpars[[1]][[1]]$CAL_bins)
lapply(MOM3@cpars,function(x)x[[1]]$CAL_bins)
lapply(MOM3@cpars,function(x)print(dim(x[[1]]$SLarray)))

saveRDS(MOM3, file = "MOM/MOM_latest.rds")


# === make a current effort MP and results ============================================

library(openMSE)
library(dplyr)

# --- Current effort MMP -------------------------------------------------------

E1 <- function(x, DataList, reps = 1, ...) {
  np <- length(DataList)
  nf <- length(DataList[[1]])
  RecList <- lapply(1:np, function(p) replicate(nf, new("Rec")))
  for(p in 1:np) { 
    for(f in 1:nf) { 
      RecList[[p]][[f]]@Effort <- 1
    }
  }
  return(RecList)
}
class(E1) = "MMP"



# - 2 - Projection ---------------------------------------------------------------

Hist = readRDS('MOM/multiHist_stitch_100sim.rds')
MMSE = ProjectMOM(Hist, MPs = 'E1', checkMPs = FALSE)
saveRDS(MMSE, 'MOM/MMSE_E1.rda')












#MOM <- readRDS(file = "MOM/MOM_stitch_100sim.rds")

setup(8)
multiHist <- multiMSE(MOM, Hist = TRUE, checkMPs = FALSE, parallel = TRUE)
sfStop()
saveRDS(multiHist, file = "MOM/multiHist_stitch_100sim.rds")
#multiHist <- readRDS(file = "MOM/multiHist_stitch_100sim.rds")