| 1 | ;;--------------------------------------- |
|---|
| 2 | ; |
|---|
| 3 | ; fill here plot parameters then type: |
|---|
| 4 | ; |
|---|
| 5 | ; @trace.idl |
|---|
| 6 | ; @report.idl |
|---|
| 7 | ; @movie.idl |
|---|
| 8 | ; |
|---|
| 9 | ;;--------------------------------------- |
|---|
| 10 | |
|---|
| 11 | ;----------------- |
|---|
| 12 | ; initialization |
|---|
| 13 | ;----------------- |
|---|
| 14 | winds=0. |
|---|
| 15 | range=0. |
|---|
| 16 | topo=0. |
|---|
| 17 | colors=32 |
|---|
| 18 | backup_data='no' ;backup_data='yes' |
|---|
| 19 | already_data='no' ;already_data='yes' |
|---|
| 20 | folder='./WHERE_RESULTS_ARE/' |
|---|
| 21 | |
|---|
| 22 | |
|---|
| 23 | ;----------------- |
|---|
| 24 | ; check this first |
|---|
| 25 | ;----------------- |
|---|
| 26 | num=6 |
|---|
| 27 | ; |
|---|
| 28 | ;coord='height' |
|---|
| 29 | coord='model_level' |
|---|
| 30 | |
|---|
| 31 | ;---------------------------- |
|---|
| 32 | ; general plot settings |
|---|
| 33 | ;---------------------------- |
|---|
| 34 | plot='map' & level=-1 |
|---|
| 35 | plot='meridional' & nlon=-1 ;bug a corriger: 0 |
|---|
| 36 | ;plot='zonal' & nlat=-1 ;bug a corriger: 0 |
|---|
| 37 | |
|---|
| 38 | ;------------------------------------------- |
|---|
| 39 | ; easy access to already defined plots |
|---|
| 40 | ;------------------------------------------- |
|---|
| 41 | ;@map_sea_level.idl |
|---|
| 42 | @map_winds_only.idl |
|---|
| 43 | ;@map_ertel_vorticity.idl |
|---|
| 44 | ;@map_wind_stress.idl |
|---|
| 45 | |
|---|
| 46 | ;------------------------ |
|---|
| 47 | ; user-defined settings |
|---|
| 48 | ;------------------------ |
|---|
| 49 | ; keywords : plot, winds, field1, field2, range, topo, colors |
|---|
| 50 | ;------------------------ |
|---|
| 51 | |
|---|
| 52 | ;field1='tk' |
|---|
| 53 | ;nlon=100 |
|---|
| 54 | ;range=[200,250] |
|---|
| 55 | |
|---|
| 56 | |
|---|
| 57 | |
|---|
| 58 | |
|---|
| 59 | title_user='' ;to avoid a warning later |
|---|
| 60 | what_I_plot=0. ;to avoid a warning later |
|---|
| 61 | missing_value=0. ;to avoid a warning later |
|---|
| 62 | space=0. |
|---|
| 63 | ;;*************************************************** |
|---|
| 64 | ;; below this line do not modify number of lines |
|---|
| 65 | ;; because this will be copied ... with command tail |
|---|
| 66 | ;; comment a line if you do not need it :) |
|---|
| 67 | ;;*************************************************** |
|---|
| 68 | ;;------------------------ |
|---|
| 69 | ;; modify plot parameters |
|---|
| 70 | |
|---|
| 71 | ;; SECTION: |
|---|
| 72 | title_user, title_axis (2) |
|---|
| 73 | xtitle='Horizontal coordinate' |
|---|
| 74 | ytitle='Altitude' |
|---|
| 75 | endif else begin |
|---|
| 76 | xtitle=title_axis(0) |
|---|
| 77 | ytitle=title_axis(1) |
|---|
| 78 | endelse |
|---|
| 79 | |
|---|
| 80 | if (n_elements(space) eq 0) then begin |
|---|
| 81 | space=findgen(n_elements(what_I_plot(*,0))) |
|---|
| 82 | endif |
|---|
| 83 | if (n_elements(altitude) eq 0) then begin |
|---|
| 84 | altitude=findgen(n_elements(what_I_plot(0,*))) |
|---|
| 85 | endif |
|---|
| 86 | |
|---|
| 87 | if (n_elements(minfield_init) eq 0) then minfield_init=0. |
|---|
| 88 | if (n_elements(maxfield_init) eq 0) then maxfield_init=0. |
|---|
| 89 | if (n_elements(minspace) eq 0) then minspace=0. |
|---|
| 90 | if (n_elements(maxspace) eq 0) then maxspace=0. |
|---|
| 91 | if (n_elements(minalt) eq 0) then minalt=min(altitude) |
|---|
| 92 | if (n_elements(maxalt) eq 0) then maxalt=max(altitude) |
|---|
| 93 | if (n_elements(colors) eq 0) then colors=32 |
|---|
| 94 | if (n_elements(title) eq 0) then title='' |
|---|
| 95 | if (n_elements(pal) eq 0) then pal=33 |
|---|
| 96 | if (n_elements(format) eq 0) then format='(F6.2)' |
|---|
| 97 | if (format eq '') then format='(F6.2)' |
|---|
| 98 | |
|---|
| 99 | if ((n_elements(overvector_x) eq 0) or (n_elements(overvector_y) eq 0)) then begin |
|---|
| 100 | overvector_x=0. |
|---|
| 101 | overvector_y=0. |
|---|
| 102 | endif |
|---|
| 103 | |
|---|
| 104 | ; |
|---|
| 105 | ; dilatation factor for winds |
|---|
| 106 | ; |
|---|
| 107 | if (n_elements(overvector_x)*n_elements(overvector_y) ne 1) then space=space/factor |
|---|
| 108 | |
|---|
| 109 | |
|---|
| 110 | |
|---|
| 111 | |
|---|
| 112 | |
|---|
| 113 | |
|---|
| 114 | ;;------------------------ |
|---|
| 115 | ;; |
|---|
| 116 | space=space |
|---|
| 117 | ;minspace=158.5 & maxspace=160. |
|---|
| 118 | minalt=-3 & maxalt=7 |
|---|
| 119 | ;factor=15. |
|---|
| 120 | ;windex=20. |
|---|
| 121 | ;format='(F6.3)' |
|---|
| 122 | ;title_user='Sea-level surface pressure'+title_user |
|---|
| 123 | pal=33 |
|---|
| 124 | ;w=where(abs(what_I_plot) ge 0.03) |
|---|
| 125 | ;if (w[0] ne -1) then what_I_plot[w]=0.03 |
|---|
| 126 | ;;**** END **** |
|---|
