how to know endpoint of each released particle

[ZhenhuaHao]

Hi I am using Flexpart-WRF to simulate the trajectory of seed dispersion. So I want to know the trajectory of each seed released and their deposition position for each particles. so I have two questions:

1 In order to draw th trajectory of each seed released, is there any way that I can have the ID number for each particles I released so I can match them through different files of hourly dumped particles?

2 Is there any way that I can have the end point for each particle I released includes the ones that I lost before the simulation finished.

Thank you for addressing the possiblities and modifications on subrountines regarding to these two questions.

Best regards,
Iris

[pesei]

You have to set the

:! ipout            0 no particle dump, 1 every output time, 3 only at end  *

to 1. If you need, you can also modify partouput.f90

[ZhenhuaHao]

Yes, I did set ipout to 1 and it produce the particle location when it's still in the air. But when I start to lose particel, it will not leave as NAN in the next file, but be replaced by other particels, so I cant identify which particle I lost.

can you help sovle this issue?

Thank you.

cheers,
Iris

[ZhenhuaHao]

Hi

Can you please don't close this ticket as I still need help regarding to get particles landing locations. I've tried to fiddle with partouput.f90 file but failed.

Can you let me know if there's any of commented out codes in the subrountines can help generating landing position?

Thanks

[pesei]

Replying to ZhenhuaHao:

Yes, I did set ipout to 1 and it produce the particle location when it's still in the air. But when I start to lose particel, it will not leave as NAN in the next file, but be replaced by other particels, so I cant identify which particle I lost.

OK, I understand the problem. Maybe you need also the mquasilag set to 1. I have to admit that I haven't worked with this feature myself. Please try and check.

[ZhenhuaHao]

Thank you for your reply. Can you point me to the file that I can set the value of mquasilag?

Thanks

[pesei]

Replying to ZhenhuaHao:

Thank you for your reply. Can you point me to the file that I can set the value of mquasilag?

Good question. This feature is present in the regular FLEXPART model, and it seems that it has not been fully integrated into the WRF version. While the code in releaseparticles (which is the decisive part) is preserved, it is missing in reading the input. My suggestion would be to add a line

read(unitcommand,*) mquasilag

to readinput.f90 and a corresponding line in your input file, giving a value of 1.

[ZhenhuaHao]

Thank you so much for helping out!