Opened 5 years ago

Last modified 12 months ago

#216 assigned Defect

how to know endpoint of each released particle

Reported by: ZhenhuaHao Owned by: igpis
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 (9)

comment:1 Changed 5 years 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 readinput.f90

Version 0, edited 5 years ago by pesei (next)

comment:2 Changed 5 years 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 5 years 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 5 years 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 5 years 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 5 years 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 5 years 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.

comment:8 Changed 5 years ago by ZhenhuaHao

Thank you so much for helping out!

comment:9 Changed 12 months ago by pesei

  • Owner changed from pesei to igpis
  • Status changed from accepted to assigned
  • Type changed from Support to Defect
Note: See TracTickets for help on using tickets.
hosted by ZAMG