wrf-flexpart输出数据的处理讨论
MeteoInfoLab脚本示例：数据投影-FLEXPART
FLEXPART应用
WRF-FLEXPART安装和基本运行
FLEXPART-WRF
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
European Integrated project on Aerosol Cloud Climate and Air Quality Interactions (EUCAARI) is a project within EU's Sixth Framework Program.
FLEXPART Products
The Lagrangian particle dispersion model FLEXPART
Investigating Source Regions of Persistent
Organic Pollutants with Lagrangian Modelling
工具
reflexible on github
readthedocs.io
FLEXPART postprocessing tools
QUICKLOOK
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:
脚本
partoutput.py
clustering.f90
plot_tracjectory.ncl
plot_trajs.ncl