Hi,all I am trying to use qpack to retrieve the amsub measurements which Ta polarization is the form of Tv*cos(sensor_los)^2+Th*sin(sensor_los)^2,I try to hack the art_oem.m with yf and J as follows:
%--- Run ARTS -------------------------------------------------------- if nargout == 3 do_j = 1; cfile = R.cfile_yj; else do_j = 0; cfile = R.cfile_y; end % arts( cfile ); % y2 = xmlLoad( fullfile( R.workfolder, 'y.xml' ) ); I=y2(1:2:end); Q=y2(2:2:end); tv=I+Q; th=I-Q; y=tv*cos((180-theta)*pi/180)*cos((180-theta)*pi/180)+th*sin((180-theta)*pi/180)*sin((180-theta)*pi/180); % if do_j % Load Jacobian J2 = xmlLoad( fullfile( R.workfolder, 'jacobian.xml' ) ); IJ=J2(1:2:end,:); QJ=J2(2:2:end,:); VJ=IJ+QJ; HJ=IJ-QJ; J=VJ*cos((180-theta)*pi/180)*cos((180-theta)*pi/180)+HJ*sin((180-theta)*pi/180)*sin((180-theta)*pi/180); But I fail,How to solve the trouble like that? Thanks very much Huang Yingzhu
_______________________________________________ qpack mailing list [email protected] https://www.sat.ltu.se/mailman/listinfo/qpack
