|
Revision 495
(checked in by cbc, 11 years ago)
|
Initial import of Stark code.
|
| Line | |
|---|
| 1 |
newplot='y'; |
|---|
| 2 |
while newplot == 'y', |
|---|
| 3 |
% |
|---|
| 4 |
% Load data from files |
|---|
| 5 |
% |
|---|
| 6 |
load_adrdep2kmg |
|---|
| 7 |
ls *.s2r |
|---|
| 8 |
filename=input('Enter the name of .s2r file: ','s'); |
|---|
| 9 |
[s2r,gname]=read_s2r(filename); |
|---|
| 10 |
scalar=s2r(:,2); |
|---|
| 11 |
% |
|---|
| 12 |
% Set contour levels and make plot |
|---|
| 13 |
% |
|---|
| 14 |
newcont='y'; |
|---|
| 15 |
while newcont == 'y', |
|---|
| 16 |
scrange(scalar) |
|---|
| 17 |
cint=input('Enter the contour interval: '); |
|---|
| 18 |
cmin=cint*ceil(min(scalar)/cint); |
|---|
| 19 |
cmax=cint*floor(max(scalar)/cint); |
|---|
| 20 |
clear cval |
|---|
| 21 |
i=1; |
|---|
| 22 |
cval(i)=cmin; |
|---|
| 23 |
while cval(i) < cmax |
|---|
| 24 |
i=i+1; |
|---|
| 25 |
cval(i)=cval(i-1)+cint; |
|---|
| 26 |
end |
|---|
| 27 |
cval |
|---|
| 28 |
[cs,h]=lcontour3(in,x,y,scalar,cval); |
|---|
| 29 |
hlabel=extclabel(cs); |
|---|
| 30 |
delete(h) |
|---|
| 31 |
title(filename) |
|---|
| 32 |
zoom on |
|---|
| 33 |
newcont=input('New contour interval? (y/n): ','s'); |
|---|
| 34 |
if newcont == 'y', |
|---|
| 35 |
delete(hlabel); |
|---|
| 36 |
end |
|---|
| 37 |
end |
|---|
| 38 |
newplot=input('New plot? (y/n): ','s'); |
|---|
| 39 |
end |
|---|
