;This is the batch file  for simap.pro
;Intended for stitching of L1B 1km direct broadcast granules,
;or any other L1B 1km.
;Can be reused for all other resolutions or L2. Just give the right parameters
;names and channels.



.full_reset_session
.compile simap.pro

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


;====Input and Output path for source and output data== 
inpath="/var/scratch/"
outpath="/var/scratch"
;======================================================



;===This is your file selection:======
hdfname="filelist"
;=====================================




;===This is your parameters and channel selections:===
   sel[0].name="EV_250_Aggr1km_RefSB"
   sel[0].chans=['1']

   sel[1].name="EV_500_Aggr1km_RefSB"
   sel[1].chans=['4','3']

;======================================================




;====================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=3	;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=0	;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=0	;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=15.
	latmax=50.
	lonmin=-95.	;westboundcoordinate
	lonmax=-70.	;eastboundcoordinate

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


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

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

	;IMAGE INTENSITY
	;default ranges for RGB planes:
	intens=FLTARR(2,3)

	intens[*,0]=[5,90]	;red
	intens[*,1]=[15,90]	;green
	intens[*,2]=[40,110]	;blue


	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=800  ;Request in pxl on the size of map 
	NYM=1000 ;to be produced. 


	ctab=0		;IDL color table

	backgnd=200	;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