From 708fdfa8b8fcb71b66f754907c71761d08f67b00 Mon Sep 17 00:00:00 2001 From: Min-Yang Lee Date: Thu, 30 Nov 2023 11:06:18 -0500 Subject: [PATCH] add lengths of a, b1, b2 by mode to catch summary. --- .../processing/catch_summaries.txt | 286 ++++++++++++++---- 1 file changed, 229 insertions(+), 57 deletions(-) diff --git a/stata_code/data_extraction_processing/processing/catch_summaries.txt b/stata_code/data_extraction_processing/processing/catch_summaries.txt index 3b0481d..a3df36c 100644 --- a/stata_code/data_extraction_processing/processing/catch_summaries.txt +++ b/stata_code/data_extraction_processing/processing/catch_summaries.txt @@ -1,6 +1,9 @@ <> - - -Effort statistics -=============================================================== - -We use Wave 1-5 data if it is available. Before Wave 5 data is available, we fill in with Wave 5 from the previous year. As always, we use November and December of <> in place of November and December of <>. And we will use Jan-April, <> in place of Jan-April, <>. - -BLAST inputs: - -1. Annual Length-frequency for cod and haddock. -2. Monthly Catch-per-trip for cod and haddock. -3. Total trips. We will use FY 2023 data, but use Waves 1 and 2 from 2023 in place of the Waves 1 and 2 of 2024 (as usual). There's no trips in Wave 6. - -For haddock, we use the survey weights to compute all catch and effort statistics (lengths and counts for a, b1, and b2). - -For cod, we use the survey weights for all catch and effort statistics *except* one. We use unweighted data to compute the length distribution of the B2s. - ~~~~ <> /* If you dyndoc this right after running your MRIP data processing code, you can comment this out.*/ @@ -52,10 +39,124 @@ local my_common2 "haddock" /*global this_working_year = 2023 */ global last_year = $this_working_year - 1 global two_years = $this_working_year - 2 + +/****************************/ +/****************************/ +/* Haddock and Cod rec ACLs in mt*/ +/****************************/ +/****************************/ + +global haddock_mort2022=666 + +global hadd_recACL2022=3634 +global cod_recACL2022=192 + + +global hadd_recACL2023=610 +global cod_recACL2023=192 + + +global mrip_vintage "2023_11_08" + +global hadd_recACL2024=0 +global cod_recACL2024=0 + <> ~~~~ -Here are monthly number of trips in Calendar Years <> through <> + + + + + +2021 Regulations +============================================================== + +Haddock: + +* 15 fish, 17" +* Open May1-Feb 28 and April 1- April 30 + +Cod: + +* 1 fish, 21" +* Open Sept 15-Sept 30 and April 1-April 14 for private anglers +* Open Sept 8-October 7 and April 1-April 14 for the for-hire anglers + +2022 Regulations +============================================================== +The recreational sub-ACL for GoM cod was <> mt for FY 2022. The recreational sub-ACL for GOM haddock was <> mt. Regulation changes from previous year in **bold**. + +Haddock: + +* **20 fish**, 17" +* Open May1-Feb 28 and April 1- April 30 + + +Cod: + +* 1 fish, **22"** +* Open Sept 1-October 7 and April 1-April 14 for private anglers +* Open Sept 1-October 7 and April 1-April 14 for the for-hire anglers + + + +2023 Regulations +============================================================== +The recreational sub-ACL for GoM cod was <> mt for FY 2023. The recreational sub-ACL for GOM haddock was <> mt. Regulation changes from previous year in **bold**. + + +Haddock: + +* For-Hire sector: **15** fish with an **18"** minimum size. +* Private Anglers: **10** fish with a 17" minimum size. +* Open May1-Feb 28 and April 1- April 30. Closed for the month of March. + +Cod: + +* 1 fish, 22" +* Open Sept 1- **October 31** for all anglers. + + +We should expect to see fewer trips in April 2024 compared to 2023 (cod not open). The May through August of 2023 should be roughly comparable to the 2022 time period, decreases in haddock regulations should reduce trips a bit. We expect to see more trips in October 2023 compared to October 2022 because cod is open for the entire month, instead of just for a week. + +Calibrating the model is tricky. I should run 1 set of models for the For-Hire sector and another set for the Private sector. I need to cross this with a set of models that is open for april and a set of models that is closed for April. An easier thing to do is just to match the April of 2023 trips and then double the number of choice occasions. + + + + +Effort statistics +=============================================================== + +We use data from Waves 1-4 (January-October). Wave 5 data is not available until mid December. We fill in with Wave 5 from the previous year, just for prototyping purposes. As always, we use November and December of <> in place of November and December of <>. And we will use Jan-April, <> in place of Jan-April, <>. + +BLAST inputs: + +1. Monthly Length-frequency for cod and haddock. +2. Monthly Catch-per-trip for cod and haddock. +3. Total trips. + +For haddock, we use the survey weights to compute all catch and effort statistics (lengths and counts for a, b1, and b2). + +For cod, we use the survey weights for all catch and effort statistics *except* one. We use unweighted data to compute the length distribution of the B2s. This is consistent with how the stock assessment is done. + + + +Here are monthly number of trips in Calendar Years <> through <>. Prior to mid-December, we fill in with Wave 5 from the previous year. ~~~~ <> @@ -117,16 +218,13 @@ collapse (sum) dtrip, by(fishing_year) list, abbreviate(16) <> ~~~~ -Data from FY <>. - - Trips by Mode =============================================================== -The regulations for Haddock were different based on mode +In FY 2023, the regulations for Haddock were different based on mode Haddock: @@ -168,6 +266,9 @@ list year month mode dtrip if mode=="ForHire", abbreviate(16) sepby(year) Trips by Disaggregated Mode =============================================================== +All Private and Shore trips are "Private". All Party and Charter are ForHire. Costs are different for each, so I need to disaggregate further. + + ~~~~ <> @@ -191,66 +292,70 @@ replace mode="Charter" if mode=="CH" replace mode="Head" if mode=="HE" replace mode="Private" if mode=="PR" replace mode="Shore" if mode=="SH" +<> +~~~~ -di "Private" +Monthly Trips for Private Vessels +~~~~ +<> list year month mode dtrip if mode=="Private", abbreviate(16) sepby(year) +<> +~~~~ -di "Headboat" -list year month mode dtrip if mode=="Head", abbreviate(16) sepby(year) +Monthly Trips for Shore +~~~~ +<> +list year month mode dtrip if mode=="Shore", abbreviate(16) sepby(year) +<> +~~~~ -di "Charter" +Monthly Trips for Headboats +~~~~ +<> +list year month mode dtrip if mode=="Head", abbreviate(16) sepby(year) +<> +~~~~ + + +Monthly Trips for Charters +~~~~ +<> list year month mode dtrip if mode=="Charter", abbreviate(16) sepby(year) +<> +~~~~ -di "Shore" -list year month mode dtrip if mode=="Shore", abbreviate(16) sepby(year) + + +All Trips by year +~~~~ +<> collapse (sum) dtrip, by(year mode) bysort year: egen t_trip=total(dtrip) gen pct=dtrip/t_trip * 100 format pct %04.2fc -sort year pct -di "Total" +gsort year - pct list year mode dtrip pct , abbreviate(16) sepby(year) +<> +~~~~ - - +All Trips, 2022 and 2023 +~~~~ +<> +keep if inlist(year,2022, 2023) collapse (sum) dtrip, by(mode) egen t_trip=total(dtrip) gen pct=dtrip/t_trip * 100 format pct %04.2fc -sort pct -di "Total" +gsort - pct list mode dtrip pct , abbreviate(16) - - - - <> ~~~~ - - - - - - - - - - - - - - - - - - - Catch Statistics for Cod =============================================================== @@ -600,6 +705,73 @@ list year month ab1mt b2mt b2dead_mt ab1_count b2_count ab1_lbs_per_fish b2_lbs_ ~~~~ +Length distributions for haddock landed. Just the Fish in the 15-19" range. There were no B2's measured on Private boats, so this assumes the B2's on private boats is similar to the B2s on ForHire and Charter vessels. + + + +<> + +use "${data_main}/MRIP_${vintage_string}/monthly/haddock_ab1_counts_${this_working_year}.dta", clear + + +merge 1:1 year month l_in_bin using "${data_main}/MRIP_${vintage_string}/monthly/haddock_b2_counts_${this_working_year}.dta" +drop _merge + + +replace ab1_count=0 if ab1_count==. +replace b2_count=0 if b2_count==. + +drop b2_UWcount + +gen total=ab1_count+b2_count +gen pct_kept=ab1_count/total*100 + +format ab1 b2 total %8.0fc + +format pct_kept %4.1fc +rename l_in_bin inches +rename ab1 landed +<> + +~~~~ +<> +list year month inches total landed pct_kept if inches>=15 & inches<=18, sepby(month) +<> +~~~~ + +Private size limit in FY 2023 is 17". MRIP data rounds down to the (a 16.9" fish is reported as 16"). Starting in May, 2023 nothing under 17" should be retained. + +There were no B2's measured on Private boats, so I've set the length distribution of the Private B2's to empty. + + +<> + +global mode_vintage 2023_11_30 + +global mode_dir "${data_main}/MRIP_$mode_vintage/mode" + +use "${mode_dir}/haddock_size_mode_class_${working_year}.dta", clear +keep if year==$working_year +rename lngcat inches +sort mode year month inches + + +replace b2=. if mode=="Private" +format ab1 b2 %8.0fc +<> + +~~~~ +<> +list month mode inches ab1_count b2 if inches>=15 & inches<=19 & mode=="Private", sepby(month) +<> +~~~~ + +For-Hire size limit in FY 2023 is 18". MRIP data rounds down to the (a 17.9" fish is reported as 17"). Starting in May, 2023 nothing under 18" should be retained. +~~~~ +<> +list month mode inches ab1_count b2 if inches>=15 & inches<=19 & mode=="ForHire", sepby(month) +<> +~~~~