;This is the batch file  for simap.pro
;The sample is for Ocean L2 Group 1 gallery



.full_reset_session
.compile simap

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


inpath="/ftp/data/modis/samples/oceans/L2/"
outpath="/var/scratch"



hdfname="MODOCL2A.A2000338.1810.002.2001033042045.hdf"



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

intrpl=1	;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=1	;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=1		;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=1		;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=lonmin+ (lonmax-lonmin)/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 data planes:
	intens=FLTARR(2,20)

	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=800 ;to be produced. 


	ctab=33		;IDL color table

	backgnd=70	;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])
;================================================================


sel[0:4].name=['chlor_fluor_ht','chlor_fluor_base','chlor_fluor_effic',$
	'phycoeryth_conc','phycou_conc']
sel[0:4].chans=['1']

zlog=0
intens[*,0]=[0.,0.5]	;chlor_fluor_ht
intens[*,1]=[0.,5.]	;chlor_fluor_base
intens[*,2]=[0,0.4]	;chlor_fluor_effic
intens[*,3]=[0,1]	;phycoeryth_conc
intens[*,4]=[0,0.1]	;phycou_conc


simap,hdfname




sel[0:7].name=['CZCS_pigment','chlor_MODIS','pigment_c1_total',$
	'susp_solids_conc','cocco_pigmnt_conc','calcite_conc','K_490','cocco_conc_detach']
sel[0:7].chans=['1']

zlog=1
intens[*,0]=[.1,10.] 	;CZCS
intens[*,1]=[.1,10.]	;chlor_MODIS
intens[*,2]=[.1,10.]	;pigment_c1
intens[*,3]=[1.,1000.]	;susp_solids
intens[*,4]=[.1,10]	;cocco_pigmnt
intens[*,5]=[.1,10]	;calcite_conc
intens[*,6]=[0.01,10]	;K_490
intens[*,7]=[10.^9,10.^11];cocco_conc_detach


simap,hdfname