A comprehensive climate model assessment from multiple dimensions is critical for model selection to reduce uncertainties. Here, we evaluated the performance of seven models involved in the Coupled Model Intercomparison Project-6 (CMIP6) by comparing the simulated meteorological variables at the near-surface (sea ice cover and sea surface temperature) and pressure levels (air temperature, specific humidity, zonal wind, meridional wind, and geopotential height at 850 hPa, 500 hPa, and 200 hPa) to those using ERA5 reanalysis data for the 1995-2014 period from the perspectives of climatological and interannual variability. Then, a comprehensive rating index approach was applied to rank the models. The results show that the CMIP6 models could mostly reproduce the spatial variability of the climatology. However, there were also systematic biases. The sea ice cover and sea surface temperature exhibited noticeable biases of approximately -13.2% and 0.6 °C, respectively. Additionally, all models underestimated air temperature and geopotential height, while overestimating specific humidity in the middle troposphere and zonal and meridional wind speeds in the upper troposphere. Regarding interannual variability, the CMIP6 models performed well for the variables at a pressure level of 850 hPa and with sea ice cover. Taken together, all variables showed that NorESM2-MM exhibited good performance in terms of climatology, while MPI-ESM1-2-LR exhibited good interannual variability. Overall, in regard to the comprehensive rate index, MPI-ESM1-2-LR performed the best among the seven models. This study provides valuable scientific references for selecting the available CMIP6 models as lateral boundary conditions towards dynamical downscaling over Asian-Pacific area.