des_coupe_z.C
/*
* Copyright (c) 1999-2001 Eric Gourgoulhon
*
* This file is part of LORENE.
*
* LORENE 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.
*
* LORENE 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 LORENE; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*/
char des_coupe_z_C[] = "$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Graphics/des_coupe_z.C,v 1.5 2014/10/13 08:53:22 j_novak Exp $" ;
/*
* $Id: des_coupe_z.C,v 1.5 2014/10/13 08:53:22 j_novak Exp $
* $Log: des_coupe_z.C,v $
* Revision 1.5 2014/10/13 08:53:22 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.4 2014/10/06 15:16:04 j_novak
* Modified #include directives to use c++ syntax.
*
* Revision 1.3 2008/08/19 06:42:00 j_novak
* Minor modifications to avoid warnings with gcc 4.3. Most of them concern
* cast-type operations, and constant strings that must be defined as const char*
*
* Revision 1.2 2004/03/25 10:29:25 j_novak
* All LORENE's units are now defined in the namespace Unites (in file unites.h).
*
* Revision 1.1.1.1 2001/11/20 15:19:29 e_gourgoulhon
* LORENE
*
* Revision 1.8 2000/02/12 11:19:20 eric
* Ajout de la version avec determination automatique des bornes de la fenetre
* graphique.
*
* Revision 1.7 2000/02/11 18:44:13 eric
* Ajout de l'argument draw_bound.
*
* Revision 1.6 2000/02/11 16:53:36 eric
* Utilisation des coordonnees cartesiennes abolues (X,Y,Z) et non plus
* des coordonnees relatives (x,y,z).
*
* Revision 1.5 1999/12/28 09:02:38 eric
* *** empty log message ***
*
* Revision 1.4 1999/12/27 12:18:51 eric
* Ajout du dessin des frontieres de domains.
*
* Revision 1.3 1999/12/24 13:01:02 eric
* On appelle desormais la routine des_surface_y pour le dessin de la
* surface.
*
* Revision 1.2 1999/12/23 16:16:16 eric
* Ajout du dessin de la surface (argument defsurf).
*
* Revision 1.1 1999/12/09 16:38:18 eric
* Initial revision
*
*
* $Header: /cvsroot/Lorene/C++/Source/Non_class_members/Graphics/des_coupe_z.C,v 1.5 2014/10/13 08:53:22 j_novak Exp $
*
*/
// Header C
#include <cmath>
// Header Lorene
#include "cmp.h"
#include "graphique.h"
#include "param.h"
#include "utilitaires.h"
#include "unites.h"
namespace Lorene {
//******************************************************************************
void des_coupe_z(const Cmp& uu, double z0, int nzdes, const char* title,
const Cmp* defsurf, double zoom, bool draw_bound,
int ncour, int nx, int ny) {
const Map& mp = *(uu.get_mp()) ;
double a1 = mp.val_r(nzdes-1, 1., M_PI/2., 0.) ;
double a2 = mp.val_r(nzdes-1, 1., M_PI/2., M_PI/2.) ;
double a3 = mp.val_r(nzdes-1, 1., M_PI/2., M_PI) ;
double ray = mp.val_r(nzdes-1, 1., 0., 0.) ;
ray = ( a1 > ray ) ? a1 : ray ;
ray = ( a2 > ray ) ? a2 : ray ;
ray = ( a3 > ray ) ? a3 : ray ;
ray *= zoom ;
double x_min = mp.get_ori_x() - ray ;
double x_max = mp.get_ori_x() + ray ;
double y_min = mp.get_ori_y() - ray ;
double y_max = mp.get_ori_y() + ray ;
des_coupe_z(uu, z0, x_min, x_max, y_min, y_max, title, defsurf, draw_bound,
ncour, nx, ny) ;
}
//******************************************************************************
void des_coupe_z(const Cmp& uu, double z0, double x_min, double x_max,
double y_min, double y_max, const char* title, const Cmp* defsurf,
bool draw_bound, int ncour, int nx, int ny) {
using namespace Unites ;
const Map& mp = *(uu.get_mp()) ;
// Plot of isocontours
// -------------------
float* uutab = new float[ny*nx] ;
double hy = (y_max - y_min) / double(ny-1) ;
double hx = (x_max - x_min) / double(nx-1) ;
for (int j=0; j<ny; j++) {
double y = y_min + hy * j ;
for (int i=0; i<nx; i++) {
double x = x_min + hx * i ;
// Computation of (r,theta,phi) :
double r, theta, phi ;
mp.convert_absolute(x, y, z0, r, theta, phi) ;
uutab[nx*j+i] = float(uu.val_point(r, theta, phi)) ;
}
}
float ymin1 = float(y_min / km) ;
float ymax1 = float(y_max / km) ;
float xmin1 = float(x_min / km) ;
float xmax1 = float(x_max / km) ;
const char* nomy = "y [km]" ;
const char* nomx = "x [km]" ;
if (title == 0x0) {
title = "" ;
}
const char* device = 0x0 ;
int newgraph = ( (defsurf != 0x0) || draw_bound ) ? 1 : 3 ;
des_equipot(uutab, nx, ny, xmin1, xmax1, ymin1, ymax1, ncour, nomx, nomy,
title, device, newgraph) ;
delete [] uutab ;
// Plot of the surface
// -------------------
if (defsurf != 0x0) {
assert(defsurf->get_mp() == uu.get_mp()) ;
newgraph = draw_bound ? 0 : 2 ;
des_surface_z(*defsurf, z0, device, newgraph) ;
} // End of the surface drawing
// Plot of the domains outer boundaries
// ------------------------------------
if (draw_bound) {
int ndom = mp.get_mg()->get_nzone() ; // total number of domains
for (int l=0; l<ndom-1; l++) { // loop on the domains (except the
// last one)
newgraph = (l == ndom-2) ? 2 : 0 ;
des_domaine_z(mp, l, z0, device, newgraph) ;
}
}
}
}
Lorene::Cmp
Component of a tensorial field *** DEPRECATED : use class Scalar instead ***.
Definition cmp.h:446
Lorene
Lorene prototypes.
Definition app_hor.h:64
Lorene::y
Coord y
y coordinate centered on the grid
Definition map.h:727
Lorene::Map
Map(const Mg3d &)
Constructor from a multi-domain 3D grid.
Definition map.C:139
Lorene::phi
Coord phi
coordinate centered on the grid
Definition map.h:720
Lorene::x
Coord x
x coordinate centered on the grid
Definition map.h:726
Lorene::r
Coord r
r coordinate centered on the grid
Definition map.h:718