diary4cvs_C3PO_yp

 C3PO - diary output for CVS error tracking in C3PO examples 

 Diary output for CVS error tracking in C3PO examples

 INPUT:

    - id_wave: id of the C3PO example
    - dkepath: dkepath structure
    - waves: list of waves

 OUTPUT: none

 By Yves Peysson (CEA-DRFC, yves.peysson@cea.fr) and Joan Decker (CEA-DRFC, joan.decker@cea.fr)

0001 function [] = diary4cvs_C3PO_yp(id_wave,dkepath,waves)
0002 % C3PO - diary output for CVS error tracking in C3PO examples
0003 %
0004 % Diary output for CVS error tracking in C3PO examples
0005 %
0006 % INPUT:
0007 %
0008 %    - id_wave: id of the C3PO example
0009 %    - dkepath: dkepath structure
0010 %    - waves: list of waves
0011 %
0012 % OUTPUT: none
0013 %
0014 % By Yves Peysson (CEA-DRFC, yves.peysson@cea.fr) and Joan Decker (CEA-DRFC, joan.decker@cea.fr)
0015 %
0016 diaryname = ['res_',id_wave,'.log'];
0017 [LUKEversion,matver] = LUKEversion_jd;
0018 %
0019 if exist(diaryname,'file'),
0020     delete(diaryname)
0021 end
0022 %
0023 diary(diaryname)
0024 %
0025 disp(['Simulation id: ',id_wave]);
0026 disp(['Simulation date: ',timeid_jd(timeid_jd(clock))]);
0027 disp(['User id: ',userid_jd]);
0028 disp(['LUKE Version: ',LUKEversion]);
0029 disp(['MATLAB Version: ',matver]);
0030 disp(['MACHINE: ',dkepath.hostname]);
0031 %
0032 for ii = 1:7
0033     for jj = 1:4,
0034         format_disp{ii,jj} = '%17.15f';
0035     end
0036 end
0037 %
0038 if ~iscell(waves), waves = {waves};end
0039 %
0040 for iw = 1:length(waves),
0041     disp(['---------------------- ',waves{iw}.id,' ----------------------'])
0042     %
0043     jj = 1;
0044     %
0045     for ir = [1,10,100,1000],
0046         if length(waves{iw}.rays{1}.sx) >= ir,
0047             disp(['sx(',int2str(ir),')    = ',num2str(waves{iw}.rays{1}.sx(ir),format_disp{1,jj})])
0048             disp(['sy(',int2str(ir),')    = ',num2str(waves{iw}.rays{1}.sy(ir),format_disp{2,jj})])
0049             disp(['sz(',int2str(ir),')    = ',num2str(waves{iw}.rays{1}.sz(ir),format_disp{3,jj})])
0050             disp(['ss(',int2str(ir),')    = ',num2str(waves{iw}.rays{1}.ss(ir),format_disp{4,jj})])
0051             disp(['skrho(',int2str(ir),') = ',num2str(waves{iw}.rays{1}.skrho(ir),format_disp{5,jj})])
0052             disp(['sm(',int2str(ir),')    = ',num2str(waves{iw}.rays{1}.sm(ir),format_disp{6,jj})])
0053             disp(['sn(',int2str(ir),')    = ',num2str(waves{iw}.rays{1}.sn(ir),format_disp{7,jj})])
0054             if jj < 4,disp('----------------------');end
0055         end
0056         jj = jj + 1;
0057     end
0058 end
0059 %
0060 diary off