;This is the batch file  for simap.pro
;The sample is for Atmosphere L2 gallery, showing how to 
;map one parameter that is a three-dimensional array.  



.full_reset_session
.compile simap

sel={name:"",chans:STRARR(38),numch:0}		
sel=REPLICATE(sel,20)


inpath="/ftp/data/modis/samples/atmosphere/L2/"
outpath="/ftp/data/3/"



hdfname="MOD06_L2.A2000110.0220.002.2000207233006.hdf"

sel[0].name=['Brightness_Temperature']
sel[0].chans=['1','2','3','4','5','6','7']




;====================FLAGS DEFINITIONS============================

intrpl=0	;Set to 1 activate filling
		;up of holes in the images. This may slow 
		;the process, but will clean the images.
		;Set to 0 to deactivate.


extgeo=0	;set to 1 to use EXTERNAL geolocations
		;set to 0 to use INTERNAL geolocations
		
igeo=1		;set to 1 to interpolate geolocations
		;It must be 1 when using internal geolocations
		;and L1B MODIS files

sdsprint=0	;Flag to list SDS
		;0=no list
		;1=short list, SDS names only, no attributes
		;2=long list, SDS names, attributes and their
		;values


True_col=1	;Flag to generate image after all procedures
		;are completed.
		;Set to 3 to save RGB image composite from the
		;first 3 channels selected.
		;Set to 1 to save grey/color scale separate
		;images of all selected channels.

flcnt=0		;set to 1 to fill continents
		;set to 0 not to fill continents

automap=1	;Set to  1 to automaticaly
		;set map limits. 
		; Set to 0 to use the default
		; map limits. 	
		

autointen=0	;Set to 1 to automatically set image intensity.
		;Set to 0 to use default image intensity
		;ranges "intens"
		

zlog=0		;Set to different than 0 to show image
		;in log scale

gallery=1	;Flag to generate galleries of image files.
		;0=no gallery, will do global map
		;1=gif or jpeg gallery of gray/color scale,
		;of true color, RGB, image files. If 
		;True_col=1 then grey/color scale.
		;If True_col=3 then RGB
	

debug=-1	;debug mode: -1 will print to monitor
		;             4 will print to log file
;============================================================




;============PARAMETERS DEFINITIONS===========================

	;---------SUBSET-------
	xstart=""	;leave as NULL string for full size
	xend=""		;or give numerical value

	ystart=""
	yend=""


	;---------Default map settings-----------------
	latmin=21.
	latmax=31.
	lonmin=-118.	;westboundcoordinate
	lonmax=-106.	;eastboundcoordinate

	map_limit=[latmin,lonmin,latmax,lonmax]
	center_lat=0;latmin+(latmax-latmin)/2.
        center_lon=atan((sin(lonmin*!dtor)+sin(lonmax*!dtor)),$
                        (cos(lonmin*!dtor)+cos(lonmax*!dtor)))*!radeg


	mappos=[.05,.08,.95,.97]	;map position and size,
					; relative to the graphic window

	proj=1	;1 - cylindrical
		;2 - stereo
	;--------------------------------------------

	;IMAGE INTENSITY
	;default ranges for data planes:
	intens=FLTARR(2,20)
	intens[0,*]=150.	;Brightness_Temp
	intens[1,*]=300.

	slice=3	;For 3D arrays, this is the dimension
		;to fix.  In IDL, dimensions are
		;e.g. array[x,y,z]. Set slice=1 or 2 or 3 to fix
		; x or y or z dimension, correspondingly.
		;For most SDS's slice=3, i.e. they are stacked
		;over the z-dimension. However, for example MODIS 
		;Atmosphere L2 Quality_Assurance needs slice=1 !!!
		;I.e. it is stacked over the x-dimension


	NXM=700  ;Request in pxl on the size of map 
	NYM=700 ;to be produced. 


	ctab=33		;IDL color table

	backgnd=255	;background inside map
			;and ouside the granule(s)

	;this will be the size of the graphic window to
	;hold the map:
	NX=NXM/(mappos[2]-mappos[0])
	NY=NYM/(mappos[3]-mappos[1])
;================================================================





simap,hdfname