Re: [arts-users] writeOut/readIn error
Thank you for helping out, Oliver! Best wishes, Rita On Wed, 15 Jan 2020, 07:58 Oliver Lemke, wrote: > Hi Rita, > > First of all thank you very much for including a working test case. :-) > > There is indeed a bug in the ReadXML routine for the CovarianceMatrix. > Instead of overwriting the workspace variable with the contents from the > file, it was appending the blocks to the existing covmat_sx, ending up with > 6 blocks after reading instead of 3. > > I will commit a fix soon. > > As a workaround in your version, you can delete the covmat_sx WSV before > reading it: > > WriteXML( "ascii", covmat_sx, "covmat_sx1.xml") > Delete(covmat_sx) > ReadXML( covmat_sx, "covmat_sx1.xml") > > Thanks a lot for reporting this. > > Cheers, > Oliver > > > > On 14. Jan 2020, at 21:17, Rita Kajtar wrote: > > > > Hello! > > > > > > In one of my tests with TestOEM.arts if I add: > > > > WriteXML( "ascii", covmat_sx, "covmat_sx1.xml") > > ReadXML( covmat_sx, "covmat_sx1.xml") > > > > right after the retrieval quantities definition, but before > retrievalDefClose (file attached) in TestOEM.arts, I get the message > included below. > > I have also tried considering a binary format instead of ascii and I get > the same error message. > > > > Is it possible to have a bug in the Write/Read methods? > > > > > > > > > > Best regards, > > Rita Kajtar > > ___ > > arts_users.mi mailing list > > arts_users.mi@lists.uni-hamburg.de > > https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi > > ___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
[arts-users] writeOut/readIn error
Hello! In one of my tests with TestOEM.arts if I add: WriteXML( "ascii", covmat_sx, "covmat_sx1.xml") ReadXML( covmat_sx, "covmat_sx1.xml") right after the retrieval quantities definition, but before retrievalDefClose (file attached) in TestOEM.arts, I get the message included below. I have also tried considering a binary format instead of ascii and I get the same error message. Is it possible to have a bug in the Write/Read methods? [image: Screen Shot 2020-01-14 at 9.02.03 PM.png] Best regards, Rita Kajtar TestOEM.arts Description: Binary data ___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
[arts-users] atm.scenario & cov-mat issues
Hello! I have a few questions, the first ones related to the issues I encounter when I try to use a subarctic atmospheric scenario with user-defined pressure and retrieval pressure grids (*), and the next questions related to what precisely does the covariance matrix that my ARTS file produces contain and how can one actually extract certain blocks out of it once one knows what it contains; the documentation on covmat_sx does not provide too much information on this (**) (*) ### - I am running an adapted version on TestOEM.arts and while the pressure grid and the pressure retrieval grid I set work well with the tropical atm. background that TestOEM.arts comes with, I get the following error when I try to use a subarctic-winter scenario: Run-time error in method: AtmFieldsCalc There is a problem with the grids for the following interpolation: Raw field to p_grid The minimum of the new grid must be inside the original grid. (We allow a bit of extrapolation, but not so much). Minimum of original grid: -2.3969 (0.091) Minimum allowed value for new grid: -2.7956 (0.0610782) Actual minimum of new grid: -2.81347 (0.0599964) - If I then replace the pressure grid with a 'shorter' one, so that I make sure I fall within the minimum allowed value of the grid, I hit an assertion: - xaStandard Assertion failed: (og_min <= tng), function gridpos, file /Users/mypath/arts/src/interpolation.cc, line 345. Abort trap: 6 (**) ## - Why for a simulation containing three retrieval species. i.e. O3, freq.shift and polyfit (no winds considered this time), the retrieval covariance matrix (Sx) contains six retrieval errors? and what are those representing? The terminal output shows a 17 x 17 sized covariance matrix, which corresponds to the 15 elem O3 array + 1 elem of freq-shift + 1 elem of polyfit. (covmat_sx attached) - I can import the generated covariance matrix in my Python code, but how can I extract the different blocks out of it for further use, i.e plotting? With best regards, Rita Kajtar <> ___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
[arts-users] O3 Jacobians
Hello! Since it is hard to trust that the e-mails sent from my LTU account are reaching everybody in the ARTS users list, I will, from now on, use my Gmail account. Thank you to those who helped out with my first problem. Below is another short question I addressed to the community yesterday, but I am not sure who received it or if anyone received it at all. If someone did, my apologies for spamming again. This time my question is related to the calculation of Jacobians for O3. I am considering a no winds, no specific atmospheric background case and the jacobianAddAbsSpecies in my control file is set as it follows: *jacobianAddAbsSpecies( g1=p_grid, g2=gas1_jac_lat_grid, g3=gas1_jac_lon_grid, species="O3", method="analytical", unit="logrel", )* When looking at the change in Bt in my plots, in what units the change of O3 is being represented? Best regards, Rita Kajtar ___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
Re: [arts-users] Test e-mail
This is another test. Sorry for spamming. /Rita On Wed, Apr 3, 2019 at 10:01 PM Rita Edit Kajtar wrote: > Hello! > > > > This is just to check whether e-mails sent from LTU go through :) > > > > /Rita > ___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi