source: trunk/xxx_oneplot.pro @ 2

Last change on this file since 2 was 2, checked in by pinsard, 18 years ago

initial import from /usr/work/fvi/OPA/geomag/

File size: 1.6 KB
Line 
1pro xxx_oneplot
2@common
3noerase = 1
4portrait = 0
5; beginning of reading the field to draw
6; field number: 1
7; declaration de la structure extra:
8extra=create_struct('NOTHING',0)
9; declaration de la structure readparameters:
10readparameters=create_struct('FUNCLEC_NAME','read_ncdf','JPIDTA',182,'JPJDTA',149,'JPKDTA',31,'IXMINDTA',0,'IXMAXDTA',181,'IYMINDTA',0,'IYMAXDTA',148,'IZMINDTA',0,'IZMAXDTA',30)
11; on change de lecture?:
12nothing=changeread(readparameters)
13; declaration de la structure meshparameters:
14meshparameters=create_struct('FILENAME','/usr/work/sur/fvi/OPA/geomag/INTERP2/initorca_reg.pro','FILETYPE','batch file','JPIGLO',182,'JPJGLO',149,'JPKGLO',31,'JPI',180,'JPJ',148,'JPK',31,'IXMINMESH',1,'IXMAXMESH',180,'IYMINMESH',0,'IYMAXMESH',147,'IZMINMESH',0,'IZMAXMESH',30,'KEY_SHIFT',0,'KEY_PERIODIQUE',1,'TRIANGULATION',1,'BOUNDARY',[80.0000,439.000])
15; on change de grille?:
16nothing=changegrid(meshparameters)
17; lecture du champ
18if keyword_set(date1in) then date1 = date1in else date1 = 10115
19if keyword_set(date2in) then date2 = date2in else date2 = 10115
20field = litchamp(read_ncdf('Diver', date1, date2,'',parent=1,boite=[80, 440, -79, 90, 0, 5] $
21, filename='/usr/work/sur/fvi/OPA/ORCA2/DivBustar_5d_1993_reg.nc' $
22, _extra=mixstru(ex, extra)))
23field = field
24field = create_struct({tab:field, grille:vargrid, unite:varunit, experience:varexp, nom:varname})
25; end of reading the field to draw
26;
27extra = create_struct('NOTHING',0,'MIN',0.00000,'MAX',0.00000,'INTER',0.00000,'LCT',39)
28plt, field $
29, boite=[80, 440, -79, 90, 0, 5], typein='plt', petit=[1,1,1] $
30, _extra= mixstru(ex, extra), NOERASE=noerase
31return
32end
Note: See TracBrowser for help on using the repository browser.