makedisp_simparam_jd

0001 function makedisp_simparam_jd(xpsin,xB,xTe,xne,xNpar,omega,knr,nrt_nr,nrt_fr,kflag,GNperp,nn_nr,tol_ht,nd_ht);
0002 %
0003 if nargin < 14
0004     nd_ht = 11;
0005 end
0006 if nargin < 13
0007     tol_ht = 0.2;
0008 end
0009 %
0010 nr = length(xpsin) - 1;%-1 because the first point is read separately in R2D2
0011 %
0012 rtparam = [knr,nrt_nr,nrt_fr,tol_ht,nd_ht];
0013 %
0014 GNperpRe = real(GNperp);
0015 GNperpIm = imag(GNperp);
0016 %
0017 bt0 = max(xB);
0018 Te0 = max(xTe);
0019 ne0 = max(xne);
0020 param = [bt0*1e4,Te0*1e3,ne0*1e-6,omega/2/pi];
0021 data = [xpsin;xTe/Te0;xne/ne0;xB/bt0;xNpar];
0022 %
0023 s_rparam = smat2f_jd(sprintf('%d',nr));
0024 s_root = smat2f_jd(sprintf('%d,%d,%d,%d,%d',rtparam));
0025 s_nterms = smat2f_jd(sprintf('%d',nn_nr));
0026 s_kflag = smat2f_jd(sprintf('%d',kflag));
0027 s_rootRe = smat2f_jd(sprintf('%d ',GNperpRe));
0028 s_rootIm = smat2f_jd(sprintf('%d ',GNperpIm));
0029 s_param = smat2f_jd(sprintf('%d,%d,%d,%d',param));
0030 s_data = smat2f_jd(sprintf('%d,%d,%d,%d,%d\n',data));
0031 
0032 savestr = ['simparam.m'];
0033 delete(savestr);
0034 
0035 fid = fopen(savestr,'w');
0036 fprintf(fid,'%s\n',s_rparam);
0037 fprintf(fid,'%s\n',s_root);
0038 fprintf(fid,'%s\n',s_nterms);
0039 fprintf(fid,'%s\n',s_kflag);
0040 fprintf(fid,'%s\n',s_rootRe);
0041 fprintf(fid,'%s\n',s_rootIm);
0042 fprintf(fid,'%s\n',s_param);
0043 fprintf(fid,'%s',s_data);
0044 fclose(fid);