Hi Witali, The functions in the README require 2 inputs. You need to give them an ArrayOfQuantumIdentifier, qid, and a Vector, gs. The xml-data/spec/zeeman/ contains these files. Load O2-66.g.xml into a Vector. Load O2-66.qid.xml into an ArrayOfQuantumIdentifier. You now call either of the two provided functions with these variables and the data will be loaded into your line-list. To be sure that you have read them in correctly, and that you have done the right thing, save the line-list and look at your lines.
They should look something like this: 52021404431.109 1.83322270858126e-21 2.82226897398507e-20 61 63 5.647e-10 0.00213411290322581 -0.0645374497084855 10059.6531675302 0.72 0 0 0 0.72 0 0 10266.7636739206 0.72 0 0 0 0.72 0 0 30 31 31 31 52542394856.9423 4.08755046918299e-21 2.47600265761908e-20 57 59 5.82e-10 0.00241711264367816 -0.0689937941931067 10355.5253195164 0.72 0 0 0 0.72 0 0 10296.3508891192 0.72 0 0 0 0.72 0 0 28 29 29 29 53066912640.8365 8.58174744782781e-21 2.15218139163761e-20 53 55 5.994e-10 0.00276410582010582 -0.0741096869563658 10651.3974715026 0.72 0 0 0 0.72 0 0 10355.5253195164 0.72 0 0 0 0.72 0 0 26 27 27 27 53595757329.277 1.69616777391079e-20 1.8508572209221e-20 49 51 6.173e-10 0.00319595076923077 -0.0800431763685832 10947.2696234888 0.72 0 0 0 0.72 0 0 10414.6997499136 0.72 0 0 0 0.72 0 0 24 25 25 25 The numbers of interest are the two before the first 10k number. Note that you need to have the quantum numbers defined in your line catalog to make this work. In ARTS 2.2 this was optional, I think, but there is an additional and strong J dependency on the Zeeman effect numbers as we showed in a paper from last year: "Updated Zeeman effect splitting coefficients for molecular oxygen in planetary applications". We had to change the method to read these or the code would simply be wrong. Note that the lines above are not going to be directly usable by you. You need to add line mixing. There are other methods in ARTS for that but generally your original line database must be tuned for this to begin with or your forward calculations will make no sense anyways. With hope, //Richard Den mån 16 nov. 2020 kl 18:18 skrev <witali.kroc...@iap.unibe.ch>: > > Hi everyone > > I am trying to initialize the Zeeman effect in ARTS 2.4 using the updated > qpack version. > The molecule of interest is oxygen, the frequency range is around > 52-53GHz (the rotational emission lines around 52GHz and 53 GHz). > > I included the method "propmat clearskyAddZeeman" in the propmat clearsky > Agenda. > > However I strugle with setting the quantum numbers and constants. There is > a Zeeman file in the xml Data folder. For me it seems that the required > constants are stored there. The README tell me to use the methods: > > - abs_lines_per_speciesSetZeemanCoefficients or > - abs_linesSetZeemanCoefficients > > I checked the built in documentary, but still not understand how to use > these methods. > For example I don't know in which workspace variable the information > should be feeded in. > > In Arts 2.2 I used readXML to read the file "Zeeman_constants.xml" (was > included in xml data) and feeded the data in the Variable WSMS_AT_START. > > Furhtermore there was an Zeeman demo in the older atmlab version. That was > very helpfull and seems to be missed now. > > Can anyone help me to set this constants? Just a few example lines would > help. > > Best regards > Witali > > > _______________________________________________ > 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