src/Environment/fgwind.cxx

   1 // fgwind.cxx -- routines to create variable winds
   2 //
   3 // Written by Torsten Dreyer
   4 //
   5 // Copyright (C) 2009 Torsten Dreyer - Torsten _at_ t3r*de
   6 //
   7 // This program is free software; you can redistribute it and/or
   8 // modify it under the terms of the GNU General Public License as
   9 // published by the Free Software Foundation; either version 2 of the
  10 // License, or (at your option) any later version.
  11 //
  12 // This program is distributed in the hope that it will be useful, but
  13 // WITHOUT ANY WARRANTY; without even the implied warranty of
  14 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  15 // General Public License for more details.
  16 //
  17 // You should have received a copy of the GNU General Public License
  18 // along with this program; if not, write to the Free Software
  19 // Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.
  20 //
  21 // $Id$
  22 #include "fgwind.hxx"
  23 #include <math.h>
  24 #include <stdio.h>
  25
  26 FGWindModulator::FGWindModulator() :
  27         direction_offset_norm(0.0),
  28         magnitude_factor_norm(1.0)
  29 {
  30 }
  31
  32 FGWindModulator::~FGWindModulator()
  33 {
  34 }
  35
  36 FGBasicWindModulator::FGBasicWindModulator() :
  37         elapsed(0.0),
  38         direction_period(17),
  39         speed_period(1)
  40 {
  41 }
  42
  43 FGBasicWindModulator::~FGBasicWindModulator()
  44 {
  45 }
  46
  47 void FGBasicWindModulator::update( double dt)
  48 {
  49         elapsed += dt;
  50         double t = elapsed/direction_period;
  51
  52         direction_offset_norm = (sin(t)*sin(2*t)+sin(t/3)) / 1.75;
  53
  54         t = elapsed/speed_period;
  55         magnitude_factor_norm = sin(t)* sin(5*direction_offset_norm*direction_offset_norm);;
  56         magnitude_factor_norm = magnitude_factor_norm < 0 ? 0 : magnitude_factor_norm;
  57 }