Hi,all I am trying to use qpack to retrieve the AMSUB measurements which Ta polarization is quasi-Vertical pol with the form of Tv*cos(sensor_los)^2+Th*sin(sensor_los)^2.
How to set the parameters of Q structure when I retrieve the atmosphere state of the quasi-Vertical polarization measurement provide? 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; But I failed. If somebody can give me some advise? Thanks very much and best regards Huang Yingzhu
_______________________________________________ qpack mailing list [email protected] https://www.sat.ltu.se/mailman/listinfo/qpack
