1 | PRO icecontour, field, hi, hs, nlevels, colors, levels, ct, z, time, nz, nt, $ |
---|
2 | tmin, tmax, zmin, zmax, title |
---|
3 | |
---|
4 | ; |
---|
5 | ;----------------------------------------- |
---|
6 | ; Make the ice thickness, snow depth plot |
---|
7 | ;----------------------------------------- |
---|
8 | ; |
---|
9 | LOADCT, 0 |
---|
10 | PLOT, [tmin, tmax], [zmin,zmax], /NODATA, $ |
---|
11 | title = title, ytitle = "z[m]", CHARSIZE=1, XSTYLE = 1, YSTYLE = 1;, $ |
---|
12 | OPLOT, time , - hi |
---|
13 | OPLOT, time , hs |
---|
14 | ; |
---|
15 | ;------------ |
---|
16 | ; Dummy axes |
---|
17 | ;------------ |
---|
18 | ; time |
---|
19 | time_matrix = FLTARR(nz+2, nt) |
---|
20 | FOR index = 0, nz+1 DO BEGIN |
---|
21 | time_matrix[index,*] = time[*] |
---|
22 | ENDFOR |
---|
23 | |
---|
24 | ; depth |
---|
25 | z_matrix = FLTARR(nz+2, nt) |
---|
26 | z_matrix[0,*] = 0.0 |
---|
27 | FOR index = 1, nz DO BEGIN |
---|
28 | z_matrix[index,*] = z[index-1,*] |
---|
29 | ENDFOR |
---|
30 | z_matrix[nz+1,*] = hi |
---|
31 | |
---|
32 | ; field to plot |
---|
33 | field_matrix = FLTARR(nz+2,nt) |
---|
34 | field_matrix[0,*] = field[0,*] |
---|
35 | FOR index = 1, nz DO BEGIN |
---|
36 | field_matrix[index,*] = field[index-1,*] |
---|
37 | ENDFOR |
---|
38 | field_matrix[nz+1,*] = field[nz-1,*] |
---|
39 | |
---|
40 | ;print, ' time_matrix : ', time_matrix[*,10] |
---|
41 | ;print, ' z_matrix : ', z_matrix[*,10] |
---|
42 | ;print, ' field_matrix : ', field_matrix[*,10] |
---|
43 | |
---|
44 | ; |
---|
45 | ;---------- |
---|
46 | ; Contours |
---|
47 | ;---------- |
---|
48 | ; colors |
---|
49 | LOADCT, ct |
---|
50 | color00 = MAX(colors) |
---|
51 | color99 = MIN(colors) |
---|
52 | |
---|
53 | CONTOUR, field_matrix, time_matrix, - z_matrix, /FILL, /OVERPLOT, C_COLORS = colors, levels = levels |
---|
54 | |
---|
55 | ;CONTOUR, s_i, time_matrix, - z_ip, /OVERPLOT, levels = levels, $ |
---|
56 | ;C_LABELS=[ 0, 1, 0, 0, 1, 1, 0, 0, 0, 1, 0, 0, 1], $ |
---|
57 | ;C_ANNOTATION=['0','1','','','4','5','','','','9','','','12'] |
---|
58 | |
---|
59 | end |
---|