1 // dome.cxx -- model sky with an upside down "bowl"
3 // Written by Curtis Olson, started December 1997.
4 // SSG-ified by Curtis Olson, February 2000.
6 // Copyright (C) 1997-2000 Curtis L. Olson - curt@flightgear.org
8 // This library is free software; you can redistribute it and/or
9 // modify it under the terms of the GNU Library General Public
10 // License as published by the Free Software Foundation; either
11 // version 2 of the License, or (at your option) any later version.
13 // This library is distributed in the hope that it will be useful,
14 // but WITHOUT ANY WARRANTY; without even the implied warranty of
15 // MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU
16 // Library General Public License for more details.
18 // You should have received a copy of the GNU Library General Public
19 // License along with this library; if not, write to the
20 // Free Software Foundation, Inc., 59 Temple Place - Suite 330,
21 // Boston, MA 02111-1307, USA.
27 # include <simgear_config.h>
37 #include <simgear/xgl/xgl.h>
41 #include <simgear/debug/logstream.hxx>
47 # pragma global_optimizer off
51 // in meters of course
52 #define CENTER_ELEV 25000.0
54 #define UPPER_RADIUS 50000.0
55 #define UPPER_ELEV 20000.0
57 #define MIDDLE_RADIUS 70000.0
58 #define MIDDLE_ELEV 8000.0
60 #define LOWER_RADIUS 80000.0
61 #define LOWER_ELEV 0.0
63 #define BOTTOM_RADIUS 50000.0
64 #define BOTTOM_ELEV -2000.0
67 // Set up dome rendering callbacks
68 static int sgSkyDomePreDraw( ssgEntity *e ) {
69 /* cout << endl << "Dome Pre Draw" << endl << "----------------"
72 ssgLeaf *f = (ssgLeaf *)e;
73 if ( f -> hasState () ) f->getState()->apply() ;
75 glPushAttrib( GL_DEPTH_BUFFER_BIT | GL_FOG_BIT );
76 // cout << "push error = " << glGetError() << endl;
78 glDisable( GL_DEPTH_TEST );
84 static int sgSkyDomePostDraw( ssgEntity *e ) {
85 /* cout << endl << "Dome Post Draw" << endl << "----------------"
89 // cout << "pop error = " << glGetError() << endl;
96 SGSkyDome::SGSkyDome( void ) {
101 SGSkyDome::~SGSkyDome( void ) {
105 // initialize the sky object and connect it into our scene graph
106 ssgBranch * SGSkyDome::build( ) {
113 dome_state = new ssgSimpleState();
114 dome_state->setShadeModel( GL_SMOOTH );
115 dome_state->disable( GL_LIGHTING );
116 dome_state->disable( GL_CULL_FACE );
117 dome_state->disable( GL_TEXTURE_2D );
118 dome_state->enable( GL_COLOR_MATERIAL );
119 dome_state->setColourMaterial( GL_AMBIENT_AND_DIFFUSE );
120 dome_state->setMaterial( GL_EMISSION, 0, 0, 0, 1 );
121 dome_state->setMaterial( GL_SPECULAR, 0, 0, 0, 1 );
122 dome_state->disable( GL_BLEND );
123 dome_state->disable( GL_ALPHA_TEST );
126 center_disk_vl = new ssgVertexArray( 14 );
127 center_disk_cl = new ssgColourArray( 14 );
129 upper_ring_vl = new ssgVertexArray( 26 );
130 upper_ring_cl = new ssgColourArray( 26 );
132 middle_ring_vl = new ssgVertexArray( 26 );
133 middle_ring_cl = new ssgColourArray( 26 );
135 lower_ring_vl = new ssgVertexArray( 26 );
136 lower_ring_cl = new ssgColourArray( 26 );
138 // initially seed to all blue
139 sgSetVec4( color, 0.0, 0.0, 1.0, 1.0 );
141 // generate the raw vertex data
142 sgVec3 center_vertex;
143 sgVec3 upper_vertex[12];
144 sgVec3 middle_vertex[12];
145 sgVec3 lower_vertex[12];
146 sgVec3 bottom_vertex[12];
148 sgSetVec3( center_vertex, 0.0, 0.0, CENTER_ELEV );
150 for ( i = 0; i < 12; i++ ) {
151 theta = (i * 30.0) * SGD_DEGREES_TO_RADIANS;
153 sgSetVec3( upper_vertex[i],
154 cos(theta) * UPPER_RADIUS,
155 sin(theta) * UPPER_RADIUS,
158 sgSetVec3( middle_vertex[i],
159 cos((double)theta) * MIDDLE_RADIUS,
160 sin((double)theta) * MIDDLE_RADIUS,
163 sgSetVec3( lower_vertex[i],
164 cos((double)theta) * LOWER_RADIUS,
165 sin((double)theta) * LOWER_RADIUS,
168 sgSetVec3( bottom_vertex[i],
169 cos((double)theta) * BOTTOM_RADIUS,
170 sin((double)theta) * BOTTOM_RADIUS,
174 // generate the center disk vertex/color arrays
175 center_disk_vl->add( center_vertex );
176 center_disk_cl->add( color );
177 for ( i = 11; i >= 0; i-- ) {
178 center_disk_vl->add( upper_vertex[i] );
179 center_disk_cl->add( color );
181 center_disk_vl->add( upper_vertex[11] );
182 center_disk_cl->add( color );
184 // generate the upper ring
185 for ( i = 0; i < 12; i++ ) {
186 upper_ring_vl->add( middle_vertex[i] );
187 upper_ring_cl->add( color );
189 upper_ring_vl->add( upper_vertex[i] );
190 upper_ring_cl->add( color );
192 upper_ring_vl->add( middle_vertex[0] );
193 upper_ring_cl->add( color );
195 upper_ring_vl->add( upper_vertex[0] );
196 upper_ring_cl->add( color );
198 // generate middle ring
199 for ( i = 0; i < 12; i++ ) {
200 middle_ring_vl->add( lower_vertex[i] );
201 middle_ring_cl->add( color );
203 middle_ring_vl->add( middle_vertex[i] );
204 middle_ring_cl->add( color );
206 middle_ring_vl->add( lower_vertex[0] );
207 middle_ring_cl->add( color );
209 middle_ring_vl->add( middle_vertex[0] );
210 middle_ring_cl->add( color );
212 // generate lower ring
213 for ( i = 0; i < 12; i++ ) {
214 lower_ring_vl->add( bottom_vertex[i] );
215 lower_ring_cl->add( color );
217 lower_ring_vl->add( lower_vertex[i] );
218 lower_ring_cl->add( color );
220 lower_ring_vl->add( bottom_vertex[0] );
221 lower_ring_cl->add( color );
223 lower_ring_vl->add( lower_vertex[0] );
224 lower_ring_cl->add( color );
226 // force a repaint of the sky colors with ugly defaults
228 sgSetVec4( fog_color, 1.0, 1.0, 1.0, 1.0 );
229 repaint( color, fog_color, 0.0, 5000.0 );
231 // build the ssg scene graph sub tree for the sky and connected
232 // into the provide scene graph branch
233 ssgVtxTable *center_disk, *upper_ring, *middle_ring, *lower_ring;
235 center_disk = new ssgVtxTable( GL_TRIANGLE_FAN,
236 center_disk_vl, NULL, NULL, center_disk_cl );
238 upper_ring = new ssgVtxTable( GL_TRIANGLE_STRIP,
239 upper_ring_vl, NULL, NULL, upper_ring_cl );
241 middle_ring = new ssgVtxTable( GL_TRIANGLE_STRIP,
242 middle_ring_vl, NULL, NULL, middle_ring_cl );
244 lower_ring = new ssgVtxTable( GL_TRIANGLE_STRIP,
245 lower_ring_vl, NULL, NULL, lower_ring_cl );
247 center_disk->setState( dome_state );
248 upper_ring->setState( dome_state );
249 middle_ring->setState( dome_state );
250 lower_ring->setState( dome_state );
252 dome_transform = new ssgTransform;
253 dome_transform->addKid( center_disk );
254 dome_transform->addKid( upper_ring );
255 dome_transform->addKid( middle_ring );
256 dome_transform->addKid( lower_ring );
258 // not entirely satisfying. We are depending here that the first
259 // thing we add to a parent is the first drawn
260 center_disk->setCallback( SSG_CALLBACK_PREDRAW, sgSkyDomePreDraw );
261 center_disk->setCallback( SSG_CALLBACK_POSTDRAW, sgSkyDomePostDraw );
263 upper_ring->setCallback( SSG_CALLBACK_PREDRAW, sgSkyDomePreDraw );
264 upper_ring->setCallback( SSG_CALLBACK_POSTDRAW, sgSkyDomePostDraw );
266 middle_ring->setCallback( SSG_CALLBACK_PREDRAW, sgSkyDomePreDraw );
267 middle_ring->setCallback( SSG_CALLBACK_POSTDRAW, sgSkyDomePostDraw );
269 lower_ring->setCallback( SSG_CALLBACK_PREDRAW, sgSkyDomePreDraw );
270 lower_ring->setCallback( SSG_CALLBACK_POSTDRAW, sgSkyDomePostDraw );
272 return dome_transform;
276 // repaint the sky colors based on current value of sun_angle, sky,
277 // and fog colors. This updates the color arrays for ssgVtxTable.
278 // sun angle in degrees relative to verticle
279 // 0 degrees = high noon
280 // 90 degrees = sun rise/set
281 // 180 degrees = darkest midnight
282 bool SGSkyDome::repaint( sgVec4 sky_color, sgVec4 fog_color, double sun_angle,
286 sgVec3 outer_param, outer_amt, outer_diff;
287 sgVec3 middle_param, middle_amt, middle_diff;
290 // Check for sunrise/sunset condition
291 if ( (sun_angle > 80.0) && (sun_angle < 100.0) ) {
293 sgSetVec3( outer_param,
294 (10.0 - fabs(90.0 - sun_angle)) / 20.0,
295 (10.0 - fabs(90.0 - sun_angle)) / 40.0,
296 -(10.0 - fabs(90.0 - sun_angle)) / 30.0 );
298 sgSetVec3( middle_param,
299 (10.0 - fabs(90.0 - sun_angle)) / 40.0,
300 (10.0 - fabs(90.0 - sun_angle)) / 80.0,
303 sgScaleVec3( outer_diff, outer_param, 1.0 / 6.0 );
305 sgScaleVec3( middle_diff, middle_param, 1.0 / 6.0 );
307 sgSetVec3( outer_param, 0.0, 0.0, 0.0 );
308 sgSetVec3( middle_param, 0.0, 0.0, 0.0 );
310 sgSetVec3( outer_diff, 0.0, 0.0, 0.0 );
311 sgSetVec3( middle_diff, 0.0, 0.0, 0.0 );
313 // printf(" outer_red_param = %.2f outer_red_diff = %.2f\n",
314 // outer_red_param, outer_red_diff);
316 // calculate transition colors between sky and fog
317 sgCopyVec3( outer_amt, outer_param );
318 sgCopyVec3( middle_amt, middle_param );
321 // First, recalulate the basic colors
325 sgVec4 upper_color[12];
326 sgVec4 middle_color[12];
327 sgVec4 lower_color[12];
328 sgVec4 bottom_color[12];
332 if ( vis < 3000.0 ) {
333 vis_factor = (vis - 1000.0) / 2000.0;
334 if ( vis_factor < 0.0 ) {
341 for ( j = 0; j < 3; j++ ) {
342 diff = sky_color[j] - fog_color[j];
343 center_color[j] = sky_color[j] - diff * ( 1.0 - vis_factor );
345 center_color[3] = 1.0;
347 for ( i = 0; i < 6; i++ ) {
348 for ( j = 0; j < 3; j++ ) {
349 diff = sky_color[j] - fog_color[j];
351 // printf("sky = %.2f fog = %.2f diff = %.2f\n",
352 // l->sky_color[j], l->fog_color[j], diff);
354 upper_color[i][j] = sky_color[j] - diff * ( 1.0 - vis_factor * 0.7);
355 middle_color[i][j] = sky_color[j] - diff * ( 1.0 - vis_factor * 0.1)
357 lower_color[i][j] = fog_color[j] + outer_amt[j];
359 if ( upper_color[i][j] > 1.0 ) { upper_color[i][j] = 1.0; }
360 if ( upper_color[i][j] < 0.0 ) { upper_color[i][j] = 0.0; }
361 if ( middle_color[i][j] > 1.0 ) { middle_color[i][j] = 1.0; }
362 if ( middle_color[i][j] < 0.0 ) { middle_color[i][j] = 0.0; }
363 if ( lower_color[i][j] > 1.0 ) { lower_color[i][j] = 1.0; }
364 if ( lower_color[i][j] < 0.0 ) { lower_color[i][j] = 0.0; }
366 upper_color[i][3] = middle_color[i][3] = lower_color[i][3] = 1.0;
368 for ( j = 0; j < 3; j++ ) {
369 outer_amt[j] -= outer_diff[j];
370 middle_amt[j] -= middle_diff[j];
374 printf("upper_color[%d] = %.2f %.2f %.2f %.2f\n", i, upper_color[i][0],
375 upper_color[i][1], upper_color[i][2], upper_color[i][3]);
376 printf("middle_color[%d] = %.2f %.2f %.2f %.2f\n", i,
377 middle_color[i][0], middle_color[i][1], middle_color[i][2],
379 printf("lower_color[%d] = %.2f %.2f %.2f %.2f\n", i,
380 lower_color[i][0], lower_color[i][1], lower_color[i][2],
385 sgSetVec3( outer_amt, 0.0, 0.0, 0.0 );
386 sgSetVec3( middle_amt, 0.0, 0.0, 0.0 );
388 for ( i = 6; i < 12; i++ ) {
389 for ( j = 0; j < 3; j++ ) {
390 diff = sky_color[j] - fog_color[j];
392 // printf("sky = %.2f fog = %.2f diff = %.2f\n",
393 // sky_color[j], fog_color[j], diff);
395 upper_color[i][j] = sky_color[j] - diff * ( 1.0 - vis_factor * 0.7);
396 middle_color[i][j] = sky_color[j] - diff * ( 1.0 - vis_factor * 0.1)
398 lower_color[i][j] = fog_color[j] + outer_amt[j];
400 if ( upper_color[i][j] > 1.0 ) { upper_color[i][j] = 1.0; }
401 if ( upper_color[i][j] < 0.0 ) { upper_color[i][j] = 0.0; }
402 if ( middle_color[i][j] > 1.0 ) { middle_color[i][j] = 1.0; }
403 if ( middle_color[i][j] < 0.0 ) { middle_color[i][j] = 0.0; }
404 if ( lower_color[i][j] > 1.0 ) { lower_color[i][j] = 1.0; }
405 if ( lower_color[i][j] < 0.0 ) { lower_color[i][j] = 0.0; }
407 upper_color[i][3] = middle_color[i][3] = lower_color[i][3] = 1.0;
409 for ( j = 0; j < 3; j++ ) {
410 outer_amt[j] += outer_diff[j];
411 middle_amt[j] += middle_diff[j];
415 printf("upper_color[%d] = %.2f %.2f %.2f %.2f\n", i, upper_color[i][0],
416 upper_color[i][1], upper_color[i][2], upper_color[i][3]);
417 printf("middle_color[%d] = %.2f %.2f %.2f %.2f\n", i,
418 middle_color[i][0], middle_color[i][1], middle_color[i][2],
420 printf("lower_color[%d] = %.2f %.2f %.2f %.2f\n", i,
421 lower_color[i][0], lower_color[i][1], lower_color[i][2],
426 for ( i = 0; i < 12; i++ ) {
427 sgCopyVec4( bottom_color[i], fog_color );
431 // Second, assign the basic colors to the object color arrays
437 // update the center disk color arrays
439 slot = center_disk_cl->get( counter++ );
441 // sgSetVec4( red, 1.0, 0.0, 0.0, 1.0 );
442 sgCopyVec4( slot, center_color );
443 for ( i = 11; i >= 0; i-- ) {
444 slot = center_disk_cl->get( counter++ );
445 sgCopyVec4( slot, upper_color[i] );
447 slot = center_disk_cl->get( counter++ );
448 sgCopyVec4( slot, upper_color[11] );
450 // generate the upper ring
452 for ( i = 0; i < 12; i++ ) {
453 slot = upper_ring_cl->get( counter++ );
454 sgCopyVec4( slot, middle_color[i] );
456 slot = upper_ring_cl->get( counter++ );
457 sgCopyVec4( slot, upper_color[i] );
459 slot = upper_ring_cl->get( counter++ );
460 sgCopyVec4( slot, middle_color[0] );
462 slot = upper_ring_cl->get( counter++ );
463 sgCopyVec4( slot, upper_color[0] );
465 // generate middle ring
467 for ( i = 0; i < 12; i++ ) {
468 slot = middle_ring_cl->get( counter++ );
469 sgCopyVec4( slot, lower_color[i] );
471 slot = middle_ring_cl->get( counter++ );
472 sgCopyVec4( slot, middle_color[i] );
474 slot = middle_ring_cl->get( counter++ );
475 sgCopyVec4( slot, lower_color[0] );
477 slot = middle_ring_cl->get( counter++ );
478 sgCopyVec4( slot, middle_color[0] );
480 // generate lower ring
482 for ( i = 0; i < 12; i++ ) {
483 slot = lower_ring_cl->get( counter++ );
484 sgCopyVec4( slot, bottom_color[i] );
486 slot = lower_ring_cl->get( counter++ );
487 sgCopyVec4( slot, lower_color[i] );
489 slot = lower_ring_cl->get( counter++ );
490 sgCopyVec4( slot, bottom_color[0] );
492 slot = lower_ring_cl->get( counter++ );
493 sgCopyVec4( slot, lower_color[0] );
499 // reposition the sky at the specified origin and orientation
500 // lon specifies a rotation about the Z axis
501 // lat specifies a rotation about the new Y axis
502 // spin specifies a rotation about the new Z axis (and orients the
503 // sunrise/set effects
504 bool SGSkyDome::reposition( sgVec3 p, double lon, double lat, double spin ) {
505 sgMat4 T, LON, LAT, SPIN;
508 // Translate to view position
509 // Point3D zero_elev = current_view.get_cur_zero_elev();
510 // xglTranslatef( zero_elev.x(), zero_elev.y(), zero_elev.z() );
511 sgMakeTransMat4( T, p );
513 // printf(" Translated to %.2f %.2f %.2f\n",
514 // zero_elev.x, zero_elev.y, zero_elev.z );
516 // Rotate to proper orientation
517 // printf(" lon = %.2f lat = %.2f\n",
518 // lon * SGD_RADIANS_TO_DEGREES,
519 // lat * SGD_RADIANS_TO_DEGREES);
520 // xglRotatef( lon * SGD_RADIANS_TO_DEGREES, 0.0, 0.0, 1.0 );
521 sgSetVec3( axis, 0.0, 0.0, 1.0 );
522 sgMakeRotMat4( LON, lon * SGD_RADIANS_TO_DEGREES, axis );
524 // xglRotatef( 90.0 - f->get_Latitude() * SGD_RADIANS_TO_DEGREES,
526 sgSetVec3( axis, 0.0, 1.0, 0.0 );
527 sgMakeRotMat4( LAT, 90.0 - lat * SGD_RADIANS_TO_DEGREES, axis );
529 // xglRotatef( l->sun_rotation * SGD_RADIANS_TO_DEGREES, 0.0, 0.0, 1.0 );
530 sgSetVec3( axis, 0.0, 0.0, 1.0 );
531 sgMakeRotMat4( SPIN, spin * SGD_RADIANS_TO_DEGREES, axis );
535 sgCopyMat4( TRANSFORM, T );
536 sgPreMultMat4( TRANSFORM, LON );
537 sgPreMultMat4( TRANSFORM, LAT );
538 sgPreMultMat4( TRANSFORM, SPIN );
541 sgSetCoord( &skypos, TRANSFORM );
543 dome_transform->setTransform( &skypos );