将MM5输出文件转化，运行run.mcip出现的错误

shepherdmeng

想将MM5生成的MMOUT_DOMAIN文件用MCIP转化成CM AQ可输入的格式，运行run.mcip时出现如下错误 ******************************************************************************************************************************************** *** SUBROUTINE: GETVERSION
***   ERROR READING MM5 FILE ON UNIT   10
***   IOSTAT =   39
**********************************************************************

     *** ERROR ABORT in subroutine GETVERSION
     ABNORMAL TERMINATION IN GETVERSION
     Date and time  8:00:00   Dec. 24, 2004   (2004359:080000)


Error running mcip
下面是run.mcip设置
source ../config.cmaq


set APPL       = 041224
set CoordName  = LAMCON  # 16-character maximum
set GridName   = vgrd  # 16-character maximum



#set HomePath   = /model3/CMAQv5.0/data/mcip/data/output
set DataPath   = /model3/CMAQv5.0/data/mcip
set InMetDir   = /model3/mm5/MM5/Run
set InTerDir   = /model3/mm5/TERRAIN
#set InSatDir   = $DataPath/goes
set OutDir     = $DataPath/data/output
set ProgDir    = $cwd/src
set WorkDir    = $OutDir


#-----------------------------------------------------------------------
# Set name(s) of input meteorology file(s)
#
#   File name(s) must be set inside parentheses since "InMetFiles" is
#   a C-shell script array.  Multiple file names should be space-
#   delimited.  Additional lines can be used when separated by a
#   back-slash (\) continuation marker.  The file names can be as
#   they appear on your system; MCIP will link the files in by a
#   Fortran unit number and the explicit name via a namelist.  The
#   files must be listed in chronological order.  The maximum number
#   of input meteorology files must be less than or equal to the number
#   in MAX_MM in file_mod.F (default is 300).
#
#   Example:
#     set InMetFiles = ( $InMetDir/MMOUT_DOMAIN2.time1 \
#                        $InMetDir/MMOUT_DOMAIN2.time2 )
#
#-----------------------------------------------------------------------


set InMetFiles = ( $InMetDir/MMOUT_DOMAIN1 )


set IfTer      = "T"
set InTerFile  = $InTerDir/TERRAIN_DOMAIN1


#set InSatFiles = ( )


#-----------------------------------------------------------------------
# Set user control options.
#
#   LPV:     0 = Do not compute and output potential vorticity
#            1 = Compute and output potential vorticity
#
#   LWOUT:   0 = Do not output vertical velocity
#            1 = Output vertical velocity
#
#   LUVCOUT: 0 = Do not output u- and v-component winds on C-grid
#            1 = Output u- and v-component winds on C-grid
#
#   LSAT:    0 = No satellite input is available (default)
#            1 = GOES observed cloud info replaces model-derived input
#-----------------------------------------------------------------------


set LPV     = 0
set LWOUT   = 0
set LUVCOUT = 1
set LSAT    = 0


#-----------------------------------------------------------------------
# Set run start and end date.  (YYYY-MO-DD-HH:MI:SS.SSSS)
#   MCIP_START:  First date and time to be output [UTC]
#   MCIP_END:    Last date and time to be output  [UTC]
#   INTVL:       Frequency of output [minutes]
#-----------------------------------------------------------------------


set MCIP_START = 2004-12-24-08:00:00.0000  # [UTC]
set MCIP_END   = 2004-12-28-06:00:00.0000  # [UTC]


set INTVL      = 60 # [min]


#-----------------------------------------------------------------------
# Set CTM layers.  Should be in descending order starting at 1 and
# ending with 0.  There is currently a maximum of 100 layers allowed.
# To use all of the layers from the input meteorology without
# collapsing (or explicitly specifying), set CTMLAYS = -1.0.
#-----------------------------------------------------------------------


#set CTMLAYS = "-1.0"
set CTMLAYS = "1.000, 0.990, 0.980, 0.970, 0.960, 0.800, 0.700, 0.600, \
               0.500,  0.300, 0.200, 0.100,0.000"

#set CTMLAYS = "1.000, 0.995, 0.990, 0.980, 0.960, 0.940, 0.910, 0.860, \
#               0.800, 0.740, 0.650, 0.550, 0.400, 0.200, 0.000"


#-----------------------------------------------------------------------
# Determine whether or not static output (GRID) files will be created.
#-----------------------------------------------------------------------


set MKGRID = T


#-----------------------------------------------------------------------
# Set number of meteorology "boundary" points to remove on each of four
# horizontal sides of MCIP domain.  This affects the output MCIP domain
# dimensions by reducing meteorology domain by 2*BTRIM + 2*NTHIK + 1,
# where NTHIK is the lateral boundary thickness (in BDY files), and the
# extra point reflects conversion from grid points (dot points) to grid
# cells (cross points).  Setting BTRIM = 0 will use maximum of input
# meteorology.  To remove MM5 lateral boundaries, set BTRIM = 5.
#
# *** If windowing a specific subset domain of input meteorology, set
#     BTRIM = -1, and BTRIM will be ignored in favor of specific window
#     information in X0, Y0, NCOLS, and NROWS.
#-----------------------------------------------------------------------


set BTRIM = -1


#-----------------------------------------------------------------------
# Define MCIP subset domain.  (Only used if BTRIM = -1.  Otherwise,
# the following variables will be set automatically from BTRIM and
# size of input meteorology fields.)
#   X0:     X-coordinate of lower-left corner of full MCIP "X" domain
#           (including MCIP lateral boundary) based on input MM5 domain.
#           X0 refers to the east-west dimension.  Minimum value is 1.
#   Y0:     Y-coordinate of lower-left corner of full MCIP "X" domain
#           (including MCIP lateral boundary) based on input MM5 domain.
#           Y0 refers to the north-south dimension.  Minimum value is 1.
#   NCOLS:  Number of columns in output MCIP domain (excluding MCIP
#           lateral boundaries).
#   NROWS:  Number of rows in output MCIP domain (excluding MCIP
#           lateral boundaries).
#-----------------------------------------------------------------------


set X0    =   6
set Y0    =   6
set NCOLS = 279
set NROWS = 240


#-----------------------------------------------------------------------
# Set coordinates for cell for diagnostic prints on output domain.
# If coordinate is set to 0, domain center cell will be used.
#-----------------------------------------------------------------------


set LPRT_COL = 0
set LPRT_ROW = 0


#-----------------------------------------------------------------------
# Optional:  Set WRF Lambert conformal reference latitude.
#            (Handy for matching WRF grids to existing MM5 grids.)
#            If not set, MCIP will use average of two true latitudes.
# To "unset" this variable, set the script variable to "-999.0".
# Alternatively, if the script variable is removed here, remove it
# from the setting of the namelist (toward the end of the script).
#-----------------------------------------------------------------------


set WRF_LC_REF_LAT = -999.0


#=======================================================================
#=======================================================================
# Set up and run MCIP.
#   Should not need to change anything below here.
#=======================================================================
#=======================================================================


set PROG = mcip


date
*** SUBROUTINE: GETVERSION
***   ERROR READING MM5 FILE ON UNIT   10
***   IOSTAT =   39
**********************************************************************

     *** ERROR ABORT in subroutine GETVERSION
     ABNORMAL TERMINATION IN GETVERSION
     Date and time  8:00:00   Dec. 24, 2004   (2004359:080000)


Error running mcip
Meteorological Numerical Model Union of China (MNMUC) 中国气象数值模式联盟% Y( u$ w: b" x9 W0 N# W+ Q