#!/bin/csh -x

set platform = "MODIS/Terra"

set input_data_directory = "/archive/MOD08/2002/014/"

set short_filename = "MOD08_D3.A2002014.003.2002020031427.hdf"

set image_date = "14Jan2002"

set output_data_directory = "/products/MOD08_D3/_BROWSE_FIXEDSCALE/2002_01_14/CURRENT/LL/"



cat > global_variables_shortlist.dat <<eof
Cloud_Phase_Infrared_Histogram_Counts
0.0
800.0
eof



cat > global_variables.dat <<eof
Scattering_Angle_Mean
0.0
180.0
Scattering_Angle_Standard_Deviation
0.0
60.0
Scattering_Angle_Minimum
0.0
180.0
Scattering_Angle_Maximum
0.0
180.0
Scattering_Angle_Pixel_Counts
0.0
200.0
Optical_Depth_Land_And_Ocean_Mean
0.0
0.8
Optical_Depth_Land_And_Ocean_Standard_Deviation
0.0
0.4
Optical_Depth_Land_And_Ocean_Minimum
0.0
0.8
Optical_Depth_Land_And_Ocean_Maximum
0.0
0.8
Optical_Depth_Land_And_Ocean_Pixel_Counts
0.0
200.0
Continental_Optical_Depth_Land_Mean
0.0
0.8
Continental_Optical_Depth_Land_Standard_Deviation
0.0
0.4
Continental_Optical_Depth_Land_Minimum
0.0
0.8
Continental_Optical_Depth_Land_Maximum
0.0
0.8
Continental_Optical_Depth_Land_QA_Mean
0.0
0.8
Continental_Optical_Depth_Land_QA_Standard_Deviation
0.0
0.4
Corrected_Optical_Depth_Smoke_Mean
0.0
0.8
Corrected_Optical_Depth_Smoke_Standard_Deviation
0.0
0.4
Corrected_Optical_Depth_Smoke_Minimum
0.0
0.8
Corrected_Optical_Depth_Smoke_Maximum
0.0
0.8
Corrected_Optical_Depth_Smoke_QA_Mean
0.0
0.8
Corrected_Optical_Depth_Smoke_QA_Standard_Deviation
0.0
0.4
Corrected_Optical_Depth_Sulfate_Mean
0.0
0.8
Corrected_Optical_Depth_Sulfate_Standard_Deviation
0.0
0.4
Corrected_Optical_Depth_Sulfate_Minimum
0.0
0.8
Corrected_Optical_Depth_Sulfate_Maximum
0.0
0.8
Corrected_Optical_Depth_Sulfate_QA_Mean
0.0
0.8
Corrected_Optical_Depth_Sulfate_QA_Standard_Deviation
0.0
0.4
Corrected_Optical_Depth_Dust_Mean
0.0
0.8
Corrected_Optical_Depth_Dust_Standard_Deviation
0.0
0.4
Corrected_Optical_Depth_Dust_Minimum
0.0
0.8
Corrected_Optical_Depth_Dust_Maximum
0.0
0.8
Corrected_Optical_Depth_Dust_QA_Mean
0.0
0.8
Corrected_Optical_Depth_Dust_QA_Standard_Deviation
0.0
0.4
Mass_Concentration_Land_Mean
0.0
20.0
Mass_Concentration_Land_Standard_Deviation
0.0
6.66
Mass_Concentration_Land_Minimum
0.0
20.0
Mass_Concentration_Land_Maximum
0.0
20.0
Mass_Concentration_Land_QA_Mean
0.0
20.0
Mass_Concentration_Land_QA_Standard_Deviation
0.0
6.66
Angstrom_Exponent_Land_Mean
-0.5
3.0
Angstrom_Exponent_Land_Standard_Deviation
0.0
1.0
Angstrom_Exponent_Land_Minimum
-0.5
3.0
Angstrom_Exponent_Land_Maximum
-0.5
3.0
Angstrom_Exponent_Land_QA_Mean
-0.5
3.0
Angstrom_Exponent_Land_QA_Standard_Deviation
-0.5
3.0
Transmitted_Flux_Land_Mean
0.5
1.0
Transmitted_Flux_Land_Standard_Deviation
0.0
0.166
Transmitted_Flux_Land_Minimum
0.5
1.0
Transmitted_Flux_Land_Maximum
0.5
1.0
Transmitted_Flux_Land_QA_Mean
0.5
1.0
Transmitted_Flux_Land_QA_Standard_Deviation
0.0
0.166
Reflected_Flux_Land_Mean
0.0
0.5
Reflected_Flux_Land_Standard_Deviation
0.0
0.166
Reflected_Flux_Land_Minimum
0.0
0.5
Reflected_Flux_Land_Maximum
0.0
0.5
Reflected_Flux_Land_QA_Mean
0.0
0.5
Reflected_Flux_Land_QA_Standard_Deviation
0.0
0.166
Mean_Reflectance_Land_All_QA47_Mean
0.0
0.5
Mean_Reflectance_Land_All_QA47_Standard_Deviation
0.0
0.166
Mean_Reflectance_Land_All_QA47_Minimum
0.0
0.5
Mean_Reflectance_Land_All_QA47_Maximum
0.0
0.5
Mean_Reflectance_Land_All_QA47_QA_Mean
0.0
0.5
Mean_Reflectance_Land_All_QA47_QA_Standard_Deviation
0.0
0.166
Mean_Reflectance_Land_All_QA66_Mean
0.0
0.5
Mean_Reflectance_Land_All_QA66_Standard_Deviation
0.0
0.166
Mean_Reflectance_Land_All_QA66_Minimum
0.0
0.5
Mean_Reflectance_Land_All_QA66_Maximum
0.0
0.5
Mean_Reflectance_Land_All_QA66_QA_Mean
0.0
0.5
Mean_Reflectance_Land_All_QA66_QA_Standard_Deviation
0.0
0.166
Path_Radiance_Land_QA47_Mean
0.0
0.25
Path_Radiance_Land_QA47_Standard_Deviation
0.0
0.08
Path_Radiance_Land_QA47_Minimum
0.0
0.25
Path_Radiance_Land_QA47_Maximum
0.0
0.25
Path_Radiance_Land_QA47_QA_Mean
0.0
0.25
Path_Radiance_Land_QA47_QA_Standard_Deviation
0.0
0.08
Path_Radiance_Land_QA66_Mean
0.0
0.25
Path_Radiance_Land_QA66_Standard_Deviation
0.0
0.08
Path_Radiance_Land_QA66_Minimum
0.0
0.25
Path_Radiance_Land_QA66_Maximum
0.0
0.25
Path_Radiance_Land_QA66_QA_Mean
0.0
0.25
Path_Radiance_Land_QA66_QA_Standard_Deviation
0.0
0.08
Critical_Reflectance_Land_QA47_Mean
0.0
1.0
Critical_Reflectance_Land_QA47_Standard_Deviation
0.0
0.3
Critical_Reflectance_Land_QA47_Minimum
0.0
1.0
Critical_Reflectance_Land_QA47_Maximum
0.0
1.0
Critical_Reflectance_Land_QA47_QA_Mean
0.0
1.0
Critical_Reflectance_Land_QA47_QA_Standard_Deviation
0.0
0.3
Critical_Reflectance_Land_QA66_Mean
0.0
1.0
Critical_Reflectance_Land_QA66_Standard_Deviation
0.0
0.3
Critical_Reflectance_Land_QA66_Minimum
0.0
1.0
Critical_Reflectance_Land_QA66_Maximum
0.0
1.0
Critical_Reflectance_Land_QA66_QA_Mean
0.0
1.0
Critical_Reflectance_Land_QA66_QA_Standard_Deviation
0.0
0.3
Mass_Concentration_Ocean_Mean
0.0
3.0
Mass_Concentration_Ocean_Standard_Deviation
0.0
1.0
Mass_Concentration_Ocean_Minimum
0.0
3.0
Mass_Concentration_Ocean_Maximum
0.0
3.0
Mass_Concentration_Ocean_QA_Mean
0.0
3.0
Mass_Concentration_Ocean_QA_Standard_Deviation
0.0
1.0
Effective_Optical_Depth_Average_Ocean_Mean
0.0
0.8
Effective_Optical_Depth_Average_Ocean_Standard_Deviation
0.0
0.4
Effective_Optical_Depth_Average_Ocean_Minimum
0.0
0.8
Effective_Optical_Depth_Average_Ocean_Maximum
0.0
0.8
Effective_Optical_Depth_Average_Ocean_QA_Mean
0.0
0.8
Effective_Optical_Depth_Average_Ocean_QA_Standard_Deviation
0.0
0.4
Optical_Depth_Ratio_Small_Ocean_Mean
0.0
1.0
Optical_Depth_Ratio_Small_Ocean_Standard_Deviation
0.0
0.3
Optical_Depth_Ratio_Small_Ocean_Minimum
0.0
1.0
Optical_Depth_Ratio_Small_Ocean_Maximum
0.0
1.0
Optical_Depth_Ratio_Small_Ocean_QA_Mean
0.0
1.0
Optical_Depth_Ratio_Small_Ocean_QA_Standard_Deviation
0.0
0.3
Effective_Radius_Ocean_Mean
0.0
3.0
Effective_Radius_Ocean_Standard_Deviation
0.0
1.0
Effective_Radius_Ocean_Minimum
0.0
3.0
Effective_Radius_Ocean_Maximum
0.0
3.0
Effective_Radius_Ocean_QA_Mean
0.0
3.0
Effective_Radius_Ocean_QA_Standard_Deviation
0.0
1.0
Cloud_Condensation_Nuclei_Ocean_Mean
0.0
2.0
Cloud_Condensation_Nuclei_Ocean_Standard_Deviation
0.0
1.0
Cloud_Condensation_Nuclei_Ocean_Minimum
0.0
2.0
Cloud_Condensation_Nuclei_Ocean_Maximum
0.0
2.0
Cloud_Condensation_Nuclei_Ocean_QA_Mean
0.0
2.0
Cloud_Condensation_Nuclei_Ocean_QA_Standard_Deviation
0.0
1.0
Angstrom_Exponent_1_Ocean_Mean
-0.5
3.0
Angstrom_Exponent_1_Ocean_Standard_Deviation
0.0
1.0
Angstrom_Exponent_1_Ocean_Minimum
-0.5
3.0
Angstrom_Exponent_1_Ocean_Maximum
-0.5
3.0
Angstrom_Exponent_1_Ocean_QA_Mean
-0.5
3.0
Angstrom_Exponent_1_Ocean_QA_Standard_Deviation
0.0
1.0
Angstrom_Exponent_2_Ocean_Mean
-0.5
3.0
Angstrom_Exponent_2_Ocean_Standard_Deviation
0.0
1.0
Angstrom_Exponent_2_Ocean_Minimum
-0.5
3.0
Angstrom_Exponent_2_Ocean_Maximum
-0.5
3.0
Angstrom_Exponent_2_Ocean_QA_Mean
-0.5
3.0
Angstrom_Exponent_2_Ocean_QA_Standard_Deviation
0.0
1.0
Asymmetry_Factor_Average_Ocean_Mean
0.0
1.0
Asymmetry_Factor_Average_Ocean_Standard_Deviation
0.0
0.3
Asymmetry_Factor_Average_Ocean_Minimum
0.0
1.0
Asymmetry_Factor_Average_Ocean_Maximum
0.0
1.0
Asymmetry_Factor_Average_Ocean_QA_Mean
0.0
1.0
Asymmetry_Factor_Average_Ocean_QA_Standard_Deviation
0.0
0.3
Backscattering_Ratio_Average_Ocean_Mean
0.0
0.5
Backscattering_Ratio_Average_Ocean_Standard_Deviation
0.0
0.2
Backscattering_Ratio_Average_Ocean_Minimum
0.0
0.5
Backscattering_Ratio_Average_Ocean_Maximum
0.0
0.5
Backscattering_Ratio_Average_Ocean_QA_Mean
0.0
0.5
Backscattering_Ratio_Average_Ocean_QA_Standard_Deviation
0.0
0.2
Transmitted_Flux_Average_Ocean_Mean
0.5
1.0
Transmitted_Flux_Average_Ocean_Standard_Deviation
0.0
0.166
Transmitted_Flux_Average_Ocean_Minimum
0.5
1.0
Transmitted_Flux_Average_Ocean_Maximum
0.5
1.0
Transmitted_Flux_Average_Ocean_QA_Mean
0.5
1.0
Transmitted_Flux_Average_Ocean_QA_Standard_Deviation
0.0
0.166
Reflected_Flux_Average_Ocean_Mean
0.0
0.5
Reflected_Flux_Average_Ocean_Standard_Deviation
0.0
0.166
Reflected_Flux_Average_Ocean_Minimum
0.0
0.5
Reflected_Flux_Average_Ocean_Maximum
0.0
0.5
Reflected_Flux_Average_Ocean_QA_Mean
0.0
0.5
Reflected_Flux_Average_Ocean_QA_Standard_Deviation
0.0
0.166
Optical_Depth_by_models_ocean_Mean
0.0
0.3
Optical_Depth_by_models_ocean_Standard_Deviation
0.0
0.1
Optical_Depth_by_models_ocean_Minimum
0.0
0.3
Optical_Depth_by_models_ocean_Maximum
0.0
0.3
Optical_Depth_by_models_ocean_QA_Mean
0.0
0.3
Optical_Depth_by_models_ocean_QA_Standard_Deviation
0.0
0.1
Water_Vapor_Near_Infrared_Clear_Mean
0.0
7.5
Water_Vapor_Near_Infrared_Clear_Standard_Deviation
0.0
3.0
Water_Vapor_Near_Infrared_Clear_Minimum
0.0
7.5
Water_Vapor_Near_Infrared_Clear_Maximum
0.0
7.5
Water_Vapor_Near_Infrared_Clear_QA_Mean
0.0
7.5
Water_Vapor_Near_Infrared_Clear_QA_Standard_Deviation
0.0
3.0
Water_Vapor_Near_Infrared_Cloud_Mean
0.0
7.5
Water_Vapor_Near_Infrared_Cloud_Standard_Deviation
0.0
3.0
Water_Vapor_Near_Infrared_Cloud_Minimum
0.0
7.5
Water_Vapor_Near_Infrared_Cloud_Maximum
0.0
7.5
Water_Vapor_Near_Infrared_Cloud_QA_Mean
0.0
7.5
Water_Vapor_Near_Infrared_Cloud_QA_Standard_Deviation
0.0
3.0
Cirrus_Reflectance_Mean
0.0
0.3
Cirrus_Reflectance_Standard_Deviation
0.0
0.2
Cirrus_Reflectance_Minimum
0.0
0.1
Cirrus_Reflectance_Maximum
0.0
1.0
Cirrus_Reflectance_Pixel_Counts
0.0
1000.0
Cirrus_Fraction_SWIR
0.0
1.0
Cirrus_Fraction_SWIR_Pixel_Counts
0.0
600.0
Cloud_Top_Temperature_Mean
200.0
300.0
Cloud_Top_Temperature_Standard_Deviation
0.0
50.0
Cloud_Top_Temperature_Minimum
200.0
300.0
Cloud_Top_Temperature_Maximum
200.0
300.0
Cloud_Top_Temperature_QA_Mean
200.0
300.0
Cloud_Top_Temperature_QA_Standard_Deviation
0.0
50.0
Cloud_Top_Temperature_Day_Mean
200.0
300.0
Cloud_Top_Temperature_Day_Standard_Deviation
0.0
50.0
Cloud_Top_Temperature_Day_Minimum
200.0
300.0
Cloud_Top_Temperature_Day_Maximum
200.0
300.0
Cloud_Top_Temperature_Day_QA_Mean
200.0
300.0
Cloud_Top_Temperature_Day_QA_Standard_Deviation
0.0
50.0
Cloud_Top_Temperature_Night_Mean
200.0
300.0
Cloud_Top_Temperature_Night_Standard_Deviation
0.0
50.0
Cloud_Top_Temperature_Night_Minimum
200.0
300.0
Cloud_Top_Temperature_Night_Maximum
200.0
300.0
Cloud_Top_Temperature_Night_QA_Mean
200.0
300.0
Cloud_Top_Temperature_Night_QA_Standard_Deviation
0.0
50.0
Cloud_Top_Pressure_Mean
300.0
1000.0
Cloud_Top_Pressure_Standard_Deviation
0.0
300.0
Cloud_Top_Pressure_Minimum
300.0
1000.0
Cloud_Top_Pressure_Maximum
300.0
1000.0
Cloud_Top_Pressure_QA_Mean
300.0
1000.0
Cloud_Top_Pressure_QA_Standard_Deviation
0.0
300.0
Cloud_Top_Pressure_Day_Mean
300.0
1000.0
Cloud_Top_Pressure_Day_Standard_Deviation
0.0
300.0
Cloud_Top_Pressure_Day_Minimum
300.0
1000.0
Cloud_Top_Pressure_Day_Maximum
300.0
1000.0
Cloud_Top_Pressure_Day_QA_Mean
300.0
1000.0
Cloud_Top_Pressure_Day_QA_Standard_Deviation
0.0
300.0
Cloud_Top_Pressure_Night_Mean
300.0
1000.0
Cloud_Top_Pressure_Night_Standard_Deviation
0.0
300.0
Cloud_Top_Pressure_Night_Minimum
300.0
1000.0
Cloud_Top_Pressure_Night_Maximum
300.0
1000.0
Cloud_Top_Pressure_Night_QA_Mean
300.0
1000.0
Cloud_Top_Pressure_Night_QA_Standard_Deviation
0.0
300.0
Cloud_Effective_Emissivity_Mean
0.0
1.0
Cloud_Effective_Emissivity_Standard_Deviation
0.0
0.3
Cloud_Effective_Emissivity_Minimum
0.0
1.0
Cloud_Effective_Emissivity_Maximum
0.0
1.0
Cloud_Effective_Emissivity_QA_Mean
0.0
1.0
Cloud_Effective_Emissivity_QA_Standard_Deviation
0.0
0.3
Cloud_Effective_Emissivity_Day_Mean
0.0
1.0
Cloud_Effective_Emissivity_Day_Standard_Deviation
0.0
0.3
Cloud_Effective_Emissivity_Day_Minimum
0.0
1.0
Cloud_Effective_Emissivity_Day_Maximum
0.0
1.0
Cloud_Effective_Emissivity_Day_QA_Mean
0.0
1.0
Cloud_Effective_Emissivity_Day_QA_Standard_Deviation
0.0
0.3
Cloud_Effective_Emissivity_Night_Mean
0.0
1.0
Cloud_Effective_Emissivity_Night_Standard_Deviation
0.0
0.3
Cloud_Effective_Emissivity_Night_Minimum
0.0
1.0
Cloud_Effective_Emissivity_Night_Maximum
0.0
1.0
Cloud_Effective_Emissivity_Night_QA_Mean
0.0
1.0
Cloud_Effective_Emissivity_Night_QA_Standard_Deviation
0.0
0.3
Cloud_Fraction_Infrared_Mean
0.0
1.0
Cloud_Fraction_Infrared_Standard_Deviation
0.0
0.3
Cloud_Fraction_Infrared_Minimum
0.0
1.0
Cloud_Fraction_Infrared_Maximum
0.0
1.0
Cloud_Fraction_Infrared_QA_Mean
0.0
1.0
Cloud_Fraction_Infrared_QA_Standard_Deviation
0.0
0.3
Cloud_Fraction_Infrared_Day_Mean
0.0
1.0
Cloud_Fraction_Infrared_Day_Standard_Deviation
0.0
0.3
Cloud_Fraction_Infrared_Day_Minimum
0.0
1.0
Cloud_Fraction_Infrared_Day_Maximum
0.0
1.0
Cloud_Fraction_Infrared_Day_QA_Mean
0.0
1.0
Cloud_Fraction_Infrared_Day_QA_Standard_Deviation
0.0
0.3
Cloud_Fraction_Infrared_Night_Mean
0.0
1.0
Cloud_Fraction_Infrared_Night_Standard_Deviation
0.0
0.3
Cloud_Fraction_Infrared_Night_Minimum
0.0
1.0
Cloud_Fraction_Infrared_Night_Maximum
0.0
1.0
Cloud_Fraction_Infrared_Night_QA_Mean
0.0
1.0
Cloud_Fraction_Infrared_Night_QA_Standard_Deviation
0.0
0.3
Cirrus_Fraction_Infrared
0.0
1.0
Cirrus_Fraction_Infrared_Pixel_Counts
0.0
1200.0
High_Cloud_Fraction_Infrared
0.0
1.0
High_Cloud_Fraction_Infrared_Pixel_Counts
0.0
1200.0
Cloud_Phase_Infrared_Histogram_Counts
0.0
800.0
Cloud_Phase_Infrared_Day_Histogram_Counts
0.0
800.0
Cloud_Phase_Infrared_Night_Histogram_Counts
0.0
800.0
Cloud_Optical_Thickness_Combined_Mean
0.0
20.0
Cloud_Optical_Thickness_Combined_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Combined_Minimum
0.0
10.0
Cloud_Optical_Thickness_Combined_Maximum
0.0
100.0
Cloud_Optical_Thickness_Combined_QA_Mean
0.0
20.0
Cloud_Optical_Thickness_Combined_QA_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Water_Mean
0.0
20.0
Cloud_Optical_Thickness_Water_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Water_Minimum
0.0
10.0
Cloud_Optical_Thickness_Water_Maximum
0.0
100.0
Cloud_Optical_Thickness_Water_QA_Mean
0.0
20.0
Cloud_Optical_Thickness_Water_QA_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Ice_Mean
0.0
20.0
Cloud_Optical_Thickness_Ice_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Ice_Minimum
0.0
10.0
Cloud_Optical_Thickness_Ice_Maximum
0.0
100.0
Cloud_Optical_Thickness_Ice_QA_Mean
0.0
20.0
Cloud_Optical_Thickness_Ice_QA_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Undetermined_Mean
0.0
20.0
Cloud_Optical_Thickness_Undetermined_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Undetermined_Minimum
0.0
10.0
Cloud_Optical_Thickness_Undetermined_Maximum
0.0
100.0
Cloud_Optical_Thickness_Undetermined_QA_Mean
0.0
20.0
Cloud_Optical_Thickness_Undetermined_QA_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Mixed_Phase_Mean
0.0
20.0
Cloud_Optical_Thickness_Mixed_Phase_Standard_Deviation
0.0
30.0
Cloud_Optical_Thickness_Mixed_Phase_Minimum
0.0
10.0
Cloud_Optical_Thickness_Mixed_Phase_Maximum
0.0
100.0
Cloud_Optical_Thickness_Mixed_Phase_QA_Mean
0.0
20.0
Cloud_Optical_Thickness_Mixed_Phase_QA_Standard_Deviation
0.0
30.0
Cloud_Effective_Radius_Combined_Mean
4.0
40.0
Cloud_Effective_Radius_Combined_Standard_Deviation
0.0
20.0
Cloud_Effective_Radius_Combined_Minimum
4.0
40.0
Cloud_Effective_Radius_Combined_Maximum
4.0
40.0
Cloud_Effective_Radius_Combined_QA_Mean
4.0
40.0
Cloud_Effective_Radius_Combined_QA_Standard_Deviation
0.0
20.0
Cloud_Effective_Radius_Water_Mean
2.0
30.0
Cloud_Effective_Radius_Water_Standard_Deviation
0.0
10.0
Cloud_Effective_Radius_Water_Minimum
2.0
30.0
Cloud_Effective_Radius_Water_Maximum
2.0
30.0
Cloud_Effective_Radius_Water_QA_Mean
2.0
30.0
Cloud_Effective_Radius_Water_QA_Standard_Deviation
0.0
10.0
Cloud_Effective_Radius_Ice_Mean
6.0
50.0
Cloud_Effective_Radius_Ice_Standard_Deviation
0.0
20.0
Cloud_Effective_Radius_Ice_Minimum
6.0
50.0
Cloud_Effective_Radius_Ice_Maximum
6.0
50.0
Cloud_Effective_Radius_Ice_QA_Mean
6.0
50.0
Cloud_Effective_Radius_Ice_QA_Standard_Deviation
0.0
20.0
Cloud_Effective_Radius_Undetermined_Mean
2.0
30.0
Cloud_Effective_Radius_Undetermined_Standard_Deviation
0.0
10.0
Cloud_Effective_Radius_Undetermined_Minimum
2.0
30.0
Cloud_Effective_Radius_Undetermined_Maximum
2.0
30.0
Cloud_Effective_Radius_Undetermined_QA_Mean
2.0
30.0
Cloud_Effective_Radius_Undetermined_QA_Standard_Deviation
0.0
10.0
Cloud_Effective_Radius_Mixed_Phase_Mean
4.0
40.0
Cloud_Effective_Radius_Mixed_Phase_Standard_Deviation
0.0
20.0
Cloud_Effective_Radius_Mixed_Phase_Minimum
4.0
40.0
Cloud_Effective_Radius_Mixed_Phase_Maximum
4.0
40.0
Cloud_Effective_Radius_Mixed_Phase_QA_Mean
4.0
40.0
Cloud_Effective_Radius_Mixed_Phase_QA_Standard_Deviation
0.0
20.0
Cloud_Fraction_Water
0.0
1.0
Cloud_Fraction_Water_Pixel_Counts
0.0
1000.0
Cloud_Fraction_Ice
0.0
1.0
Cloud_Fraction_Ice_Pixel_Counts
0.0
1000.0
Cloud_Fraction_Undetermined
0.0
1.0
Cloud_Fraction_Undetermined_Pixel_Counts
0.0
1000.0
Cloud_Fraction_Mixed_Phase
0.0
1.0
Cloud_Fraction_Mixed_Phase_Pixel_Counts
0.0
1000.0
Cloud_Fraction_Total
0.0
1.0
Cloud_Fraction_Total_Pixel_Counts
0.0
1000.0
Cloud_Fraction_SWIR_Water
0.0
1.0
Cloud_Fraction_SWIR_Water_Pixel_Counts
0.0
1000.0
Cloud_Fraction_SWIR_Ice
0.0
1.0
Cloud_Fraction_SWIR_Ice_Pixel_Counts
0.0
1000.0
Cloud_Fraction_SWIR_Total
0.0
1.0
Cloud_Fraction_SWIR_Total_Pixel_Counts
0.0
1000.0
Water_Path_Combined_Mean
0.0
400.0
Water_Path_Combined_Standard_Deviation
0.0
130.0
Water_Path_Combined_Minimum
0.0
400.0
Water_Path_Combined_Maximum
0.0
400.0
Water_Path_Combined_QA_Mean
0.0
400.0
Water_Path_Combined_QA_Standard_Deviation
0.0
130.0
Water_Path_Liquid_Mean
0.0
400.0
Water_Path_Liquid_Standard_Deviation
0.0
130.0
Water_Path_Liquid_Minimum
0.0
400.0
Water_Path_Liquid_Maximum
0.0
400.0
Water_Path_Liquid_QA_Mean
0.0
400.0
Water_Path_Liquid_QA_Standard_Deviation
0.0
130.0
Water_Path_Ice_Mean
0.0
400.0
Water_Path_Ice_Standard_Deviation
0.0
130.0
Water_Path_Ice_Minimum
0.0
400.0
Water_Path_Ice_Maximum
0.0
400.0
Water_Path_Ice_QA_Mean
0.0
400.0
Water_Path_Ice_QA_Standard_Deviation
0.0
130.0
Water_Path_Undetermined_Mean
0.0
400.0
Water_Path_Undetermined_Standard_Deviation
0.0
130.0
Water_Path_Undetermined_Minimum
0.0
400.0
Water_Path_Undetermined_Maximum
0.0
400.0
Water_Path_Undetermined_QA_Mean
0.0
400.0
Water_Path_Undetermined_QA_Standard_Deviation
0.0
130.0
Water_Path_Mixed_Phase_Mean
0.0
400.0
Water_Path_Mixed_Phase_Standard_Deviation
0.0
130.0
Water_Path_Mixed_Phase_Minimum
0.0
400.0
Water_Path_Mixed_Phase_Maximum
0.0
400.0
Water_Path_Mixed_Phase_QA_Mean
0.0
400.0
Water_Path_Mixed_Phase_QA_Standard_Deviation
0.0
130.0
Total_Ozone_Mean
100.0
500.0
Total_Ozone_Standard_Deviation
0.0
200.0
Total_Ozone_Minimum
100.0
500.0
Total_Ozone_Maximum
100.0
500.0
Total_Ozone_QA_Mean
100.0
500.0
Total_Ozone_QA_Standard_Deviation
0.0
200.0
Total_Totals_Mean
0.0
80.0
Total_Totals_Standard_Deviation
0.0
25.0
Total_Totals_Minimum
0.0
80.0
Total_Totals_Maximum
0.0
80.0
Total_Totals_QA_Mean
0.0
80.0
Total_Totals_QA_Standard_Deviation
0.0
25.0
Lifted_Index_Mean
-20.0
40.0
Lifted_Index_Standard_Deviation
0.0
20.0
Lifted_Index_Minimum
-20.0
40.0
Lifted_Index_Maximum
-20.0
40.0
Lifted_Index_QA_Mean
-20.0
40.0
Lifted_Index_QA_Standard_Deviation
0.0
20.0
K_Index_Mean
-5.0
65.0
K_Index_Standard_Deviation
0.0
25.0
K_Index_Minimum
-5.0
65.0
K_Index_Maximum
-5.0
65.0
K_Index_QA_Mean
-5.0
65.0
K_Index_QA_Standard_Deviation
0.0
25.0
Atmospheric_Water_Vapor_Mean
0.0
7.5
Atmospheric_Water_Vapor_Standard_Deviation
0.0
3.0
Atmospheric_Water_Vapor_Minimum
0.0
7.5
Atmospheric_Water_Vapor_Maximum
0.0
7.5
Atmospheric_Water_Vapor_QA_Mean
0.0
7.5
Atmospheric_Water_Vapor_QA_Standard_Deviation
0.0
3.0
Atmospheric_Water_Vapor_Low_Mean
0.0
2.5
Atmospheric_Water_Vapor_Low_Standard_Deviation
0.0
1.0
Atmospheric_Water_Vapor_Low_Minimum
0.0
2.5
Atmospheric_Water_Vapor_Low_Maximum
0.0
2.5
Atmospheric_Water_Vapor_Low_QA_Mean
0.0
2.5
Atmospheric_Water_Vapor_Low_QA_Standard_Deviation
0.0
1.0
Atmospheric_Water_Vapor_High_Mean
0.0
2.5
Atmospheric_Water_Vapor_High_Standard_Deviation
0.0
1.0
Atmospheric_Water_Vapor_High_Minimum
0.0
2.5
Atmospheric_Water_Vapor_High_Maximum
0.0
2.5
Atmospheric_Water_Vapor_High_QA_Mean
0.0
2.5
Atmospheric_Water_Vapor_High_QA_Standard_Deviation
0.0
1.0
eof


## test ## 
## cp global_variables_shortlist.dat global_variables.dat


cat > make_gif.pro <<eof
pro make_gif



;;;;;;;;;;;;;;;;
; set up display
;;;;;;;;;;;;;;;;
;device, pseudo=8
;window, /free



num_colors = !D.N_Colors
thisDevice = !D.Name
help,/device








; *************************************************
; ***** Need this for display to TV Screen!!! *****
; *************************************************

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Open Graphics Window
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
window, XSize=600, YSize=210, /free

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Setup Hardware Fonts
;   these will only work on display to X (hardware)
;   you cannot use hardware fonts to Z (software)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; use hardware fonts
!p.font=0

; regular 10pt helvetica font
Device, Set_Font='-adobe-helvetica-medium-r-normal--10-100-75-75-p-56-iso8859-1'

; *************************************************
; ***** END *****
; *************************************************






;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Set up Z-Graphics Buffer
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;Set_Plot,'Z',/Copy

;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Use vector fonts (Z or X)
;;;;;;;;;;;;;;;;;;;;;;;;;;;
;Device, Set_Colors=256, Set_Resolution=[600,210]

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; END
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;






help,/device




print,"FIXED SCALE PROCESSING!"


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Set Color Table Assignments
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

print,' # of colors = ', !D.N_Colors

red=bytarr(256)
grn=bytarr(256)
blu=bytarr(256)









red=bytarr(256)
grn=bytarr(256)
blu=bytarr(256)








;;;;;;;;;;;;;;;;;;;;
;  FROM ALLEN CHU
;;;;;;;;;;;;;;;;;;;;
;  color table and load them into IDL for displaying color images. After
;  compiling the procedure, users need to run the
;  procedure with three parameters. The resulting three parameters
;  contain the red, green, and blue color tables. Finally,  the procedure
;  loads the color tables using TVLCT, Param1, Param2, Param3
;

LUTR=BYTARR(256)
LUTG=BYTARR(256)
LUTB=BYTARR(256)


;;;;;;;;;;;;;;;;;;
; red


;173 to 0
FOR INDEX = 0,50 DO BEGIN
 DYDX = 173.0 - ((173.0*INDEX)/50.0)
 LUTR(INDEX) = BYTE(DYDX)
END

LUTR(51:113)=0

FOR I=114,145 DO BEGIN
 DX=FLOAT(I)-113.
 DYDX=(255.-0.)/(146.-113.)
 LUTR(I)=BYTE((DYDX*DX)+0.)
END

LUTR(146:230)=255

;255 to 180
FOR I=231,253 DO BEGIN
 DYDX=255.0 - 75.0*(I-231.0)/(22.0)
 LUTR(I)=BYTE((DYDX))
END

;;;;;;;;;;;;;;;;;;
; green

LUTG(0:36)=0

FOR I=37,129 DO BEGIN
 DX=FLOAT(I)-36.
 DYDX=(255.-0.)/(130.-36.)
 LUTG(I)=BYTE((DYDX*DX)+0.)
END

LUTG(130:154)=255

FOR I=155,250 DO BEGIN
 DX=FLOAT(I)-154.
 DYDX=(0.-255.)/(251.-154.)
 LUTG(I)=BYTE((DYDX*DX)+255.)
END

LUTG(251:255)=0

;;;;;;;;;;;;;;;;;;
; blue

;Pauls Mod (goes to indigo at 0)
;170 to 255
FOR I=0,60 DO BEGIN
 DYDX = 255.0 - (85.0*(60-I)/60.0)
 LUTB(I)=BYTE(DYDX)
END

LUTB(61)=255

FOR I=62,64 DO BEGIN
 DX=FLOAT(I)-61.
 DYDX=(238.-255.)/(65.-61.)
 LUTB(I)=BYTE((DYDX*DX)+255.)

END

LUTB(65:77)=238

FOR I=78,124 DO BEGIN
 DX=FLOAT(I)-77.
 DYDX=(0.-238.)/(125.-77.)
 LUTB(I)=BYTE((DYDX*DX)+238.)
END

LUTB(125:255)=0


goto, jump4577









; load rainbow Color Table
; note rainbow table has the best contrast with
; the black fill color
; but is too green heavy!
;loadct,34
; load EOS-B (white+rainbow) Color Table
;loadct,27
; load rainbow+white Color Table
;loadct,39
; load EOS-A (black+rainbow) Color Table
;;;;;loadct,26

;get the EOS-A rgb color assignments
tvlct,r,g,b,/Get

help,r,g,b

print,' r(0:5) =',r(0:5)
print,' g(0:5) =',g(0:5)
print,' b(0:5) =',b(0:5)

; EOS-A goes Black/Black/Gray/Purple in 0/1/2/3

;for i=0,252 do begin
;red(i)=r(i+3)
;grn(i)=g(i+3)
;blu(i)=b(i+3)
;endfor
;assign 253 to 255 colors
;red(253)=r(255)
;grn(253)=g(255)
;blu(253)=b(255)

; trunc was to 180 ... now its 129
;Rescale EOS-A so that there is less red
for i=0,255 do begin
index = fix(129.0*(i/255.0)+3.0)
red(i)=r(index)
grn(i)=g(index)
blu(i)=b(index)
endfor



jump4577:


for i=0,255 do begin
red(i)=LUTR(i)
grn(i)=LUTG(i)
blu(i)=LUTB(i)
endfor


;assign 254 to white
red(254)=255
grn(254)=255
blu(254)=255


;assign 255 to black
red(255)=000
grn(255)=000
blu(255)=000






print,'modified red gun table (240:255) = ',red(240:255)
print,'modified grn gun table (240:255) = ',grn(240:255)
print,'modified blu gun table (240:255) = ',blu(240:255)

;reload modified color table
tvlct,red,grn,blu




;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;Set scale size of title and annotations (1.0 is default)
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
sizenum_title=0.8
sizenum=0.75
sizenum_footer=0.75


;a test gausian dist plot to test tv command
;tv,bytscl(dist(500))

; start loop over all variables using unit 19

openr,19,"global_variables.dat"
while not eof(19) do begin
variable = ' ' & readf,19,variable,rmin,rmax

long_filename = "$input_data_directory" + "$short_filename"
iret=MOD08_Reader(long_filename,variable,data=array,FillValue=FillValue,unit=unit)




asize = size(array)

print,"  "
print,"  "
print," ******************** "
print,"  "


print,variable
help,array
print,'fill=',FillValue

print,"  "
print," ******************** "
print,"  "
print,"  "

;number of dimensions in array
dim=asize(0)

;size of 1st, 2nd, 3rd, 4th dimension
xdim=asize(1)
ydim=asize(2)
ndim=asize(3)
mdim=asize(4)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Begin Building Master Image
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

;create new master image
new=bytarr(420,210)

;set master image background color to white
new(0:419,0:209)=254


;;;;;;;;;;;;;;;;;;;;;;;;
; Build a Color Bar
;;;;;;;;;;;;;;;;;;;;;;;;

; build colorbar (0 to 253 only)
bar=fltarr(5,180)
colorbar=bytarr(5,180)
for i=0,179 do begin
realval=253.0*(float(i)/179.0)
bar(0:4,i)=realval
endfor
colorbar=byte(bar)
print,"barmin = ",min(colorbar),"barmax = ",max(colorbar)
print,"barbot = ",colorbar(0,0),"bartop = ",colorbar(0,179)




;;;;;;;;;;;;;;;;;;;;
; Process Image Data
;;;;;;;;;;;;;;;;;;;;

byte_array = bytarr(360,180)


if ( dim eq 2 ) then begin

flipped=rotate(array,7)

;check for all missing data values (no data)
if (MIN(flipped) eq FillValue and MAX(flipped) eq FillValue) then begin
minvalue = 0.0
maxvalue = 0.0
goto, JUMP1
endif

;compute min and max data values
if (MIN(flipped) eq FillValue) then begin
minvalue = MIN(flipped[WHERE(flipped ne FillValue)])
minvalue = rmin
endif

if (MIN(flipped) ne FillValue) then begin
minvalue = MIN(flipped)
minvalue = rmin
endif

maxvalue = MAX(flipped)
maxvalue = rmax

print, "datamin = ", minvalue, "datamax = ", maxvalue

;scale data from 0 to 253
byte_array=bytscl(flipped,min=minvalue,max=maxvalue,top=253)

;print the scaled byte min and max
print,"bytemin (0) = ",min(byte_array),"bytemax (253) = ",max(byte_array)

JUMP1:

;set fillvalue to 255 (black)
if (MIN(flipped) eq FillValue) then begin
print,'setting fillvalue to 255'
fill_array=bytarr(360,180)
fill_array(0:359,0:179)=255
byte_array[WHERE(flipped eq FillValue)]=fill_array[WHERE(flipped eq FillValue)]
endif


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Assign Lat-Lon Lines to White
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

;Lat-Lon lines assigned to White
llcolor=254
;
; Longitude Lines
;for i=1,11 do begin
;byte_array(i*30-1,0:179)=llcolor
;endfor
;
; Latitude Lines
;for i=1,5 do begin
;byte_array(0:359,i*30-1)=llcolor
;endfor





;;;;;;;;;;;;;;;;;;;;;;;;;;
; LATLON PROJECTION
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; define relative coordinates of LL and UR corners 
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;   XMap(02:361) = 360  
; XImage(00:419) = 420
; X   0.0=pixel(0)   1.0=pixel(419)
; X(LL)=2/419   X(UR)=361/419
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; YCoords 0.0=pixel0 1.0=pixel209 
;   YMap(14:193) = 180
; YImage(00:209) = 210
; Y(LL)=14/209   Y(UR)=193/209
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;if position is set to [0,0,1,1] it produces an image that is 420x210
;set position to make lat-lon 360x180
position=[0.004773269689737,0.06698564593301,0.8615751789976,0.9234449760766]
limit=[-90.0,-180.0,90.0,180.0]
Map_Set,0.0,0.0,0.0,POSITION=position,LIMIT=limit,/isotropic,/noborder


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Now Create LatLons and Coastlines On TvScreen Only
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;tv,byte_array
;Map_Continents, color=254, mlinethick=1
;Map_Grid, color=254, glinestyle=0, glinethick=linethick
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Now Create LatLons and Coastlines To Put In Final Array
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
templines=bytarr(360,180)
fillcolor = 0
templines(0:359,0:179)=fillcolor
tv,templines

 Map_Continents, color=254, mlinethick=1
;;;;; Map_Grid, color=254, glinestyle=0, glinethick=1, latdel=30, londel=30

print,'read tv screen'
templines = tvrd(80,14,360,180)
print,'done reading tv screen!'

print,'min templines=',min(templines),'max templines=',max(templines)

;print,'redisplay saved templines'
;tv,templines







;;;;;;;;;;;;;;;;;;;
; TEXT ANNOTATIONS
;;;;;;;;;;;;;;;;;;;

nnn=minvalue
dd1=minvalue+(maxvalue-minvalue)*1.0/4.0
dd2=minvalue+(maxvalue-minvalue)*2.0/4.0
dd3=minvalue+(maxvalue-minvalue)*3.0/4.0
xxx=maxvalue
print,nnn,dd1,dd2,dd3,xxx

fillcolor = 000

;Create Temporary Array for Text

temptext=bytarr(600,210)
temptext(0:599,0:209)=fillcolor

;Display image on screen
TV,temptext

textcolor = 255

;set vertical text baseline (0 is at bottom of image)
base=0

;Write Annotations

; sds name
XYOuts,005,base+198,variable,/Device,Color=textcolor,CharSize=sizenum_title

; date
;XYOuts,320,base+198,'$image_date',/Device,Color=textcolor,CharSize=sizenum_title
XYOuts,357,base+198,'$image_date',/Device,Color=textcolor,CharSize=sizenum

; color bar annotation
qxxx=strtrim(xxx,2)
qdd3=strtrim(dd3,2)
qdd2=strtrim(dd2,2)
qdd1=strtrim(dd1,2)
qnnn=strtrim(nnn,2)
XYOuts,372,base+187,qxxx,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+143,qdd3,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+99,qdd2,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+55,qdd1,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+11,qnnn,/Device,Color=textcolor,CharSize=sizenum

; data/platform
XYOuts,005,base+2,"$platform",/Device,Color=textcolor,CharSize=sizenum_footer

; hdf filename
XYOuts,100,base+2,"$short_filename",/Device,Color=textcolor,CharSize=sizenum_footer


; add asterisk for optimized (fixed) scale
XYOuts,340,base+3,'*',/Device,Color=textcolor,CharSize=sizenum


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Fix incorrect units specified for Atmospheric_Water_Vapor
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
if ( variable eq 'Atmospheric_Water_Vapor_Low_Mean' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_Low_Standard_Deviation' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_Low_Minimum' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_Low_Maximum' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_Low_QA_Mean' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_Low_QA_Standard_Deviation' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_High_Mean' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_High_Standard_Deviation' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_High_Minimum' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_High_Maximum' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_High_QA_Mean' ) then begin
unit='cm'
endif
if ( variable eq 'Atmospheric_Water_Vapor_High_QA_Standard_Deviation' ) then begin
unit='cm'
endif
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; End Fix
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;



; units
qunit=strtrim(unit,2)
print,'qunit=',qunit
XYOuts,345,base+2,qunit,/Device,Color=textcolor,CharSize=sizenum




;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; NOTE: If you load the whole color image (and bar) and
; then use tvrd to grab the annotation ...
; it screws up the user specified color table!
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;










;Read the screen to get text annotation
print,'read tv screen'
temptext = tvrd(0,0,600,210)
print,'done reading tv screen!'

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Place Text Into Master Image & Write to GIF
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

mintmp = min(temptext)
maxtmp = max(temptext)

print,'min temptext=',mintmp,'max temptext=',maxtmp

; SET BACKGROUND TO WHITE!!!
temptext[where(temptext eq mintmp)]=254

; SET TEXT TO BLACK!!
temptext[where(temptext eq maxtmp)]=255

; note new image is 420,210
 
print,'test1'
;ADD ANNOTATION
new(0,0)=temptext(0:419,0:209)

print,'test2'
;ADD IMAGE
new(2,14)=byte_array

print, 'test3'
;ADD COLOR BAR
new(365,14)=colorbar

;ADD LATLON + COASTLINES
temp2lines=bytarr(420,210)
temp2lines(0:419,0:209)=0
temp2lines(2,14)=templines

print, 'max temp2lines=', max(temp2lines)


; SET LATLON and COAST TO WHITE!!
maxlines = max(temp2lines)
new[where(temp2lines eq maxlines)]=254


;WRITE IMAGE
print,'write image'


; write_gif,variable+".gif",new,red,grn,blu

image_filename="$output_data_directory"+variable+".png"

write_image,image_filename,'PNG',new,red,grn,blu


endif















if ( dim eq 3 ) then begin

print,"third dimension = ",ndim

for j=0,ndim-1 do begin

flipped=rotate(array[*,*,j],7)

;check for all missing data values (no data)
if (MIN(flipped) eq FillValue and MAX(flipped) eq FillValue) then begin
minvalue = 0.0
maxvalue = 0.0
goto, JUMP2
endif

;compute min and max data values
if (MIN(flipped) eq FillValue) then begin
minvalue = MIN(flipped[WHERE(flipped ne FillValue)])
minvalue = rmin
endif

if (MIN(flipped) ne FillValue) then begin
minvalue = MIN(flipped)
minvalue = rmin
endif

maxvalue = MAX(flipped)
maxvalue = rmax

print, "datamin = ", minvalue, "datamax = ", maxvalue

;scale data from 0 to 253

byte_array=bytscl(flipped,min=minvalue,max=maxvalue,top=253)

;print the scaled byte min and max
print,"bytemin (0) = ",min(byte_array),"bytemax (253) = ",max(byte_array)

JUMP2:

;set fillvalue to 255 (black)
if (MIN(flipped) eq FillValue) then begin
print,'setting fillvalue to 255'
fill_array=bytarr(360,180)
fill_array(0:359,0:179)=255
byte_array[WHERE(flipped eq FillValue)]=fill_array[WHERE(flipped eq FillValue)]
endif

llcolor=254

;
; Longitude Lines
;for i=1,11 do begin
;byte_array(i*30-1,0:179)=llcolor
;endfor
;
; Latitude Lines
;for i=1,5 do begin
;byte_array(0:359,i*30-1)=llcolor
;endfor






;;;;;;;;;;;;;;;;;;;;;;;;;;
; LATLON PROJECTION
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; define relative coordinates of LL and UR corners 
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;   XMap(02:361) = 360  
; XImage(00:419) = 420
; X   0.0=pixel(0)   1.0=pixel(419)
; X(LL)=2/419   X(UR)=361/419
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; YCoords 0.0=pixel0 1.0=pixel209 
;   YMap(14:193) = 180
; YImage(00:209) = 210
; Y(LL)=14/209   Y(UR)=193/209
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;if position is set to [0,0,1,1] it produces an image that is 420x210
;set position to make lat-lon 360x180
position=[0.004773269689737,0.06698564593301,0.8615751789976,0.9234449760766]
limit=[-90.0,-180.0,90.0,180.0]
Map_Set,0.0,0.0,0.0,POSITION=position,LIMIT=limit,/isotropic,/noborder


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Now Create LatLons and Coastlines On TvScreen Only
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;tv,byte_array
;Map_Continents, color=254, mlinethick=1
;Map_Grid, color=254, glinestyle=0, glinethick=linethick
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Now Create LatLons and Coastlines To Put In Final Array
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
templines=bytarr(360,180)
fillcolor = 0
templines(0:359,0:179)=fillcolor
tv,templines

 Map_Continents, color=254, mlinethick=1
;;;;;;;;;; Map_Grid, color=254, glinestyle=0, glinethick=1, latdel=30, londel=30

print,'read tv screen'
templines = tvrd(80,14,360,180)
print,'done reading tv screen!'

print,'min templines=',min(templines),'max templines=',max(templines)

;print,'redisplay saved templines'
;tv,templines







;;;;;;;;;;;;;;;;;;;
; TEXT ANNOTATIONS
;;;;;;;;;;;;;;;;;;;

nnn=minvalue
dd1=minvalue+(maxvalue-minvalue)*1.0/4.0
dd2=minvalue+(maxvalue-minvalue)*2.0/4.0
dd3=minvalue+(maxvalue-minvalue)*3.0/4.0
xxx=maxvalue
print,nnn,dd1,dd2,dd3,xxx

fillcolor = 000

;Display image on screen
temptext=bytarr(600,210)
temptext(0:599,0:209)=fillcolor
TV,temptext

textcolor = 255

;set vertical text baseline (0 is at bottom of image)
base=0

;Write Annotations
; sds name
XYOuts,005,base+198,variable,/Device,Color=textcolor,CharSize=sizenum_title

; date
;XYOuts,320,base+198,'$image_date',/Device,Color=textcolor,CharSize=sizenum_title
XYOuts,357,base+198,'$image_date',/Device,Color=textcolor,CharSize=sizenum

; color bar annotation
qxxx=strtrim(xxx,2)
qdd3=strtrim(dd3,2)
qdd2=strtrim(dd2,2)
qdd1=strtrim(dd1,2)
qnnn=strtrim(nnn,2)
XYOuts,372,base+187,qxxx,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+143,qdd3,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+99,qdd2,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+55,qdd1,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+11,qnnn,/Device,Color=textcolor,CharSize=sizenum

XYOuts,005,base+2,"$platform",/Device,Color=textcolor,CharSize=sizenum_footer

XYOuts,100,base+2,"$short_filename",/Device,Color=textcolor,CharSize=sizenum_footer

XYOuts,340,base+3,'*',/Device,Color=textcolor,CharSize=sizenum

qunit=strtrim(unit,2)
XYOuts,345,base+2,qunit,/Device,Color=textcolor,Charsize=sizenum

;Read the screen to get annotation
temptext = tvrd(0,0,600,210)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Place Text Into Master Image & Write to GIF
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

mintmp = min(temptext)
maxtmp = max(temptext)

print,'min temptext=',mintmp,'max temptext=',maxtmp

; SET BACKGROUND TO WHITE!!!
temptext[where(temptext eq mintmp)]=254

; SET TEXT TO BLACK!!
temptext[where(temptext eq maxtmp)]=255

; note new image is 420,210

;ADD ANNOTATION
new(0,0)=temptext(0:419,0:209)

;ADD IMAGE
new(2,14)=byte_array

;ADD COLOR BAR
new(365,14)=colorbar

;ADD LATLON + COASTLINES
temp2lines=bytarr(420,210)
temp2lines(0:419,0:209)=0
temp2lines(2,14)=templines


; SET LATLON and COAST TO WHITE!!
maxlines = max(temp2lines)
new[where(temp2lines eq maxlines)]=254


image_filename="$output_data_directory"+variable+".cata-"+strtrim(string(j),2)+".png"
write_image,image_filename,'PNG',new,red,grn,blu


endfor
endif



if ( dim eq 4 ) then begin
print,"third dimension = ",ndim
print,"fourth dimension = ",mdim
for j=0,ndim-1 do begin
for k=0,mdim-1 do begin

flipped=rotate(array[*,*,j,k],7)

;check for all missing data values (no data)
if (MIN(flipped) eq FillValue and MAX(flipped) eq FillValue) then begin
minvalue = 0.0
maxvalue = 0.0
goto, JUMP3
endif

;compute min and max data values
if (MIN(flipped) eq FillValue) then begin
minvalue = MIN(flipped[WHERE(flipped ne FillValue)])
minvalue = rmin
endif

if (MIN(flipped) ne FillValue) then begin
minvalue = MIN(flipped)
minvalue = rmin
endif

maxvalue = MAX(flipped)
maxvalue = rmax

print, "datamin = ", minvalue, "datamax = ", maxvalue

;scale data from 0 to 253
byte_array=bytscl(flipped,min=minvalue,max=maxvalue,top=253)

;print the scaled byte min and max
print,"bytemin (0) = ",min(byte_array),"bytemax (253) = ",max(byte_array)

JUMP3:

;set fillvalue to 255 (black)
if (MIN(flipped) eq FillValue) then begin
print,'setting fillvalue to 255'
fill_array=bytarr(360,180)
fill_array(0:359,0:179)=255
byte_array[WHERE(flipped eq FillValue)]=fill_array[WHERE(flipped eq FillValue)]
endif

llcolor=254

; Longitude Lines
;for i=1,11 do begin
;byte_array(i*30-1,0:179)=llcolor
;endfor

; Latitude Lines
;for i=1,5 do begin
;byte_array(0:359,i*30-1)=llcolor
;endfor









;;;;;;;;;;;;;;;;;;;;;;;;;;
; LATLON PROJECTION
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; define relative coordinates of LL and UR corners 
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;   XMap(02:361) = 360  
; XImage(00:419) = 420
; X   0.0=pixel(0)   1.0=pixel(419)
; X(LL)=2/419   X(UR)=361/419
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; YCoords 0.0=pixel0 1.0=pixel209 
;   YMap(14:193) = 180
; YImage(00:209) = 210
; Y(LL)=14/209   Y(UR)=193/209
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;if position is set to [0,0,1,1] it produces an image that is 420x210
;set position to make lat-lon 360x180
position=[0.004773269689737,0.06698564593301,0.8615751789976,0.9234449760766]
limit=[-90.0,-180.0,90.0,180.0]
Map_Set,0.0,0.0,0.0,POSITION=position,LIMIT=limit,/isotropic,/noborder


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Now Create LatLons and Coastlines On TvScreen Only
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
;tv,byte_array
;Map_Continents, color=254, mlinethick=1
;Map_Grid, color=254, glinestyle=0, glinethick=linethick
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Now Create LatLons and Coastlines To Put In Final Array
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
templines=bytarr(360,180)
fillcolor = 0
templines(0:359,0:179)=fillcolor
tv,templines

 Map_Continents, color=254, mlinethick=1
;;;;;;;;;; Map_Grid, color=254, glinestyle=0, glinethick=1, latdel=30, londel=30

print,'read tv screen'
templines = tvrd(80,14,360,180)
print,'done reading tv screen!'

print,'min templines=',min(templines),'max templines=',max(templines)

;print,'redisplay saved templines'
;tv,templines









;;;;;;;;;;;;;;;;;;;
; TEXT ANNOTATIONS
;;;;;;;;;;;;;;;;;;;

nnn=minvalue
dd1=minvalue+(maxvalue-minvalue)*1.0/4.0
dd2=minvalue+(maxvalue-minvalue)*2.0/4.0
dd3=minvalue+(maxvalue-minvalue)*3.0/4.0
xxx=maxvalue
print,nnn,dd1,dd2,dd3,xxx

fillcolor = 000

;Display image on screen
temptext=bytarr(600,210)
temptext(0:599,0:209)=fillcolor
TV,temptext

textcolor = 255

;set vertical text baseline (0 is at bottom of image)
base=0

;Write Annotations
; sds name
XYOuts,005,base+198,variable,/Device,Color=textcolor,CharSize=sizenum_title

; date
;XYOuts,320,base+198,'$image_date',/Device,Color=textcolor,CharSize=sizenum_title
XYOuts,357,base+198,'$image_date',/Device,Color=textcolor,CharSize=sizenum

; color bar annotation
qxxx=strtrim(xxx,2)
qdd3=strtrim(dd3,2)
qdd2=strtrim(dd2,2)
qdd1=strtrim(dd1,2)
qnnn=strtrim(nnn,2)
XYOuts,372,base+187,qxxx,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+143,qdd3,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+99,qdd2,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+55,qdd1,/Device,Color=textcolor,CharSize=sizenum
XYOuts,372,base+11,qnnn,/Device,Color=textcolor,CharSize=sizenum

; data/platform
XYOuts,005,base+2,"$platform",/Device,Color=textcolor,CharSize=sizenum_footer

; hdf filename
XYOuts,100,base+2,"$short_filename",/Device,Color=textcolor,CharSize=sizenum_footer

XYOuts,340,base+3,'*',/Device,Color=textcolor,CharSize=sizenum

; units
qunit=strtrim(unit,2)
XYOuts,345,base+2,qunit,/Device,Color=textcolor,CharSize=sizenum

;Read the screen to get annotation
temptext = tvrd(0,0,600,210)

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Place Text Into Master Image & Write to GIF
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

mintmp = min(temptext)
maxtmp = max(temptext)

print,'min temptext=',mintmp,'max temptext=',maxtmp

; SET BACKGROUND TO WHITE!!!
temptext[where(temptext eq mintmp)]=254

; SET TEXT TO BLACK!!
temptext[where(temptext eq maxtmp)]=255

; note new image is 420,210

;ADD ANNOTATION
new(0,0)=temptext(0:419,0:209)

;ADD IMAGE
new(2,14)=byte_array

;ADD COLOR BAR
new(365,14)=colorbar

;ADD LATLON + COASTLINES
temp2lines=bytarr(420,210)
temp2lines(0:419,0:209)=0
temp2lines(2,14)=templines


; SET LATLON and COAST TO WHITE!!
maxlines = max(temp2lines)
new[where(temp2lines eq maxlines)]=254


image_filename="$output_data_directory"+variable+".catb-"+strtrim(string(j),2)+".cata-"+strtrim(string(k),2)+".png"

write_image,image_filename,'PNG',new,red,grn,blu


endfor
endfor
endif

endwhile
close,19

end

eof



# setenv IDL_DEVICE "PS"
timex idl <<eof
.rnew mod08_reader.pro
.rnew make_gif.pro
make_gif
exit
eof

#test#  xv -expand 2 -wait 1 *.gif