| 1 | PRO plt_def |
|---|
| 2 | ; |
|---|
| 3 | ; Defaults plot attributes |
|---|
| 4 | ; |
|---|
| 5 | @com_eg |
|---|
| 6 | @common |
|---|
| 7 | |
|---|
| 8 | ; boxes used in horizontal/vertical display |
|---|
| 9 | ; ----------------------------------------- |
|---|
| 10 | ; |
|---|
| 11 | ; horizontal domain |
|---|
| 12 | |
|---|
| 13 | ; box_h = [20,380,-30,30] |
|---|
| 14 | box_h = [120,290,-20,20] |
|---|
| 15 | ; box_h = [10,180,-20,20] |
|---|
| 16 | box_h = [120,290,-30,30] |
|---|
| 17 | ; box_h = [-180,180,-90,90] |
|---|
| 18 | ; box_h = [20,380,-90,90] |
|---|
| 19 | ; box_h = [20,150,-30,60] |
|---|
| 20 | ; box_h = [120,290,-20,20] |
|---|
| 21 | box_h = [120,290,-30,30] |
|---|
| 22 | ; box_h = [100,300,-40,40] |
|---|
| 23 | ; box_h = [120,290,0,60] |
|---|
| 24 | ; box_h = [40,120,-30,30] |
|---|
| 25 | ; box_h = [-180,180,-90,90] |
|---|
| 26 | box_h = [20,380,-90,90] |
|---|
| 27 | ; box_h = [20,380,-60,60] |
|---|
| 28 | ; box_h = [-20,340,-90,90] |
|---|
| 29 | ; box_h = [160,210,-5,5] ; nino 4 |
|---|
| 30 | ; box_h = [210,270,-5,5] ; nino 3 |
|---|
| 31 | ; box_h = [100,300,-30,30] |
|---|
| 32 | |
|---|
| 33 | ; working grid longitudes |
|---|
| 34 | |
|---|
| 35 | glamboundary_box = [20, 380] |
|---|
| 36 | ; glamboundary_box = [-20, 340] |
|---|
| 37 | |
|---|
| 38 | ; time domain |
|---|
| 39 | ; |
|---|
| 40 | time_domain = ['199609', '199802']; 'auto' or ['date1','date2'] |
|---|
| 41 | ; time_domain = 'auto'; 'auto' or ['date1','date2'] |
|---|
| 42 | |
|---|
| 43 | common_time = 0 ; set to 1 to have common time axis (using first one in t_ plots) |
|---|
| 44 | |
|---|
| 45 | ; vertical domain |
|---|
| 46 | |
|---|
| 47 | depth_z = 1000 |
|---|
| 48 | zoom_z = 300 |
|---|
| 49 | |
|---|
| 50 | pres_min = 10 |
|---|
| 51 | pres_max = 500 |
|---|
| 52 | |
|---|
| 53 | msf_mean = 0 |
|---|
| 54 | |
|---|
| 55 | ; vertical average for 3D fields |
|---|
| 56 | |
|---|
| 57 | vert_type = '0' ; 'z' for depth/altitude or 'level' or '0' for nothing |
|---|
| 58 | vert_mean = [0, 0] ; [depth1,depth2] or [level1,level2] in C notation i.e 0 to jpk-1 |
|---|
| 59 | ; vert_mean = [90000, 100000] |
|---|
| 60 | |
|---|
| 61 | ; density domain (sigma) + delta sigma |
|---|
| 62 | |
|---|
| 63 | sig_min = 20. |
|---|
| 64 | sig_max = 28. |
|---|
| 65 | sig_del = 0.2 |
|---|
| 66 | |
|---|
| 67 | sig_bowl = 0 ; 0/1 bowl overlay on density binned plots |
|---|
| 68 | |
|---|
| 69 | ; Statistics |
|---|
| 70 | ; ----------- |
|---|
| 71 | |
|---|
| 72 | ; max date in spectrums (in days) |
|---|
| 73 | ; spectrum window (in days) |
|---|
| 74 | ; tukey window (in years) |
|---|
| 75 | ; bootstrap window to compute error bar for 1d time serie (in units of time array) |
|---|
| 76 | |
|---|
| 77 | max_spec = 7*360 |
|---|
| 78 | spec_win = long(20)*long(360) |
|---|
| 79 | tukey = 20 |
|---|
| 80 | boot_win = [24, 36, 48] |
|---|
| 81 | |
|---|
| 82 | ; cross-correlation and normalisation: |
|---|
| 83 | c_normal = 0 ; normalise data serie |
|---|
| 84 | c_correl = 1 ; when overlay of 1d curves : 0/1 |
|---|
| 85 | lag_correl = 12L ; bound for lagged-correlation [-lag_correl,..,0,...,lag_correl] |
|---|
| 86 | |
|---|
| 87 | ; standard deviation |
|---|
| 88 | stddev_diff = 0 ; 0/1 (if = 1 makes a difference between the |
|---|
| 89 | ; std dev for a given month and the std dev of the whole time serie) |
|---|
| 90 | ; Active with the make_stddev macro for a given month (ex : @@sosstdev_01 for |
|---|
| 91 | ; sst std dev only for January) |
|---|
| 92 | ; ensembles |
|---|
| 93 | ensbl_mean = 0 ; option for pltt 't' ensembles |
|---|
| 94 | ; 0 : just members |
|---|
| 95 | ; 1 : members + ensemble mean |
|---|
| 96 | ; 2 : min/max envelop of members + ensemble mean |
|---|
| 97 | |
|---|
| 98 | ; axis |
|---|
| 99 | ; ---- |
|---|
| 100 | ; lat_axis : latitudinal axis : 'reg'ular or 'sin'us latitude |
|---|
| 101 | lat_axis = 'reg' |
|---|
| 102 | |
|---|
| 103 | ; min /max |
|---|
| 104 | ; -------- |
|---|
| 105 | ; free minmax in 1D-plots (yes) or take fld_glo_mmx.def value (no) [yes] |
|---|
| 106 | free_1d_minmax = 'no' |
|---|
| 107 | ; free_1d_minmax = 'yes' |
|---|
| 108 | |
|---|
| 109 | ; Hovmoellers |
|---|
| 110 | ; ----------- |
|---|
| 111 | ; trend_typ : 0 no trend |
|---|
| 112 | ; 1 remove initial value to serie (trend) |
|---|
| 113 | ; 2 remove previous value to current (drift) |
|---|
| 114 | ; 3[<n>] remove serie to mean of <n> final values (inverse trend) |
|---|
| 115 | ; 4[<n>] remove [n]-average running mean of serie (anomaly) |
|---|
| 116 | ; for monthly serie use <n>=12 to remove mean |
|---|
| 117 | ; seasonal cycle |
|---|
| 118 | ; this value can be field-overriden by specifying @t<trend_typ> after |
|---|
| 119 | ; timave in post-it line |
|---|
| 120 | ; 6[<n>] time integral of field activation. <n> is the |
|---|
| 121 | ; number of previous time steps to integrate from |
|---|
| 122 | ; (default is from beginning of time serie) |
|---|
| 123 | ; |
|---|
| 124 | ; |
|---|
| 125 | ; field_int : field time integral (0/1) |
|---|
| 126 | ; nb_cycles : number of repeated cycles for seasonal time series |
|---|
| 127 | ; def_stride : default stride in hovmoeller sampling (see domain_boxes.def) |
|---|
| 128 | |
|---|
| 129 | trend_typ = '0' |
|---|
| 130 | field_int = 0 |
|---|
| 131 | nb_cycles = 2 |
|---|
| 132 | def_stride = 1 |
|---|
| 133 | asciidir = homedir+'/ENSO_works/data/SST/' |
|---|
| 134 | ; |
|---|
| 135 | ; calendar type |
|---|
| 136 | ; 0: 365 days/year |
|---|
| 137 | ; 1: gregorian |
|---|
| 138 | ; N: N days per month |
|---|
| 139 | ; |
|---|
| 140 | calendar_type = 1 |
|---|
| 141 | ; |
|---|
| 142 | ; Graphic options |
|---|
| 143 | ; ---------------- |
|---|
| 144 | ; shadind : (fill) 1/0 |
|---|
| 145 | ; pal_type: 'bw', 'col', '2dom' (2dom = read 2 domains in isolignes - ISOSP) |
|---|
| 146 | ; grey_shade: % of black for 2dom case |
|---|
| 147 | ; grey_shade_1/2: % of black for 2dom case : difference plots (-/+) |
|---|
| 148 | ; col_palette: 'yes' or 'no' (put colorbar or not) |
|---|
| 149 | ; |
|---|
| 150 | shading = 1 |
|---|
| 151 | pal_type = 'col' |
|---|
| 152 | grey_shade = 20 |
|---|
| 153 | grey_shade_1 = 5 |
|---|
| 154 | grey_shade_2 = 20 |
|---|
| 155 | col_palette = 'no' |
|---|
| 156 | cont_fill = 1 ; continent fill for atmosphere grids (1/0) |
|---|
| 157 | cont_real = 1 ; draw real continents (0) |
|---|
| 158 | ; = 0 (no real continents drawn) |
|---|
| 159 | ; = 1 (real continents drawn instead of the mask) |
|---|
| 160 | ; = 2 (mask + real continents drawn) |
|---|
| 161 | vector_sample = 2 ; vector sampling (1) |
|---|
| 162 | nino_plot = 1 ; draw nino 1,2,3,4 boxes on xy plots |
|---|
| 163 | ;; |
|---|
| 164 | ; line_thick = [2,2,2,2,2,4,4,4,4,4] ; 1D plot line thickness (for each overlay curves) |
|---|
| 165 | ; line_style = [1,1,1,1,1,3,3,3,3,3] ; 1D plot line style (for each overlay curves) |
|---|
| 166 | ; line_color = [1,2,3,4,5,1,2,3,4,5] ; 1D plot line color (for each overlay curves) |
|---|
| 167 | line_thick = [2,2,2,4,4,4,4,4,4,6,6,6] ; 1D plot line thickness (for each overlay curves) |
|---|
| 168 | line_style = [1,1,1,2,2,2,3,3,3,1,1,1] ; 1D plot line style (for each overlay curves) |
|---|
| 169 | line_color = [1,2,4,1,2,4,1,2,4,1,2,4] ; 1D plot line color (for each overlay curves) |
|---|
| 170 | |
|---|
| 171 | ensbl_color_member = 1 ; 0 = same color(from line_color(overlay#)) + attributes below, |
|---|
| 172 | ; 1 = use line attributes above |
|---|
| 173 | ensbl_thick_member = 1 |
|---|
| 174 | ensbl_thick_mean = 7 |
|---|
| 175 | ensbl_style_member = 1 |
|---|
| 176 | ensbl_style_mean = 1 |
|---|
| 177 | |
|---|
| 178 | ; 1=black; 2=red; 3=green; 4=blue; 5=black; 6=purple |
|---|
| 179 | ; 43 ---- 90 ---- 160 ---- 190 ---- 210 ---- 254 |
|---|
| 180 | ; blue turq green yell orange red |
|---|
| 181 | |
|---|
| 182 | ;; |
|---|
| 183 | ;; Scatter/bining plot variables |
|---|
| 184 | ;; |
|---|
| 185 | symbol_families = '4x3' ; 4x3 number of colors in time symbol plots (modulo) |
|---|
| 186 | ; use <n> or <n>x<m> (n same colors, modulo nxm) |
|---|
| 187 | symbol_style = [1,1,1,1,2,2,1,1,1,1,1,1] ; Symbol style (for each overlay symbol) |
|---|
| 188 | symbol_color = [1,2,3,4,2,3,7,8,9,10,11,12] ; Symbol line color (for each overlay symbol) |
|---|
| 189 | ; symbol_color = [1,1,1,1,1,3,7,8,9,10,11,12] ; Symbol line color (for each overlay symbol) |
|---|
| 190 | symbol_size = 1 ; Symbol size (1) |
|---|
| 191 | mean_sc_only = 0 ; 0 = default yfx (no SC plotted) |
|---|
| 192 | ; 1 = add mean SC curve |
|---|
| 193 | ; 2 = Only plot mean SC in yfx when symbol_families='4x3' and hotyp=t |
|---|
| 194 | ; (3) = Only plot SC of std dev when symbol_families='4x3' and hotyp=t |
|---|
| 195 | ; (require 1m@t412) |
|---|
| 196 | ; 4 = Only plot y=f(x) per month + errorbar when symbol_families='12x1' |
|---|
| 197 | ; |
|---|
| 198 | linfit_sep = 0.0 ; separation for two domains linear fit computation (on x axis) |
|---|
| 199 | ; set to -99999 to desactivate |
|---|
| 200 | linfit_map = '' ; ''/'p'/'m' to compute fit for all/before/after linfit_sep (maps only) |
|---|
| 201 | corr_thres = 0.1 ; correlation threshold to plot regression maps (between 0.1 and 0.2) |
|---|
| 202 | ; bining interval (one number) or |
|---|
| 203 | ; binning array [b1,b2,b3,b4,...] |
|---|
| 204 | ; bin_interval1 = [-60, -40, -20, 0, 20, 40] |
|---|
| 205 | bin_interval1 = [22, 23.5, 25, 25.5, 26, 26.5, 27, 27.5, 28, 28.5, 29, 29.5, 30, 31] |
|---|
| 206 | ; bin_interval = [-.1, -.08, -.06, -.04, -.02, 0, .02, .04, .06, .1] ; taux |
|---|
| 207 | ; bin_interval1 = [0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10] |
|---|
| 208 | ; bin_interval1 = [-10, -9, -8, -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10] |
|---|
| 209 | ; bin_interval2 = [-10, -9, -8, -7, -6, -5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10] |
|---|
| 210 | bin_interval2 = [-5, -4, -3, -2, -1, 0, 1, 2, 3, 4, 5] |
|---|
| 211 | |
|---|
| 212 | large_domain = 1 ; remove large domain average in binning |
|---|
| 213 | |
|---|
| 214 | ; Graphic keywords for axis, etc... |
|---|
| 215 | ; |
|---|
| 216 | contour_options = ',/nocontour,cell_fill=2' ;;;;;; include this line to delete contour labels |
|---|
| 217 | ; contour_options = ',cell_fill=2' ;;;;;; include this line to include contour label |
|---|
| 218 | contour_options = '' ;;;;;; include this line to include contour label |
|---|
| 219 | look = ',xthick=2,ythick=2,zthick=2,sepdate='' ''' |
|---|
| 220 | |
|---|
| 221 | ; use of page space and titles |
|---|
| 222 | |
|---|
| 223 | windowsize_scale = 1.0 ; size of screen display (visu option 'v') |
|---|
| 224 | fill_space = 1 ; fill space on plot (0/1) (/rempli option) ? |
|---|
| 225 | marge_option = 'marge=[0,0,2,2]' ; (use [0, 0, -2, -2] to fill up even more space |
|---|
| 226 | title_type = 'TS' ; default type of titles |
|---|
| 227 | ; 'T' for title only, 'S' for subtitle only |
|---|
| 228 | ; 'TS' for both |
|---|
| 229 | ; 'off' for no titles |
|---|
| 230 | xchartxt = 1.0 ; x axis font size (1.0) |
|---|
| 231 | ychartxt = 1.0 ; y axis font size (1.0) |
|---|
| 232 | |
|---|
| 233 | |
|---|
| 234 | default_txt_format = 'E' ; default text in legend for 1D plots (E)xperiment(V)ariable[L]ong name[U]nits(B)ox |
|---|
| 235 | |
|---|
| 236 | ; Multiple windows (0/1) turn on for one window per plot |
|---|
| 237 | |
|---|
| 238 | multi_win = 0 |
|---|
| 239 | |
|---|
| 240 | ; Save PostScript/pdf/gif files |
|---|
| 241 | |
|---|
| 242 | save_ps = 2 ; set to 2 to convert to pdf |
|---|
| 243 | file_naming = 'prompt' ; 'auto' or 'prompt' |
|---|
| 244 | |
|---|
| 245 | ; Averages |
|---|
| 246 | ; --------- |
|---|
| 247 | ; land-sea mask use in atmospheric boxes 0 : take all |
|---|
| 248 | ; 1 : ocean points |
|---|
| 249 | ; 2 : land points |
|---|
| 250 | ; |
|---|
| 251 | atmos_msk = 1 |
|---|
| 252 | |
|---|
| 253 | ; |
|---|
| 254 | ; Data domain (OPA only) |
|---|
| 255 | ; ----------- |
|---|
| 256 | v_config_oce = 'L30' ; L46 for Drakkar ORCA05/ORCA2 |
|---|
| 257 | orca_mask_version = 'V3' ;V1(ST7/CT runs), V2 (? runs), V3 (CH, Dan runs, CM4) |
|---|
| 258 | ; Zoomed domains ORCA2 (indexes are within array [0:181,0:148]) |
|---|
| 259 | ; global <file>_grid_<grd>.nc 27118 pts |
|---|
| 260 | ; equator <file>_grid_<grd>_eq.nc (j=73) 182 pts |
|---|
| 261 | ; equator_band <file>_grid_<grd>_eqb.nc (j=64,82 -> [5S,5N]) 2912 pts = 10% |
|---|
| 262 | ; pacific <file>_grid_<grd>_pac.nc (j=49,97 -> [30S,30N], |
|---|
| 263 | ; i=16,106 -> [110E,70W]) 4320 pts = 16% |
|---|
| 264 | ; |
|---|
| 265 | ; pacific_eq <file>_grid_<grd>_paceq.nc (j=64,82 -> [5S,5N], |
|---|
| 266 | ; i=16,106 -> [110E,70W]) 1620 pts = 6% |
|---|
| 267 | ; zonal <file>_grid_<grd>_diaznl.nc |
|---|
| 268 | ; |
|---|
| 269 | ; Zoom ORCA05 (indexes are within array [0:721,0:510]) |
|---|
| 270 | ; pacific i=54,434 & j=183,312 |
|---|
| 271 | ; |
|---|
| 272 | ; data_domain = 'glosea' |
|---|
| 273 | ; data_domain = 'equator' |
|---|
| 274 | data_domain = 'equator_band' |
|---|
| 275 | ; data_domain = 'pacific_eq' |
|---|
| 276 | data_domain = 'zonal' |
|---|
| 277 | data_domain = 'pacific' |
|---|
| 278 | data_domain = 'global' |
|---|
| 279 | |
|---|
| 280 | ; grids list (IPCC atmos regular) |
|---|
| 281 | nc_grids_list = ['ct42', 'mt42', 'mt106', 'nt85', 'giss4x5', 'giss4x3', 'cgcm3_1', 'mk3', 'gfdl2.5x2', 'fgoals1', 'inmcm3', 'cgcm232', 'nt42', 'n96', 'mt63', 'lmdzl', 'lmdzh', 'lmdz142', 'ncpt62', 'ert106', 'it30', 'bt42', 'reg1', 'hadcm3t', 'reg2.5'] |
|---|
| 282 | |
|---|
| 283 | ; machine type ('x' or 'WIN') |
|---|
| 284 | dev_type='x' |
|---|
| 285 | |
|---|
| 286 | ; debug mode |
|---|
| 287 | debug_w = 1 |
|---|
| 288 | |
|---|
| 289 | END |
|---|
