Opened 2 weeks ago

Last modified 4 hours ago

#216 accepted Support

how to know endpoint of each released particle

Reported by: ZhenhuaHao Owned by: pesei
Priority: minor Milestone: FLEXPART_WRF_3.4_FPbase_9
Component: FP input data Version: FLEXPART-WRF
Keywords: endpoint, Particle ID number Cc:

Description

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

Change History (7)

comment:1 Changed 2 weeks ago by pesei

  • Owner set to pesei
  • Status changed from new to accepted

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

Last edited 2 weeks ago by pesei (previous) (diff)

comment:2 Changed 2 weeks ago by pesei

  • Component changed from FP coding/compilation to FP input data
  • Priority changed from critical to minor
  • Type changed from New feature to Support

comment:3 follow-up: Changed 2 weeks ago by 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

comment:4 Changed 11 days ago by 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

comment:5 in reply to: ↑ 3 Changed 3 days ago by 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.

comment:6 follow-up: Changed 36 hours ago by ZhenhuaHao

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

Thanks

comment:7 in reply to: ↑ 6 Changed 4 hours ago by 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.

Note: See TracTickets for help on using tickets.
hosted by ZAMG