LORENE
des_coupe_y.C
1 /*
2  * Copyright (c) 1999-2001 Eric Gourgoulhon
3  *
4  * This file is part of LORENE.
5  *
6  * LORENE is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License as published by
8  * the Free Software Foundation; either version 2 of the License, or
9  * (at your option) any later version.
10  *
11  * LORENE is distributed in the hope that it will be useful,
12  * but WITHOUT ANY WARRANTY; without even the implied warranty of
13  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  * GNU General Public License for more details.
15  *
16  * You should have received a copy of the GNU General Public License
17  * along with LORENE; if not, write to the Free Software
18  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
19  *
20  */
21 
22 
23 char des_coupe_y_C[] = "$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Graphics/des_coupe_y.C,v 1.6 2014/10/13 08:53:22 j_novak Exp $" ;
24 
25 /*
26  * $Id: des_coupe_y.C,v 1.6 2014/10/13 08:53:22 j_novak Exp $
27  * $Log: des_coupe_y.C,v $
28  * Revision 1.6 2014/10/13 08:53:22 j_novak
29  * Lorene classes and functions now belong to the namespace Lorene.
30  *
31  * Revision 1.5 2014/10/06 15:16:04 j_novak
32  * Modified #include directives to use c++ syntax.
33  *
34  * Revision 1.4 2008/08/19 06:42:00 j_novak
35  * Minor modifications to avoid warnings with gcc 4.3. Most of them concern
36  * cast-type operations, and constant strings that must be defined as const char*
37  *
38  * Revision 1.3 2004/03/25 10:29:24 j_novak
39  * All LORENE's units are now defined in the namespace Unites (in file unites.h).
40  *
41  * Revision 1.2 2002/09/06 15:18:52 e_gourgoulhon
42  * Changement du nom de la variable "hz" en "hza"
43  * pour assurer la compatibilite avec le compilateur xlC_r
44  * sur IBM Regatta (le preprocesseur de ce compilateur remplace
45  * "hz" par "100").
46  *
47  * Revision 1.1.1.1 2001/11/20 15:19:29 e_gourgoulhon
48  * LORENE
49  *
50  * Revision 1.8 2000/02/12 11:18:59 eric
51  * Ajout de la version avec determination automatique des bornes de la fenetre
52  * graphique.
53  *
54  * Revision 1.7 2000/02/11 18:44:06 eric
55  * Ajout de l'argument draw_bound.
56  *
57  * Revision 1.6 2000/02/11 16:53:25 eric
58  * Utilisation des coordonnees cartesiennes abolues (X,Y,Z) et non plus
59  * des coordonnees relatives (x,y,z).
60  *
61  * Revision 1.5 1999/12/28 09:02:30 eric
62  * *** empty log message ***
63  *
64  * Revision 1.4 1999/12/27 12:18:43 eric
65  * Ajout du dessin des frontieres de domains.
66  *
67  * Revision 1.3 1999/12/24 13:00:46 eric
68  * On appelle desormais la routine des_surface_y pour le dessin de la
69  * surface.
70  *
71  * Revision 1.2 1999/12/23 16:16:26 eric
72  * Ajout du dessin de la surface (argument defsurf).
73  *
74  * Revision 1.1 1999/12/09 16:38:10 eric
75  * Initial revision
76  *
77  *
78  * $Header: /cvsroot/Lorene/C++/Source/Non_class_members/Graphics/des_coupe_y.C,v 1.6 2014/10/13 08:53:22 j_novak Exp $
79  *
80  */
81 
82 // Header C
83 #include <cmath>
84 
85 // Header Lorene
86 #include "cmp.h"
87 #include "graphique.h"
88 #include "param.h"
89 #include "utilitaires.h"
90 #include "unites.h"
91 
92 namespace Lorene {
93 //******************************************************************************
94 
95 void des_coupe_y(const Cmp& uu, double y0, int nzdes, const char* title,
96  const Cmp* defsurf, double zoom, bool draw_bound,
97  int ncour, int nx, int nz) {
98 
99  const Map& mp = *(uu.get_mp()) ;
100 
101  double a1 = mp.val_r(nzdes-1, 1., M_PI/2., 0.) ;
102  double a2 = mp.val_r(nzdes-1, 1., M_PI/2., M_PI/2.) ;
103  double a3 = mp.val_r(nzdes-1, 1., M_PI/2., M_PI) ;
104  double ray = mp.val_r(nzdes-1, 1., 0., 0.) ;
105 
106  ray = ( a1 > ray ) ? a1 : ray ;
107  ray = ( a2 > ray ) ? a2 : ray ;
108  ray = ( a3 > ray ) ? a3 : ray ;
109 
110  ray *= zoom ;
111 
112  double x_min = mp.get_ori_x() - ray ;
113  double x_max = mp.get_ori_x() + ray ;
114  double z_min = mp.get_ori_z() - ray ;
115  double z_max = mp.get_ori_z() + ray ;
116 
117  des_coupe_y(uu, y0, x_min, x_max, z_min, z_max, title, defsurf, draw_bound,
118  ncour, nx, nz) ;
119 
120 }
121 
122 //******************************************************************************
123 
124 void des_coupe_y(const Cmp& uu, double y0, double x_min, double x_max,
125  double z_min, double z_max, const char* title, const Cmp* defsurf,
126  bool draw_bound, int ncour, int nx, int nz) {
127 
128  using namespace Unites ;
129 
130  const Map& mp = *(uu.get_mp()) ;
131 
132  // Plot of isocontours
133  // -------------------
134 
135  float* uutab = new float[nx*nz] ;
136 
137  double hx = (x_max - x_min) / double(nx-1) ;
138  double hza = (z_max - z_min) / double(nz-1) ;
139 
140 
141 
142  for (int j=0; j<nz; j++) {
143 
144  double z = z_min + hza * j ;
145 
146  for (int i=0; i<nx; i++) {
147 
148  double x = x_min + hx * i ;
149 
150  // Computation of (r,theta,phi) :
151  double r, theta, phi ;
152  mp.convert_absolute(x, y0, z, r, theta, phi) ;
153 
154  uutab[nx*j+i] = float(uu.val_point(r, theta, phi)) ;
155  }
156  }
157 
158  float xmin1 = float(x_min / km) ;
159  float xmax1 = float(x_max / km) ;
160  float zmin1 = float(z_min / km) ;
161  float zmax1 = float(z_max / km) ;
162 
163  const char* nomx = "x [km]" ;
164  const char* nomz = "z [km]" ;
165 
166  if (title == 0x0) {
167  title = "" ;
168  }
169 
170 
171  const char* device = 0x0 ;
172  int newgraph = ( (defsurf != 0x0) || draw_bound ) ? 1 : 3 ;
173 
174  des_equipot(uutab, nx, nz, xmin1, xmax1, zmin1, zmax1, ncour, nomx, nomz,
175  title, device, newgraph) ;
176 
177  // Plot of the surface
178  // -------------------
179 
180  if (defsurf != 0x0) {
181 
182  assert(defsurf->get_mp() == uu.get_mp()) ;
183 
184  newgraph = draw_bound ? 0 : 2 ;
185 
186  des_surface_y(*defsurf, y0, device, newgraph) ;
187 
188  } // End of the surface drawing
189 
190  delete [] uutab ;
191 
192  // Plot of the domains outer boundaries
193  // ------------------------------------
194 
195  if (draw_bound) {
196 
197  int ndom = mp.get_mg()->get_nzone() ; // total number of domains
198 
199  for (int l=0; l<ndom-1; l++) { // loop on the domains (except the
200  // last one)
201 
202  newgraph = (l == ndom-2) ? 2 : 0 ;
203 
204  des_domaine_y(mp, l, y0, device, newgraph) ;
205  }
206  }
207 
208 
209 }
210 }
void des_surface_y(const Scalar &defsurf, double y0, const char *device=0x0, int newgraph=3, double x_min=-1, double x_max=1, double z_min=-1, double z_max=1, const char *nomx=0x0, const char *nomz=0x0, const char *title=0x0, int nxpage=1, int nypage=1)
Basic routine for drawing a stellar surface in a plane Y=constant.
void des_coupe_y(const Scalar &uu, double y0, int nzdes, const char *title=0x0, const Scalar *defsurf=0x0, double zoom=1.2, bool draw_bound=true, int ncour=15, int nx=100, int nz=100)
Draws isocontour lines of a Scalar in a plane Y=constant.
Lorene prototypes.
Definition: app_hor.h:64
Standard units of space, time and mass.
void des_equipot(float *uutab, int nx, int ny, float xmin, float xmax, float ymin, float ymax, int ncour, const char *nomx, const char *nomy, const char *title, const char *device, int newgraph, int nxpage, int nypage)
Basic routine for drawing isocontours.
Definition: des_equipot.C:74
void des_domaine_y(const Map &mp, int l0, double y0, const char *device=0x0, int newgraph=3, double x_min=-1, double x_max=1, double z_min=-1, double z_max=1, const char *nomx=0x0, const char *nomz=0x0, const char *title=0x0, int nxpage=1, int nypage=1)
Basic routine for drawing the outer boundary of a given domain in a plane Y=constant.