Astro/sky.cxx

   1 /**************************************************************************
   2  * sky.c -- model sky with an upside down "bowl"
   3  *
   4  * Written by Curtis Olson, started December 1997.
   5  *
   6  * Copyright (C) 1997  Curtis L. Olson  - curt@infoplane.com
   7  *
   8  * This program is free software; you can redistribute it and/or
   9  * modify it under the terms of the GNU General Public License as
  10  * published by the Free Software Foundation; either version 2 of the
  11  * License, or (at your option) any later version.
  12  *
  13  * This program is distributed in the hope that it will be useful, but
  14  * WITHOUT ANY WARRANTY; without even the implied warranty of
  15  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the GNU
  16  * General Public License for more details.
  17  *
  18  * You should have received a copy of the GNU General Public License
  19  * along with this program; if not, write to the Free Software
  20  * Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  21  *
  22  * $Id$
  23  * (Log is kept at end of this file)
  24  **************************************************************************/
  25
  26
  27 #ifdef HAVE_CONFIG_H
  28 #  include <config.h>
  29 #endif
  30
  31 #ifdef HAVE_WINDOWS_H
  32 #  include <windows.h>
  33 #endif
  34
  35 #include <math.h>
  36
  37 #include <GL/glut.h>
  38 #include <XGL/xgl.h>
  39
  40 #include <Aircraft/aircraft.h>
  41 #include <Debug/fg_debug.h>
  42 #include <Flight/flight.h>
  43 #include <Include/fg_constants.h>
  44 #include <Main/views.hxx>
  45 #include <Math/fg_random.h>
  46 #include <Time/event.hxx>
  47 #include <Time/fg_time.hxx>
  48
  49 #include "sky.hxx"
  50
  51
  52 /*
  53 #include <Include/general.h>
  54 */
  55
  56 /* in meters of course */
  57 #define CENTER_ELEV   25000.0
  58
  59 #define INNER_RADIUS  50000.0
  60 #define INNER_ELEV    20000.0
  61
  62 #define MIDDLE_RADIUS 70000.0
  63 #define MIDDLE_ELEV    8000.0
  64
  65 #define OUTER_RADIUS  80000.0
  66 #define OUTER_ELEV        0.0
  67
  68 #define BOTTOM_RADIUS 50000.0
  69 #define BOTTOM_ELEV   -2000.0
  70
  71
  72 static float inner_vertex[12][3];
  73 static float middle_vertex[12][3];
  74 static float outer_vertex[12][3];
  75 static float bottom_vertex[12][3];
  76
  77 static float inner_color[12][4];
  78 static float middle_color[12][4];
  79 static float outer_color[12][4];
  80
  81
  82 /* Calculate the sky structure vertices */
  83 void fgSkyVerticesInit( void ) {
  84     float theta;
  85     int i;
  86
  87     fgPrintf(FG_ASTRO, FG_INFO, "  Generating the sky dome vertices.\n");
  88
  89     for ( i = 0; i < 12; i++ ) {
  90         theta = (i * 30.0) * DEG_TO_RAD;
  91
  92         inner_vertex[i][0] = cos(theta) * INNER_RADIUS;
  93         inner_vertex[i][1] = sin(theta) * INNER_RADIUS;
  94         inner_vertex[i][2] = INNER_ELEV;
  95
  96         /* printf("    %.2f %.2f\n", cos(theta) * INNER_RADIUS,
  97                sin(theta) * INNER_RADIUS); */
  98
  99         middle_vertex[i][0] = cos((double)theta) * MIDDLE_RADIUS;
 100         middle_vertex[i][1] = sin((double)theta) * MIDDLE_RADIUS;
 101         middle_vertex[i][2] = MIDDLE_ELEV;
 102
 103         outer_vertex[i][0] = cos((double)theta) * OUTER_RADIUS;
 104         outer_vertex[i][1] = sin((double)theta) * OUTER_RADIUS;
 105         outer_vertex[i][2] = OUTER_ELEV;
 106
 107         bottom_vertex[i][0] = cos((double)theta) * BOTTOM_RADIUS;
 108         bottom_vertex[i][1] = sin((double)theta) * BOTTOM_RADIUS;
 109         bottom_vertex[i][2] = BOTTOM_ELEV;
 110     }
 111 }
 112
 113
 114 /* (Re)calculate the sky colors at each vertex */
 115 void fgSkyColorsInit( void ) {
 116     fgLIGHT *l;
 117     double sun_angle, diff;
 118     double outer_param[3], outer_amt[3], outer_diff[3];
 119     double middle_param[3], middle_amt[3], middle_diff[3];
 120     int i, j;
 121
 122     l = &cur_light_params;
 123
 124     fgPrintf( FG_ASTRO, FG_INFO,
 125               "  Generating the sky colors for each vertex.\n" );
 126
 127     /* setup for the possibility of sunset effects */
 128     sun_angle = l->sun_angle * RAD_TO_DEG;
 129     // fgPrintf( FG_ASTRO, FG_INFO,
 130     //           "  Sun angle in degrees = %.2f\n", sun_angle);
 131
 132     if ( (sun_angle > 80.0) && (sun_angle < 100.0) ) {
 133         /* 0.0 - 0.4 */
 134         outer_param[0] = (10.0 - fabs(90.0 - sun_angle)) / 25.0;
 135         outer_param[1] = (10.0 - fabs(90.0 - sun_angle)) / 35.0;
 136         outer_param[2] = 0.0;
 137
 138         middle_param[0] = (10.0 - fabs(90.0 - sun_angle)) / 20.0;
 139         middle_param[1] = (10.0 - fabs(90.0 - sun_angle)) / 40.0;
 140         middle_param[2] = 0.0;
 141
 142         outer_diff[0] = outer_param[0] / 6.0;
 143         outer_diff[1] = outer_param[1] / 6.0;
 144         outer_diff[2] = outer_param[2] / 6.0;
 145
 146         middle_diff[0] = middle_param[0] / 6.0;
 147         middle_diff[1] = middle_param[1] / 6.0;
 148         middle_diff[2] = middle_param[2] / 6.0;
 149     } else {
 150         outer_param[0] = outer_param[1] = outer_param[2] = 0.0;
 151         middle_param[0] = middle_param[1] = middle_param[2] = 0.0;
 152
 153         outer_diff[0] = outer_diff[1] = outer_diff[2] = 0.0;
 154         middle_diff[0] = middle_diff[1] = middle_diff[2] = 0.0;
 155     }
 156     /* printf("  outer_red_param = %.2f  outer_red_diff = %.2f\n",
 157            outer_red_param, outer_red_diff); */
 158
 159     /* calculate transition colors between sky and fog */
 160     for ( j = 0; j < 3; j++ ) {
 161         outer_amt[j] = outer_param[j];
 162         middle_amt[j] = middle_param[j];
 163     }
 164
 165     for ( i = 0; i < 6; i++ ) {
 166         for ( j = 0; j < 3; j++ ) {
 167             diff = l->sky_color[j] - l->fog_color[j];
 168
 169             /* printf("sky = %.2f  fog = %.2f  diff = %.2f\n",
 170                    l->sky_color[j], l->fog_color[j], diff); */
 171
 172             inner_color[i][j] = l->sky_color[j] - diff * 0.3;
 173             middle_color[i][j] = l->sky_color[j] - diff * 0.9 + middle_amt[j];
 174             outer_color[i][j] = l->fog_color[j] + outer_amt[j];
 175
 176             if ( inner_color[i][j] > 1.00 ) { inner_color[i][j] = 1.00; }
 177             if ( inner_color[i][j] < 0.10 ) { inner_color[i][j] = 0.10; }
 178             if ( middle_color[i][j] > 1.00 ) { middle_color[i][j] = 1.00; }
 179             if ( middle_color[i][j] < 0.10 ) { middle_color[i][j] = 0.10; }
 180             if ( outer_color[i][j] > 1.00 ) { outer_color[i][j] = 1.00; }
 181             if ( outer_color[i][j] < 0.10 ) { outer_color[i][j] = 0.10; }
 182         }
 183         inner_color[i][3] = middle_color[i][3] = outer_color[i][3] =
 184             l->sky_color[3];
 185
 186         for ( j = 0; j < 3; j++ ) {
 187             outer_amt[j] -= outer_diff[j];
 188             middle_amt[j] -= middle_diff[j];
 189         }
 190
 191         /*
 192         printf("inner_color[%d] = %.2f %.2f %.2f %.2f\n", i, inner_color[i][0],
 193                inner_color[i][1], inner_color[i][2], inner_color[i][3]);
 194         printf("middle_color[%d] = %.2f %.2f %.2f %.2f\n", i,
 195                middle_color[i][0], middle_color[i][1], middle_color[i][2],
 196                middle_color[i][3]);
 197         printf("outer_color[%d] = %.2f %.2f %.2f %.2f\n", i,
 198                outer_color[i][0], outer_color[i][1], outer_color[i][2],
 199                outer_color[i][3]);
 200         */
 201     }
 202
 203     for ( j = 0; j < 3; j++ ) {
 204         outer_amt[j] = 0.0;
 205         middle_amt[j] = 0.0;
 206     }
 207
 208     for ( i = 6; i < 12; i++ ) {
 209
 210         for ( j = 0; j < 3; j++ ) {
 211             diff = l->sky_color[j] - l->fog_color[j];
 212
 213             /* printf("sky = %.2f  fog = %.2f  diff = %.2f\n",
 214                    l->sky_color[j], l->fog_color[j], diff); */
 215
 216             inner_color[i][j] = l->sky_color[j] - diff * 0.3;
 217             middle_color[i][j] = l->sky_color[j] - diff * 0.9 + middle_amt[j];
 218             outer_color[i][j] = l->fog_color[j] + outer_amt[j];
 219
 220             if ( inner_color[i][j] > 1.00 ) { inner_color[i][j] = 1.00; }
 221             if ( inner_color[i][j] < 0.10 ) { inner_color[i][j] = 0.10; }
 222             if ( middle_color[i][j] > 1.00 ) { middle_color[i][j] = 1.00; }
 223             if ( middle_color[i][j] < 0.10 ) { middle_color[i][j] = 0.10; }
 224             if ( outer_color[i][j] > 1.00 ) { outer_color[i][j] = 1.00; }
 225             if ( outer_color[i][j] < 0.15 ) { outer_color[i][j] = 0.15; }
 226         }
 227         inner_color[i][3] = middle_color[i][3] = outer_color[i][3] =
 228             l->sky_color[3];
 229
 230         for ( j = 0; j < 3; j++ ) {
 231             outer_amt[j] += outer_diff[j];
 232             middle_amt[j] += middle_diff[j];
 233         }
 234
 235         /*
 236         printf("inner_color[%d] = %.2f %.2f %.2f %.2f\n", i, inner_color[i][0],
 237                inner_color[i][1], inner_color[i][2], inner_color[i][3]);
 238         printf("middle_color[%d] = %.2f %.2f %.2f %.2f\n", i,
 239                middle_color[i][0], middle_color[i][1], middle_color[i][2],
 240                middle_color[i][3]);
 241         printf("outer_color[%d] = %.2f %.2f %.2f %.2f\n", i,
 242                outer_color[i][0], outer_color[i][1], outer_color[i][2],
 243                outer_color[i][3]);
 244         */
 245     }
 246 }
 247
 248
 249 /* Initialize the sky structure and colors */
 250 void fgSkyInit( void ) {
 251     fgPrintf(FG_ASTRO, FG_INFO, "Initializing the sky\n");
 252
 253     fgSkyVerticesInit();
 254
 255     /* regester fgSkyColorsInit() as an event to be run periodically */
 256     global_events.Register( "fgSkyColorsInit()", fgSkyColorsInit,
 257                             fgEVENT::FG_EVENT_READY, 30000);
 258 }
 259
 260
 261 /* Draw the Sky */
 262 void fgSkyRender( void ) {
 263     fgFLIGHT *f;
 264     fgLIGHT *l;
 265     fgVIEW *v;
 266     float inner_color[4];
 267     float middle_color[4];
 268     float outer_color[4];
 269     double diff;
 270     int i;
 271
 272     f = current_aircraft.flight;
 273     l = &cur_light_params;
 274     v = &current_view;
 275
 276     /* printf("Rendering the sky.\n"); */
 277
 278     // calculate the proper colors
 279     for ( i = 0; i < 3; i++ ) {
 280         diff = l->sky_color[i] - l->adj_fog_color[i];
 281
 282         /* printf("sky = %.2f  fog = %.2f  diff = %.2f\n",
 283            l->sky_color[j], l->adj_fog_color[j], diff); */
 284
 285         inner_color[i] = l->sky_color[i] - diff * 0.3;
 286         middle_color[i] = l->sky_color[i] - diff * 0.9;
 287         outer_color[i] = l->adj_fog_color[i];
 288     }
 289     inner_color[3] = middle_color[3] = outer_color[3] = l->adj_fog_color[3];
 290
 291     xglPushMatrix();
 292
 293     /* Translate to view position */
 294     xglTranslatef( v->cur_zero_elev.x, v->cur_zero_elev.y, v->cur_zero_elev.z );
 295     /* printf("  Translated to %.2f %.2f %.2f\n",
 296            v->cur_zero_elev.x, v->cur_zero_elev.y, v->cur_zero_elev.z ); */
 297
 298     /* Rotate to proper orientation */
 299     /* printf("  lon = %.2f  lat = %.2f\n", FG_Longitude * RAD_TO_DEG,
 300            FG_Latitude * RAD_TO_DEG); */
 301     xglRotatef( FG_Longitude * RAD_TO_DEG, 0.0, 0.0, 1.0 );
 302     xglRotatef( 90.0 - FG_Latitude * RAD_TO_DEG, 0.0, 1.0, 0.0 );
 303     xglRotatef( l->sun_rotation * RAD_TO_DEG, 0.0, 0.0, 1.0 );
 304
 305     /* Draw inner/center section of sky*/
 306     xglBegin( GL_TRIANGLE_FAN );
 307     xglColor4fv(l->sky_color);
 308     xglVertex3f(0.0, 0.0, CENTER_ELEV);
 309     for ( i = 0; i < 12; i++ ) {
 310         xglColor4fv( inner_color );
 311         xglVertex3fv( inner_vertex[i] );
 312     }
 313     xglColor4fv( inner_color );
 314     xglVertex3fv( inner_vertex[0] );
 315     xglEnd();
 316
 317     // Draw the middle ring
 318     xglBegin( GL_TRIANGLE_STRIP );
 319     for ( i = 0; i < 12; i++ ) {
 320         xglColor4fv( middle_color );
 321         /* printf("middle_color[%d] = %.2f %.2f %.2f %.2f\n", i,
 322                middle_color[i][0], middle_color[i][1], middle_color[i][2],
 323                middle_color[i][3]); */
 324         // xglColor4f(1.0, 0.0, 0.0, 1.0);
 325         xglVertex3fv( middle_vertex[i] );
 326         xglColor4fv( inner_color );
 327         /* printf("inner_color[%d] = %.2f %.2f %.2f %.2f\n", i,
 328                inner_color[i][0], inner_color[i][1], inner_color[i][2],
 329                inner_color[i][3]); */
 330         // xglColor4f(0.0, 0.0, 1.0, 1.0);
 331         xglVertex3fv( inner_vertex[i] );
 332     }
 333     xglColor4fv( middle_color );
 334     // xglColor4f(1.0, 0.0, 0.0, 1.0);
 335     xglVertex3fv( middle_vertex[0] );
 336     xglColor4fv( inner_color );
 337     // xglColor4f(0.0, 0.0, 1.0, 1.0);
 338     xglVertex3fv( inner_vertex[0] );
 339     xglEnd();
 340
 341     /* Draw the outer ring */
 342     xglBegin( GL_TRIANGLE_STRIP );
 343     for ( i = 0; i < 12; i++ ) {
 344         xglColor4fv( outer_color );
 345         xglVertex3fv( outer_vertex[i] );
 346         xglColor4fv( middle_color );
 347         xglVertex3fv( middle_vertex[i] );
 348     }
 349     xglColor4fv( outer_color );
 350     xglVertex3fv( outer_vertex[0] );
 351     xglColor4fv( middle_color );
 352     xglVertex3fv( middle_vertex[0] );
 353     xglEnd();
 354
 355     /* Draw the bottom skirt */
 356     xglBegin( GL_TRIANGLE_STRIP );
 357     xglColor4fv( outer_color );
 358     for ( i = 0; i < 12; i++ ) {
 359         xglVertex3fv( bottom_vertex[i] );
 360         xglVertex3fv( outer_vertex[i] );
 361     }
 362     xglVertex3fv( bottom_vertex[0] );
 363     xglVertex3fv( outer_vertex[0] );
 364     xglEnd();
 365
 366     xglPopMatrix();
 367 }
 368
 369
 370 /* $Log$
 371 /* Revision 1.10  1998/08/29 13:07:16  curt
 372 /* Rewrite of event manager thanks to Bernie Bright.
 373 /*
 374  * Revision 1.9  1998/08/22 01:18:59  curt
 375  * Minor tweaks to avoid using unitialized memory.
 376  *
 377  * Revision 1.8  1998/08/12 21:40:44  curt
 378  * Sky now tracks adjusted fog color so it blends well with terrain.
 379  *
 380  * Revision 1.7  1998/07/22 21:39:21  curt
 381  * Lower skirt tracks adjusted fog color, not fog color.
 382  *
 383  * Revision 1.6  1998/05/23 14:07:14  curt
 384  * Use new C++ events class.
 385  *
 386  * Revision 1.5  1998/04/28 01:19:02  curt
 387  * Type-ified fgTIME and fgVIEW
 388  *
 389  * Revision 1.4  1998/04/26 05:10:01  curt
 390  * "struct fgLIGHT" -> "fgLIGHT" because fgLIGHT is typedef'd.
 391  *
 392  * Revision 1.3  1998/04/25 22:06:25  curt
 393  * Edited cvs log messages in source files ... bad bad bad!
 394  *
 395  * Revision 1.2  1998/04/24 00:45:03  curt
 396  * Wrapped "#include <config.h>" in "#ifdef HAVE_CONFIG_H"
 397  * Fixed a bug when generating sky colors.
 398  *
 399  * Revision 1.1  1998/04/22 13:21:32  curt
 400  * C++ - ifing the code a bit.
 401  *
 402  * Revision 1.9  1998/04/03 21:52:50  curt
 403  * Converting to Gnu autoconf system.
 404  *
 405  * Revision 1.8  1998/03/09 22:47:25  curt
 406  * Incorporated Durk's updates.
 407  *
 408  * Revision 1.7  1998/02/19 13:05:49  curt
 409  * Incorporated some HUD tweaks from Michelle America.
 410  * Tweaked the sky's sunset/rise colors.
 411  * Other misc. tweaks.
 412  *
 413  * Revision 1.6  1998/02/07 15:29:32  curt
 414  * Incorporated HUD changes and struct/typedef changes from Charlie Hotchkiss
 415  * <chotchkiss@namg.us.anritsu.com>
 416  *
 417  * Revision 1.5  1998/01/27 00:47:48  curt
 418  * Incorporated Paul Bleisch's <pbleisch@acm.org> new debug message
 419  * system and commandline/config file processing code.
 420  *
 421  * Revision 1.4  1998/01/26 15:54:28  curt
 422  * Added a "skirt" to try to help hide gaps between scenery and sky.  This will
 423  * have to be revisited in the future.
 424  *
 425  * Revision 1.3  1998/01/19 19:26:59  curt
 426  * Merged in make system changes from Bob Kuehne <rpk@sgi.com>
 427  * This should simplify things tremendously.
 428  *
 429  * Revision 1.2  1998/01/19 18:40:17  curt
 430  * Tons of little changes to clean up the code and to remove fatal errors
 431  * when building with the c++ compiler.
 432  *
 433  * Revision 1.1  1998/01/07 03:16:19  curt
 434  * Moved from .../Src/Scenery/ to .../Src/Astro/
 435  *
 436  * Revision 1.11  1997/12/30 22:22:38  curt
 437  * Further integration of event manager.
 438  *
 439  * Revision 1.10  1997/12/30 20:47:53  curt
 440  * Integrated new event manager with subsystem initializations.
 441  *
 442  * Revision 1.9  1997/12/30 13:06:57  curt
 443  * A couple lighting tweaks ...
 444  *
 445  * Revision 1.8  1997/12/23 04:58:38  curt
 446  * Tweaked the sky coloring a bit to build in structures to allow finer rgb
 447  * control.
 448  *
 449  * Revision 1.7  1997/12/22 23:45:48  curt
 450  * First stab at sunset/sunrise sky glow effects.
 451  *
 452  * Revision 1.6  1997/12/22 04:14:34  curt
 453  * Aligned sky with sun so dusk/dawn effects can be correct relative to the sun.
 454  *
 455  * Revision 1.5  1997/12/19 23:34:59  curt
 456  * Lot's of tweaking with sky rendering and lighting.
 457  *
 458  * Revision 1.4  1997/12/19 16:45:02  curt
 459  * Working on scene rendering order and options.
 460  *
 461  * Revision 1.3  1997/12/18 23:32:36  curt
 462  * First stab at sky dome actually starting to look reasonable. :-)
 463  *
 464  * Revision 1.2  1997/12/18 04:07:03  curt
 465  * Worked on properly translating and positioning the sky dome.
 466  *
 467  * Revision 1.1  1997/12/17 23:14:30  curt
 468  * Initial revision.
 469  * Begin work on rendering the sky. (Rather than just using a clear screen.)
 470  *
 471  */