Opened 4 years ago
Last modified 4 years ago
#233 accepted Defect
Error: NetCDF: Invalid dimension ID or name
Reported by: | lchong | Owned by: | pasko |
---|---|---|---|
Priority: | major | Milestone: | |
Component: | FP coding/compilation | Version: | FLEXPART-WRF |
Keywords: | Cc: |
Description
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
Change History (4)
comment:1 Changed 4 years ago by lchong
comment:2 Changed 4 years ago by pasko
- Owner set to pasko
- Status changed from new to accepted
comment:3 Changed 4 years ago by 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)
comment:4 Changed 4 years ago by pasko
- Component changed from FP other to FP coding/compilation
Note: See
TracTickets for help on using
tickets.
The backward simulation is ok, but the forward simulation has this problem