| 1 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 2 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 3 | filename = '/d5/aslmd/WRF_Mars/EN_COURS/OMG_cstTI/geo_em.d03.nc' |
|---|
| 4 | save_ps = 'meridiani_ti' |
|---|
| 5 | field1 = 'THERMAL_INERTIA' ;; omettre trace vitesse horizontale |
|---|
| 6 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 7 | format = '(I0)' |
|---|
| 8 | colors = 32 |
|---|
| 9 | pal = 4 ;33 ;4 ;22 ;33 ;39 ;6 ;11 ;0 ;11 ;6 ;0 ;6 ;39 ;19 |
|---|
| 10 | title_user = 'Apparent nighttime thermal inertia (J m!U-2!N s!U-0.5!N K!U-1!N)' |
|---|
| 11 | title_axis = ['East longitude','North latitude'] |
|---|
| 12 | poscb=0.96 |
|---|
| 13 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 14 | ;;; FILL LIMITS ;;; |
|---|
| 15 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 16 | ;minfield_init= 10. |
|---|
| 17 | ;maxfield_init= 110. |
|---|
| 18 | minfield_init= 30. ;50. |
|---|
| 19 | maxfield_init= 230. ;250. |
|---|
| 20 | minfield_init= 50. |
|---|
| 21 | maxfield_init= 250. |
|---|
| 22 | ;minfield_init=1.4 |
|---|
| 23 | ;maxfield_init=2.4 |
|---|
| 24 | ;what_I_plot = alog(what_I_plot)/alog(10) |
|---|
| 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. + 0.2*findgen(100) |
|---|
| 38 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 39 | ;;; AXIS ;;; |
|---|
| 40 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 41 | intervalx=0.5 |
|---|
| 42 | intervaly=0.5 |
|---|
| 43 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 44 | ;;; MAP LIMITS ;;; |
|---|
| 45 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 46 | windowx=[14.0,16.0] |
|---|
| 47 | windowy=[10.0,12.0] |
|---|
| 48 | |
|---|
| 49 | |
|---|
| 50 | filename = '/donnees/aslmd/MODELES/LMD_MM_MARS/g95_32_single/WPS/geo_em.d01.nc' |
|---|
| 51 | filename = './geo_em.d01.nc' |
|---|
| 52 | ;filename = '/donnees/aslmd/MODELES/LMD_MM_MARS/mpi_64/WPS/geo_em.d01.nc' |
|---|
| 53 | filename = '/d5/aslmd/LMD_MM_MARS_SIMUS/OM/OM6_TI85/geo_em.d01.nc' |
|---|
| 54 | save_ps = 'OM6_TI85_ti' |
|---|
| 55 | lev = -10. + 2.*findgen(20) |
|---|
| 56 | intervalx = 2. |
|---|
| 57 | intervaly = 1. |
|---|
| 58 | windowx = [-146.,-126.] |
|---|
| 59 | windowy = [10.,28.] |
|---|
| 60 | windowx = [-146.,-128.] |
|---|
| 61 | windowy = [12.,26.] |
|---|
| 62 | windowx = [-146.,-126.] |
|---|
| 63 | windowy = [11.,27.] |
|---|
| 64 | |
|---|
| 65 | |
|---|
| 66 | title='Slope angle (deg)' |
|---|
| 67 | format='(I0)' |
|---|
| 68 | pal=16 |
|---|
| 69 | |
|---|
