#! /usr/bin/gawk -f
# Last edited on 2009-12-05 14:57:16 by stolfi

# GAWK library meant to be included in other GAWk programs with "-f"
# Defines parameters for the first WP growth model (2009-11-29),
# used in the IC-UNICAMP Technical Report 09-45.

function set_model_params()
{
  # Phase 1
  bph1 = -14.0;     # Initial day.
  eph1 = 1862;      # Final day.

  Sph1 = 0.00       # Reference day.
  Kph1 = 1050;      # Multiplier.
  Rph1 = +0.00217;  # Growth rate (day^{-1}).
  Aph1 = 0.00;      # Amplitude of seasonal factor.
  Tph1 = 365.25;    # Period of seasonal factor (days).
  Dph1 =  0.50;     # Delay of seasonal factor (cycles).
  
  # Phase 2
  bph2 = 1910;      # Initial day.
  eph2 = 5000;      # Final day.

  Sph2 = 1826;      # Reference day (= 2006-01-01).
  Kph2 = 57670;     # Multiplier.
  Rph2 = -0.000407; # Growth rate (day^{-1}).
  Aph2 = 0.17;      # Amplitude of seasonal factor.
  Tph2 = 182.625;   # Period of seasonal factor (days).
  Dph2 = 0.15;      # Delay of seasonal factor (cycles).
  
  # Transition
  wblur = 365.25;   # Half-width of transition zone (days).
}