Calculate initial N// et PLH for C3PO from simplified or full power spectrum (from SWAN coupling code for example) by Y.Peysson <yves.peysson@cea.fr> (CEA-DRFC) and J. Decker <joan.decker@cea.fr> (CEA-DRFC)
0001 function [Plhtot,Npar0] = lhspectrum_cronos(directivity,signe,cons,lob) 0002 % 0003 % Calculate initial N// et PLH for C3PO from simplified or full power spectrum (from SWAN coupling code for example) 0004 % 0005 % 0006 %by Y.Peysson <yves.peysson@cea.fr> (CEA-DRFC) and J. Decker <joan.decker@cea.fr> (CEA-DRFC) 0007 % 0008 if strcmp(lob,'Yes') 0009 Plhtot = abs(cons)*[directivity,1-directivity];%initial power in the ray (W) (main lobe) [n_antenna,m_lobe] 0010 Npar0 = [angle(cons),-6.0];%Initial Npar value (main lobe) [n_antenna,m_lobe] 0011 else 0012 Plhtot = abs(cons)*directivity;%initial power in the ray (W) (main lobe) [n_antenna,m_lobe] 0013 Npar0 = angle(cons);%Initial Npar value (main lobe) [n_antenna,m_lobe] 0014 end 0015 % 0016 % For CRONOS, N//> 0 means always co-current, whatever the toroidal Bfield or the direction of Ip 0017 % 0018 Npar0 = -signe.ip*signe.b0*Npar0;%Initial Npar value (main lobe) [n_antenna,m_lobe] 0019 % 0020