1 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
2 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
3 | filename = '/d5/aslmd/LMD_MM_MARS_SIMUS/OM/OM6_TI85/geo_em.d01.nc' |
---|
4 | save_ps = 'olympus_slope' |
---|
5 | field1 = 'THERMAL_INERTIA' ;; omettre trace pentes |
---|
6 | save_ps = 'olympus_ti' |
---|
7 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
8 | format = '(I0)' |
---|
9 | colors = 32 ;64 ;128 ;32 |
---|
10 | pal = 6 ;39 ;4 ;33 ;4 ;33 ;4 ;22 ;33 ;39 ;6 ;11 ;0 ;11 ;6 ;0 ;6 ;39 ;19 |
---|
11 | pal = 4 |
---|
12 | title_user = 'Slope angle (deg)' |
---|
13 | title_user = 'Apparent nighttime thermal inertia (J m!U-2!N s!U-0.5!N K!U-1!N)' |
---|
14 | title_axis = ['East longitude','North latitude'] |
---|
15 | poscb=0.80 |
---|
16 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
17 | ;what_I_plot = what_I_plot / 1000. |
---|
18 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
19 | ;;; FILL LIMITS ;;; |
---|
20 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
21 | minfield_init= 0. |
---|
22 | maxfield_init= 25. |
---|
23 | minfield_init= 50. |
---|
24 | maxfield_init= 250. |
---|
25 | missing_value=1.e30 |
---|
26 | lim_max = maxfield_init & w=where((what_I_plot ge lim_max) and (what_I_plot le 1e9)) & if (w[0] ne -1) then what_I_plot[w]=lim_max |
---|
27 | lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min |
---|
28 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
29 | ;;; WINDS ;;; |
---|
30 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
31 | windex=20. ;; default: 20. |
---|
32 | stride=3. ;; default: 5. |
---|
33 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
34 | ;;; CONTOUR ;;; |
---|
35 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
36 | overcontour=overcontour/1000. |
---|
37 | lev=-10. + 2.*findgen(20) |
---|
38 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
39 | ;;; AXIS ;;; |
---|
40 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
41 | intervalx=2.0 |
---|
42 | intervaly=1.0 |
---|
43 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
44 | ;;; MAP LIMITS ;;; |
---|
45 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
---|
46 | windowx = [-146.,-126.] |
---|
47 | windowy = [11.,27.] |
---|