# HDFLook script: California Case, Terra MODIS 2005273.1920
# MOD06 Effective Radius

# Startup
verbose
clear_data

# Input file
set_input_hdf_file          "MOD06_L2.A2005273.1920.005.2006229083418.hdf"
select_SDS                  SDSname="Cloud_Effective_Radius"
set_SDS_enhance             min=5 max=20
set_color_scale_name        "palette_brewer.txt"

# Projection
set_projection_to_geometry  ProjectionTo=LINEAR PIXELSIZEXTO=500 PIXELSIZEYTO=500
set_image_background        Color=Black

# Labels
add_text                    x=1 y=0 ForegroundColor="Black"  BackgroundColor="White" \
                            Border=0 FontName="Verdana.ttf" SizeFactor=2 Rotation=0 \
                            Text=" Terra MODIS 2005273 19:20 UTC Effective Radius "
#add_text                    x=1 y=19 ForegroundColor="Black"  BackgroundColor="White" \
#                            Border=0 FontName="Verdana.ttf" SizeFactor=1.5 Rotation=0 \
#                            Text=" Effective Radius (microns) "
#add_scale_to_image          x=40 y=100 width=800 height=30 Rotation=90.0 \
#                            TickColor=White ForegroundColor="Black" BackgroundColor="White" \
#                            TextFormat="%.1f" FontName="Verdana.ttf" SizeFactor=1.2 \
#                            UpLeftSide=Yes NoTicks=No

# Grids
add_maps                    Coastlines=Yes Boundaries=Yes CoastlineColor="Yellow" BoundaryColor="Yellow"
add_grids                   Color=White LatitudeMin=-90  LatitudeMax=90 LatitudeStep=2 \
                            LongitudeMin=-180 LongitudeMax=180 LongitudeStep=2
add_latitude_labels         ForegroundColor="Black" BackgroundColor="Grey0" \
                            Border=0 FontName="Verdana.ttf" \
                            SizeFactor=1.5 Rotation=0. Format="%.0f" \
                            LatitudeMin=-90  LatitudeMax=90 LatitudeStep=2
add_longitude_labels        ForegroundColor="Black" BackgroundColor="White" \
                            Border=0 FontName="Verdana.ttf" \
                            SizeFactor=1.5 Rotation=90. Format="%.0f" \
                            LongitudeMin=-180 LongitudeMax=180 LongitudeStep=2

# Subset output file
create_MODIS_SDS_GEOTIFF_image FileName="MOD06_L2.A2005273.1920.reff.subset.tif" Scaling="y=ax+b" \
  LatitudeMin=30.0 LatitudeMax=37.0 LongitudeMin=-128.5 LongitudeMax=-121.5