Error: NetCDF: Invalid dimension ID or name

[lchong]

Dear all,
I'm trying to run flexpart-wrf-3.3 on LINUX (INTEL) but I try running it, it crashes with
WARNING: No (or incomplete) flux data contained in WRF output file

wrfout_d01_2015-04-27_18:00:00

WARNING: No (or incomplete) flux data contained in WRF output file
wrfout_d01_2015-04-27_19:00:00

WARNING: No (or incomplete) flux data contained in WRF output file
wrfout_d01_2015-04-27_20:00:00

readwind 0.2240000
calcpar 0.1655000
verttran 0.1810000
Error:
NetCDF: Invalid dimension ID or name

following is my setup files

=====================FORMER COMMAND FILE=====================

     1                LDIRECT:          1 for forward simulation, -1 for backward simulation
    20150427 180000  YYYYMMDD HHMISS   beginning date of simulation
    20150504 180000  YYYYMMDD HHMISS   ending date of simulation
    3600             SSSSS  (int)      output every SSSSS seconds
    3600             SSSSS  (int)      time average of output (in SSSSS seconds)
    180              SSSSS  (int)      sampling rate of output (in SSSSS seconds)
    999999999        SSSSS  (int)      time constant for particle splitting (in seconds)
    180              SSSSS  (int)      synchronisation interval of flexpart (in seconds)
    10.              CTL    (real)     factor by which time step must be smaller than tl
    10               IFINE  (int)      decrease of time step for vertical motion by factor ifine
    5                IOUT              1 concentration, 2 mixing ratio, 3 both, 4 plume traject, 5=1+4
    1                IPOUT             particle dump: 0 no, 1 every output interval, 2 only at end
    0                LSUBGRID          subgrid terrain effect parameterization: 1 yes, 0 no
    0                LCONVECTION       convection: 3 yes, 0 no
    3600.            DT_CONV  (real)   time interval to call convection, seconds
    0                LAGESPECTRA       age spectra: 1 yes, 0 no
    0                IPIN              continue simulation with dumped particle data: 1 yes, 0 no
    1                IFLUX             calculate fluxes: 1 yes, 0 no
    1                IOUTPUTFOREACHREL CREATE AN OUPUT FILE FOR EACH RELEASE LOCATION: 1 YES, 0 NO
    0                MDOMAINFILL       domain-filling trajectory option: 1 yes, 0 no, 2 strat. o3 tracer
    1                IND_SOURCE        1=mass unit , 2=mass mixing ratio unit
    2                IND_RECEPTOR      1=mass unit , 2=mass mixing ratio unit
    0                NESTED_OUTPUT     shall nested output be used? 1 yes, 0 no
    0                LINIT_COND   INITIAL COND. FOR BW RUNS: 0=NO,1=MASS UNIT,2=MASS MIXING RATIO UNIT
    1                TURB_OPTION       0=no turbulence; 1=diagnosed as in flexpart_ecmwf; 2 and 3=from tke.
    1                LU_OPTION         0=old landuse (IGBP.dat); 1=landuse from WRF
    1                CBL SCHEME        0=no, 1=yes. works if TURB_OPTION=1
    0                SFC_OPTION        0=default computation of u*, hflux, pblh, 1=from wrf
    0                WIND_OPTION       0=snapshot winds, 1=mean winds,2=snapshot eta-dot,-1=w based on divergence
    0                TIME_OPTION       1=correction of time validity for time-average wind,  0=no need
    1                OUTGRID_COORD     0=wrf grid(meters), 1=regular lat/lon grid
    1                RELEASE_COORD     0=wrf grid(meters), 1=regular lat/lon grid
    2                IOUTTYPE          0=default binary, 1=ascii (for particle dump only),2=netcdf
    1                NCTIMEREC (int)   Time frames per output file, only used for netcdf
    0                VERBOSE           VERBOSE MODE,0=minimum, 100=maximum

[lchong]

The backward simulation is ok, but the forward simulation has this problem

[pasko]

There is a typo/bug in

netcdf_output_mod.F90

which causes this error.

This patch should fix it

--- a/netcdf_output_mod.F90
+++ b/netcdf_output_mod.F90
@@ -621,7 +621,7 @@ module netcdf_output_mod
   !   write(*,*) ncdimsid2
       if ((ldirect.eq.1).and.(maxpointspec_act.gt.1)) then
 #ifdef NETCDF4_OUTPUT
-        ncret = NF90_DEF_VAR(ncid,'DRYDEP',NF90_FLOAT,ncdimsid2,ncddvid,&
+        ncret = NF90_DEF_VAR(ncid,'DRYDEP',NF90_FLOAT,ncdimsid22,ncddvid,&
             shuffle = shuffle, deflate_level=deflate_level)
 #else
         ncret = NF90_DEF_VAR(ncid,'DRYDEP',NF90_FLOAT,ncdimsid22,ncddvid)