/*
 * EasyWave - A realtime tsunami simulation program with GPU support.
 * Copyright (C) 2014  Andrey Babeyko, Johannes Spazier
 * GFZ German Research Centre for Geosciences (http://www.gfz-potsdam.de)
 *
 * Parts of this program (especially the GPU extension) were developed
 * within the context of the following publicly funded project:
 * - TRIDEC, EU 7th Framework Programme, Grant Agreement 258723
 *   (http://www.tridec-online.eu)
 *
 * This program is free software: you can redistribute it and/or modify
 * it under the terms of the GNU Affero 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 Affero General Public License for more details.
 *
 * You should have received a copy of the GNU Affero General Public License
 * along with this program.  If not, see <http://www.gnu.org/licenses/>.
 */

#ifndef OKADAFAULT_H
#define OKADAFAULT_H

// Position of the fault reference point (to which lon/lat/depth are related)
#define FLT_POS_C  0      // center of the fault, Okada's coordinates: (L/2;W/2)
#define FLT_POS_MT 1      // middle of the top edge: (L/2;W)
#define FLT_POS_BT 2      // beginning of the top edge: (0;W)
#define FLT_POS_BB 3      // beginning of the bottom edge: (0;0)
#define FLT_POS_MB 4      // middle of the bottom edge: (L/2;0)

#define FLT_ERR_DATA      1
#define FLT_ERR_ZTOP      3
#define FLT_ERR_INTERNAL  4
#define FLT_ERR_STRIKE    6

class cOkadaFault
{

protected:

  int checked;
  int adjust_depth;
  double sind;
  double cosd;
  double sins;
  double coss;
  double tand;
  double coslat;
  double zbot;
  double wp;
  double dslip;
  double sslip;

  double mw2m0();

public:

  int refpos;                 // 0- one of POS_XX (see above definitions)
  double mw;
  double slip;
  double lon,lat,depth;
  double strike;
  double dip;
  double rake;
  double length,width;
  double mu;

  cOkadaFault();
  ~cOkadaFault();
  int read( char *faultparam );
  int check();
  double getMw();
  double getM0();
  double getZtop();
  int global2local( double glon, double glat, double& lx, double& ly );
  int local2global( double lx, double ly, double& glon, double& glat );
  int getDeformArea( double& lonmin, double& lonmax, double& latmin, double& latmax );
  int calculate( double lon, double lat, double& uz, double& ulon, double &ulat );
  int calculate( double lon, double lat, double& uz );
};

int okada( double L,double W,double D,double sinD,double cosD,double U1,double U2,double x,double y,int flag_xy, double *Ux,double *Uy,double *Uz );

#endif // OKADAFAULT_H