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