Custom Query (210 matches)

Dear colleagues, It seems that deposition units in output NETCDF files of FLEXPART-WRF are defined erroneously (pg·m-2) while It should be (ng·m-2)

This problem was first noted because of extremely high underestimation of deposition following release scenario as compared to results of simulations of other models. E.g. dry deposition calculated by the RODOS system for the same site, same release scenario and using the same meteorological data was on average by about 3 orders of magnitude higher as compared to FLEXPART results.

The source of potential error was identified by analysis of the FLEXPART code. The variable drygrid2 which is written to output netcdf file is defined in subroutine nc_write_output (source file netcdf_output_mod) as:

drygrid2(ix,jy,kp,nage)=1.e12*drygrid(ix,jy)/area(ix,jy)

Thus converstion to pico- units is done by multiplication on 1E12 factor here. Further inspection of code showed that the variable drygrid is in case of 1 age class and single specie simply equal to the variable drygridunc2, while the latter is defined in time loop when running model in the subroutine drydepokernel:

drygridunc2(ix,jy,ks,kp,nunc,nage)= &

drygridunc2(ix,jy,ks,kp,nunc,nage)+deposit(ks)

Regarding the variable deposit the comment in the same subroutine states that:

! deposit amount (kg) to be deposited Thus final deposition should be in pKg·m-2=ng·m-2.

Attached is file with formulas used for conversion from kg to Bq.

I have tried to run flexpart-wrf with different version(from 3.0 to 3.2). However, I found the valuse of potential voricity equal zeros, while I set particle dump for every output interval. Is any bugs in the programes? thanks for your help in advance.

BTW: There seems a bug in convmix_kfeta.f90 programe, if set variable "maxpar" not equal "numpart", which could make the value of "maxval(igrid)" too big and programe stop at the line of "too much x and y grid modify convmix_kfeta.f". This bug could be solved by add another loop behind the line 154.

do ipart=1,maxpart igrid(ipart)=-1 enddo

Hello,

I used flexpart-wrf6.2 for INDIVIDUAL track.But I can not found the MQUASILAG (command file 21 opt )in the new flexwrf3.x. I also find that the out put of particles for individual path program and "MQUASILAG=1" was commented in the timemanager.f90. For this reason, I think the first column in the particle output file was the release area only. I don't know where should I find the id of each particles? Thanks!