Opened 3 weeks ago

Last modified 3 weeks ago

#268 accepted Support

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

Download all attachments as: .zip

Change History (5)

Changed 3 weeks ago by sprasad

here iam attachingthe txt file getting while running for trajectory

comment:1 Changed 3 weeks ago by pesei

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

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

  • Milestone Long-term tasks deleted

comment:4 Changed 3 weeks 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.

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