| Dear RegCMers, I have an error each time I run the postproc.x so I didn't get any result file. The message and the error are: " ifort -tpp7 -O3 -cm -w -w90 -w95 -Vaxlib -convert big_endian -o postproc postproc.o netcdf.o -L../Commons/env/liblinux/ -lnetcdf -ludunits -Vaxlib -I../Commons/env/include ENTER THE TYPE OF REGCM OUTPUT TO BE PROCESSED: ICBC (0=no; 1=yes) 0 ATM (0=no; 1=yes) 1 SRF (0=no; 1=yes) 0 RAD (0=no; 1=yes) 0 CHE (0=no; 1=yes) 0 SUB (0=no; 1=yes) 0 ******************************************************* IDATE0 = 2000050100 IDATE1 = 2000050100 IDATE2 = 2000060100 ******************************************************* $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ POSTPROCESSING OF 3D DATA WILL BE PERFORMED ON SIGMA LEVELS $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ OUTPUT WILL BE WRITTEN IN NETCDF 4-BYTE FORMAT $$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$$ OUTPUT FILE: ATM1..nc mdate0,ibltyp,icup,imoist,iboudy,ni,nj,nk,ds= 40 40 18 1189765120 1108567654 1092846879 1108567654 1092846879 -2.3087023E+20 sigf= -71.62398 1.2684776E-26 -1.5200503E+14 9.1119407E-33 -3.3707561E-11 -71.62359 89214.41 -1.0760833E+08 -2.3337048E-17 -3.3707204E-11 -5.0126368E-05 -7.4876450E-36 -1.5200069E+14 -1.0760630E+08 -1.0760560E+08 2.1989176E-41 3.4561065E+09 5.0672908E+07 0.7155668 pt,clat,clon,xplat,xplon,proj,dto,dtb,dtr,dtc= 4.0665847E+29 2.3183082E-41 0.0000000E+00 9.1835496E-41 1.1548941E-40 Bp 0.0000000E+00 0.0000000E+00 0.0000000E+00 0.0000000E+00 Grid Dimensions DO NOT MATCH nyf= 40 nxf= 40 nz= 18 ni= 1092846879 nj= 1108567654 nk= 1092846879 Also check ibyte in postproc.param: ibyte= 4 BAD DIMENSIONS (SUBROUTINE RDHEAD) " As a result, I have no result files, so please help me fixing this error! I thank you very much for interest yours SASSI The following files are postproc.in, postproc.param and postproc1.param: postproc.in: 2000050100, ! idate0 = First date in File (yymmddhh) 2000050100, ! idate1 = Start date for averaging and re-writing 2000060100, ! idate2 = End date for averaging and re-writing 2, ! iotype: 1=I*2 NetCDF; 2=r*4 NetCDF; 3=grads; 4=vis5d .false., ! plev .false., ! usgs .true., ! Write out header? .false., ! Write out all RegCM data b/twn idate1 & idate2? .false., ! Average RegCM data b/twn idate1 & idate2? .false., ! Diurnali avg of RegCM data b/twn idate1 & idate2? .true., ! Continually average ATM data b/twn idate1 & idate2? 1., ! No. Days for continual averaging (-1=monthly;1=daily;5=5day) '200005AVE', ! postproc output filename (not including type & ext) '../Input', ! ICBC directory 'output', ! RegCM Output directory 'DOMAIN.INFO', ! Domain Info Filename (from terrain) 'OUT_HEAD', ! Header File name (from RegCM) '2000050100', ! 1stRegCM Output File Extension postproc.param: cccc SET DOMAIN DIMENSIONS cccc ny = number of north-south points cccc nx = number of east-west points cccc nz = number of vertical levels cccc SET MODEL OUTPUT INTERVALS cccc dtbc = ibdyfrq = ICBC output interval (hrs) cccc dtout = atmfrq = atmospheric output interval (hrs) cccc dtbat = srffrq = surface output interval (hrs) cccc dtrad = radfrq = radiation output interval (hrs) cccc DIRECT ACCESS BINARY SPECIFICATION cccc ibyte = 4 for PGI, IFC; 1 for SUN, SGI, DEC, IBM integer nxf,nyf,nz,nxs,nys parameter (nyf = 40) parameter (nxf = 40) parameter (nz = 18) parameter (nxs = 1) parameter (nys = 1) integer ibyte parameter (ibyte = 4) real dtbc,dtout,dtbat,dtrad,dtche,dtsub parameter (dtbc = 6.00) parameter (dtout = 6.00) parameter (dtbat = 3.00) parameter (dtrad = 6.00) parameter (dtche = 6.00) parameter (dtsub = 3.00) and postproc1.param: integer nx,ny,nx1,ny1,nb1,nxb,nyb & , ni3d,ni2d,ni3d2,ni2d2,nbc3d,nbc2d,nitot & , no3d,no2d,no3d2,no2d2,nout3d,nout2d,notot & , nbat,nbat2,nsub,nsub2,nr2d,nr3d,nrtot & , nhrbc,nhrout,nhrbat,nhrsub,nhrrad,nhrche & , ntrac,nc2d,nc3d,nctot,nrtot1,nctot1 & , nxsf,nysf,nxsb,nysb,npl parameter (ny=nyf-2, nx=nxf-2) parameter (nysf=nyf*nys,nxsf=nxf*nxs) parameter (nysb=ny*nys,nxsb=nx*nxs) parameter (ny1=ny,nx1=nx) parameter (nb1=1,nyb=ny,nxb=nx) parameter (nhrbc=24.001/dtbc,nhrout=24.001/dtout & , nhrbat=24.001/dtbat,nhrsub=24.001/dtsub & , nhrrad=24.001/dtrad,nhrche=24.001/dtche) c ******************* SET NUMBER OF OUTPUT FIELDS ****************** c c **** ni3d = number of 3D ICBC atmospheric fields **** c c **** ni2d = number of 2D ICBC atmospheric fields **** c c **** no3d = number of 3D atmospheric fields (from model) **** c c **** no2d = number of 2D atmospheric fields (from model) **** c c **** nbat = number of surface fields (from model) **** c c **** nsub = number of subbats-surface fields (from model) **** c c **** nr3d = number of 3D radiation fields (from model) **** c c **** nr2d = number of 2D radiation fields (from model) **** c c **** **** c c **** The other fields (ni3d2, ni2d2, no3d2, no2d2, nbat2 and **** c c **** nsub2) should only be changed if additional fields are **** c c **** computed from the model output. Examples of additional **** c c **** computed fields are geopotential height, moist static **** c c **** energy, and relative humidity. **** c c ****************************************************************** c parameter (ni3d=4, ni2d=2) ! 3D and 2D variables in ICBC files parameter (no3d=6, no2d=5) ! 3D and 2D variables in ATM files parameter (nbat=27) ! 2d variables in SRF files parameter (nsub=16) ! 2d variables in SUB files parameter (nr3d=4,nr2d=10) ! 2D and 3D variables in RAD files. This includes PS from ATM parameter (ntrac = 3) ! Number of CHEM tracers parameter (nc3d=ntrac*1+3, nc2d=ntrac*7+3) ! 3D and 2D varaibles in CHEM files parameter (ni3d2=7, ni2d2=2) ! new 3D and 2D varaibles calculated in PostProc from ICBC parameter (no3d2=7, no2d2=2) ! new 3D and 2D varaibles calculated in PostProc from ATM parameter (nbat2=nbat+2) ! total SRF varaibles (nbat) + new calculated parameter (nsub2=nsub+4) ! total SUB varaibles (nsub) + new calculated parameter (nbc3d=ni3d+ni3d2, nbc2d=ni2d+ni2d2 & , nitot=nbc3d+nbc2d) ! Total 3D, 2D and 3D+2D varaibles for ICBC data parameter (nout3d=no3d+no3d2, nout2d=no2d+no2d2 & , notot=nout3d+nout2d) ! Total 3D, 2D ad 3D+2D varaibles for ATM data parameter (nrtot1=nr3d+nr2d, nrtot=nrtot1-1) ! Total RAD varaibles with and without PS parameter (nctot1=nc3d+nc2d, nctot=nctot1-1) ! Total CHEM variables with and without PS parameter (npl=11) ! Number of Pressure levels real plev(npl),plevr(npl) ! Pressure levels in normal and reverse order data plev/1000.,925.,850.,700.,500.,400.,300.,250.,200. & ,150.,100./ data plevr/100.,150.,200.,250.,300.,400.,500.,700.,850. & ,925., 1000./ c ****************************************************************** c c **** nfmax is the maximum number of files outputed by the **** c c **** model simulation. **** c c ****************************************************************** c integer nfmax parameter (nfmax=999) c ********************* NAME OF THE INPUT FILE ********************* c c **** This file ("plist") contains the information on the **** c c **** dates and filenames to processed and the type of **** c c **** processing to be performed. **** c c **** **** c c **** If "plist" is modified, the postprocessing code does not **** c c **** need to be recompile. **** c c ****************************************************************** c character plist*70, ulist*70 parameter(plist = 'postproc.in') parameter(ulist = 'user.in') |