| 1 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 2 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 3 | ntime = 7 ;; le vrai indice IDL |
|---|
| 4 | nlevel = 2 |
|---|
| 5 | filename = '/tmp15/aslmd/polar_61pts/wrfout_d01_2024-03-04_06:00:00_zabg' |
|---|
| 6 | ;filename = '/tmp15/aslmd/polar_61pts/wrfout_d01_2024-03-02_06:00:00_zabg' |
|---|
| 7 | save_ps = 'mesoscale_polar_winds_friction' |
|---|
| 8 | ;coord2d = 'true' |
|---|
| 9 | field1 = 'USTM' ;; omettre trace vitesse horizontale |
|---|
| 10 | ;field1 = 'HGT' ;; shaded |
|---|
| 11 | ;field2 = 'XLAT' ;; contour |
|---|
| 12 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 13 | format = '(F5.2)' |
|---|
| 14 | colors = 32 |
|---|
| 15 | pal = 22 ;33 ;39 ;6 ;11 ;0 ;11 ;6 ;0 ;6 ;39 ;19 |
|---|
| 16 | title_user = '' |
|---|
| 17 | title_axis = ['grid points','grid points'] |
|---|
| 18 | title_user = 'Friction velocity 50m ABG (m s!U-1!N) + Horizontal wind vectors' |
|---|
| 19 | ;title_user = 'Topography (m) + Horizontal wind vectors 50m ABG' |
|---|
| 20 | subtitle_user = 'LMD Mars Mesoscale / UTC = 01:00pm / Ls = 61.5!Uo!N / dx = 25km / single domain' |
|---|
| 21 | poscb=0.85 |
|---|
| 22 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 23 | ;;; FILL LIMITS ;;; |
|---|
| 24 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 25 | minfield_init= +00.1 |
|---|
| 26 | maxfield_init= +00.6 |
|---|
| 27 | ;format = '(I0)' |
|---|
| 28 | ;minfield_init= -6500. |
|---|
| 29 | ;maxfield_init= -1500. |
|---|
| 30 | ;pal = 19 |
|---|
| 31 | missing_value=1.e30 |
|---|
| 32 | 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 |
|---|
| 33 | lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min |
|---|
| 34 | ;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 |
|---|
| 35 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 36 | ;;; WINDS ;;; |
|---|
| 37 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 38 | windex=12. ;; default: 20. |
|---|
| 39 | stride=2. ;; default: 5. |
|---|
| 40 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 41 | ;;; CONTOUR ;;; |
|---|
| 42 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 43 | ;overcontour=overcontour-mean(overcontour) |
|---|
| 44 | overcontour=overcontour/1000. |
|---|
| 45 | lev=-10. + 0.2*findgen(100) |
|---|
| 46 | lev=-10. + 0.5*findgen(50) |
|---|
| 47 | ;lev=-10. + 0.1*findgen(200) |
|---|
| 48 | ;lev=0 |
|---|
| 49 | overcontour=-overcontour |
|---|
| 50 | lev=0.2*findgen(200) |
|---|
| 51 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 52 | ;;; AXIS ;;; |
|---|
| 53 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 54 | ;isotropic='false' ;; default: true |
|---|
| 55 | ;intervalx=0.5 |
|---|
| 56 | ;intervaly=0.5 |
|---|
| 57 | intervalx=5. |
|---|
| 58 | intervaly=5. |
|---|
| 59 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 60 | ;;; MAP LIMITS ;;; |
|---|
| 61 | ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; |
|---|
| 62 | ;windowx=[-51.0,-41.0] |
|---|
| 63 | ;windowy=[78.0,88.0] |
|---|
| 64 | windowx=[5.,55.] |
|---|
| 65 | windowy=[5.,55.] |
|---|
