Opened 4 years ago

Closed 4 years ago

#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:

Description

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 4 years ago.
here iam attachingthe txt file getting while running for trajectory

Download all attachments as: .zip

Change History (7)

Changed 4 years ago by sprasad

here iam attachingthe txt file getting while running for trajectory

comment:1 Changed 4 years ago by pesei

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

comment:2 Changed 4 years 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', &
      form='formatted',err=998)

      if (ldirect.eq.1) then
      write(unitouttraj,'(i8,1x,i6,1x,a)') ibdate,ibtime,'FLEXWRF  V3.0'
      else
      write(unitouttraj,'(i8,1x,i6,1x,a)') iedate,ietime,'FLEXWRF  V3.0'
      endif
      write(unitouttraj,*) method,lsubgrid,lconvection
      write(unitouttraj,*) numpoint
      do i=1,numpoint
       if (outgrid_option .eq. 0) then
         xp1=xpoint1(i)*dx+xmet0
         yp1=ypoint1(i)*dy+ymet0
         xp2=xpoint2(i)*dx+xmet0
         yp2=ypoint2(i)*dy+ymet0
        endif
       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 )
        endif
!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), &
        xp1,yp1,xp2,yp2,zpoint1(i),zpoint2(i),kindz(i),npart(i)
        write(unitouttraj,'(a20)') compoint(i)(1:20)
       enddo
101     format(2i5,4f11.5,2f11.3,2i5)

comment:3 Changed 4 years ago by pesei

  • Milestone Long-term tasks deleted

comment:4 Changed 4 years ago by pesei

and plumetraj.f90

! Write out results in trajectory data file
!******************************************

      write(unitouttraj,'(i5,1x,i8,1x,2f9.4,1x,4f8.1,1x,f8.2,1x,4f8.1,1x,3f6.1,&
           &5(1x,2f9.3,1x,f7.0,1x,f6.1,1x,f8.1))')&
           &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), &
           k=1,ncluster)

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

comment:5 Changed 4 years ago by pesei

Can we close the ticket?

comment:6 Changed 4 years 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