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