source: trunk/mesoscale/LMD_MM_MARS/SRC/ARWpost/idl/goodscale.pro @ 13

Last change on this file since 13 was 11, checked in by aslmd, 14 years ago

spiga@svn-planeto:ajoute le modele meso-echelle martien

File size: 1.3 KB
Line 
1;;Fixe l'echelle pour une tab2diable de type (x,y,z,t) ou (x,y,t) ou (z,t)
2        ;;Possibilite de restreindre a une zone et des heures donnees
3
4pro goodscale, $
5    tab2d, $
6    mini, $
7    maxi, $
8    coeff=coeff_sigma
9       
10print,'Yeah !'
11
12;
13; input: the 2D array which is ready to be plotted
14;
15
16
17;if (n_elements(coeff_sigma) eq 0) then coeff_sigma=3.3 ;; vincent
18if (n_elements(coeff_sigma) eq 0) then coeff_sigma=2.5 ;; aymeric
19
20;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;     
21       
22;;calcul de la moyenne
23        N=n_elements(tab2d)
24        moy=(1/float(N))*total(tab2d)   
25
26;;calcul de la tab2diance
27        moyarray = make_array(N, 1, /integer, value = 1)       
28        varia=(1/float(N-1))*total( ( tab2d - (moyarray * moy) )^2 )
29       
30;;calcul de sigma
31        sigma = sqrt( varia )
32       
33;;echelle
34        mini = moy - float(coeff_sigma) * sigma & stratmin = 'sigma'
35        maxi = moy + float(coeff_sigma) * sigma & stratmax = 'sigma'
36
37;;au cas ou
38if (mini lt min(tab2d)) then begin
39        mini=min(tab2d)
40        stratmin = 'real min'
41endif
42if (maxi gt max(tab2d)) then begin
43        maxi=max(tab2d)
44        stratmax = 'real max'
45endif
46       
47       
48;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;     
49;;affichage a l'ecran :
50print,'min :'+STRTRIM(mini,2)+' '+stratmin     
51print,'max :'+STRTRIM(maxi,2)+' '+stratmax
52
53
54end
Note: See TracBrowser for help on using the repository browser.