calc_external_forcing.cpp

Go to the documentation of this file.

 
/* ==============================================================================================================================
 
   This file is part of CoastalME, the Coastal Modelling Environment.
 
   CoastalME is free software; you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation; either version 3 of the License, or (at your option) any later version.
 
   This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
 
   You should have received a copy of the GNU General Public License along with this program; if not, write to the Free Software Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
 
==============================================================================================================================*/
#include <iostream>
using std::cerr;
using std::endl;
using std::cout;
 
#include "cme.h"
#include "simulation.h"
 
//===============================================================================================================================
//===============================================================================================================================
int CSimulation::nCalcExternalForcing(void)
{
   // Increment SWL (note that increment may be zero)
   m_dAccumulatedSeaLevelChange += m_dDeltaSWLPerTimestep;
 
   // This-iteration mean SWL includes only long-term SWL change
   m_dThisIterMeanSWL = m_dInitialMeanSWL + m_dAccumulatedSeaLevelChange;
 
   int nSize = static_cast<int>(m_VdTideData.size());
 
   if (nSize == 0)
   {
      // No tide data
      m_dThisIterSWL = m_dThisIterMeanSWL;
   }
 
   else
   {
      // We have tide data
      static int snTideDataCount = 0;
 
      // Wrap the tide data, i.e. start again with the first record if we do not have enough
      if (snTideDataCount > nSize - 1)
         snTideDataCount = 0;
 
      // This-iteration SWL includes both tidal change and long-term SWL change
      m_dThisIterSWL = m_dInitialMeanSWL + m_VdTideData[snTideDataCount];
 
      snTideDataCount++;
   }
 
   // Update min and max still water levels
   m_dMaxSWL = tMax(m_dThisIterSWL, m_dMaxSWL);
   m_dMinSWL = tMin(m_dThisIterSWL, m_dMinSWL);
 
   // Update the wave height, orientation and period for this time step and start again with the first record if we do not have enough
   if (m_bHaveWaveStationData)
   {
      static int snWaveStationDataCount = 0;
 
      if (snWaveStationDataCount > m_nDeepWaterWaveDataNumTimeSteps - 1)
      {
         // Wrap the tide data, i.e. start again with the first record if we do not have enough
         snWaveStationDataCount = 0;
      }
 
      // Do we have just a single wave station?
      if (m_bSingleDeepWaterWaveValues)
      {
         // Yes, just a single wave station
         m_dAllCellsDeepWaterWaveHeight = m_VdTSDeepWaterWaveStationHeight[snWaveStationDataCount];
         m_dAllCellsDeepWaterWaveAngle = m_VdTSDeepWaterWaveStationAngle[snWaveStationDataCount];
         m_dAllCellsDeepWaterWavePeriod = m_VdTSDeepWaterWaveStationPeriod[snWaveStationDataCount];
      }
 
      else
      {
         // More than one wave station, so update this time step's deep water wave values for use in the nInterpolateAllDeepWaterWaveValues() routine. Note that the order on the vector is determined by the points ID i.e. to ensure that stations match with time series
         int nNumberDeepWaterWaveStations = static_cast<int>(m_VnDeepWaterWaveStationID.size());
         int nTot = nNumberDeepWaterWaveStations * snWaveStationDataCount;
 
         for (int j = 0; j < nNumberDeepWaterWaveStations; j++)
         {
            m_VdThisIterDeepWaterWaveStationHeight[j] = m_VdTSDeepWaterWaveStationHeight[(m_VnDeepWaterWaveStationID[j] - 1) + nTot];
            m_VdThisIterDeepWaterWaveStationAngle[j] = m_VdTSDeepWaterWaveStationAngle[(m_VnDeepWaterWaveStationID[j] - 1) + nTot];
            m_VdThisIterDeepWaterWaveStationPeriod[j] = m_VdTSDeepWaterWaveStationPeriod[(m_VnDeepWaterWaveStationID[j] - 1) + nTot];
         }
      }
 
      snWaveStationDataCount++;
   }
 
   return RTN_OK;
}