| 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.] |
|---|