source: trunk/MESOSCALE/PLOT/SPEC/MAP/defs/map_uvt_inc.pro @ 205

Last change on this file since 205 was 91, checked in by aslmd, 14 years ago

mars: test outliers [dans initracer.F, commente] LMD_MM_MARS: modifications mineures [retrocompatibilite ancienne physique, callphys.def test pour nouvelle physique] PLOT: generalisation de la routine map_uvt pour pouvoir tracer des figures en projection polaire complete
File size: 5.5 KB
Line 
1;;;;;*************************************** FOLDER
2        folder  = '/tmp7/aslmd/'
3        coord2d = 'false'       ;; for non-regular projections  ;;CHANGER windowx and windowy
4        coord2d = 'polar'       ;; for polar projection with map_set ;; utiliser use_lt=99
5        filename        = folder + 'wrfout_d01_2024-05-02_01:00:00'
6        ini_utc         = 1           ;; cf. name of file
7        freq            = 1           ;; cf. 1 output per ... hour
8        save_ps         = 'POLAR_newphys_RAC_TAUTES'
9
10;;;;;*************************************** LVL & FLD & TIME
11        nlevel          = 0
12        use_lt          = 99               ;; cf. what user wants (99 pour tous)
13        field1          = 'TAU_ICE'        ;; comment to trace horizontal velocity
14        ;field2         = 'XLAT'           ;; contour
15        no3d            = 'true'
16        ;overvector_x    = 0               ;; comment out to get rid of vectors
17        ;overvector_y    = 0               ;; comment out to get rid of vectors
18
19;;;;;*************************************** TWEAK VAR
20        ;what_I_plot = what_I_plot * 0.
21        ;print, max(what_I_plot), min(what_I_plot)
22
23;;;;;*************************************** PLOT TITLES
24        title_user      =       'Water ice cloud optical depth at 825 cm!U-1!N'
25        title_axis      =       ['East longitude','North latitude']
26        ;subtitle_user   = 'LMD_MM'
27        ;;subtitle_user   = subtitle_user + ' / LT = '+string(use_lt,'(I0)')+'h'
28        ;;subtitle_user   = subtitle_user + ' / UTC!D0!N = 06:00am' ;; pour les series
29        ;subtitle_user   = subtitle_user + ' / Ls = 120!Uo!N'
30        ;subtitle_user   = subtitle_user + ' / dx = 10km [S]'
31
32;;;;;*************************************** COLOR TABLES
33        flag_cb         =       'true'
34        poscb           =       0.85
35        format          =       '(F5.2)'
36        colors          =       32
37        pal             =       22              ;; GOOD: 4, 18, 22, 16, 37, 33, 39, 6, 11
38
39;;;;;*************************************** FILL LIMITS
40        minfield_init   =       0.
41        maxfield_init   =       0.55
42        ndiv            =       11
43
44        ;;;;;*************************************** LIMIT TRICKS
45        ;;;;;********************* must always follow FILL LIMITS
46                missing_value=1.e30
47                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
48                lim_min = minfield_init & w=where(what_I_plot le lim_min) & if (w[0] ne -1) then what_I_plot[w]=lim_min
49                ;lim_blank = 2. & w=where(abs(what_I_plot) le lim_blank) & if (w[0] ne -1) then what_I_plot[w]=missing_value
50
51;;;;;*************************************** WINDS
52        windex          =       30.             ;; DEF: 20.
53        stride          =       3.              ;; DEF: 5.
54
55;;;;;*************************************** CONTOUR
56        overcontour     =       overcontour/1000.
57        lev             =       -10. + 0.5*findgen(80)
58
59;;;;;*************************************** AXIS
60        isotropic       =       'true'
61        intervalx       =       30.0
62        intervaly       =       05.0
63
64;;;;;*************************************** MAP LIMITS
65        windowx         =       [-180.,180.]       
66        windowy         =       [65.,90.]
67
68;;;;;***************************************;;;;;
69;;;;;***************************************;;;;;
70;;;;;***************************************;;;;;
71;;;;;***************************************;;;;;
72
73;;;;;*************************************** SETTING TIME (do not modify)
74        if (n_elements(windowx) eq 0) then windowx = 0
75        utc_to_lt       = mean(windowx) / 15.    ;; cf. longitude -- LT = UTC + utc_to_lt
76        use_utc         = use_lt - utc_to_lt
77        zentime         = floor(((24 + use_utc - ini_utc) MOD 24)/freq)        ;; TRUE IDL SUBSCRIPT...
78        if (use_lt ne 99) then ntime = zentime else ntime = 99    ;; ou commenter pour avoir tous les pas de temps
79        print, zentime, ntime
80
81;;;;;***************************************
82;;;;;***************************************
83;;;;;***************************************
84;;;;;***************************************
85
86;;;
87;;; VERTICAL SECTION LIMITS
88;;;
89;minalt=-5.             ;; grepSEC
90;maxalt=40.             ;; grepSEC
91;minspace=0.0           ;; grepSEC
92;maxspace=35.0          ;; grepSEC
93
94;;;
95;;; METRIC UNITS FOR VERTICAL SECTION
96;;;
97;factor=10.        ;; grepSEC
98;space=space*60. & spacekm='true'       ;; grepSEC
99;minspace=minspace*60./factor           ;; grepSEC
100;maxspace=maxspace*60./factor           ;; grepSEC
101;intervalx=round(intervalx*60./factor)   ;; grepSEC
102
103
104;;;
105;;; TRICKS
106;;;
107
108;; pour tracer juste les vecteurs sur un fond uni ou vide
109;pal=0                                          ;; 1/4 grepMAP
110;what_I_plot(*,*)=what_I_plot(*,*)*0.+0.2       ;; 2/4 grepMAP
111;what_I_plot(0,0)=0.                            ;; 3/4 grepMAP
112;flag_cb='false'                                        ;; 4/4 grepMAP
113
114;;; truc pour tracer juste les contours et une zone grisee de topo
115;pal=0                                           ;; 1/9 grepALL
116;w=where(abs(what_I_plot) lt missing_value)     ;; 2/9 grepALL
117;what_I_plot[w]=0.                              ;; 3/9 grepALL
118;w=where(abs(what_I_plot) gt missing_value)     ;; 4/9 grepALL
119;what_I_plot[w]=0.75                            ;; 5/9 grepALL
120;w=where(what_I_plot eq 0.)                     ;; 6/9 grepALL
121;what_I_plot[w]=missing_value                           ;; 7/9 grepALL
122;what_I_plot(0,0)=1. & what_I_plot(1,0)=0.      ;; 8/9 grepALL
123;flag_cb='false'                                ;; 9/9 grepALL
124
Note: See TracBrowser for help on using the repository browser.