[RegCNET] PostProc Error!
Smz Rcm
smzrcm at yahoo.fr
Thu Dec 3 11:40:29 CET 2009
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')
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.ictp.it/pipermail/regcnet/attachments/20091203/df86a6dd/attachment.html>
More information about the RegCNET
mailing list