Opened 5 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)
Change History (7)
Changed 5 years ago by sprasad
comment:1 Changed 5 years ago by pesei
- Owner set to pesei
- Status changed from new to accepted
comment:2 Changed 5 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 5 years ago by pesei
- Milestone Long-term tasks deleted
comment:4 Changed 5 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.
here iam attachingthe txt file getting while running for trajectory