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.

Hi all,

has anybody ran Flexpart 9.0.3 successfully with GFS 0.25 degree data in grib2 format? (e.g. from here: ​https://rda.ucar.edu/datasets/ds083.3/#!description)

After compiling with par_mod.f90 set to the larger grid integer,parameter :: nxmax=1441,nymax=721,nuvzmax=64,nwzmax=64,nzmax=64, the error with t not in [K] occurs for any input data (1deg or 0.25deg, both grib2). Hence I suspect a problem related to the grid size.

Any workaround? Thanks in advance, Martin

I am using Flexpart 10.4 and am trying to restart from the positions of the last simulation in partposit_end. But I am not sure how to modify the RELEASES file for the simulation to run. When I try to set PARTS = 0 (Total number of particles to be released), I get the error :

FLEXPART MODEL ERROR RELEASES file is corrupt. At least for one release point, there are zero

If I try to set PARTS = 1, I get :

##################################################### #### FLEXPART MODEL SUBROUTINE RELEASEPARTICLES: #### #### #### #### ERROR - TOTAL NUMBER OF PARTICLES REQUIRED #### #### EXCEEDS THE MAXIMUM ALLOWED NUMBER. REDUCE #### #### EITHER NUMBER OF PARTICLES PER RELEASE POINT#### #### OR REDUCE NUMBER OF RELEASE POINTS. #### #####################################################

I also tried to set NSPEC = 0, but then I get

#### FLEXPART MODEL ERROR IN READPARTPOSITIONS#### #### THE NUMBER OF SPECIES TO BE READ IN DOES #### #### NOT AGREE WITH CURRENT SETTINGS! ####

Any help will be appreciated, thank you!