[4] | 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 |
---|