Regarding your second issue: The negative water (or other species) values probably comes from the fact that your retrieved atmosphere contains such values. The Q.ABS_SPECIES(2).ATMDATA field is only you apriori (1st guess) atmosphere, the retrieved result might, and should, differ from this.
Negative water vapour is not necessarily a problem in itself, depending on the amount of noise in your data, but to avoid the arts error message there are two options. 1) allow for negative water vapour in your simulations by setting Q.VMR_NEGATIVE_OK = true 2) force a min/max constraint on the vmr by using the Q.ABS_SPECIES.MINMAX field. For example Q.ABS_SPECIES.MINMAX = [1e-12 1e-3]. This means that whenever your retrievals create water vapour below 1e-12 kg/m3 or above 1e-3 kg/m3 Qpack/Qarts will change this value to the min/max value before returning it to ARTS. I hope this helps Ole Martin ________________________________ From: Ole Martin Christensen Sent: Tuesday, September 09, 2014 2:40 PM To: Ford, George P.; MARTINET Pauline; [email protected] Subject: RE: [Qpack] retrieval of humidity profiles with qpack and multi viewing angles Regarding your second issue: The negative water (or other species) values probably comes from the fact that your retrieved atmosphere contains such values. The Q.ABS_SPECIES(2).ATMDATA field is only you apriori (1st guess) atmosphere, the retrieved result might, and should, differ from this. Negative water vapour is not necessarily a problem in itself, depending on the amount of noise in your data, but to avoid the arts error message there are two options. 1) allow for negative water vapour in your simulations by setting Q.VMR_NEGATIVE_OK = true 2) force a min/max constraint on the vmr by using the Q.ABS_SPECIES.MINMAX field. For example Q.ABS_SPECIES.MINMAX = [1e-12 1e-3]. This means that whenever your retrievals create water vapour below 1e-12 kg/m3 or above 1e-3 kg/m3 Qpack/Qarts will change this value to the min/max value before returning it to ARTS. I hope this helps Ole Martin ________________________________ From: [email protected] [[email protected]] on behalf of Ford, George P. [[email protected]] Sent: Tuesday, September 09, 2014 11:53 AM To: MARTINET Pauline; [email protected] Subject: Re: [Qpack] retrieval of humidity profiles with qpack and multi viewing angles Re. Your first issue, I’ve been using the qpack2_wind3d_demo which already has multiple viewing angles – you could probably modify it for your purposes, or it would at least be a good reference. It appears it appends each spectrum on to the end of the Y vector, rather than creating a new column, hence avoiding the tensor1 issue. From: MARTINET Pauline [mailto:[email protected]] Sent: 09 September 2014 10:30 To: [email protected] Subject: [Qpack] retrieval of humidity profiles with qpack and multi viewing angles Dear all, I have some issues to run humidity profile retrievals with qpack and also to run the retrieval with different viewing angle. I am using the observations from a ground based microwave radiometer with 7 channels and 8 viewing angles. I) Problem with different viewing angles Fistly, I managed to retrieve temperature profiles with qpack but only with one viewing angle. I use a matlab program similar to the qpack_demo.m file. So I modify the q_demo routine to include a second angle in Q.MBLOCK_ZA_GRID ( for example Q.MBLOCK_ZA_GRID= [0 48] instead of 0). I also create the structure Y with a function similar to the y_demo function in the qpack_demo.m file. After creating the Y structure, I copy the brightness temperature observed by the instrument in the Y(1).Y field: Y(1).Y=TB_mesure but this time the size of the array TB_mesure is : 2 * 7 instead of 1* 7 (so 2 angles and 7 channels). When I add this line, I have the error message: Error using rqre_datatype (line 68) Y(1).Y must be tensor1. Error in qpack2>checkY (line 154) rqre_datatype( Y(1).Y, @istensor1, 'Y(1).Y' ); Error in qpack2 (line 52) empty_y = checkY( Y ); I assume this comes from the fact that there are two dimensions instead of 1 due to the addition of the new viewing angle. Does anyone know how we can run a retrieval with multiple viewing angles ? By the way, does the variable Y.ZA refers to the wiewing angles like Q.MBLOCK_ZA_GRID ? (I tried to modifiy this variable to include the second angle also.) II) Problem with the humidity retrieval In my ARTS simulations, I use a humidity profile directly read from a background profile (thus it is not a profile given with the atmlab package). When I declare the water vapor species in the q_demo file: Q.ABS_SPECIES(2).TAG = {'H2O-PWR98'}; Q.ABS_SPECIES(2).RETRIEVE = false; If Q.ABS_SPECIES(2).RETRIEVE is false, I have no problem during the computation of the temperature retrieval (even if I use my own humidity profile for the simulation of the brightness temperature). Then, when I try to retrieve the humidity profile in addition to the temperature profile, I changed the Q.ABS_SPECIES(2).RETRIEVE variable to true, and I have the following error: Run-time error in controlfile: /tmp/atmlab-martinet-tp9a5f653c_5510_4107_adc5_2463807d555b/cfile_y.arts Run-time error in method: atmfields_checkedCalc All values in *vmr_field* must be >= 0. Stopping ARTS execution. Goodbye. However, I checked the vmr_field of the Q.ABS_SPECIES(2).ATMDATA.GRID and all the values look higher than 0 (see attached file) and I already used this humidity profile even if I did not retrieve the humidity without no problem in ARTS. Thus, I don't really understand why I have this error. Does anyone have on idea ? Thank you very much for your help. Best regards, Pauline ----- Météo-France ----- Dr. Pauline Martinet Chercheur CNRM/GMEI/LISA [email protected]<mailto:[email protected]> Fixe : +33 561079031 Site web: https://sites.google.com/site/martinetpauline31/ ________________________________ This message (and any attachments) is for the recipient only. NERC is subject to the Freedom of Information Act 2000 and the contents of this email and any reply you make may be disclosed by NERC unless it is exempt from release under the Act. Any material supplied to NERC may be stored in an electronic records management system.
_______________________________________________ qpack mailing list [email protected] https://www.sat.ltu.se/mailman/listinfo/qpack
