#268 closed Support (fixed)

Cant understand the variables in the trajectory.txt file

Reported by: sprasad Owned by: pesei
Priority: major Milestone:
Component: FP post-processing Version: FLEXPART-WRF
Keywords: output, post processing Cc:


I already registered a complaint that I can't understand what are that numbers in the txt file are indicating.
the text file has lots of numbers and there are no mentioning about which number is indicating lat,long etc
how can i understand that..?
please help me with some answers

Attachments (1)

trajectories.txt (2.4 KB) - added by sprasad 13 months ago.
here iam attachingthe txt file getting while running for trajectory

Download all attachments as: .zip

Change History (7)

Changed 13 months ago by sprasad

here iam attachingthe txt file getting while running for trajectory

comment:1 Changed 13 months ago by pesei

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

comment:2 Changed 13 months ago by pesei

This is duplicate with #258. Please do not open new tickets for the same issue!

Also, please understand we do not have the ressources do give everybody detailed support on simple issues. If you are a student, your supervisor should be your first address, and they should be able answer a question like this. As an exception, I have looked up in the source code file openouttraj.f90 what is written out (you could have done that as well). I think this should be self-explaining:

      open(unitouttraj,file=path(1)(1:length(1))//'trajectories.txt', &

      if (ldirect.eq.1) then
      write(unitouttraj,'(i8,1x,i6,1x,a)') ibdate,ibtime,'FLEXWRF  V3.0'
      write(unitouttraj,'(i8,1x,i6,1x,a)') iedate,ietime,'FLEXWRF  V3.0'
      write(unitouttraj,*) method,lsubgrid,lconvection
      write(unitouttraj,*) numpoint
      do i=1,numpoint
       if (outgrid_option .eq. 0) then
       if (outgrid_option .eq. 1) then
          xtmp = xpoint1(i)*dx+xmet0
          ytmp = ypoint1(i)*dy+ymet0
          call xymeter_to_ll_wrf( xtmp, ytmp, xp1, yp1 )
          xtmp = xpoint2(i)*dx+xmet0
          ytmp = ypoint2(i)*dy+ymet0
          call xymeter_to_ll_wrf( xtmp, ytmp, xp2, yp2 )
!jdf    write(unitouttraj,*) ireleasestart(i),ireleaseend(i),
!jdf +  xp1,yp1,xp2,yp2,zpoint1(i),zpoint2(i),kindz(i),npart(i)
!jdf    write(unitouttraj,'(a)') compoint(i)(1:40)
        write(unitouttraj,*) ireleasestart(i),ireleaseend(i), &
        write(unitouttraj,'(a20)') compoint(i)(1:20)
101     format(2i5,4f11.5,2f11.3,2i5)

comment:3 Changed 13 months ago by pesei

  • Milestone Long-term tasks deleted

comment:4 Changed 13 months ago by pesei

and plumetraj.f90

! Write out results in trajectory data file

           &j,itime-(ireleasestart(j)+ireleaseend(j))/2, &
           xcenter,ycenter,zcenter,topocenter,hmixcenter,tropocenter, &
           pvcenter,rmsdist,rms,zrmsdist,zrms,hmixfract,pvfract, &
           tropofract, &
           (xclust(k),yclust(k),zclust(k),fclust(k),rmsclust(k), &

The plume trajectory output is described in the FLEXPART v8 documentation, see FpDocumentation.

comment:5 Changed 12 months ago by pesei

Can we close the ticket?

comment:6 Changed 10 months ago by pesei

  • Resolution set to fixed
  • Status changed from accepted to closed

as there is no reply, I close the ticket.

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