source: trunk/UTIL/IDLplot/SPEC/MAP/defs/map_uvt_inc.pro.map @ 1306

Last change on this file since 1306 was 85, checked in by aslmd, 15 years ago

LMD_MM_MARS et LMD_LES_MARS: ajout des routines IDL pour tracer les sorties --> voir mesoscale/PLOT
File size: 6.1 KB
Line 
1;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; FOLDER
2;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; FOLDER
3        folder  = '/d5/aslmd/LMD_MM_MARS_SIMUS/OM/'
4        coord2d = 'false'       ;; for non-regular projections 
5filename        = folder + 'OM6_TI85/wrfout_d01_2024-06-43_06:00:00_zabg'
6save_ps         = 'OM6_TI85_winds'
7;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; TIME
8;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; TIME
9ini_utc         = 6     ;; cf. name of file
10freq            = 1     ;; cf. 1 output per ... hour
11utc_to_lt       = -9    ;; cf. longitude -- LT = UTC + utc_to_lt
12        use_lt  = 2     ;; cf. what user wants
13        use_utc = use_lt - utc_to_lt
14        ntime   = (use_utc - ini_utc)/freq      ;; TRUE IDL SUBSCRIPT ... ajouter un modulo...
15        print, 'CHECK ... ', filename, ntime
16;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; LVL & FLD
17;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; LVL & FLD
18nlevel          = 0
19field1          = 'WAVE'        ;; comment to trace horizontal velocity
20;no3d           = 'true'
21;overvector_x   = 0             ;; comment out to get rid of vectors
22;overvector_y   = 0             ;; comment out to get rid of vectors
23;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; TWEAK VAR
24;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; TWEAK VAR
25        ;a0 = 64.839418 & a1 = 3.1521937 & a2 = 0.0054538848 ;; Ls ~ 173 -- 2AM
26        ;;a0 = 65.052165 & a1 = 3.1228993 & a2 = 0.0053787417 ;; Ls ~ 173
27        ;;;a0 = 64.039300 & a1 = 3.1378104 & a2 = 0.0055369148 ;; Ls ~ 120
28        ;what_I_plot = - ( alog ( a1 - what_I_plot / a0 ) ) /  a2
29        ;print, max(what_I_plot), min(what_I_plot)
30        ;;overvector_x=0
31        ;;overvector_y=0
32;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
33;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
34
35
36;;;;;*************************************** PLOT TITLES
37        ;title_user     =       'Surface temperature (K) and Winds 10m ABG (m s!U-1!N)'
38        ;title_user     =       'Wind-induced apparent thermal inertia (J m!U-2!N s!U-0.5!N K!U-1!N)'
39        title_user      =       'Winds 10m ABG (m s!U-1!N)'
40        ;title_user      =       'Surface temperature (K)'
41        ;subtitle_user   =       'LT = 03:00am  /  Ls = 173!Uo!N  /  dx = 10km [single]  /  Uniform TI = 85 J m!U-2!N s!U-0.5!N K!U-1!N'
42        subtitle_user   =       'LT = 02:00am  /  Ls = 173!Uo!N  /  dx = 6km [single]  /  Uniform TI = 85 J m!U-2!N s!U-0.5!N K!U-1!N'
43        title_axis      =       ['East longitude','North latitude']
44
45;;;;;*************************************** COLOR TABLES
46        flag_cb         =       'true'
47        poscb           =       0.96
48        format          =       '(I0)'
49        ;format          =       '(F4.1)'
50        ;colors          =       128
51        colors          =       32
52        pal             =       22              ;; GOOD: 4, 18, 22, 16, 37, 33, 39, 6, 11
53        pal             =       4
54        pal             =       0
55
56;;;;;*************************************** FILL LIMITS
57        ;minfield_init   =       150.
58        ;maxfield_init   =       180.
59        minfield_init   =       -10.
60        maxfield_init   =       0.
61
62        ;;;;;*************************************** LIMIT TRICKS
63        ;;;;;********************* must always follow FILL LIMITS
64                missing_value=1.e30
65                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
66                lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min
67                ;lim_blank = 2. & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value
68
69;;;;;*************************************** WINDS
70        ;windex          =       10.            ;; DEF: 20.
71        windex          =       20.             ;; DEF: 20.
72        ;stride          =       1.             ;; DEF: 5.
73        ;stride          =       2.              ;; DEF: 5.
74        stride          =       5.              ;; DEF: 5.
75
76;;;;;*************************************** CONTOUR
77        overcontour     =       overcontour/1000.
78        ;lev             =       50. + 50.*findgen(20)
79        ;lev             =       -10. + 0.2*findgen(20)
80        lev             =       -10. + 2.*findgen(20)
81        ;lev             =       -10. + 1.*findgen(40)
82
83;;;;;*************************************** AXIS
84        isotropic       =       'true'
85        ;intervalx       =       0.5
86        intervalx       =       2.0
87        ;intervaly       =      intervalx
88        intervaly       =       1.0
89
90;;;;;*************************************** MAP LIMITS
91        ;windowx         =       [-144.,-126.]
92        ;windowy         =       [10.,26.]
93        windowx         =       [-146.,-126.]       
94        ;windowy         =      [10.,28.]
95        ;windowx         =      [-146.,-128.]
96        ;windowy         =      [12.,26.]
97        windowy         =       [11.,27.]
98
99;;;;;***************************************
100;;;;;***************************************
101;;;;;***************************************
102;;;;;***************************************
103
104;;;
105;;; VERTICAL SECTION LIMITS
106;;;
107;minalt=-5.             ;; grepSEC
108;maxalt=40.             ;; grepSEC
109;minspace=0.0           ;; grepSEC
110;maxspace=35.0          ;; grepSEC
111
112;;;
113;;; METRIC UNITS FOR VERTICAL SECTION
114;;;
115;factor=10.        ;; grepSEC
116;space=space*60. & spacekm='true'       ;; grepSEC
117;minspace=minspace*60./factor           ;; grepSEC
118;maxspace=maxspace*60./factor           ;; grepSEC
119;intervalx=round(intervalx*60./factor)   ;; grepSEC
120
121
122;;;
123;;; TRICKS
124;;;
125
126;; pour tracer juste les vecteurs sur un fond uni ou vide
127;pal=0                                          ;; 1/4 grepMAP
128;what_I_plot(*,*)=what_I_plot(*,*)*0.+0.2       ;; 2/4 grepMAP
129;what_I_plot(0,0)=0.                            ;; 3/4 grepMAP
130;flag_cb='false'                                        ;; 4/4 grepMAP
131
132;;; truc pour tracer juste les contours et une zone grisee de topo
133;pal=0                                           ;; 1/9 grepALL
134;w=where(abs(what_I_plot) lt missing_value)     ;; 2/9 grepALL
135;what_I_plot[w]=0.                              ;; 3/9 grepALL
136;w=where(abs(what_I_plot) gt missing_value)     ;; 4/9 grepALL
137;what_I_plot[w]=0.75                            ;; 5/9 grepALL
138;w=where(what_I_plot eq 0.)                     ;; 6/9 grepALL
139;what_I_plot[w]=missing_value                           ;; 7/9 grepALL
140;what_I_plot(0,0)=1. & what_I_plot(1,0)=0.      ;; 8/9 grepALL
141;flag_cb='false'                                ;; 9/9 grepALL
142
Note: See TracBrowser for help on using the repository browser.