| 1 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 2 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 3 | ;ntime = 7 ;; le vrai indice IDL |
|---|
| 4 | nlevel = 0 |
|---|
| 5 | filename = '/donnees/aslmd/MODELES/LMD_MM_MARS/OUTILS_CONVERSION/wrfout_d03_2024-06-20_06:00:00_z' |
|---|
| 6 | save_ps = 'nest_polar_winds' |
|---|
| 7 | filename = '/donnees/aslmd/MODELES/LMD_MM_MARS/OUTILS_CONVERSION/wrfout_d01_2024-06-20_06:00:00_z' |
|---|
| 8 | save_ps = 'single_polar_winds' |
|---|
| 9 | ;filename = '/tmp7/aslmd/one_more_nest/test/wrfout_d02_2024-06-19_06:00:00_zabg' |
|---|
| 10 | ;save_ps = '2nest_polar_winds' |
|---|
| 11 | |
|---|
| 12 | ;coord2d = 'true' |
|---|
| 13 | field1 = 'HGT' ;; omettre trace vitesse horizontale |
|---|
| 14 | ;field2 = 'XLAT' ;; contour |
|---|
| 15 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 16 | format = '(I0)' ;'(F5.2)' |
|---|
| 17 | colors = 32 |
|---|
| 18 | pal = 19 ;22 ;33 ;39 ;6 ;11 ;0 ;11 ;6 ;0 ;6 ;39 ;19 |
|---|
| 19 | title_user = 'Topography (m) + Horizontal wind vectors 6km ABG' |
|---|
| 20 | title_axis = ['grid points','grid points'] |
|---|
| 21 | subtitle_user = 'LMD Mars Mesoscale' |
|---|
| 22 | ;subtitle_user = subtitle_user + ' / UTC = 01:00pm' |
|---|
| 23 | ;subtitle_user = subtitle_user + ' / Ls = 61.5!Uo!N' |
|---|
| 24 | ;subtitle_user = subtitle_user + ' / dx = 10km' ;25km' |
|---|
| 25 | ;subtitle_user = subtitle_user + ' / single domain' |
|---|
| 26 | poscb=0.85 |
|---|
| 27 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 28 | ;;; FILL LIMITS ;;; |
|---|
| 29 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 30 | minfield_init= -3000. |
|---|
| 31 | maxfield_init= +3000. |
|---|
| 32 | missing_value=1.e30 |
|---|
| 33 | 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 |
|---|
| 34 | lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min |
|---|
| 35 | ;lim_blank = 0.2 & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value |
|---|
| 36 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 37 | ;;; WINDS ;;; |
|---|
| 38 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 39 | windex=25. ;; default: 20. |
|---|
| 40 | stride=5. ;; default: 5. |
|---|
| 41 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 42 | ;;; CONTOUR ;;; |
|---|
| 43 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 44 | ;overcontour=overcontour-mean(overcontour) |
|---|
| 45 | overcontour=overcontour/1000. |
|---|
| 46 | lev=-10. + 0.2*findgen(100) |
|---|
| 47 | lev=-10. + 0.5*findgen(50) |
|---|
| 48 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 49 | ;;; AXIS ;;; |
|---|
| 50 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 51 | ;isotropic='false' ;; default: true |
|---|
| 52 | intervalx=20. |
|---|
| 53 | intervaly=20. |
|---|
| 54 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 55 | ;;; MAP LIMITS ;;; |
|---|
| 56 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 57 | windowx=[0.,120.] |
|---|
| 58 | windowy=[0.,120.] |
|---|
| 59 | windowx=[40.,140.] |
|---|
| 60 | windowy=[40.,140.] |
|---|
| 61 | |
|---|
