source: trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_CTL.m @ 16

Last change on this file since 16 was 16, checked in by pinsard, 15 years ago

remove trailing blanks, split lines with multiple statements

File size: 3.1 KB
Line 
1%
2
3% ce programme prepare les fichiers d'OLR de simulations IRCAAM CNRM
4% pour le programme de filtrage
5% ncload('rlut_d.AfNQIVIV.197106-09.nc')
6% whos
7% Name        Size                  Bytes  Class     Attributes
8%
9% ans         1x4                     268  cell
10% lat        64x1                     512  double
11% lon       128x1                    1024  double
12% rlut      122x64x128            7995392  double
13% time      122x1                     976  double
14%
15
16close all;
17clear;
18
19olr=zeros(122,29,64,128);
20cd('CTL')
21ncload('rlut_d.CtIV.197106-09.nc','rlut');
22olr(:,1,:,:)=rlut; clear rlut;
23
24ncload('rlut_d.CtIV.197206-09.nc','rlut');
25olr(:,2,:,:)=rlut; clear rlut;
26
27ncload('rlut_d.CtIV.197306-09.nc','rlut');
28olr(:,3,:,:)=rlut; clear rlut;
29
30ncload('rlut_d.CtIV.197406-09.nc','rlut');
31olr(:,4,:,:)=rlut; clear rlut;
32
33ncload('rlut_d.CtIV.197506-09.nc','rlut');
34olr(:,5,:,:)=rlut; clear rlut;
35
36ncload('rlut_d.CtIV.197606-09.nc','rlut');
37olr(:,6,:,:)=rlut; clear rlut;
38
39ncload('rlut_d.CtIV.197706-09.nc','rlut');
40olr(:,7,:,:)=rlut; clear rlut;
41
42ncload('rlut_d.CtIV.197806-09.nc','rlut');
43olr(:,8,:,:)=rlut; clear rlut;
44
45ncload('rlut_d.CtIV.197906-09.nc','rlut');
46olr(:,9,:,:)=rlut; clear rlut;
47
48ncload('rlut_d.CtIV.198006-09.nc','rlut');
49olr(:,10,:,:)=rlut; clear rlut;
50
51ncload('rlut_d.CtIV.198106-09.nc','rlut');
52olr(:,11,:,:)=rlut; clear rlut;
53
54ncload('rlut_d.CtIV.198206-09.nc','rlut');
55olr(:,12,:,:)=rlut; clear rlut;
56
57ncload('rlut_d.CtIV.198306-09.nc','rlut');
58olr(:,13,:,:)=rlut; clear rlut;
59
60ncload('rlut_d.CtIV.198406-09.nc','rlut');
61olr(:,14,:,:)=rlut; clear rlut;
62
63ncload('rlut_d.CtIV.198506-09.nc','rlut');
64olr(:,15,:,:)=rlut; clear rlut;
65
66ncload('rlut_d.CtIV.198606-09.nc','rlut');
67olr(:,16,:,:)=rlut; clear rlut;
68
69ncload('rlut_d.CtIV.198706-09.nc','rlut');
70olr(:,17,:,:)=rlut; clear rlut;
71
72ncload('rlut_d.CtIV.198806-09.nc','rlut');
73olr(:,18,:,:)=rlut; clear rlut;
74
75ncload('rlut_d.CtIV.198906-09.nc','rlut');
76olr(:,19,:,:)=rlut; clear rlut;
77
78ncload('rlut_d.CtIV.199106-09.nc','rlut');
79olr(:,20,:,:)=rlut; clear rlut;
80
81ncload('rlut_d.CtIV.199206-09.nc','rlut');
82olr(:,21,:,:)=rlut; clear rlut;
83
84ncload('rlut_d.CtIV.199306-09.nc','rlut');
85olr(:,22,:,:)=rlut; clear rlut;
86
87ncload('rlut_d.CtIV.199406-09.nc','rlut');
88olr(:,23,:,:)=rlut; clear rlut;
89
90ncload('rlut_d.CtIV.199506-09.nc','rlut');
91olr(:,24,:,:)=rlut; clear rlut;
92
93ncload('rlut_d.CtIV.199606-09.nc','rlut');
94olr(:,25,:,:)=rlut; clear rlut;
95
96ncload('rlut_d.CtIV.199706-09.nc','rlut');
97olr(:,26,:,:)=rlut; clear rlut;
98
99ncload('rlut_d.CtIV.199806-09.nc','rlut');
100olr(:,27,:,:)=rlut; clear rlut;
101
102ncload('rlut_d.CtIV.199906-09.nc','rlut');
103olr(:,28,:,:)=rlut; clear rlut;
104
105ncload('rlut_d.CtIV.200006-09.nc','rlut');
106olr(:,29,:,:)=rlut; clear rlut;
107
108cd('..')
109
110% selction du domaine geographique
111olr1=zeros(122,29,64,128);
112olr1(:,:,:,1:64)=olr(:,:,:,65:128);
113olr1(:,:,:,65:128)=olr(:,:,:,1:64);
114clear olr
115olr=olr1;
116clear olr1;
117olr=olr(:,:,29:43,54:77);
118taille=(43-29+1)*(77-54+1);
119tps=122*29;
120% taille de OLR
121% 122 28 14 24
122olr=reshape(olr,3538,15,24);
123olr=permute(olr,[3 2 1]);
124olr=reshape(olr,360,3538);
125fid=fopen('olr_ctl.dat','w');
126fwrite(fid,olr,'float');
127fclose(fid);
Note: See TracBrowser for help on using the repository browser.