source: trunk/mesoscale/PLOT/SPEC/POLAR/defs/polar_inc_uv.pro.masse.2 @ 85

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