[arts-users] **SPAM** Re: catch exception in ARTS
Hi, Oliver Thank you very much for your reply! Best regards, Shaofei Wang Replied Message | From | Lemke, Oliver | | Date | 11/20/2023 18:32 | | To | Shaofei Wang | | Cc | ARTS Users List | | Subject | Re: [arts-users] catch exception in ARTS | Hi Shaofei, You can pass the robust=1 option to ybatchCalc to avoid termination on failed batch cases. Cheers, Oliver On 20. Nov 2023, at 11:20, Shaofei Wang wrote: Hi, all When i use ybatchCalc to simulate bt, some errors are triggered and then ARTS terminates. I would like to ask if there is a way to skip the wrong simulation and move on to the next simulation. I've tried converting .arts to.py and adding "try" and "except" in the .py file. However, it seems that try and except must be placed in the ybatch_calc_agenda. Therefore, i set allow_callbacks=True and ws.Touch(y, y_aux, jacobian). But it still can't work correctly. Best regards, Shaofei Wang
[arts-users] **SPAM** catch exception in ARTS
Hi, all When i use ybatchCalc to simulate bt, some errors are triggered and then ARTS terminates. I would like to ask if there is a way to skip the wrong simulation and move on to the next simulation. I've tried converting .arts to.py and adding "try" and "except" in the .py file. However, it seems that try and except must be placed in the ybatch_calc_agenda. Therefore, i set allow_callbacks=True and ws.Touch(y, y_aux, jacobian). But it still can't work correctly. Best regards, Shaofei Wang
[arts-users] **SPAM** Problem for surface_rtprop_agenda in scattering
Hi, all When I simulated the effect of cloud on the microwave bt, I encountered a error. The ARTS version is 2.5.11, stokes_dim=2. The option of surface_rtprop_agenda used is "Specular_WithPol_ReflFix_SurfTFromt_surface". I have 10 frequencies, and I entered the corresponding surfaceFlatReflectivity (Dimensions = [10, 2, 2]). However, i got the following error for both RT4 and DOIT. The error of RT4 is: - RT4Calc Could not increase nstreams sufficiently (current: 32) to satisfy scattering matrix norm at f[4]=89 GHz. Continuing with nstreams=32. Output for this frequency might be erroneous.} This run took 86.72s Run-time error in controlfile: cloud.arts Run-time error in method: RT4Calc Run-time error in agenda: surface_rtprop_agenda Run-time error in method: surfaceFlatReflectivity The number of pages in *surface_reflectivity* should match length of *f_grid* or be 1. length of *f_grid* : 1 dimension of *surface_reflectivity* : 10 and the error of DOIT is: User Error: failed Error is found at: Filename: "/home/wsf/arts/src/m_doit.cc" Function Name: "void DoitCalc(Workspace&, Tensor7&, const Index&, const Index&, const Index&, const Index&, const Index&, const Vector&, const Agenda&, const Index&, const Verbosity&)" Line Number: 2832 Column Number: 5 Please follow these instructions to correct your error: Error for f_index = 0 (1.065e+10 Hz) Run-time error in agenda: doit_mono_agenda Run-time error in method: cloudbox_field_monoIterate Run-time error in agenda: doit_rte_agenda Run-time error in method: cloudbox_fieldUpdateSeq1D Run-time error in agenda: surface_rtprop_agenda Run-time error in method: surfaceFlatReflectivity The number of pages in *surface_reflectivity* should match length of *f_grid* or be 1. length of *f_grid* : 1 dimension of *surface_reflectivity* : 10 Any one can help this? Best regards, Shaofei Wang
[arts-users] **SPAM** Re: Zeeman effect of complete oxygen models
Hi Richard, Thank you very much for your reply. Best regards, Shaofei Wang Replied Message | From | Richard Larsson | | Date | 10/8/2023 14:29 | | To | Shaofei Wang | | Cc | arts users mi | | Subject | Re: [arts-users] **SPAM** Zeeman effect of complete oxygen models | Dear Shaofei, No. If you are only in the upper atmosphere, you need to use line-by-line when considering the Zeeman Effect. The full models fail there anyways because they don't consider Doppler effects or Zeeman effects. If your beam leaves the upper atmosphere, you need to do additional work. With hope, Richard On Sun, Oct 8, 2023, 04:19 Shaofei Wang wrote: Hi all, Recently, I've been doing some work on the BT simulations for 50-60 GHz. I have a question. Can I just use complete oxygen models provided by ARTS when considering the Zeeman effect for upper atmosphere? Or do I have to use the absorption line and the required line parameters to calculate the Zeeman effect? Best regards, Shaofei Wang
[arts-users] **SPAM** Zeeman effect of complete oxygen models
Hi all, Recently, I've been doing some work on the BT simulations for 50-60 GHz. I have a question. Can I just use complete oxygen models provided by ARTS when considering the Zeeman effect for upper atmosphere? Or do I have to use the absorption line and the required line parameters to calculate the Zeeman effect? Best regards, Shaofei Wang
Re: [arts-users] Question about Fastem
Dear Patrick, Thank you very much for your reply! I'm sorry that I can't make some suggestions for modifying surfaceFASTEM for the short time, since I'm a beginner in Vector Radiation Transfer. In addition, I would like to ask whether e+r≈1 is true for U and V stokes parameters. Best regards, Shaofei Wang Replied Message | From | Patrick Eriksson | | Date | 2/23/2023 21:29 | | To | Shaofei Wang , arts_users.mi@lists.uni-hamburg.de | | Subject | Re: [arts-users] Question about Fastem | Dear Shaofei Wang, I had a quick look at the code, and yes FASTEM values matching U and V are not used. Despite I stand as author for the FASTEM methods, I can answer why we left it like this. There are several possibilities: 1. It was left to be fixed later. 2. We were unsure about how to map those values to surface_rmatrix. 3. We did not trust FASTEM for those values I think it was a combination of 2 and 3. If you have looked into this, do you have a suggestion for how to modify surfaceFASTEM to fix this? Please notice that Section 6.9.3 in ARTS theory guide deals with this. Notice especially Eq 6.79. If you can relate what FASTEM outputs for U and V to that equation, it would be easy to adopt the code. Bye, Patrick On 2023-02-23 13:28, Shaofei Wang wrote: Hi all, Fastem provided by ARTS can calculate the sea emissivities of H, V, the third Stokes component, and the fourth Stokes component. However, the source code of WSM surfaceFastem shows that the reflection coefficients of the third and the fourth stokes components are equal to (rv+rh)/2, which confuses me. Why are they not equal to the r3/r4 calculated by the Fastem? For surfaceFastem, surface_rmatrix seems to be calculated in the same way as when only rv/rh is available (if I understand correctly) Best regards, Shaofei Wang
[arts-users] Question about Fastem
Hi all, Fastem provided by ARTS can calculate the sea emissivities of H, V, the third Stokes component, and the fourth Stokes component. However, the source code of WSM surfaceFastem shows that the reflection coefficients of the third and the fourth stokes components are equal to (rv+rh)/2, which confuses me. Why are they not equal to the r3/r4 calculated by the Fastem? For surfaceFastem, surface_rmatrix seems to be calculated in the same way as when only rv/rh is available (if I understand correctly) Best regards, Shaofei Wang
Re: [arts-users] the absorption by particles
Dear Patrick, Thank you very much for your reply and reminding! Best regards, Shaofei Wang | | | | | On 10/24/2022 19:12,Patrick Eriksson wrote: Dear Shaofei Wang, rain-MPM93 works in a stand-alone fashion. That is, you can forget about ScatElementsToabs_speciesAdd. But please note that you have to be at very low frequencies (<10 GHz?) to be in a situation where neglecting scattering by rain is a fair approximation. Regards, Patrick On 2022-10-24 11:09, Shaofei Wang wrote: Hi, all I want to use ARTS-2.4.0 to calculate the the absorption by particles, not the scattering. In the arts user guide 6.5.8, it said that the ARTS method for this is ScatElementsToabs_speciesAdd. If i want to use the condensate absorption models (e.g., rain-MPM93), whether i still need the method ScatElementsToabs_speciesAdd or i just need the add "rain" to the abs_species, and read the rain profile by AtmRawRead. I've tried the latter, and it does produce different results from the clear sky, but I don't know if this setup is correct. Can anyone help this? Best regards, Shaofei Wang
[arts-users] the absorption by particles
Hi, all I want to use ARTS-2.4.0 to calculate the the absorption by particles, not the scattering. In the arts user guide 6.5.8, it said that the ARTS method for this is ScatElementsToabs_speciesAdd. If i want to use the condensate absorption models (e.g., rain-MPM93), whether i still need the method ScatElementsToabs_speciesAdd or i just need the add "rain" to the abs_species, and read the rain profile by AtmRawRead. I've tried the latter, and it does produce different results from the clear sky, but I don't know if this setup is correct. Can anyone help this? Best regards, Shaofei Wang
Re: [arts-users] Some questions about WSMS
Hi, Patrick Thank you very much for your reply. This is exactly what i need. Best Regards, Shaofei Wang At 2022-02-26 02:00:04, "Patrick Eriksson" wrote: >Hi, > >1. I think AtmFieldPRegrid is what you are looking for > >2. WriteXMLIndexed should solve your problem here > >Bye, > >Patrick > > >On 2022-02-25 08:28, Shaofei Wang wrote: >> Hi, all >> (1) Recently I was using ARTS to simulate scattering effects and I >> generated pnd_field (Tensor4) using WSM >> pnd_fieldCalcFromParticleBulkProps. The pressure grid of the IWC is not >> consistent with my p_grid. Therefore, i need to create >> ArrayOfGriddedField3 and convert pnd_field from Tensor4 to >> ArrayOfGriddedField3 , and then i can >> use pnd_fieldCalcFrompnd_field_raw. However, I didn't find any method to >> set the ArrayOfGriddedField3 (grid and data), only create. Is there any >> way to achieve this conversion? >> (2) Is there any way to achieve conversion from Index/Numeric to >> string. This would help a lot with the naming of the files I write out >> while looping. >> >> Best Regards, >> Shaofei Wang >> >> >> >> ___ >> 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] Some questions about WSMS
Hi, all (1) Recently I was using ARTS to simulate scattering effects and I generated pnd_field (Tensor4) using WSM pnd_fieldCalcFromParticleBulkProps. The pressure grid of the IWC is not consistent with my p_grid. Therefore, i need to create ArrayOfGriddedField3 and convert pnd_field from Tensor4 to ArrayOfGriddedField3 , and then i can use pnd_fieldCalcFrompnd_field_raw. However, I didn't find any method to set the ArrayOfGriddedField3 (grid and data), only create. Is there any way to achieve this conversion? (2) Is there any way to achieve conversion from Index/Numeric to string. This would help a lot with the naming of the files I write out while looping. Best Regards, Shaofei Wang___ arts_users.mi mailing list arts_users.mi@lists.uni-hamburg.de https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi