*********************************** *** Streamer REVISION HISTORY *** *********************************** Version 3.1, 13 February 2013: - Minor code changes for use with the g77 compile on OS X. The old PowerPC architecture is no longer supported on OS X, and g77 is no longer being updated, so the change was necessary. Wanted to use gfortran, but couldn't get the binary to run on different computers. Still investigating. - Rebuilt with gfortran for Windows. See Makefile.win. - Removed Sun binary executable (no Sun computer available for rebuilding) Version 3.0b8, 22 March 2006: - Fixed an error in the calculation of Legendre coefficients for the double Henyey-Greenstein function (optics.f), where g2 was not taken to the k power (thanks Elsa Cattani). - Added makefile for Mac OS X. Version 3.0b4, (beta), 19 February 2002; v3.0b1, 1 August 2001; v3.0a (alpha), 1 March 2001: Significant changes: - Incorporated new version of DISORT (v2.0) - Added eight ice particle "habit" optical property models - Revised the ice cloud optical properties for the spherical particle model, increasing the resolution of particle size at the lower end. The interpolation step was too large for small particles, which became evident in brightness temperatures around 8-9 microns. Thanks to Ralph Holz for discovering this. - Added capability to input surface bidirectional reflectance and added BRDF models - Changed cloud input format to allow for the calculation of water content based on optical and geometrical thickness - Cloud thickness can now default to the layer thickness - Added capability to input cloud optical properties, including a phase function - Switched from the single to the double Henyey-Greenstein function for ice clouds in the shortwave with more than two streams - Increased CO2 from 355 to 360 ppmv to better represent present day levels - Increased the maximum number of clouds to 50 - Changed the structure of the default input file and added a command line option to specify a defaults file - Added an interactive mode of operation and the commands $ASSIGN, $LOAD, $RUN, and $EXIT - Added the option of specifying spectral range with wavelengths, wavenumbers, or band numbers - Added one snow surface type - ***Streamer no longer sets the relative humidity within a cloud to 100%. Code has been added to do this, but it is not done by default. - Added defaults file alternative to editing 'defaults.f' and recompiling - Added GOES-8 imager, ATSR2, SEVIRI, NOAA-15 and -16 AVHRR band weights; recreated MODIS band weights with MODIS spectral response functions. (Note: band 28 in versions 3.0a1 through 3.0b2 was incorrect.) - Made the following available for output via CPRINT/PRINT: a cloud water path output variable and atmospheric temp, humid, etc. profiles - Changed a number of the cloud variable names. Switched around the dimensions of 'averad'. - Added the "-s" command line option to print input data as read. Added additional output information for input read errors. - As a result of the above changes there are a number of changes to the input format and features. See the file 'input_2vs3.txt' for details. - Lots of improvements to the web interface. See 'swirevisions.txt'. Minor changes: - Input unit variables are now replaceable - Added "/n" and "/t" for "/newline" and "/tab" in CPRINT/PRINT - Added warning for negative radiances and reset to zero. - Eliminated backspacing with $replace - Restructured the make files, adding the automatic copying of the appropriate 'd1mach*.f' file for Linux, Sun, SGI, IBM, and HP platforms. - Revised 'getcloud', 'getalb', and 'gethaze, adding rudimentary help and improving consistency -------------------------------------------------------------------------------- Version 2.6.2, 2 Nov 1999: - Made a few minor changes for compatibility with Watcom Fortran (no effect on other platforms). Version 2.6.1, 31 Oct 1999: - Fixed bug in 'readdata.f' that caused DISORT to issue a fatal error when 'tsurf' is negative with band 106 - Added allowance for blank lines in input file - Reworked the test files. Version 2.6, 26 Oct 1999: - Added Web interface and discontinued IDL widget GUI - Added an HTML version of the manuals - Added grass, dry grass, deciduous, and coniferous vegetation albedo types and replaced the generic vegetation type with the average of these - Redid the longwave ice cloud parameterization down to 5 micron Re - Added rudimentory command line help - Eliminated all traces of UVRAD (it wasn't used) - Added -999 tsurf option. Version 2.5.1, 18 Apr 1999: - Made the "dynamic" version of 'd1mach' the default, renaming the original 'd1mach_const.f' - Modified 'getcommand' adding a save command for variables needed in 'parse_replace' (reducing the dependence on compiling "static") - Built Linux executable - Rewrote the makefiles. Version 2.5, 9 Mar 1999: - Changed year input to 2- or 4-digit for Y2K compliance - Increased base CO2 amount to approximately 350 ppm (with a scaling factor) - Rewrote 'parse_replace' and added parallel loop structure - Added automatic conversion of satellite scan angles - Added simple error handling in 'readdata'. Version 2.4.2p, 23 Nov 1998: - Changed 'cldparm.f' so that the shortwave ice cloud parameterization is used for band 106 only if solar radiation is present, otherwise the longwave properties are used. Version 2.4.1p, 17 Sep 1998: - Fixed bug in 'substrmr.f' (2nd call to save_values). Version 2.4p, 8 Sep 1998: - Added mixed-phase cloud optical property parameterization. - Added warning for least accurate conditions in 2-stream shortwave fluxes. Version 2.3p, 15 Jul 1998: - Replaced Ebert and Curry shortwave ice cloud parameterization with Fu and Liou. Version 2.2p, 15 Apr 1998: - Further improved the IDL GUI - Replaced the solar flux data - Added a solar flux adjustment for Earth-Sun distance. Version 2.1.6p, 28 Sep 1997: - Added the revised IDL graphical user interface, v2.0. Version 2.1.5p, 27 Jun 1997: - Fixed a problem with 'bstart'/'bend' for radiances in the subroutine version of streamer. Version 2.1.4p, 18 Jun 1997: - Changed SOLAR and INFRD assignments for radiances in readdata.f. Version 2.1.3p, 10 Jun 1997: - Now requires values on the NCOEF, IMTHD options line - Other minor fixes. Version 2.1.2p, 2 Nov 1996: - Added the capability of $REPLACEing entire arrays with the "{}" syntax. Version 2.1.1p, 15 Oct 1996: - Added bidirectionally reflecting surface option (not completely implemented) - Removed user interface from the streamer subroutine. Version 2.1p, 28 Aug 1996: - Completed Beta testing of cloud overlap. Version 2.1pB, 10 Jun 1996: - Added output for total shortwave and net fluxes - Added heating/cooling rate, increased max line length for $REPLACE - Removed IDL interface (at least temporarily) - Modified cloud overlap to allow for 10 overlap sets and 10 overlapping clouds - Developed a streamer subroutine - Corrected bug in open water albedo routine - Removed command line options capability, removed cloud model definition from options and added it to the input data - Increased the maximum number of polar and azimuthal angles. The input file is not backward compatible. Version 2.04p, 18 Jul 1995: - Fixed bug in water cloud param for long wavelengths and small Re. - Required the albedo and emissivity input lines to have values. Version 2.03p, 20 May 1995: - Fixed bug with cloud base and cloud overlap specification. Version 2.02p, 22 Apr 1995: - Changed the number of profile levels to 100. - Added an optical depth threshold for the use of diffuse albedo under cloud. Version 2.01p, 16 Apr 1995: - Replaced the ice cloud parameterization with one using spherical particles. Version 2.00p, 15 Apr 1994: - Changed interface (added $EXCHANGE, $OPTIONS, $PRINT, $COMMENT) - Added other aerosol models, broadband albedo input - Added "missing" cloud top temperatures, cloud bottom option - Added heating rates to flux output - Added IDL interface - Added DISORT version 2 (preliminary, not complete) - Replaced water cloud parameterization. -------------------------------------------------------------------------------- Version 1.25, 23 Feb 1994: - Added trap for out of range adjusted surface albedo (if rsurf is too high). Version 1.24, 16 Feb 1994: - Fixed bug in strats related to the gas absorber amount. Important! Version 1.23, 23 Dec 1993: - Fixed bug in replace.f, reindexed the array surfalb, added code so that the RH within clouds is 100%. All changes are transparent to the user. Version 1.22, 13 Dec 1993: - Changed the computation and output of satalb & sattb - Created the READRAW utility. Version 1.21, 8 Dec 1993: - Fixed bug with writeusr in streamer.f. Version 1.2, 6 Dec 1993: - Added: user output option, standard profiles, "$replace" input option - MODIS & HIRS weights. Version 1.11, 30 Oct 1993 - More work with band 106 - Lots of code changes and some input options and data changes - Added Rayleigh scattering option. Version 1.11, 8 Oct 1993: - Changed method of specifying standard profiles and surface fractions - Modified the way band 106 is treated (radiances) Version 1.01, 6 Sep 1993: - Added the weighting function option Version 1.7B, 14 Sep 1993: - Added the cloud thickness option for pressure. Version 1.6B, 9 Sep 1993: - Fixed ozone profile conversion - Reversed hzchoice and ozchoice in readdata - Reset the range of cloud parameters - Modified the cloud overlap section - Removed the maxthick option. Version 1.5B, 25 Aug 1993: - Changed use of built-in surface albedo models. Version 1.4B, 23 Aug 1993: - Added the cloud phase option to the cloud model specs. Version 1.3B, 17 Aug 1993: - Changed NSTR option to separate short and longwave - Added cloud overlap - Added the cloud forcing option - Changed input format of profiles. Version 1.2B, 8 Aug 1993: - Redefinition of "broad" bands so that shortwave and longwave results can be given for narrower intervals. Version 1.1B, 15 Jul 1993: - Allow for running without gas abs. Version 1.0B, 1 Jul 1993: - First public Beta test version -------------------------------------------------------------------------------- Version 0.0, 1992: - Development of Streamer begun with a version of Si-Chee Tsay's "strats" program that had been modified by A. Schweiger and J. Key (modified for A. Schweiger's dissertation research at the University of Colorado). ================================================================================ A = alpha test B = beta test p = preliminary In early versions, version numbers in the hundreths place indicate minor bug corrections, e.g., 1.51 - 1.59. In more recent versions minor bug fixes are denoted by numbers following a second period, e.g., 2.1.4.