diff --git a/SS_expval.tpl b/SS_expval.tpl
index df6fb5eb..d5b7d798 100644
--- a/SS_expval.tpl
+++ b/SS_expval.tpl
@@ -666,13 +666,15 @@ FUNCTION void Get_expected_values(const int y, const int t);
               {
                 if (SzFreq_Nmeth > 0) //  have some sizefreq data
                 {
-
                   if (have_data(ALK_time, f, data_type, 0) > 0)
                   {
                     for (j = 1; j <= have_data(ALK_time, f, data_type, 0); j++) // loop all obs of this type
                     {
                       iobs = have_data(ALK_time, f, data_type, j); //  observation index
                       SzFreqMethod = SzFreq_obs_hdr(iobs, 6);
+
+    {  
+      echoinput<<" begin szfreq legacy method"<<endl;
                       SzFreqMethod_seas = nseas * (SzFreqMethod - 1) + s; // index that combines sizefreqmethod and season and used in SzFreqTrans
                       if (SzFreq_obs_hdr(iobs, 9) > 0) // first occurrence of this method at this time is with fleet = f
                       {
@@ -699,7 +701,328 @@ FUNCTION void Get_expected_values(const int y, const int t);
 
                             //  NOTE:  wt_len_low is  calculated separately for each growth pattern (GPat)
                             //  but the code below still just uses GPat=1 for calculation of the sizefreq transition matrix
+    echoinput<<" ready "<<SzFreq_units(SzFreqMethod)<< " z1 "<<z1 << " "<<z2<<endl;
+                            switch (SzFreq_units(SzFreqMethod)) // biomass vs. numbers are accumulated in the bins
+                            {
+                              case (1): // units are biomass, so accumulate body weight into the bins;  Assume that bin demarcations are also in biomass
+                              {
+                                if (SzFreq_Omit_Small(SzFreqMethod) == 1)
+                                {
+                                  while (wt_len_low(s, 1, z1 + 1) < SzFreq_bins(SzFreqMethod, 1) && z1 < z2)
+                                  {
+                                    z1++;
+                                  }
+                                } // ignore tiny fish
+                                if (z1 + 1 >= z2)
+                                {
+                                  warnstream << "max population size " << wt_len_low(s, 1, z1) << " is less than first data bin " << SzFreq_bins(SzFreqMethod, 1) << " for SzFreqMethod " << SzFreqMethod;
+                                  write_message (FATAL, 0); // EXIT!
+                                }
+
+                                if (wt_len_low(s, 1, nlength2) < SzFreq_bins(SzFreqMethod, SzFreq_Nbins(SzFreqMethod)))
+                                {
+                                  warnstream << "max population size " << wt_len_low(s, 1, nlength2) << " is less than max data bin " << SzFreq_bins(SzFreqMethod, SzFreq_Nbins(SzFreqMethod)) << " for SzFreqMethod " << SzFreqMethod;
+                                  write_message (FATAL, 0); // EXIT!
+                                }
+                                for (z = z1; z <= z2; z++)
+                                {
+                                  if (ibin == SzFreq_Nbins(SzFreqMethod))
+                                  {
+                                    SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = wt_len2(s, 1, z);
+                                  }
+                                  else
+                                  {
+                                    if (wt_len_low(s, 1, z) >= topbin)
+                                    {
+                                      ibin++;
+                                      ibinsave++;
+                                    }
+                                    if (ibin > 1)
+                                    {
+                                      botbin = SzFreq_bins2(SzFreqMethod, ibin);
+                                    }
+                                    if (ibin == SzFreq_Nbins(SzFreqMethod))
+                                    {
+                                      SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = wt_len2(s, 1, z);
+                                      topbin = 99999.;
+                                    }
+                                    else
+                                    {
+                                      topbin = SzFreq_bins2(SzFreqMethod, ibin + 1);
+                                      if (wt_len_low(s, 1, z) >= botbin && wt_len_low(s, 1, z + 1) <= topbin)
+                                      {
+                                        SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = wt_len2(s, 1, z);
+                                      }
+                                      else
+                                      {
+                                        temp = (wt_len_low(s, 1, z + 1) - topbin) / wt_len_fd(s, 1, z); // frac in pop bin above (data bin +1)
+                                        temp1 = wt_len_low(s, 1, z) + (1. - temp * 0.5) * wt_len_fd(s, 1, z); // approx body wt for these fish
+                                        temp2 = wt_len_low(s, 1, z) + (1. - temp) * 0.5 * wt_len_fd(s, 1, z); // approx body wt for  fish below
+                                        SzFreqTrans(SzFreqMethod_seas, z, ibinsave + 1) = temp * temp1;
+                                        SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = (1. - temp) * temp2;
+                                      }
+                                    }
+                                  }
+                                }
+                                if (SzFreq_scale(SzFreqMethod) == 2 && gg == gender) // convert to pounds
+                                {
+                                  SzFreqTrans(SzFreqMethod_seas) /= 0.4536;
+                                }
+                                break;
+                              } //  end of units in biomass
+                              // NOTE: even though  the transition matrix is currently in units of biomass distribution, there is no need to
+                              // normalize to sum to 1.0 here because the normalization will occur after it gets used to create SzFreq_exp
 
+                              case (2): // units are numbers
+                              {
+                                  echoinput<<" ibin "<<ibin <<" max "<<SzFreq_Nbins(SzFreqMethod)<<endl;
+                                if (SzFreq_scale(SzFreqMethod) <= 2) //  bin demarcations are in weight units (1=kg, 2=lbs), so uses wt_len to compare to bins
+                                {
+                                  if (SzFreq_Omit_Small(SzFreqMethod) == 1)
+                                  {
+                                    while (wt_len_low(s, 1, z1 + 1) < SzFreq_bins(SzFreqMethod, 1) && z1 < z2)
+                                    {
+                                      z1++;
+                                    }
+                                  } // ignore tiny fish
+                                  if (z1 + 1 >= z2)
+                                  {
+                                    warnstream << "max population size " << wt_len_low(s, 1, z1) << " is less than first data bin " << SzFreq_bins(SzFreqMethod, 1) << " for SzFreqMethod " << SzFreqMethod;
+                                    write_message (FATAL, 0); // EXIT!
+                                  }
+                                  if (wt_len_low(s, 1, nlength2) < SzFreq_bins(SzFreqMethod, SzFreq_Nbins(SzFreqMethod)))
+                                  {
+                                    warnstream << "max population size " << wt_len_low(s, 1, nlength2) << " is less than max data bin " << SzFreq_bins(SzFreqMethod, SzFreq_Nbins(SzFreqMethod)) << " for SzFreqMethod " << SzFreqMethod;
+                                    write_message (FATAL, 0); // EXIT!
+                                  }
+                                  for (z = z1; z <= z2; z++)
+                                  {
+//                                    echoinput<<"ibin "<<ibin<<" z "<<z<<" popwt_low "<<wt_len_low(s, 1, z);
+                                    if (ibin == SzFreq_Nbins(SzFreqMethod))
+                                    {
+                                      SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = 1.;
+                                    } //checkup<<" got to last ibin, so put rest of popbins here"<<endl;
+                                    else
+                                    {
+                                      if (wt_len_low(s, 1, z) >= topbin)
+                                      {
+                                        ibin++;
+                                        ibinsave++;
+                                      } 
+                                      if (ibin > 1)
+                                      {
+                                        botbin = SzFreq_bins2(SzFreqMethod, ibin);
+                                      }
+                                      if (ibin == SzFreq_Nbins(SzFreqMethod)) // checkup<<" got to last ibin, so put rest of popbins here"<<endl;
+                                      {
+                                        SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = 1.;
+                                        topbin = 99999.;
+                                      }
+                                      else
+                                      {
+                                        topbin = SzFreq_bins2(SzFreqMethod, ibin + 1);
+                                        if (wt_len_low(s, 1, z) >= botbin && wt_len_low(s, 1, z + 1) <= topbin) //checkup<<" pop inside dat, put here"<<endl;
+                                        {
+                                          SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = 1.;
+                                        }
+                                        else // checkup<<" overlap"<<endl;
+                                        {
+                                          SzFreqTrans(SzFreqMethod_seas, z, ibinsave + 1) = (wt_len_low(s, 1, z + 1) - topbin) / wt_len_fd(s, 1, z);
+                                          SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = 1. - SzFreqTrans(SzFreqMethod_seas, z, ibinsave + 1);
+                                       }
+//                                        echoinput<<" result "<<SzFreqTrans(SzFreqMethod_seas, z, ibinsave) << " "<<SzFreqTrans(SzFreqMethod_seas, z, ibinsave + 1)<<endl;
+                                     }
+                                    }
+                                  }
+                                }
+
+                                else //  bin demarcations are in length unit (3=cm, 4=inch) so uses population len_bins to compare to data bins
+                                {
+                                  if (SzFreq_Omit_Small(SzFreqMethod) == 1)
+                                  {
+                                    while (len_bins2(z1 + 1) < SzFreq_bins(SzFreqMethod, 1))
+                                    {
+                                      z1++;
+                                    }
+                                    //  echoinput<<"accumulate starting at bin: "<<z1<<endl;
+                                  } // ignore tiny fish
+
+                                  if (len_bins2(nlength2) < SzFreq_bins(SzFreqMethod, SzFreq_Nbins(SzFreqMethod)))
+                                  {
+                                    warnstream << "max population len bin " << len_bins2(nlength2) << " is less than max data bin " << SzFreq_bins(SzFreqMethod, SzFreq_Nbins(SzFreqMethod)) << " for SzFreqMethod " << SzFreqMethod;
+                                    write_message (FATAL, 0); // EXIT!
+                                  }
+
+                                  for (z = z1; z <= z2; z++)
+                                  {
+                                    if (ibin == SzFreq_Nbins(SzFreqMethod))
+                                    {
+                                      SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = 1.;
+                                    } //checkup<<" got to last ibin, so put rest of popbins here"<<endl;
+                                    else
+                                    {
+                                      if (len_bins2(z) >= topbin)
+                                      {
+                                        ibin++;
+                                        ibinsave++;
+                                      } //checkup<<" incr ibin "<<z<<" "<<ibin<<" "<<len_bins(z)<<" "<<len_bins_dat(ibin);
+                                      if (ibin > 1)
+                                      {
+                                        botbin = SzFreq_bins2(SzFreqMethod, ibin);
+                                      }
+                                      if (ibin == SzFreq_Nbins(SzFreqMethod)) // checkup<<" got to last ibin, so put rest of popbins here"<<endl;
+                                      {
+                                        SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = 1.;
+                                        topbin = 99999.;
+                                      }
+                                      else
+                                      {
+                                        topbin = SzFreq_bins2(SzFreqMethod, ibin + 1);
+                                        if (len_bins2(z) >= botbin && len_bins2(z + 1) <= topbin) //checkup<<" pop inside dat, put here"<<endl;
+                                        {
+                                          SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = 1.;
+                                        }
+                                        else // checkup<<" overlap"<<endl;
+                                        {
+                                          SzFreqTrans(SzFreqMethod_seas, z, ibinsave + 1) = (len_bins2(z + 1) - topbin) / (len_bins2(z + 1) - len_bins2(z));
+                                          SzFreqTrans(SzFreqMethod_seas, z, ibinsave) = 1. - SzFreqTrans(SzFreqMethod_seas, z, ibinsave + 1);
+                                        }
+                                      }
+                                    }
+                                  }
+                                }
+                                break;
+                              } //  end of units in numbers
+                            }
+                            if (docheckup == 1 && gg == gender)
+                              echoinput << " sizefreq trans_matrix: method/season " << SzFreqMethod << " / " << s << endl
+                                        << trans(SzFreqTrans(SzFreqMethod_seas)) << endl;
+                          } // end gender loop
+                        } //  end needing to calc the matrix because it may have changed
+                      }  //  end first occurrence
+                      switch (SzFreq_obs_hdr(iobs, 5)) // discard/retained partition
+                      {
+                        case (0):
+                        {
+                          SzFreq_exp(iobs) = trans(SzFreqTrans(SzFreqMethod_seas)) * exp_l_temp;
+                          echoinput<<" sum_old: "<<sum(exp_l_temp)<<" exP " <<sum(SzFreq_exp(iobs))<<endl;
+                          break;
+                        }
+                        case (1):
+                        {
+                          SzFreq_exp(iobs) = trans(SzFreqTrans(SzFreqMethod_seas)) * (exp_l_temp - exp_l_temp_ret);
+                          break;
+                        }
+                        case (2):
+                        {
+                          SzFreq_exp(iobs) = trans(SzFreqTrans(SzFreqMethod_seas)) * exp_l_temp_ret;
+                          break;
+                        }
+                      }
+  #ifdef DO_ONCE
+                      if (do_once == 1)
+                        echoinput << y << " " << f << " szfreq_exp_initial  " << SzFreq_exp(iobs) << endl;
+  #endif
+
+                      if (gender == 2)
+                      {
+                        k = SzFreq_obs_hdr(iobs, 8); // max bins for this method
+                        switch (SzFreq_obs_hdr(iobs, 4)) //  combine, select or each gender
+                        {
+                          case (0): // combine genders
+                          {
+                            for (ibin = 1; ibin <= k; ibin++)
+                              SzFreq_exp(iobs, ibin) += SzFreq_exp(iobs, k + ibin);
+                            SzFreq_exp(iobs)(k + 1, 2 * k) = 0.0;
+                            SzFreq_exp(iobs)(1, k) /= sum(SzFreq_exp(iobs)(1, k));
+                            if (SzFreq_mincomp(SzFreqMethod) > 0.0)
+                            {
+                              SzFreq_exp(iobs)(1, k) += SzFreq_mincomp(SzFreqMethod);
+                              SzFreq_exp(iobs)(1, k) /= sum(SzFreq_exp(iobs)(1, k));
+                            }
+                            break;
+                          }
+                          case (1): // female only
+                          {
+                            SzFreq_exp(iobs)(k + 1, 2 * k) = 0.0; //  zero out the males so will not interfere with data generation
+                            SzFreq_exp(iobs)(1, k) /= sum(SzFreq_exp(iobs)(1, k));
+                            if (SzFreq_mincomp(SzFreqMethod) > 0.0)
+                            {
+                              SzFreq_exp(iobs)(1, k) += SzFreq_mincomp(SzFreqMethod);
+                              SzFreq_exp(iobs)(1, k) /= sum(SzFreq_exp(iobs)(1, k));
+                            }
+                            break;
+                          }
+                          case (2): //   male only
+                          {
+                            ibin = SzFreq_obs_hdr(iobs, 7);
+                            SzFreq_exp(iobs)(1, ibin - 1) = 0.0; //  zero out the females so will not interfere with data generation
+                            SzFreq_exp(iobs)(ibin, k) /= sum(SzFreq_exp(iobs)(ibin, k));
+                            if (SzFreq_mincomp(SzFreqMethod) > 0.0)
+                            {
+                              SzFreq_exp(iobs)(ibin, k) += SzFreq_mincomp(SzFreqMethod);
+                              SzFreq_exp(iobs)(ibin, k) /= sum(SzFreq_exp(iobs)(ibin, k));
+                            }
+                            break;
+                          }
+                          case (3): //  each gender
+                          {
+                            SzFreq_exp(iobs) /= sum(SzFreq_exp(iobs));
+                            if (SzFreq_mincomp(SzFreqMethod) > 0.0)
+                            {
+                              SzFreq_exp(iobs) += SzFreq_mincomp(SzFreqMethod);
+                              SzFreq_exp(iobs) /= sum(SzFreq_exp(iobs));
+                            }
+                            break;
+                          }
+                        } //  end gender switch
+                      } // end have 2 genders
+                      else
+                      {
+                        k = SzFreq_obs_hdr(iobs, 8); // max bins for this method
+                        SzFreq_exp(iobs)(1, k) /= sum(SzFreq_exp(iobs)(1, k));
+                        if (SzFreq_mincomp(SzFreqMethod) > 0.0)
+                        {
+                          SzFreq_exp(iobs)(1, k) += SzFreq_mincomp(SzFreqMethod);
+                          SzFreq_exp(iobs)(1, k) /= sum(SzFreq_exp(iobs)(1, k));
+                        }
+                      }
+  #ifdef DO_ONCE
+                      if (do_once == 1)
+                        echoinput << y << " " << f << " szfreq_exp_final  " << SzFreq_exp(iobs) << endl;
+  #endif
+
+                    }  //  end szfreqtrans method
+
+                    {  //  begin rebin method
+                  echoinput<<"begin rebin method"<<endl;
+                  SzFreqMethod_seas = nseas * (SzFreqMethod - 1) + s; // index that combines sizefreqmethod and season and used in SzFreqTrans
+                      if (SzFreq_obs_hdr(iobs, 9) > 0) // first occurrence of this method at this time is with fleet = f
+                      {
+                        if (do_once == 1 || (MG_active(3) > 0 && (timevary_MG(y, 3) > 0))) // calc  matrix because wtlen parameters have changed
+                        {
+                          for (gg = 1; gg <= gender; gg++)
+                          {
+                            if (gg == 1)
+                            {
+                              z1 = 1;
+                              z2 = nlength;
+                              ibin = 0;
+                              ibinsave = 0;
+                            } // female
+                            else
+                            {
+                              z1 = nlength1;
+                              z2 = nlength2;
+                              ibin = 0;
+                              ibinsave = SzFreq_Nbins(SzFreqMethod);
+                            } // male
+                            topbin = 0.;
+                            botbin = 0.;
+
+                            //  NOTE:  wt_len_low is  calculated separately for each growth pattern (GPat)
+                            //  but the code below still just uses GPat=1 for calculation of the sizefreq transition matrix
+ /*
                             switch (SzFreq_units(SzFreqMethod)) // biomass vs. numbers are accumulated in the bins
                             {
                               case (1): // units are biomass, so accumulate body weight into the bins;  Assume that bin demarcations are also in biomass
@@ -894,10 +1217,68 @@ FUNCTION void Get_expected_values(const int y, const int t);
                             if (docheckup == 1 && gg == gender)
                               echoinput << " sizefreq trans_matrix: method/season " << SzFreqMethod << " / " << s << endl
                                         << trans(SzFreqTrans(SzFreqMethod_seas)) << endl;
-                          } // end gender loop
-                        } //  end needing to calc the matrix because it may have changed
-                      } // done calculating the SzFreqTransition matrix for this method
+  */
+    dvector src_edges(1, nlength + 1); // +1 to create a dummy bin holding upper edge of the actual last bin
+    src_edges(1, nlength) = len_bins(1, nlength); 
+    src_edges(nlength + 1) = src_edges(nlength) + (src_edges(nlength) - src_edges(nlength - 1));
+    echoinput << "src_edges "<<src_edges<<endl;
+    echoinput << "SzFreq_bins2 " << SzFreq_bins2 << endl;
+    dvar_vector dest_edges(1,SzFreq_Nbins(SzFreqMethod) + 1);  // +1 to create a dummy bin holding upper edge of the actual last bin
+    dest_edges.initialize();
+    for (ibin = 1; ibin <= SzFreq_Nbins(SzFreqMethod); ibin++) {
+      dest_edges(ibin) = SzFreq_bins(SzFreqMethod, ibin);
+    }
+    if(SzFreq_Omit_Small(SzFreqMethod) < 0.0)
+      {dest_edges(1) = 0.;} // so all small fish go into this bin
+    dest_edges(SzFreq_Nbins(SzFreqMethod) + 1) = 99999.; // dest_edges(SzFreq_Nbins(SzFreqMethod)) + (dest_edges(SzFreq_Nbins(SzFreqMethod)) - dest_edges(SzFreq_Nbins(SzFreqMethod) - 1));
+    dvar_vector dest_comp(1,SzFreq_Nbins(SzFreqMethod));  //  destination composition.  one-sex at a time
+// convert dest_edges to cm
+    switch (SzFreq_scale(SzFreqMethod)) // biomass vs. numbers are accumulated in the bins
+    {
+      case (1): // scale is kg
+      {
+        // bin boundary in inch have already been converted to cm, so no further action needed
+        // fall thru to calcs in the lb case
+      }
+      case (2): // scale is lb
+      {
+// bin boundary in lb have already been converted to kg
+// wt_len(s, gp) is potentially by season and by Gpat within sex.  Here, only s==1 and Gpat==1 will be used.  It will be sex-specific
+// wt_len_low(s, GPat)(1, nlength) = wtlen_p(GPat, 1) * pow(len_bins2(1, nlength), wtlen_p(GPat, 2));
+// W=a*L^b; so L = exp((ln(W)-ln(a))/B)
+        echoinput << "sex: "<<gg << " wtlen_p " << wtlen_p(GPat) << "  z1, z2 "<<z1<<" "<<z2<<endl;
+        dest_edges = exp((log(dest_edges)-log(wtlen_p(GPat, 1)))/wtlen_p(GPat, 2));
+        echoinput << "dest_edges_in_cm "<<dest_edges<<endl;
+        break;
+      }
+    }  // end calc of new dest_edges  .shift  shallow
+    echoinput << "src_comps "<<exp_l_temp(z1, z2)<<endl<<"  sum_src: " << sum(exp_l_temp(z1, z2))<< endl;
+    dvar_vector pass_comp = exp_l_temp(z1, z2);
+    pass_comp.shift(1);
+    echoinput<<"ready to call rebin min-max: "<<pass_comp<<endl;
+    dest_comp = rebin(src_edges, pass_comp, dest_edges);
+    echoinput << "dest_comps "<<dest_comp<<endl<<"  sum_dest: " << sum(dest_comp)<< endl;
+    } // end gender loop
+  } //  end needing to calc the matrix because it may have changed
+  } // done calculating the SzFreqTransition matrix for this method
+//  practice call to rebin
+// FUNCTION dvector rebin(const dvector& src_edges, const dvector& src_counts, const dvector& dest_edges)
+//  srcedges (0, SzFreq_Nbins(SzFreqMethod)-1) corresponds to len_bins2
+//  dest_edges are SzFreq_bins2(1, SzFreq_Nbins(SzFreqMethod))
+//  nlength2 = gender * nlength; // doubled for males
+//  vector len_bins2(1,nlength2);  //vector with lower edge of population length bins
+    
 
+//  with rebin, each sex will need to be done separately, not as an end-to-end vector
+//  with rebin, the mapping of inputs is:  each sex's range of exp_l_temp goes to scr_counts.  Fact that it is not "counts" does not matter
+//  for females:   z1 = 1; z2 = nlength;
+//  for males:   z1 = nlength1; z2 = nlength2;
+//  imatrix sexbins(1,2,1,2);  // fill with 1,nlength,nlength1,nlength2
+//  with SzFreqTrans, the conversion of len to wt has already occurred in the creation of SzFreqTrans, so the numbers based exp_l_temp gets converted to weights as SzFreqTrans is applied
+//  the logic that created SzFreqTrans needs to be converted into logic that finds the length (in cm) that corresponds to the destination bin boundaries.
+//  This allows rebin to go from cm to cm
+//  with rebin it seems better to convert the numbers to weight to create an exp_wt_temp, then apply rebin to parse those weights to the new bins, using the bin boundaries that are in cm
+//accumulate body weight into the bins and 
                       switch (SzFreq_obs_hdr(iobs, 5)) // discard/retained partition
                       {
                         case (0):
@@ -988,6 +1369,7 @@ FUNCTION void Get_expected_values(const int y, const int t);
                       if (do_once == 1)
                         echoinput << y << " " << f << " szfreq_exp_final  " << SzFreq_exp(iobs) << endl;
   #endif
+                      }  // end rebin method
                     } // end loop of obs for fleet = f
                   } //  end having some obs for this method in this fleet
                 } //  end use of generalized size freq data
diff --git a/SS_miscfxn.tpl b/SS_miscfxn.tpl
index 76964cb3..f67c059d 100644
--- a/SS_miscfxn.tpl
+++ b/SS_miscfxn.tpl
@@ -120,3 +120,62 @@ FUNCTION dvariable Comp_logL_Dirichlet(const double& Nsamp, const dvariable& dir
     logL = sum(gammln(Nsamp * obs_comp + dirichlet_Parm * exp_comp)) - sum(gammln(dirichlet_Parm * exp_comp));
     return (logL);
   }
+
+  
+FUNCTION dvar_vector rebin(const dvector& src_edges, const dvar_vector& src_counts, const dvar_vector& dest_edges)
+  {
+    /*
+    This implementation takes two vectors representing the boundaries (edges) of the source
+    and destination bins, and one vector for the source counts.
+
+    Rebins frequency data from one set of boundaries to another.
+    @param src_edges Boundaries of the original bins (size N+1).
+    @param src_counts Frequency/counts in original bins (size N).
+    @param dest_edges Boundaries of the new bins (size M+1).
+    @return Vector of rebinned frequency data (size M).
+    the original code searched all source bins for each destination bins.
+    Here the ordered characteristic of the bins allows for the search for bin (i+1) to continue from search for bin(i)
+    */
+
+    dvar_vector dest_counts(1, dest_edges.size() - 1);  // size to leave off the topbin bounary
+    dest_counts.initialize();
+    echoinput<<"in rebin: "<<endl<<"dest_edges: "<<dest_edges<<endl;
+    int src_startbin = 1;
+    for (int i = 1; i <= dest_counts.size(); i++) 
+    {
+      int src_stopflag = 0;
+      j = src_startbin;
+//      echoinput<<"dest_bin "<<i<<"  src_start "<<src_startbin<<endl;
+      while(src_stopflag < 2 && j <= src_counts.size())
+      {
+        // Calculate the overlap between [d_low, d_high] and [s_low, s_high]
+        dvariable overlap_low = dest_edges[i];
+        dvariable overlap_high = dest_edges[i + 1];
+        if( src_edges[j] > dest_edges[i]) overlap_low = src_edges[j];
+        if(src_edges[j + 1] < dest_edges[i + 1]) overlap_high = src_edges[j + 1];
+//        echoinput<<" "<<j<<" "<<overlap_low<<" "<<overlap_high;
+
+        if (overlap_low < overlap_high)
+        {
+          dest_counts[i] += src_counts[j] * ((overlap_high - overlap_low)) / (src_edges[j + 1] - src_edges[j]);
+          src_startbin = j;  // so start with this bin next time; allows for overlap
+          src_stopflag = 1;  // got at least one bin assigned
+          // Distribute source count proportionally to the overlap area
+          // echoinput<<" store "<<endl;
+
+        }
+        else if (src_stopflag = 1)
+        {  // triggers end of seach for this i
+          // echoinput<<" stop_here: "<<j<<endl;
+          src_stopflag = 2;
+        }
+        else
+        {
+          //  echoinput<<"next bin "<<endl;
+        }
+        j++;
+      }
+    }
+//    echoinput<<"finish rebin "<<dest_counts<<endl<"<"sum_dest: "<<sum(dest_counts)<<endl;
+    return (dest_counts);
+  }
\ No newline at end of file
diff --git a/SS_readdata_330.tpl b/SS_readdata_330.tpl
index b8fac5ff..c06c1838 100644
--- a/SS_readdata_330.tpl
+++ b/SS_readdata_330.tpl
@@ -3136,7 +3136,7 @@
   init_ivector SzFreq_nobs(1,SzFreq_Nmeth);
 !!if (SzFreq_Nmeth > 0) echoinput << SzFreq_nobs << " Sizefreq N obs per method" << endl;
   ivector SzFreq_Nbins_seas_g(1,SzFreq_Nmeth*nseas);   //  array dimensioner used only for the SzFreqTrans array
-  ivector SzFreq_Nbins3(1,SzFreq_Nmeth);      // doubles the Nbins if gender==2
+  ivector SzFreq_Nbins3(1,SzFreq_Nmeth);      // values will be doubled if gender==2
   int SzFreqMethod_seas;
   ivector Comp_Err_Sz(1,SzFreq_Nmeth);
   ivector Comp_Err_Sz2(1,SzFreq_Nmeth);
@@ -3196,7 +3196,6 @@
     }
   }
   // clang-format off
-  echoinput<<"here"<<endl;
  END_CALCS
 
  !!echoinput<<"bins "<<SzFreq_Nbins<<endl;
@@ -3266,9 +3265,14 @@
           SzFreq_means(k, z) = SzFreq_means(k, z - 1) + (SzFreq_bins2(k, z) - SzFreq_bins2(k, z - 1));
         }
         if (gender == 2)
+        {
           SzFreq_means(k, z + SzFreq_Nbins(k)) = SzFreq_means(k, z);
+        }
       }
+//      SzFreq_bins2(k, SzFreq_Nbins(k)) = 99999.;
+//      if (gender == 2) SzFreq_bins2(k,SzFreq_Nbins(k) + SzFreq_Nbins(k) + 1) = 99999.;
       echoinput << "Processed_SizeFreqMethod_bins for method: " << k << endl
+                << "bins2: " << SzFreq_bins2(k) << endl
                 << "low: " << SzFreq_bins(k) << endl
                 << "mean: " << SzFreq_means(k) << endl;
     }
diff --git a/SS_write_report.tpl b/SS_write_report.tpl
index 0c08d496..5aadffb6 100644
--- a/SS_write_report.tpl
+++ b/SS_write_report.tpl
@@ -3624,7 +3624,7 @@ FUNCTION void write_bigoutput()
       for (f1 = 1; f1 <= N_pred; f1++)
       {
         f = predator(f1);
-        SS2out << fleetname(f) << "_M2 comsume_Bio consume_Num";
+        SS2out << fleetname(f) << "_M2 consume_Bio consume_Num";
       }
       SS2out << endl;
       for (y = styr - 2; y <= YrMax; y++)