October 20, 1999 Readme for ERBE_S7_NAT Data Note to user: Questions concerning these data or read software should contact the following: EOSDIS Langley DAAC Science, Users and Data Services Office NASA Langley Research Center Mail Stop 157D 2 South Wright Street Hampton, Virginia 23681-2199 U.S.A. Email: support-asdc@earthdata.nasa.gov Phone: (757)864-8656 FAX: (757)864-8807 URL: http://eosweb.larc.nasa.gov The sample read software has been provided by the ERBE Data Management Team in 1994. To maintain all of the information for this sample read software package, all files have bin "zipped" together using the Info-Zip software. To read the files, first unzip them with Info-Zip software. Information about downloading and using this freeware is available at the following URL: http://www.cdrom.com/pub/infozip/ Once you have "unzipped" the sample read software file package, you will find the following files: makefile : file used to compile software reads7 : script to execute code s7.f.cray : Cray Fortran V version of code s7.f.sun : Sun Fortran V version of code In order to follow the instructions provided by the ERBE Team, you must perform the following: * Decide which version of the Fortran code you plan to use on your designated platform, either s7.f.cray or 7.f.sun. * Once you have decided which version of code to use, remove the ".cray" or ".sun" extension from the name. Rename the Fortran V code to s7.f. * The Sun Fortran 77 code is able to compile and execute the s7.f.sun version of the code. * Rename your data files. The underscores must be removed from the names in order for the code to recognize the data files. Example: cp s7_9301_2 s793012 ------------------------------ Original Document ---------------------------------- READING THE S-7 DATA FROM TAPE MEDIA To read an S-7 file from tape, you will need to use the UNIX utility dd. For example, if you have a Sun workstation, you should use: dd bs=10240 if=dev of=fn where: bs = blocksize if = input file name (or 8mm drive name) of = output file name The average length of most of these files is 28 to 29 Mbytes. The dd command above will read one entire file from the tape, including the 30-byte ERBE header, and places it in the file designated by "of". Note that the above parameters may vary depending upon the type of workstation you are using. If you have any problems with this procedure, please contact User Services at support-asdc@earthdata.nasa.gov ------------------------------------------------------------------------------- USING THE S-7 READ SOFTWARE This software is meant to be only a starting point for the user. It is expected that the user will need to modify the program in order to meet the user's specific needs. In order to get a data dump, the following files must be present: Files: s7.f - read software written in FORTRAN V makefile - script that is run to compile the software and produce the executable s7.exe reads7 - a UNIX script which requests certain information from the user and then executes the read software passing on the user-entered parameters to the program The S7 read software is a FORTRAN program. To compile the program, execute the script, makefile, by entering: make This procedure will compile the program and produce an executable called s7.exe. To execute the program and get a data dump, execute the script, reads7, by entering: reads7 The following is a sample session using data from the ERBS satellite for January 13, 1993. The information after the "=>" was entered by the user. The rest of the output was generated by the script and the program s7.exe. In this particular session, records 1 through 5 were requested for day 13. Note that there can be a maximum of 5400 records for each day. -------------------------------------------------------------------- Enter file name for S-7 =>s793012 Enter day of month (or 0 for entire month) =>13 Enter dump range: Initial =>1 Final =>5 *********************************************************** S-7 Dump Request: S-7 file name is : s793012 Day : 13 Dump range: Initial : 1 Final : 5 ************************************************************ Do you wish to continue with execution of this job?(enter y/n)=>y ------------------------------------------------------------------- If you enter n, the script will exit. If you enter y, the script will execute the program s7.exe and produce file fort.9. The output for the indicated data file above (ERBS data, January 1993, day 13, records 1 through 5) is attached below. Note: Some versions of FORTRAN appear to have difficulties with underscores being used in the file name. If you used underscores (for example, s7_9003_2) and an error occurs, try renaming the file without the underscores (example, s790032). ------------------------------------------------------------------------------ CONTENTS OF SAMPLE OUTPUT IN FORT.9 product specifications: temporal span : 93/01/01 0000 - 93/01/31 2359 spacecraft : ERBS product : S-7 vsn / file : s793012 physical header contents: subsystem indicator : 5 product code : 9 spacecraft indicator : 2 whole julian date - high : 244 whole julian date - low : 8988 fractional julian date : 5000 processed version counter : 1 year processed : 93 month processed : 7 day processed : 7 hour processed : 10 minute processed : 41 second processed : 38 record information: index content 1 Julian date 2 Julian time 3 Earth - Sun distance 4 s/c position, x 6 s/c position, y 8 s/c position, z 10 s/c velocity (xdot) 12 s/c velocity (ydot) 14 s/c velocity (zdot) 16 s/c nadir position, colatitude 18 s/c nadir position, longitude 20 sun position, colatitude 21 sun position, longitude 22 orbit number key 23 wfov, radiometric, total 27 wfov, radiometric, sw 31 mfov, radiometric, total 35 mfov, radiometric, sw 39 wfov, unfiltered, sw 43 wfov, unfiltered, lw 47 mfov, unfiltered, sw 51 mfov, unfiltered, lw 55 wfov, toa est., nf, sw 56 wfov, toa est., nf, lw 57 mfov, toa est., nf, sw 58 mfov, toa est., nf, lw 59 wfov, toa est., sf, sw 60 wfov, toa est., sf, lw 61 mfov, toa est., sf, sw 62 mfov, toa est., sf, lw 63 nonscanner, fov, colatitude 67 nonscanner, fov, longitude 71 nonscanner, op. mode, flag word 73 flag, nonscanner, toa est. 74 orbit number scale factor 75 spare scale factors offsets record xmwsf ( 1) = 1.000000000 xmwoff( 1) = 0.000000000 xmwsf ( 2) = 1000000000.000000000 xmwoff( 2) = 0.000000000 xmwsf ( 3) = 1000000000.000000000 xmwoff( 3) = 0.000000000 xmwsf ( 4) = 1.000000000 xmwoff( 4) = 0.000000000 xmwsf ( 5) = 1.000000000 xmwoff( 5) = 0.000000000 xmwsf ( 6) = 1.000000000 xmwoff( 6) = 0.000000000 xmwsf ( 7) = 1.000000000 xmwoff( 7) = 0.000000000 xmwsf ( 8) = 1.000000000 xmwoff( 8) = 0.000000000 xmwsf ( 9) = 1.000000000 xmwoff( 9) = 0.000000000 xmwsf (10) = 1.000000000 xmwoff(10) = 0.000000000 xmwsf (11) = 1.000000000 xmwoff(11) = 0.000000000 xmwsf (12) = 1.000000000 xmwoff(12) = 0.000000000 xmwsf (13) = 1.000000000 xmwoff(13) = 0.000000000 xmwsf (14) = 1.000000000 xmwoff(14) = 0.000000000 xmwsf (15) = 1.000000000 xmwoff(15) = 0.000000000 xmwsf (16) = 100.000000000 xmwoff(16) = 0.000000000 xmwsf (17) = 100.000000000 xmwoff(17) = 0.000000000 xmwsf (18) = 100.000000000 xmwoff(18) = -180.000000000 xmwsf (19) = 100.000000000 xmwoff(19) = -180.000000000 xmwsf (20) = 100.000000000 xmwoff(20) = 0.000000000 xmwsf (21) = 100.000000000 xmwoff(21) = -180.000000000 xmwsf (22) = 1.000000000 xmwoff(22) = 0.000000000 xmwsf (23) = 10.000000000 xmwoff(23) = 0.000000000 xmwsf (24) = 10.000000000 xmwoff(24) = 0.000000000 xmwsf (25) = 10.000000000 xmwoff(25) = 0.000000000 xmwsf (26) = 10.000000000 xmwoff(26) = 0.000000000 xmwsf (27) = 10.000000000 xmwoff(27) = 0.000000000 xmwsf (28) = 10.000000000 xmwoff(28) = 0.000000000 xmwsf (29) = 10.000000000 xmwoff(29) = 0.000000000 xmwsf (30) = 10.000000000 xmwoff(30) = 0.000000000 xmwsf (31) = 10.000000000 xmwoff(31) = 0.000000000 xmwsf (32) = 10.000000000 xmwoff(32) = 0.000000000 xmwsf (33) = 10.000000000 xmwoff(33) = 0.000000000 xmwsf (34) = 10.000000000 xmwoff(34) = 0.000000000 xmwsf (35) = 10.000000000 xmwoff(35) = 0.000000000 xmwsf (36) = 10.000000000 xmwoff(36) = 0.000000000 xmwsf (37) = 10.000000000 xmwoff(37) = 0.000000000 xmwsf (38) = 10.000000000 xmwoff(38) = 0.000000000 xmwsf (39) = 10.000000000 xmwoff(39) = 0.000000000 xmwsf (40) = 10.000000000 xmwoff(40) = 0.000000000 xmwsf (41) = 10.000000000 xmwoff(41) = 0.000000000 scale factors offsets record xmwsf (42) = 10.000000000 xmwoff(42) = 0.000000000 xmwsf (43) = 10.000000000 xmwoff(43) = 0.000000000 xmwsf (44) = 10.000000000 xmwoff(44) = 0.000000000 xmwsf (45) = 10.000000000 xmwoff(45) = 0.000000000 xmwsf (46) = 10.000000000 xmwoff(46) = 0.000000000 xmwsf (47) = 10.000000000 xmwoff(47) = 0.000000000 xmwsf (48) = 10.000000000 xmwoff(48) = 0.000000000 xmwsf (49) = 10.000000000 xmwoff(49) = 0.000000000 xmwsf (50) = 10.000000000 xmwoff(50) = 0.000000000 xmwsf (51) = 10.000000000 xmwoff(51) = 0.000000000 xmwsf (52) = 10.000000000 xmwoff(52) = 0.000000000 xmwsf (53) = 10.000000000 xmwoff(53) = 0.000000000 xmwsf (54) = 10.000000000 xmwoff(54) = 0.000000000 xmwsf (55) = 10.000000000 xmwoff(55) = 0.000000000 xmwsf (56) = 10.000000000 xmwoff(56) = 0.000000000 xmwsf (57) = 10.000000000 xmwoff(57) = 0.000000000 xmwsf (58) = 10.000000000 xmwoff(58) = 0.000000000 xmwsf (59) = 10.000000000 xmwoff(59) = 0.000000000 xmwsf (60) = 10.000000000 xmwoff(60) = 0.000000000 xmwsf (61) = 10.000000000 xmwoff(61) = 0.000000000 xmwsf (62) = 10.000000000 xmwoff(62) = 0.000000000 xmwsf (63) = 100.000000000 xmwoff(63) = 0.000000000 xmwsf (64) = 100.000000000 xmwoff(64) = 0.000000000 xmwsf (65) = 100.000000000 xmwoff(65) = 0.000000000 xmwsf (66) = 100.000000000 xmwoff(66) = 0.000000000 xmwsf (67) = 100.000000000 xmwoff(67) = -180.000000000 xmwsf (68) = 100.000000000 xmwoff(68) = -180.000000000 xmwsf (69) = 100.000000000 xmwoff(69) = -180.000000000 xmwsf (70) = 100.000000000 xmwoff(70) = -180.000000000 xmwsf (71) = 1.000000000 xmwoff(71) = 0.000000000 xmwsf (72) = 1.000000000 xmwoff(72) = 0.000000000 xmwsf (73) = 1.000000000 xmwoff(73) = 0.000000000 xmwsf (74) = 1.000000000 xmwoff(74) = 0.000000000 xmwsf (75) = 1.000000000 xmwoff(75) = 0.000000000 day: 13 record num = 1 2 3 4 5 Julian date xmwdt( 1) = 2449000.00 2449000.00 2449000.00 2449000.00 2449000.00 Julian time xmwdt( 2) = 0.5001388 0.5003240 0.5005091 0.5006943 0.5008795 Earth - Sun distance xmwdt( 3) = 0.9835182 0.9835182 0.9835182 0.9835182 0.9835182 s/c position, x xmwdt( 4) = 2344521.00 2261630.00 2178027.00 2093716.00 2008737.00 xmwdt( 5) = 2261631.00 2178027.00 2093717.00 2008737.00 1923096.00 s/c position, y xmwdt( 6) = 5525055.00 5507752.00 5488993.00 5468781.00 5447124.00 xmwdt( 7) = 5507752.00 5488993.00 5468780.00 5447124.00 5424031.00 s/c position, z xmwdt( 8) = -3537514.00 -3617817.00 -3697017.00 -3775103.00 -3852046.00 xmwdt( 9) = -3617817.00 -3697017.00 -3775103.00 -3852046.00 -3927829.00 s/c velocity (xdot) xmwdt(10) = -5560.00 -5605.00 -5648.00 -5689.00 -5729.00 xmwdt(11) = -5605.00 -5648.00 -5689.00 -5729.00 -5768.00 s/c velocity (ydot) xmwdt(12) = -865.00 -962.00 -1059.00 -1156.00 -1252.00 xmwdt(13) = -962.00 -1059.00 -1156.00 -1252.00 -1348.00 s/c velocity (zdot) xmwdt(14) = -5052.00 -4985.00 -4916.00 -4845.00 -4773.00 xmwdt(15) = -4985.00 -4916.00 -4845.00 -4773.00 -4699.00 day: 13 record num = 1 2 3 4 5 s/c nadir position, colatitude xmwdt(16) = 120.51 121.28 122.05 122.81 123.56 xmwdt(17) = 121.28 122.05 122.81 123.56 124.31 s/c nadir position, longitude xmwdt(18) = 67.01 67.68 68.36 69.05 69.76 xmwdt(19) = 67.68 68.36 69.05 69.76 70.48 sun position, colatitude xmwdt(20) = 111.51 111.51 111.51 111.51 111.51 sun position, longitude xmwdt(21) = 182.10 182.04 181.97 181.90 181.84 orbit number key xmwdt(22) = 32767.00 32767.00 32767.00 32767.00 32767.00 wfov, radiometric, total xmwdt(23) = 247.10 246.50 246.40 245.00 243.90 xmwdt(24) = 246.80 246.50 245.80 244.80 243.50 xmwdt(25) = 246.90 246.20 245.70 244.00 243.20 xmwdt(26) = 246.70 246.10 245.50 243.70 242.90 wfov, radiometric, sw xmwdt(27) = 14.80 16.30 18.70 20.70 23.20 xmwdt(28) = 15.00 17.00 19.00 21.00 23.70 xmwdt(29) = 15.80 17.40 19.70 21.70 24.60 xmwdt(30) = 16.00 17.90 20.20 22.30 25.60 mfov, radiometric, total xmwdt(31) = 129.80 128.00 127.60 125.20 123.30 xmwdt(32) = 129.50 128.20 127.30 124.80 122.80 xmwdt(33) = 128.90 127.40 126.90 123.70 122.00 xmwdt(34) = 128.90 127.50 126.40 123.50 121.60 day: 13 record num = 1 2 3 4 5 mfov, radiometric, sw xmwdt(35) = 1.90 3.10 3.10 2.70 3.20 xmwdt(36) = 3.20 3.40 3.20 2.30 32767.00 xmwdt(37) = 1.80 32767.00 32767.00 2.80 2.70 xmwdt(38) = 2.90 2.70 32767.00 1.90 3.00 wfov, unfiltered, sw xmwdt(39) = 14.80 16.30 18.70 20.70 23.20 xmwdt(40) = 15.00 17.00 19.00 21.00 23.70 xmwdt(41) = 15.80 17.40 19.70 21.70 24.60 xmwdt(42) = 16.00 17.90 20.20 22.30 25.60 wfov, unfiltered, lw xmwdt(43) = 232.30 230.30 227.70 224.30 220.70 xmwdt(44) = 231.80 229.50 226.80 223.80 219.90 xmwdt(45) = 231.10 228.80 226.00 222.30 218.60 xmwdt(46) = 230.70 228.20 225.30 221.40 217.30 mfov, unfiltered, sw xmwdt(47) = 0.00 0.00 0.00 0.00 0.00 xmwdt(48) = 0.00 0.00 0.00 0.00 0.00 xmwdt(49) = 0.00 0.00 0.00 0.00 2.70 xmwdt(50) = 0.00 0.00 0.00 0.00 3.00 mfov, unfiltered, lw xmwdt(51) = 129.80 128.00 127.60 125.20 123.30 xmwdt(52) = 129.50 128.20 127.30 124.80 122.80 xmwdt(53) = 128.90 127.40 126.90 123.70 119.40 xmwdt(54) = 128.90 127.50 126.40 123.50 118.60 wfov, toa est., nf, sw xmwdt(55) = 32767.00 32767.00 32767.00 32767.00 32767.00 day: 13 record num = 1 2 3 4 5 wfov, toa est., nf, lw xmwdt(56) = 32767.00 32767.00 32767.00 32767.00 32767.00 mfov, toa est., nf, sw xmwdt(57) = 32767.00 32767.00 32767.00 32767.00 32767.00 mfov, toa est., nf, lw xmwdt(58) = 32767.00 32767.00 32767.00 32767.00 32767.00 wfov, toa est., sf, sw xmwdt(59) = 32767.00 32767.00 32767.00 32767.00 32767.00 wfov, toa est., sf, lw xmwdt(60) = 275.40 275.40 268.60 268.60 259.40 mfov, toa est., sf, sw xmwdt(61) = 32767.00 32767.00 32767.00 32767.00 32767.00 mfov, toa est., sf, lw xmwdt(62) = 291.70 291.70 285.30 285.30 269.10 nonscanner, fov, colatitude xmwdt(63) = 120.57 121.34 122.10 122.86 123.61 xmwdt(64) = 120.76 121.53 122.29 123.05 123.80 xmwdt(65) = 120.95 121.72 122.48 123.23 123.99 xmwdt(66) = 121.15 121.91 122.67 123.42 124.18 nonscanner, fov, longitude xmwdt(67) = 67.04 67.71 68.39 69.09 69.79 xmwdt(68) = 67.21 67.88 68.56 69.26 69.97 xmwdt(69) = 67.37 68.05 68.74 69.44 70.15 xmwdt(70) = 67.54 68.22 68.91 69.61 70.33 day: 13 record num = 1 2 3 4 5 nonscanner, op. mode, flag word (octal) xmwdt(71) = 3710 3700 3700 3700 3700 xmwdt(72) = 3750 3750 3750 3750 3750 flag, nonscanner, toa est. (octal) xmwdt(73) = 1 0 1 0 1 orbit number scale factor xmwdt(74) = 32767.00 32767.00 32767.00 32767.00 32767.00 spare xmwdt(75) = 32767.00 32767.00 32767.00 32767.00 32767.00