source: trunk/SRC/Tests/tst_plt.pro

Last change on this file was 495, checked in by pinsard, 10 years ago

fix thanks to coding rules; typo; dupe empty lines; trailing blanks

  • Property svn:keywords set to Id
File size: 3.5 KB
Line 
1;+
2;
3; @file_comments
4;
5; @categories
6;
7; @keyword IMAGE
8;
9; @examples
10;
11; To run test with Levitus temperature on a regular 1x1 grid :
12;   IDL> @tst_initlev
13;   IDL> tst_plt
14;
15; @history
16; Sebastien Masson
17;
18; @version
19; $Id$
20;
21;-
22PRO tst_plt, IMAGE=image
23;
24  compile_opt idl2, strictarrsubs
25;
26@cm_4mesh
27;
28; figure 2: regular grid
29;
30  CASE jpi*key_stride[0]+jpj*key_stride[1] OF
31    180+148:file = 'Levitus98_1m_01_12_Temperature_Pot_ORCA2.nc'
32    360+180:file = 'Levitus98_1m_01_12_Temperature_Pot_1x1.nc'
33    720+510:file = 'Levitus98_1m_01_12_Temperature_Pot_ORCA05.nc'
34  ENDCASE
35  file = isafile(file, title = 'Where is '+file+' ?', /lookalldir)
36  IF size(file, /type) NE 7 THEN return
37;
38  dummy = ''
39  cnt = 1
40;
41  domdef
42  domdef, gdept[0], gdept[0], gridtype = 'T'
43  temp = read_ncdf('votemper', 00101, 00131, file = file)
44;
45  plt, temp, /landscape
46  if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
47  ELSE read, dummy,  prompt = 'press return for the next plot'
48  cnt = cnt + 1
49;
50  plt, temp, -2, 31, int = .5, /nocontour, format = '(i3)'
51  if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
52  ELSE read, dummy,  prompt = 'press return for the next plot'
53  cnt = cnt + 1
54;
55  plt, temp, -2, 31, int = 1, /color_c
56  if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
57  ELSE read, dummy,  prompt = 'press return for the next plot'
58  cnt = cnt + 1
59;
60  plt, temp, -2, 31, int = 1, /nofill, c_thick = 1, coast_thick = 2
61  if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
62  ELSE read, dummy,  prompt = 'press return for the next plot'
63  cnt = cnt + 1
64;
65  plt, temp, cell_fill = 1 + (jpi EQ 180) ; we must use cell_fill = 2 for ORCA2
66  if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
67  ELSE read, dummy,  prompt = 'press return for the next plot'
68  cnt = cnt + 1
69; zoom
70  IF key_onearth THEN box = [40, 375, -20, 20] $
71  ELSE box = [jpi/4, 3*jpi/4, jpj/4, 3*jpj/4]
72  plt, temp, boxzoom = box, /landscape
73  if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
74  ELSE read, dummy,  prompt = 'press return for the next plot'
75  cnt = cnt + 1
76
77; projections
78  IF key_onearth THEN BEGIN
79;
80    plt, temp, boxzoom = [20, 380, 60, 90], /stereo, map = [90, 0, 0], /portrait
81    if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
82    ELSE read, dummy,  prompt = 'press return for the next plot'
83    cnt = cnt + 1
84;
85    plt, temp, boxzoom = [20, 380, -90, 50], /ortho, map = [-90, 180, 0], /portrait
86    if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
87    ELSE read, dummy,  prompt = 'press return for the next plot'
88    cnt = cnt + 1
89;
90    plt, temp, /ortho, map = [0, 0, 21], /portrait, /realcont
91    if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png $
92    ELSE read, dummy,  prompt = 'press return for the next plot'
93    cnt = cnt + 1
94;
95  ENDIF
96
97; deep plot
98; domdef, 150, 150 ; => warning message
99  domdef, gdept[jpk/2], gdept[jpk/2], gridtype = 'T'
100  temp = read_ncdf('votemper', 00101, 00131, file = file)
101  plt, temp, realcont = 2*(key_onearth), /land
102  if keyword_set(image) then saveimage, 'tst_plt_'+image+'_'+string(cnt, format='(i2.2)')+'.png', /png
103  cnt = cnt + 1
104
105  return
106end
Note: See TracBrowser for help on using the repository browser.