forked from sergio66/UMBC_LBL
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Copy patharmante.m
37 lines (30 loc) · 1018 Bytes
/
armante.m
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
topts.HITRAN = '/asl/data/geisa/g15.by.gas/';; [w,dx] = run8(3,975,1225,'IPFILES/co2one',topts);
[w,d] = run8(3,975,1225,'IPFILES/co2one');
save armante.mat w d dx topts
subplot(211); plot(w,dx,w,d)
title('p=1.1139 atm pO3=4.0418E-03 atm T=269.429K');; ylabel('OD')
grid
subplot(212); plot(w,dx./d)
ylabel('OD geisa/OD hitran')
grid
data = [w; d; dx];
whos data
fid = fopen('armante.txt','w');
fprintf(fid,'%12.5f %12.7e %12.7e \n',data);
fclose(fid);
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
clear all
topts.HITRAN = '/asl/data/geisa/g15.by.gas/';; [w,dx] = run8(3,1200,1275,'IPFILES/co2one',topts);
[w,d] = run8(3,1200,1275,'IPFILES/co2one');
save armante2.mat w d dx topts
subplot(211); plot(w,dx,w,d)
title('p=1.1139 atm pO3=4.0418E-03 atm T=269.429K');; ylabel('OD')
grid
subplot(212); plot(w,dx./d)
ylabel('OD geisa/OD hitran')
grid
data = [w; d; dx];
whos data
fid = fopen('armante2.txt','w');
fprintf(fid,'%12.5f %12.7e %12.7e \n',data);
fclose(fid);