Re: [arts-users] Zeeman effect

2020-11-25 Thread witali.krochin
Dear Mr. Eriksson


Thank you for the answer.

Since i had my script in qpack, I just generated abs_lines in ARTS 2.4 in a 
short way to demonstrate my setup.


The code was


Arts2{
  abs_speciesSet(species=["O2-Z-66-53.0469e9-53.0869e9","O2-66"])
  VectorCreate(gs)
  ArrayOfQuantumIdentifierCreate(qid)
  ReadXML(gs, "/arts-xml-data-2.4.0/spectroscopy/zeeman/O2-66.g.xml")
  ReadXML(qid, "/arts-xml-data-2.4.0/spectroscopy/zeeman/O2-66.qid.xml")
  ReadHITRAN(abs_lines,"HITRAN2012.par")
  abs_linesSetZeemanCoefficients(abs_lines, qid, gs)
 WriteXML( "ascii", abs_lines, abs_lines.xml")

It seems to have some 'nan' entries. I compared with your example from 
18.11.2020 and an other  difference is that the last 4 numbers are missed.

52021456295.2042 1.82478855875538e-21 2.82227115907551e-20 61 63 5.622e-10 
0.00213411290322581 -0.0645374497084855 9172.03671157167 0.71 0 0 0 0.71 0 0 
9467.90886355786 0.71 0 0 0 0.71 0 0

52073500265.913 1.02264066638734e-24 5.88280179901669e-20 61 63 5.64e-10 nan 
nan 9172.03671157167 0.71 0 0 0 0.71 0 0 9467.90886355786 0.71 0 0 0 0.71 0 0

52542435628.7166 4.0664650946183e-21 2.47600404813118e-20 57 59 5.79e-10 
0.00241711264367816 -0.0689937941931067 9172.03671157167 0.71 0 0 0 0.71 0 0 
9467.90886355786 0.71 0 0 0 0.71 0 0

52594539557.917 2.25794239624139e-24 5.54036634348622e-20 57 59 5.808e-10 nan 
nan 9172.03671157167 0.71 0 0 0 0.71 0 0 9467.90886355786 0.71 0 0 0 0.71 0 0

53066952513.2334 8.54258889506482e-21 2.15218198757136e-20 53 55 5.967e-10 
0.00276410582010582 -0.0741096869563658 9172.03671157167 0.71 0 0 0 0.71 0 0 
9467.90886355786 0.71 0 0 0 0.71 0 0

53119086421.6796 4.69902633155893e-24 5.22012465293933e-20 53 55 5.982e-10 nan 
nan 9172.03671157167 0.71 0 0 0 0.71 0 0 9467.90886355786 0.71 0 0 0 0.71 0 0

53595786409.1454 1.68833606335819e-20 1.85085761821127e-20 49 51 6.144e-10 
0.00319595076923077 -0.0800431763685832 9172.03671157167 0.71 0 0 0 0.71 0 0 
9467.90886355786 0.71 0 0 0 0.71 0 0

Deos line mixing affect this part? I planed to care about that as soon as the 
quantum numbers are fixed.

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


Re: [arts-users] Zeeman effect

2020-11-25 Thread Richard Larsson
Hi Witali,

The error means that your lines' quantum numbers are defined in a way the
ARTS Zeeman module cannot understand (we need F or J to do these
computations).  If you print the line catalog, you will probably see
several "0/0" in there for the lines you have defined as using Zeeman
effect.  Is this correct?

If it is correct, the steps to fix this is to limit where you activate the
Zeeman effect, and to define the quantum numbers there correctly.  In your
"abs_species", however you have defined it, you should limit the Zeeman
species to just the few lines you wish to compute it for.  This is often
done by simply adding "-FLOW-FHIGH" at the end of a species, and then have
the species without the "-Z-" defined after the limited species in a
different entry.  As an example:

abs_speciesSet(species=["O2-Z-66-50e9-55e9", "O2-66"])

Now the abs_lines_per_species variable that lbl_checkedCalc takes should
have an entry with only the few lines between 50 and 55 GHz.  All of these
should already have their quantum numbers defined correctly in all standard
ways to define them.  If you still encounter this problem, you need to
manually add the quantum numbers because your original catalog is not
well-understood...

(Note that you must ensure that you have line mixing defined in your line
catalog for O2 16-16 or the above tags will anyways just produce nonsense
around the lines.)

If it is not correct that you see many 0/0's, please tell me what it looks
like and I will try to figure out the problem.

With hope,
//Richard

Den tis 24 nov. 2020 kl 12:43 skrev :

>
> Hi everyone,
>
>
> I am still working on the Zeeman effect in ARTS 2.4.
>
> I have set the quantum constants with "abs_linesSetZeemanCoefficients",
> but still get the error message:
>
>
> "Run-time error in method: lbl_checkedCalc
>
>   Bad upper state F(s) or J(s)."
>
> What can that mean?
>
> Best regards
> Witali Krochin
>
>
>
> ___
> 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