gadget-framework / gadgetutils

gadgetutils is a set of useful utilities for developing and maintaining gadget3 models
GNU General Public License v3.0
0 stars 0 forks source link

difference between input and reported catch from g3_fit #36

Open vbartolino opened 3 days ago

vbartolino commented 3 days ago

I do not seem to find a match between input catch and reported catch from g3_fit, which I can't explain. Am I missing something?

Output below is from g3_GoB/rin_ven_her_g3/gadget3/04.01-setup

> hunt.catch.rin %>% filter(fleet=="hunt.rin") %>% filter(year %in% 2019)
# A tibble: 4 × 5
   year  step area  fleet    number
  <int> <int> <chr> <chr>     <dbl>
1  2019     1 1     hunt.rin      0
2  2019     2 1     hunt.rin    378
3  2019     3 1     hunt.rin     81
4  2019     4 1     hunt.rin     81
> fit$fleet.info %>% filter(fleet=="hunt.rin") %>% filter(year %in% 2019)
# A tibble: 4 × 7
   year  step area  fleet    harv.bio amount harv.rate
  <int> <int> <chr> <chr>       <dbl>  <dbl>     <dbl>
1  2019     1 1     hunt.rin        0    0         NaN
2  2019     2 1     hunt.rin       NA  126.         NA
3  2019     3 1     hunt.rin       NA   27.0        NA
4  2019     4 1     hunt.rin       NA   27.0        NA
> fit$res.by.year %>% filter(stock %in% c("rin_imm","rin_mat","rin_allM")) %>% filter(year %in% 2019)
# A tibble: 3 × 11
  stock     year area   step total.number total.biomass harv.biomass catch num.catch        F recruitment
  <chr>    <int> <chr> <int>        <dbl>         <dbl>        <dbl> <dbl>     <dbl>    <dbl>       <dbl>
1 rin_allM  2019 1         1        5999.       334047.            0  92.3     1.76  0.000263        891.
2 rin_imm   2019 1         1        1961.        67531.            0  18.7     0.632 0               891.
3 rin_mat   2019 1         1        3859.       291595.            0  69.0     1.08  0.000263         NA
willbutler42 commented 2 days ago

@vbartolino are you using a specific branch of g3? Cannot run the setup with master

lentinj commented 2 days ago

@willbutler42 probably wip-projection-hockeyfleet, which still isn't merged (sorry!), I'm nearly there though.

willbutler42 commented 2 days ago

@lentinj Thanks!

vbartolino commented 2 days ago

sorry, I've been away the whole day (meeting with fishermen...). Yes, wip-projection-hockeyfleet