# Last edited on 2009-11-28 20:43:19 by stolfi

  split("", Fph1); # Month-specific factor for phase 1
  for (mo = 1; mo <= 12; mo++) { Fph1[mo] = 1.000; }
  split("", Fph2); # Month-specific factor for phase 2
  for (mo = 1; mo <= 12; mo++) { Fph2[mo] = 1.000; }


function model_peak(mr,mpk,inpk,   dmr)
{
  # Models a peak centered in month {mpk} with values {inpk[mr-mpk]}
  dmr = mr - mpk;
  if (dmr in inpk)
    { return inpk[dmr]; }
  else
    { return 0; }
}