Imposing photoionization rates

[Valentin Mauerhofer]

Hello,

Yet another post from me :-)

Is it possible, instead of giving a spectrum divided in N bins, to directly impose the photoionization rate on each species ? 

Thank you again for all your help,

All the best,

Valentin

[Tommaso Grassi]

Yes,

just give use the rate coefficient directly e.g.

 @format:idx,R,P,P,rate

 1,C,C+,E,3.1e-10

for the unshielded case.

Numbers are from here http://kida.obs.u-bordeaux1.fr/reaction/611/C_+_Photon.html?filter=Both

If you want to provide the visual extinction you can do

 @commons:user_av

 @format:idx,R,P,P,rate

 1,C,C+,E,3.1e-10*exp(-3.33*user_av)

and then do

 call krome_set_user_Av(your_Av)

before calling krome, where your_Av is the av you are using in the given cell.

Numbers are again from here http://kida.obs.u-bordeaux1.fr/reaction/611/C_+_Photon.html?filter=Both

cheers,

-tg

--
You received this message because you are subscribed to the Google Groups "KROMEusers" group.
To unsubscribe from this group and stop receiving emails from it, send an email to kromeusers+...@googlegroups.com.
To post to this group, send email to krome...@googlegroups.com.
Visit this group at https://groups.google.com/group/kromeusers.
To view this discussion on the web visit https://groups.google.com/d/msgid/kromeusers/fa8407ff-46af-4b3f-bf03-b28a8d97637b%40googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

[Valentin Mauerhofer]

Thanks, just to be sure that I have explained well my question, I would like to call some functions like :

call krome_set_photorate_SI(1d-10)

call krome_set_photorate_SIj(1d-11)

call krome_set_photorate_SIjj(1d-12) etc.

Is it doable with this user_av method ?

And, when I use 

 @commons:user_av

 @format:idx,R,P,P,rate

 1,C,C+,E,3.1e-10*exp(-3.33*user_av)

I get an error when I launch krome (./krome .....) :

WARNING: wrong format for reaction 1

@commons:user_av

What is happening ?

Cheers,

Valentin

[Tommaso Grassi]

HI,

my bad there is a typo, is @common: not @commons:

with this change if you copy the lines in network.dat or whatever file you like, using

 ./krome -n network.dat -noSinkCheck -noRecCheck

should work.

You can use @common for any reaction as

###

@common:user_phrate_SI
@common:user_phrate_SIj
@format:idx,R,P,P,rate
1,Si,Si+,E,user_phrate_SI

2,Si+,Si++,E,user_phrate_SIj

(you can easily prepare a script for doing this).

Then call

call krome_set_user_phrate_SI(1d-12) etc.

-tg

To view this discussion on the web visit https://groups.google.com/d/msgid/kromeusers/07b51172-1e64-4803-9ff6-e0dd6767f9db%40googlegroups.com.