source: trunk/MESOSCALE_DEV/PLOT/RESERVE/goodscale.pro @ 273

;;Fixe l'echelle pour une tab2diable de type (x,y,z,t) ou (x,y,t) ou (z,t)
        ;;Possibilite de restreindre a une zone et des heures donnees

pro goodscale, $
    tab2d, $
    mini, $
    maxi, $
    coeff=coeff_sigma
        
print,'Yeah !'

;
; input: the 2D array which is ready to be plotted
;


;if (n_elements(coeff_sigma) eq 0) then coeff_sigma=3.3 ;; vincent
if (n_elements(coeff_sigma) eq 0) then coeff_sigma=2.5 ;; aymeric

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;     
        
;;calcul de la moyenne
        N=n_elements(tab2d)
        moy=(1/float(N))*total(tab2d)   

;;calcul de la tab2diance
        moyarray = make_array(N, 1, /integer, value = 1)        
        varia=(1/float(N-1))*total( ( tab2d - (moyarray * moy) )^2 )
        
;;calcul de sigma
        sigma = sqrt( varia )
        
;;echelle
        mini = moy - float(coeff_sigma) * sigma & stratmin = 'sigma'
        maxi = moy + float(coeff_sigma) * sigma & stratmax = 'sigma'

;;au cas ou 
if (mini lt min(tab2d)) then begin
        mini=min(tab2d)
        stratmin = 'real min'
endif
if (maxi gt max(tab2d)) then begin
        maxi=max(tab2d)
        stratmax = 'real max'
endif
        
        
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;     
;;affichage a l'ecran :
print,'min :'+STRTRIM(mini,2)+' '+stratmin      
print,'max :'+STRTRIM(maxi,2)+' '+stratmax


end