[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