bousinessq_operations.h

波浪数值模拟

/*
 * Copyright (c) 2001-2005 Falk Feddersen
 *
 * This program 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 2 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., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
 *
 */

/* ---  bousinessq_operations.h  */

#include "field.h"


/* Calculates field2D operators such as  U_xx, V_yy etc.  */



void b_calc_uxx(field2D *UXX, const field2D *U); /* Returns d^2U/dx^2 on a U grid */

void b_calc_huxx(field2D *HUXX, const field2D *U, const field2D *HBX); /* Returns d^2(hU)/dy^2 on a U grid */

void b_calc_vxy(field2D *VXY, const field2D *V); /* Returns d^2(V)/dxdy on a U grid */

/* Returns d^2(V)/dxdy on a U grid */
/* Note that V = (NX+1,NY),  VXY = (NX,NY)  - U points! */

void b_calc_hvxy(field2D *HVXY, const field2D *V, const field2D *HBY);

/* ----- V POINT ROUTINES -------- */

void b_calc_vyy(field2D *VYY, const field2D *V);   /* Returns d^2V/dx^2 on a V grid */

void b_calc_hvyy(field2D *HVYY, const field2D *V, const field2D *HBY);   /* Returns d^2(hV)/dx^2 on a V grid */

void b_calc_uxy(field2D *UXY, const field2D *U);  /* Returns d^2(U)/dxdy on a V grid */

void b_calc_huxy(field2D *HUXY, const field2D *U, const field2D *HBX);  /* Returns d^2(hU)/dxdy on a V grid */

/* Calculates  UXPVY = delta_x(U) + \delta_y(V) - Can also be used for  \delta_x(hU) etc. 
    -- U is on a U grid, V is on a VGRID,  UXPVY is on an ETA GRID, ie size (NX-1,NY)            */

void b_calc_ux_plus_vy(field2D *UXPVY, const field2D *U, const field2D *V, double idx, double idy);