- Timestamp:
- 02/02/09 17:55:24 (15 years ago)
- Location:
- trunk/src/SIMULS_IRCAAM
- Files:
-
- 6 edited
- 1 copied
Legend:
- Unmodified
- Added
- Removed
-
trunk/src/SIMULS_IRCAAM/carte_eof_df.m
r52 r53 30 30 % 31 31 % Following line read ``eof_arpege_AfNQIVIV_10-30d.mat`` and 32 % ``rlut_d. afNQIVIV.197106-09.nc``,32 % ``rlut_d.AfNQIVIV.197106-09.nc``, 33 33 % compute ++, 34 34 % plot ++, … … 65 65 % - fplod 2009-02-02T13:45:57Z aedon.locean-ipsl.upmc.fr (Darwin) 66 66 % 67 % * created from carte_eof_ircaam_1030.m .mto replace it as well as67 % * created from carte_eof_ircaam_1030.m to replace it as well as 68 68 % carte_eof_ircaam_1030_as.m, carte_eof_ircaam_1030_ctl.m, 69 69 % carte_eof_ircaam_1030_ctl_clim.m, carte_eof_ircaam_1030_tr.m, -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_AFR.m
r29 r53 1 1 % 2 2 %+ 3 % 4 % EVOLUTIONS 5 % =========== 6 % 7 % $Id$ 8 % 9 % - fplod 2009-02-02T15:39:22Z aedon.locean-ipsl.upmc.fr (Darwin) 10 % 11 % * replaced by prepare_olr_filtre_simulation.m 12 % 13 %- 3 14 % ce programme prepare les fichiers d'OLR de simulations IRCAAM CNRM 4 15 % pour le programme de filtrage -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_AS.m
r29 r53 1 1 % 2 2 %+ 3 % 4 % EVOLUTIONS 5 % =========== 6 % 7 % $Id$ 8 % 9 % - fplod 2009-02-02T15:39:22Z aedon.locean-ipsl.upmc.fr (Darwin) 10 % 11 % * replaced by prepare_olr_filtre_simulation.m 12 % 13 %- 3 14 % ce programme prepare les fichiers d'OLR de simulations IRCAAM CNRM 4 15 % pour le programme de filtrage -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_CTL.m
r29 r53 1 1 % 2 2 %+ 3 % 4 % EVOLUTIONS 5 % =========== 6 % 7 % $Id$ 8 % 9 % - fplod 2009-02-02T15:39:22Z aedon.locean-ipsl.upmc.fr (Darwin) 10 % 11 % * replaced by prepare_olr_filtre_simulation.m 12 % 13 %- 3 14 % ce programme prepare les fichiers d'OLR de simulations IRCAAM CNRM 4 15 % pour le programme de filtrage -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_CTL_clim.m
r29 r53 1 1 % 2 2 %+ 3 % 4 % EVOLUTIONS 5 % =========== 6 % 7 % $Id$ 8 % 9 % - fplod 2009-02-02T15:39:22Z aedon.locean-ipsl.upmc.fr (Darwin) 10 % 11 % * replaced by prepare_olr_filtre_simulation.m 12 % 13 %- 3 14 % ce programme prepare les fichiers d'OLR de simulations IRCAAM CNRM 4 15 % pour le programme de filtrage -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_TR.m
r29 r53 1 1 % 2 2 %+ 3 % 4 % EVOLUTIONS 5 % =========== 6 % 7 % $Id$ 8 % 9 % - fplod 2009-02-02T15:39:22Z aedon.locean-ipsl.upmc.fr (Darwin) 10 % 11 % * replaced by prepare_olr_filtre_simulation.m 12 % 13 %- 3 14 % ce programme prepare les fichiers d'OLR de simulations IRCAAM CNRM 4 15 % pour le programme de filtrage -
trunk/src/SIMULS_IRCAAM/prepare_olr_filtre_simulation.m
r29 r53 1 function prepare_olr_filtre_simulation(ircaam_dataset,simulation) 2 % prepare_olr_filtre_simulation ++explication ircaam_dataset,simulation 3 %+ 4 % module 5 % ====== 1 6 % 2 3 7 % ce programme prepare les fichiers d'OLR de simulations IRCAAM CNRM 4 8 % pour le programme de filtrage 9 % 10 % DESCRIPTION 11 % =========== 12 % 13 % From rlut_d.\ *simulation*\ .197106-09.nc, 14 % ``prepare_olr_filtre_simulation`` compute ++. 15 % 16 % ``prepare_olr_filtre_simulation`` write 17 % olr_nofiltre_arpege_ \ *simulation*\ .dat. 18 % 19 % CAUTIONS 20 % ======== 21 % 22 % can't work with octave because of ncload. 23 % 24 % EXAMPLES 25 % ======== 26 % 27 % Following line read from ``rlut_d.AfNQIVIV.197106-09.nc`` to 28 % ``rlut_d.AfNQIVIV.200006-09.nc``, 29 % compute ++, 30 % write ``olr_nofiltre_arpege_AfNQIVIV.dat`` 31 % :: 32 % 33 % >> prepare_olr_filtre_simulation('arpege','AfNQIVIV') 34 % 35 % SEE ALSO 36 % ======== 37 % 38 % carte_eof_df.m_ 39 % 40 % .. _carte_eof_df.m : carte_eof_df.m.html 41 % 42 % TODO 43 % ==== 44 % 45 % improve description 46 % 47 % improve file pb 48 % 49 % use return function 50 % 51 % EVOLUTIONS 52 % ========== 53 % 54 % $Id$ 55 % 56 % - fplod 2009-02-02T13:45:57Z aedon.locean-ipsl.upmc.fr (Darwin) 57 % 58 % * created from prepare_olr_filtre_AFR.m to replace it as well as 59 % prepare_olr_filtre_AS.m, prepare_olr_filtre_CTL.m, 60 % prepare_olr_filtre_CTL_clim.m and prepare_olr_filtre_TR.m 61 % 62 % Dataset ``arpege`` will be forced in output file names. 63 % 64 % Simulation name will be add to output file names (see parameters). 65 % 66 % 1990 is added to years set. 67 %- 5 68 6 % initialisation 69 global IRCAAM_ID; 70 global IRCAAM_OD; 71 72 if nargin==2 73 % ircaam_dataset must be equal to arpege 74 if ~strcmp(ircaam_dataset,'arpege') 75 disp(['ircaam_dataset = ', ircaam_dataset]); 76 error('eee : ircaam_dataset must be equal to arpege'); 77 end 78 if strcmp(ircaam_dataset,'arpege') 79 list_simulation={'AfNQIVIV','TrNQIVIV','AsNQIVIV','CtIV','CtCl'}; 80 end 81 switch simulation 82 case list_simulation 83 % ok 84 otherwise 85 disp(['simulation = ', simulation]); 86 disp(['list_simulation = ', mat2str(list_simulation)]); 87 error('eee : simulation must be equal to one of list_simulation'); 88 end 89 clear list_simulation; 90 else 91 usage='usage : prepare_olr_filtre_simulation(ircaam_dataset,simulation)'; 92 disp(usage); 93 error('eee : wrong arguments numbers') 94 end 95 7 96 close all; 8 clear;9 97 10 olr=zeros(122,29,64,128); 98 iyearmin=1971; 99 iyearmax=2000; 100 nbyear=iyearmax-iyearmin+1; 11 101 12 ncload('rlut_d.AfNQIVIV.197106-09.nc','rlut'); 13 olr(:,1,:,:)=rlut; clear rlut; 102 olr=zeros(122,nbyear,64,128); 14 103 15 ncload('rlut_d.AfNQIVIV.197206-09.nc','rlut'); 16 olr(:,2,:,:)=rlut; clear rlut; 104 iiyear=0; 105 for iyear=iyearmin:iyearmax 106 cyear=num2str(iyear); 107 fullfilename=[IRCAAM_ID,'rlut_d.',simulation,'.',cyear,'06-09.nc']; 108 ncload(fullfilename,'rlut'); 109 clear fullfilename; 110 iiyear=iiyear+1; 111 olr(:,iiyear,:,:)=rlut; 112 clear rlut; 113 end; 17 114 18 ncload('rlut_d.AfNQIVIV.197306-09.nc','rlut'); 19 olr(:,3,:,:)=rlut; clear rlut; 20 21 ncload('rlut_d.AfNQIVIV.197406-09.nc','rlut'); 22 olr(:,4,:,:)=rlut; clear rlut; 23 24 ncload('rlut_d.AfNQIVIV.197506-09.nc','rlut'); 25 olr(:,5,:,:)=rlut; clear rlut; 26 27 ncload('rlut_d.AfNQIVIV.197606-09.nc','rlut'); 28 olr(:,6,:,:)=rlut; clear rlut; 29 30 ncload('rlut_d.AfNQIVIV.197706-09.nc','rlut'); 31 olr(:,7,:,:)=rlut; clear rlut; 32 33 ncload('rlut_d.AfNQIVIV.197806-09.nc','rlut'); 34 olr(:,8,:,:)=rlut; clear rlut; 35 36 ncload('rlut_d.AfNQIVIV.197906-09.nc','rlut'); 37 olr(:,9,:,:)=rlut; clear rlut; 38 39 ncload('rlut_d.AfNQIVIV.198006-09.nc','rlut'); 40 olr(:,10,:,:)=rlut; clear rlut; 41 42 ncload('rlut_d.AfNQIVIV.198106-09.nc','rlut'); 43 olr(:,11,:,:)=rlut; clear rlut; 44 45 ncload('rlut_d.AfNQIVIV.198206-09.nc','rlut'); 46 olr(:,12,:,:)=rlut; clear rlut; 47 48 ncload('rlut_d.AfNQIVIV.198306-09.nc','rlut'); 49 olr(:,13,:,:)=rlut; clear rlut; 50 51 ncload('rlut_d.AfNQIVIV.198406-09.nc','rlut'); 52 olr(:,14,:,:)=rlut; clear rlut; 53 54 ncload('rlut_d.AfNQIVIV.198506-09.nc','rlut'); 55 olr(:,15,:,:)=rlut; clear rlut; 56 57 ncload('rlut_d.AfNQIVIV.198606-09.nc','rlut'); 58 olr(:,16,:,:)=rlut; clear rlut; 59 60 ncload('rlut_d.AfNQIVIV.198706-09.nc','rlut'); 61 olr(:,17,:,:)=rlut; clear rlut; 62 63 ncload('rlut_d.AfNQIVIV.198806-09.nc','rlut'); 64 olr(:,18,:,:)=rlut; clear rlut; 65 66 ncload('rlut_d.AfNQIVIV.198906-09.nc','rlut'); 67 olr(:,19,:,:)=rlut; clear rlut; 68 69 ncload('rlut_d.AfNQIVIV.199106-09.nc','rlut'); 70 olr(:,20,:,:)=rlut; clear rlut; 71 72 ncload('rlut_d.AfNQIVIV.199206-09.nc','rlut'); 73 olr(:,21,:,:)=rlut; clear rlut; 74 75 ncload('rlut_d.AfNQIVIV.199306-09.nc','rlut'); 76 olr(:,22,:,:)=rlut; clear rlut; 77 78 ncload('rlut_d.AfNQIVIV.199406-09.nc','rlut'); 79 olr(:,23,:,:)=rlut; clear rlut; 80 81 ncload('rlut_d.AfNQIVIV.199506-09.nc','rlut'); 82 olr(:,24,:,:)=rlut; clear rlut; 83 84 ncload('rlut_d.AfNQIVIV.199606-09.nc','rlut'); 85 olr(:,25,:,:)=rlut; clear rlut; 86 87 ncload('rlut_d.AfNQIVIV.199706-09.nc','rlut'); 88 olr(:,26,:,:)=rlut; clear rlut; 89 90 ncload('rlut_d.AfNQIVIV.199806-09.nc','rlut'); 91 olr(:,27,:,:)=rlut; clear rlut; 92 93 ncload('rlut_d.AfNQIVIV.199906-09.nc','rlut'); 94 olr(:,28,:,:)=rlut; clear rlut; 95 96 ncload('rlut_d.AfNQIVIV.200006-09.nc','rlut'); 97 olr(:,29,:,:)=rlut; clear rlut; 98 99 % selction du domaine geographique 100 olr1=zeros(122,29,64,128); 115 % selection du domaine geographique 116 olr1=zeros(122,nbyear,64,128); 101 117 olr1(:,:,:,1:64)=olr(:,:,:,65:128); 102 118 olr1(:,:,:,65:128)=olr(:,:,:,1:64); … … 105 121 clear olr1; 106 122 olr=olr(:,:,29:43,54:77); 107 taille=(43-29+1)*(77-54+1); 108 tps=122*29; 109 % taille de OLR 110 % 122 28 14 24 111 olr=reshape(olr,3538,15,24); 123 tps=size(olr,1)*size(olr,2); 124 olr=reshape(olr,tps,15,24); 112 125 olr=permute(olr,[3 2 1]); 113 olr=reshape(olr,360,3538); 114 for a=1:360; 115 clear vaude vaude1; 116 vaude=olr(a,:); 117 vaude=vaude'; 118 vaude=mean(vaude); 119 olr(a,:)=olr(a,:)-vaude; 120 end; 121 fid=fopen('coco.dat','w'); 126 olr=reshape(olr,360,tps); 127 128 switch simulation 129 case {'CtIV','CtCl'} 130 % nothing 131 otherwise 132 for a=1:360; 133 clear vaude; 134 vaude=olr(a,:); 135 vaude=vaude'; 136 vaude=mean(vaude); 137 olr(a,:)=olr(a,:)-vaude; 138 end; 139 end 140 141 fullfilename=[IRCAAM_OD,'olr_nofiltre_',ircaam_dataset,'_', simulation, '.dat']; 142 fid=fopen(fullfilename,'w'); 143 clear fullfilename; 122 144 fwrite(fid,olr,'float'); 123 145 fclose(fid);
Note: See TracChangeset
for help on using the changeset viewer.