Hi Patrick,
Thanks for the suggestions.
1.)
=> I tried to check if Dmean was set to 0 value anywhere in the Bulkprop_fields
when using "psdModifiedGammaMassXmean()" . The only place it is '0' is at the
border of the cloudbox and outside the cloudbox.
=> This is how "pnd_fieldCalcFromParticleBulkProps()" expects to be set as
given in the documentation.
=> When I try making the border values a positive non-zero value. It gives me
this error.
Exception: Call to ARTS WSM pnd_fieldCalcFromParticleBulkProps failed with
error: *particle_bulkprop_field* allowed to contain non-zero values only inside
the cloudbox.
2.)
=> So, this looks to me like a bug that if I make borders zero valued then
"psdModifiedGammaMassXmean()" would have a problem, else if I make it non zero
then
"pnd_fieldCalcFromParticleBulkProps()" will give me problems.
=> Am I thinking this right?
3.)
=> Would you recommend me to change the source code in "PSD.CC" and change line
321 "if (ext_pars[1] <= 0) " and remove the "=" sign and build ARTS.
Would this be a fix to the problem?
Thanks,
Renish Thomas
________________________________
From: Patrick Eriksson <patrick.eriks...@chalmers.se>
Sent: Sunday, December 6, 2020 8:25 AM
To: Thomas,Renish <renish.tho...@colostate.edu>
Cc: arts users mi <arts_users.mi@lists.uni-hamburg.de>
Subject: Re: Fwd: Re: [arts-users] ARTS internally generated PND fields
Renish,
Let me first clarify that temperature is not meant to be part of
particle_bulkprop_field.
pnd_fieldCalcFromParticleBulkProps interpolates t_field to set
pnd_agenda_input_t. As this
is done inside a method, the workspace version of pnd_agenda_input_t is left
unset.
Note that the unit of hydrometeor mass contents is kg/m3. It seems that you set
"Mass_Content" to 1 kg/m3. A more realistic value would be 1 g/m3, or lower.
Regarding the use of psdModifiedGammaMassXmean, I don't spot any obvious error.
It could
be a bug. This PSD has not been used a lot, or maybe not all. And you are
likely the first
to combine MC and this PSD system.
I checked the code and even Xmean == 0 will cause the error. Could it be that
the PSD is
called for a point where your Dmean is zero?
To check, can you try to run MC with one core, and add
Print(pnd_agenda_input,0)
just before calling psdModifiedGammaMassXmean (inside of the agenda).
Bye,
Patrick
On 2020-12-05 19:38, Thomas,Renish wrote:
> Hi Patrick / ARTS users,
>
> 1) I am trying to run full calculations using the MC_General() method, after
> setting the
> PSD fields internally.
>
> 2) =>The full calculations do run without errors when I use
> [psdModifiedGammaMass].
> => I see the effects of increase in mass content in the simulation, which
> means the
> bulkprop_fields are read into "ws.pnd_agenda_input".
> =>The only thing is that I don't see is the effect of temperature when after
> the run. Even
> when I remove the name "Temperature" from "ws.particle_bulkprop_names" it runs
> without errors. This makes me wonder if the "pnd_agenda_input_t" is set.
> => And, when I try extracting the "ws.pnd_agenda_input_t.value" after the
> successful run.
> It says that the variable is uninitialized.
>
> 3) When I run with [psdModifiedGammaMassXmean],It gives me the following
> error.
>
> Exception: Call to ARTS WSM pnd_fieldCalcFromParticleBulkProps failed with
> error:
> Run-timeerror in agenda: pnd_agenda_array
>
> Run-time error in method: psdModifiedGammaMassXmean
> Negative mean sizefound.
> This is not allowed.
>
>
>
> 4) I have attached the sub-function python control file for setting the
> scattering fields
> before running MC_General(). Any help would be great.
>
>
> Thanks,
> Renish
>
> ------------------------------------------------------------------------------------------
> *From:* Patrick Eriksson <patrick.eriks...@chalmers.se>
> *Sent:* Saturday, December 5, 2020 1:55 AM
> *To:* Thomas,Renish <renish.tho...@colostate.edu>
> *Cc:* arts users mi <arts_users.mi@lists.uni-hamburg.de>
> *Subject:* Re: Fwd: Re: [arts-users] ARTS internally generated PND fields
> Renish,
>
> Are you doing full calculations, running eg DISORT? If yes, then setting the
> input to
> pnd_agenda should be automatic.
>
> Or are you trying to use the PSD functions separately?
>
> Can you clarify this before I try to answer? And maybe also send me your
> control file.
>
> Bye,
>
> Patrick
>
>
>> *Från:* "Thomas,Renish" <renish.tho...@colostate.edu>
>> *Skickat:* 4 december 2020 16:40:25 CET
>> *Till:* Patrick Eriksson <patrick.eriks...@chalmers.se>,
>> "arts_users.mi@lists.uni-hamburg.de"
>> <arts_users...@mailman.rrz.uni-hamburg.de>
>> *Ämne:* Re: [arts-users] ARTS internally generated PND fields
>>
>> Hi Patrick/ Everyone,
>>
>> Thanks, for the inputs .
>> I was able to clear some of the errors.
>>
>> I have a couple more questions. It will be great to have some input on these.
>>
>> 1.) How do I set "pnd_agenda_input_t"?. Is the name "Temperature" in bulk
>> prop fields
>> names enough to set this, or should I set this explicitly?
>>
>> 2.) Can the particle temperature be different from the atmosphere
>> temperature around it?
>>
>> 3.) I have no problems with running psdModifiedGammaMass, while running it.
>> But while
>> running psdModifiedGammaMassXmean, I get this error.
>>
>> Exception: Call to ARTS WSM pnd_fieldCalcFromParticleBulkProps failed with
>> error: Run-time
>> error in agenda: pnd_agenda_array
>> Run-time error in method: psdModifiedGammaMassXmean
>> Negative mean sizefound.
>> This is not allowed.
>>
>> My inputs are
>> Dmean = 1e-3.
>> n0 = -999
>> mu = 2
>> la = -999
>> ga = 1
>>
>> Could I know why this error occours ?
>>
>> Cheers,
>> Renish
>> ------------------------------------------------------------------------------------------
>> *From:* Patrick Eriksson <patrick.eriks...@chalmers.se>
>> *Sent:* Thursday, November 19, 2020 11:47 AM
>> *To:* Thomas,Renish <renish.tho...@colostate.edu>;
>> arts_users.mi@lists.uni-hamburg.de
>> <arts_users...@mailman.rrz.uni-hamburg.de>
>> *Subject:* Re: [arts-users] ARTS internally generated PND fields
>> Hi Renish,
>>
>> Yes, this part is a complex and could need some more examples. However,
>> this part will be modified in a big overhaul of the handling of
>> scattering properties that Simon is working on and this makes us all a
>> bit reluctant to add examples right now.
>>
>>
>>> 1.) So, the "Mass_content" , "Xmean" and "Temperature" are the only
>>> fields that I need to provide in the particle_bulkprop_field?
>>
>> This depends on the PSDs you are using. In TestScatSolvers.arts this
>> field holds IWC and RWC, which are sufficient for the selected PSDs, but
>> if you want to use psdModifiedGammaMassXmean, you need to include data
>> on the mean size.
>>
>>
>>> 2.) Do I need to set the pnd_agenda_input manually , or does it take
>>> values from the Bulkprop_fields?
>>
>> No, but you need to link the naming in particle_bulkprop_field to the
>> pnd_agenda_input. In your case it could look this
>>
>> ArrayOfStringSet( pnd_agenda_input_names, [ "IWC", "Dmean" ]
>>
>>
>>> 3.) Are there other examples for this method that I can refer to?
>>
>> Which method? (But the answer is likely anyhow no)
>>
>>> 4.) Also, when I run the example "TestScatSolvers.arts" modified to use
>>> PSD from psdModifiedGammaMassXmean(). I get this error
>>>
>>> Exception: Call to ARTS WSM pnd_fieldCalcFromParticleBulkProps failed
>>> with error: *scat_data* and *scat_species* are inconsistent in size.
>>
>> The scattering data in TestScatSolvers.arts cover both IWC and RWC. if
>> you have changed and just work with one hydrometeor type, you need to
>> modify the scattering data accordingly.
>>
>> I hope this was of some help,
>>
>> Patrick
_______________________________________________
arts_users.mi mailing list
arts_users.mi@lists.uni-hamburg.de
https://mailman.rrz.uni-hamburg.de/mailman/listinfo/arts_users.mi
