search

jjung11 / Mshpy

Simplified Python3 magnetosheath model
MIT License
1 stars 0 forks source link

soft x-ray emissivity graph #1

Open hiscolor opened 5 months ago

hiscolor commented 5 months ago

Hello! I would like to reproduce the soft x-ray emissivity graph in the article about SMILE at a specific location, can you briefly explain the steps involved? Thank you!

jjung11 commented 5 months ago

Hi! I apologize, but the module for generating the soft X-ray emissivity plot has not been added to this repository yet. I plan to work on incorporating that module, but initially, I anticipated doing so in a few weeks due to other commitments. Do you have a specific timeline in mind for this addition?

hiscolor commented 5 months ago

Thank you for your reply. I don't have a timeline, but I just want to try your research results as soon as possible, hoping to help me with the data of my thesis. Thank you again!

jjung11 commented 5 months ago

Now it is added, though it is a bit messy right now, as I'm not accustomed to how to distribute code through Github and so. Do 'pip install Mshpy==0.0.56' In your Python session, do from Mshpy.Msh_xray import mshpy_xray

Before running the code, you should set in which solar wind condition the simulation to be done. Solar wind condition should be pre-exist as 'SW_cond.txt' as in current directory. Format of the file should be following: Year DOY HR MN Bx By Bz Vx Vy Vz n Pd Ma Mm, as in OMNIweb e.g. 1967 1 0 00 2 -2 -5 -400 0 0 10 2 10 6

now do something like: mshpy_xray('test',0,'test','jel',30,75,105,1,180,210,1,0,0,1)

Brief explanation: mshpy_xray(RUN, TI, Dout, model_bs, dsc, the1, the2, dthe, phi1, phi2, dphi, XRAYPS, MSHPY, ALL): RUN is just name of your simulation. TI is just time number. It may be useful if you want time series simulation. Dout is the output directory. model_bs means bow shock model. Mshpy support Jelinek (jel) and Jerab (jer) dsc means spacecraft distance from the earth in Re. the1/the2 means field of view range of s/c in terms of latitudinal direction (in spacecraft coordintes. It is not the same as GSE) dthe is resolution in latitudinal direction. Same for phi1/phi2/dphi, longitudinal direction. XRAYPS is 1 if you want calculate data. MSHPY is 1 if you want to plot it. ALL is 1 if you want to do both.

Hope this helps and let me know if you are stuck with something.