source: trunk/MESOSCALE/PLOT/SPEC/POLAR/defs/polar_inc_uv.pro.bak @ 134

Last change on this file since 134 was 85, checked in by aslmd, 14 years ago

LMD_MM_MARS et LMD_LES_MARS: ajout des routines IDL pour tracer les sorties --> voir mesoscale/PLOT
File size: 2.8 KB
Line 
1;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
2;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
3;ntime    = 7  ;; le vrai indice IDL
4nlevel   = 0
5filename = '/donnees/aslmd/MODELES/LMD_MM_MARS/OUTILS_CONVERSION/wrfout_d03_2024-06-20_06:00:00_zabg'
6filename = '/donnees/aslmd/MODELES/LMD_MM_MARS/OUTILS_CONVERSION/wrfout_d03_2024-06-20_06:00:00_z'
7filename = '/donnees/aslmd/MODELES/LMD_MM_MARS/OUTILS_CONVERSION/wrfout_d03_2024-06-20_06:00:00_z' & nlevel=1
8save_ps  = 'nest_polar_winds'
9;filename = '/donnees/aslmd/MODELES/LMD_MM_MARS/OUTILS_CONVERSION/wrfout_d01_2024-06-20_06:00:00_z'
10;save_ps  = 'single_polar_winds'
11;filename = '/tmp7/aslmd/one_more_nest/test/wrfout_d02_2024-06-19_06:00:00_zabg'
12;save_ps  = '2nest_polar_winds'
13
14;coord2d  = 'false'  ;; si on veut grid point coordinates ;;CHANGER windowx and windowy
15;field1   = 'HGT'  ;; omettre trace vitesse horizontale
16;field2  = 'XLAT' ;; contour
17;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
18format        = '(I0)' ;'(F4.1)' ;'(I0)'
19colors        =  32
20pal           =  22 ;19 ;22 ;33 ;39 ;6 ;11 ;0 ;11 ;6 ;0 ;6 ;39 ;19
21title_user    = 'Horizontal wind (m s!U-1!N) 20m AGL'
22title_user    = 'Horizontal wind (m s!U-1!N) 6km AMR'
23title_user    = 'Horizontal wind (m s!U-1!N) 20km AMR'
24subtitle_user   = 'LMD Mars Mesoscale'
25subtitle_user   = subtitle_user + ' / UTC!D0!N = 06:00am'
26subtitle_user   = subtitle_user + ' / Ls = 160!Uo!N'
27subtitle_user   = subtitle_user + ' / dx = 12km'
28subtitle_user   = subtitle_user + ' / nest 3'
29poscb=0.95
30;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
31;;; FILL LIMITS                                ;;;
32;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
33minfield_init= 30. ;0.
34maxfield_init= 60. ;30.
35missing_value=1.e30
36lim_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
37lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min
38;;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
39;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
40;;; WINDS                                      ;;;
41;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
42windex=50. ;25.         ;; default: 20.
43stride=5.       ;; default: 5.
44;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
45;;; CONTOUR                                    ;;;
46;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
47;overcontour=overcontour-mean(overcontour)     
48overcontour=overcontour/1000.
49lev=-10. + 0.5*findgen(50)
50;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
51;;; AXIS                                       ;;;
52;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
53;isotropic='false'       ;; default: true
54intervalx=2.
55intervaly=2.
56;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
57;;; MAP LIMITS                                 ;;;
58;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
59windowx=[-42.,-26.]
60windowy=[-34.,-20.]
61
Note: See TracBrowser for help on using the repository browser.