wrf-flexpart输出数据的处理讨论
MeteoInfoLab脚本示例：数据投影-FLEXPART
FLEXPART应用
partposition文件：
第一张很简单的partposition文件所出的，笨办法直接导入excel用GIS出图，非常简单方便~~NCL，GRADS也都可以实现，偶只是觉得GIS最简单方便就用啦..
第二张密度，是GIS做的点密度处理~同样，其他也可以实现，就知道有童鞋用fortran计算了格点上的点个数来做的，grads出图的。。
第三张轨迹，其实就是partposition文件各个点连接在一起，注意时间先后，然后做平滑。。依然GIS实现的。。。
注意的是，这个轨迹和flexpart的trajectory文件是不同的。
问题：
WIND_OPTION选择0=snapshot winds的原因，我觉得如果WIND_OPTION选择1的话应该会有移动
2.trajectories文件有效位数！
./plumetraj.f90
write(unitouttraj,'(i8,1x,i6,1x,a)') ibdate,ibtime,'FLEXWRF V3.0'
write(unitouttraj,*) method,lsubgrid,lconvection
write(unitouttraj,*) numpoint
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)
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)
ncluster=5
3.** par_mod.f90 格点范围**
integer,parameter :: nxmax=334,nymax=334,nuvzmax=61,nwzmax=61,nzmax=61
4. partposition 输出格式
partoutput.f90
! Open output file and write the output
!**************************************
if (ipout.eq.1) then
if (iouttype.eq.0 .or. iouttype .eq.2) &
open(unitpartout,file=path(1)(1:length(1))//'partposit_'//adate// &
atime,form='unformatted')
if (iouttype.eq.1) &
open(unitpartout,file=path(1)(1:length(1))//'partposit_'//adate// &
atime,form='formatted')
else
if (iouttype.eq.0 .or. iouttype .eq.2) &
open(unitpartout,file=path(1)(1:length(1))//'partposit_end', &
form='unformatted')
if (iouttype.eq.1) &
open(unitpartout,file=path(1)(1:length(1))//'partposit_end', &
form='formatted')
endif
forrtl: severe (256): unformatted I/O to unit open for formatted transfers, unit 93, file /data/c03n02/cliu/study/flexpart/scripts/output/2013_07_31_05/partposit_20130731040000
readpartpositions.f90: open(unitpartin,file=path(1)(1:length(1))//'partposit_end', &
readpartpositions.f90: open(unitpartin,file=path(1)(1:length(1))//'partposit_end', &
partoutput.f90: open(unitpartout,file=path(1)(1:length(1))//'partposit_'//adate// &
partoutput.f90: open(unitpartout,file=path(1)(1:length(1))//'partposit_'//adate// &
partoutput.f90: open(unitpartout,file=path(1)(1:length(1))//'partposit_end', &
partoutput.f90: open(unitpartout,file=path(1)(1:length(1))//'partposit_end', &
readinput.f90: open(unitpartin,file=path(1)(1:length(1))//'partposit_end', &
readinput.f90: open(unitpartin,file=path(1)(1:length(1))//'partposit_end', &
在线课程、项目、文件
Using FLEXPART with Metview
The Lagrangian particle dispersion model FLEXPART
Investigating Source Regions of Persistent
Organic Pollutants with Lagrangian Modelling
工具
NOT UPDATED FOR UIO
- Generating web pages using NCARG routines FLEXPART Model/NCARGRoutinesAndMakeWebpages
- Plot your data with Flexplorer FLEXPART Model/FlexplorerManual
- Miscellaneous Scripts FLEXPART Model/MiscScripts
- Introduction to Backward Runs FLEXPART Model/BackwardRunsIntro
- Link your project output to the transport site FLEXPART Model/AddingProjectOutput
- Link your project output to the Interactive Page FLEXPART Model/FlexpartInteractive
- PYTHON For FLEXPART - See the main Python Page
- MATLAB For FLEXPART - See the main Matlab Page and therein:
