Re: [arts-users] Fwd: Re: ARTS internally generated PND fields

2020-12-06 Thread Thomas,Renish
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 
Sent: Sunday, December 6, 2020 10:59 AM
To: Thomas,Renish 
Cc: arts users mi 
Subject: Re: Fwd: Re: [arts-users] ARTS internally generated PND fields

Ranish,

> 1.) So, if I must set the scattering particle temperature to an arbitrary 
> value after
> calling "pnd_fieldCalcFromParticleBulkProps", do I just override and assign
> 'pnd_agenda_input_t".

You got me wrong. You don't need to do anything. pnd_agenda_input_t is set 
automatically
where needed.


> 2.) Can the particle temperature be set different than the ambient 
> temperature?

Hardly. You can in principle try to "hack" pnd_agenda to set it to some other 
value. But
why do you want to do it?

And please note that the psdModifiedGamma PSD-s do not even have a temperature 
dependence.
So in this case I really don't see any reason.

The temperature also comes into play when interpolating scat_data, but that is 
handled by
another part of ARTS. And over there, there is no option beside using the local 
temperature.


> 3.) Do you recommend me to use internally generated PSD for MC or set it 
> externally?

What you have started should work. So continue on this track.

(The alternative is to set pnd_field already in Python. Then you can ignore
particle_bulkprop_field ad pnd_agenda etc.

I can also mention that the PSD system in ARTS v2.5 soon will change. Simon is 
making a
complete revision of the handling of particles.)

Bye,

Patrick
___
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] Fwd: Re: ARTS internally generated PND fields

2020-12-06 Thread Thomas,Renish
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 
Sent: Sunday, December 6, 2020 8:25 AM
To: Thomas,Renish 
Cc: arts users mi 
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 
> *Sent:* Saturday, December 5, 2020 1:55 AM
> *To:* Thomas,Renish 
> *Cc:* arts users mi 
> *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" 
>> *Skickat:* 4 december 2020 16:40:25 CET
>> *Till:* Patrick Eriksson ,
>> "arts_users.mi@lists.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 

Re: [arts-users] Fwd: Re: ARTS internally generated PND fields

2020-12-06 Thread Patrick Eriksson

Ranish,

1.) So, if I must set the scattering particle temperature to an arbitrary value after 
calling "pnd_fieldCalcFromParticleBulkProps", do I just override and assign 
'pnd_agenda_input_t".


You got me wrong. You don't need to do anything. pnd_agenda_input_t is set automatically 
where needed.




2.) Can the particle temperature be set different than the ambient temperature?


Hardly. You can in principle try to "hack" pnd_agenda to set it to some other value. But 
why do you want to do it?


And please note that the psdModifiedGamma PSD-s do not even have a temperature dependence. 
So in this case I really don't see any reason.


The temperature also comes into play when interpolating scat_data, but that is handled by 
another part of ARTS. And over there, there is no option beside using the local temperature.




3.) Do you recommend me to use internally generated PSD for MC or set it 
externally?


What you have started should work. So continue on this track.

(The alternative is to set pnd_field already in Python. Then you can ignore 
particle_bulkprop_field ad pnd_agenda etc.


I can also mention that the PSD system in ARTS v2.5 soon will change. Simon is making a 
complete revision of the handling of particles.)


Bye,

Patrick
___
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] Fwd: Re: ARTS internally generated PND fields

2020-12-06 Thread Thomas,Renish
Hi Patrick/All,

Thanks, and that explains a lot.

1.) So, if I must set the scattering particle temperature to an arbitrary value 
after calling "pnd_fieldCalcFromParticleBulkProps", do I just override and 
assign 'pnd_agenda_input_t".

2.) Can the particle temperature be set different than the ambient temperature?

3.) Do you recommend me to use internally generated PSD for MC or set it 
externally?

Thanks,
Renish Thomas

From: Patrick Eriksson 
Sent: Sunday, December 6, 2020 8:25 AM
To: Thomas,Renish 
Cc: arts users mi 
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 
> *Sent:* Saturday, December 5, 2020 1:55 AM
> *To:* Thomas,Renish 
> *Cc:* arts users mi 
> *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" 
>> *Skickat:* 4 december 2020 16:40:25 CET
>> *Till:* Patrick Eriksson ,
>> "arts_users.mi@lists.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 
>> *Sent:* Thursday, November 19, 2020 11:47 AM
>> *To:* Thomas,Renish ; 
>> arts_users.mi@lists.uni-hamburg.de
>> 
>> *Subject:* Re: 

Re: [arts-users] Fwd: Re: ARTS internally generated PND fields

2020-12-06 Thread Patrick Eriksson

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 
*Sent:* Saturday, December 5, 2020 1:55 AM
*To:* Thomas,Renish 
*Cc:* arts users mi 
*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" 
*Skickat:* 4 december 2020 16:40:25 CET
*Till:* Patrick Eriksson , 
"arts_users.mi@lists.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 
*Sent:* Thursday, November 19, 2020 11:47 AM
*To:* Thomas,Renish ; 
arts_users.mi@lists.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

Re: [arts-users] Fwd: Re: ARTS internally generated PND fields

2020-12-05 Thread Thomas,Renish
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-time error 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 
Sent: Saturday, December 5, 2020 1:55 AM
To: Thomas,Renish 
Cc: arts users mi 
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" 
> *Skickat:* 4 december 2020 16:40:25 CET
> *Till:* Patrick Eriksson ,
> "arts_users.mi@lists.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 
> *Sent:* Thursday, November 19, 2020 11:47 AM
> *To:* Thomas,Renish ; 
> arts_users.mi@lists.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


Re: [arts-users] Fwd: Re: ARTS internally generated PND fields

2020-12-05 Thread Patrick Eriksson

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" 
*Skickat:* 4 december 2020 16:40:25 CET
*Till:* Patrick Eriksson , 
"arts_users.mi@lists.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 
*Sent:* Thursday, November 19, 2020 11:47 AM
*To:* Thomas,Renish ; arts_users.mi@lists.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