Changeset 11 for trunk/private


Ignore:
Timestamp:
08/13/09 18:10:58 (15 years ago)
Author:
pinsard
Message:

make_lut.m : header

Location:
trunk/private
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/private/make_base_image.m

    r10 r11  
    1 %save Thetas Thetav Dphi Rho670 Rho765 Rho865 of one day. 
     1%MAKE_BASE_IMAGE save Thetas Thetav Dphi Rho670 Rho765 Rho865 of one day. 
    22 
    33%+ 

  • trunk/private/make_lut.m

    r2 r11  
     1%MAKE_LUT save ++ 
     2%+ 
     3% module 
     4% ====== 
     5% 
     6% ++ 
     7% 
     8% DESCRIPTION 
     9% =========== 
     10% 
     11% From ../LUT/lut_junge_nir.dat, 
     12% ``make_lut`` produce a ASCII data file ../LUT/LUT_NIR.dat 
     13% containing ++ 
     14% 
     15% EXAMPLES 
     16% ========== 
     17% 
     18% :: 
     19% 
     20% >> addpath('../toolbox') 
     21% >> tp_mlp_aerosols_startup 
     22% >> make_lut 
     23% 
     24% SEE ALSO 
     25% ======== 
     26% 
     27% tp_mlp_aerosols_startup.m_ 
     28% 
     29% .. _tp_mlp_aerosols_startup.m : tp_mlp_aerosols_startup.m.html 
     30% 
     31% TODO 
     32% ==== 
     33% 
     34% improve description 
     35% 
     36% avoid overriding existing output file 
     37% 
     38% make it launchable from everywhere (path and IO directories) 
     39% 
     40% EVOLUTIONS 
     41% ========== 
     42% 
     43% $Id$ 
     44% 
     45% - fplod 2009-08-13T16:05:03Z aedon.locean-ipsl.upmc.fr (Darwin) 
     46% 
     47%    * add header 
     48%    * heaven if no modification is needed by octave now, I use 
     49%      tp_mlp_aerosols_startup in the example to point out IO directories 
     50%      and path future modifications. 
     51%       
     52%- 
     53 
    154infile='../LUT/lut_junge_nir.dat'; 
    255outfile='../LUT/LUT_NIR.dat'; 
    3  
    456 
    557fid=fopen(infile,'r'); 
    … …  
    1870for js=1:length(theta) 
    1971  for jv=1:length(theta) 
    20 Itheta=union(Itheta,find(lut(:,3)==theta(js)&lut(:,4)==theta(jv))); 
     72    Itheta=union(Itheta,find(lut(:,3)==theta(js)&lut(:,4)==theta(jv))); 
    2173  end 
    22    
    2374end 
    2475 
    … …  
    2879 
    2980fid=fopen(outfile,'w'); 
    30 fprintf(fid,'%d %f %f %f %f %f %f\n',[type' ; lut(:,[6 3 4 5 7 ... 
    31                     8])']); 
     81fprintf(fid,'%d %f %f %f %f %f %f\n',[type' ; lut(:,[6 3 4 5 7 ... 8])']); 
    3282fclose(fid); 
Note: See TracChangeset for help on using the changeset viewer.