Skip to content

Instantly share code, notes, and snippets.

@lizadams
Created April 11, 2019 19:09
Show Gist options
  • Save lizadams/4d5a3c76d28621f30485f398e491c7ff to your computer and use it in GitHub Desktop.
Save lizadams/4d5a3c76d28621f30485f398e491c7ff to your computer and use it in GitHub Desktop.
mcipv4.5 run script error
[lizadams@dogwood-login2 scripts]$ ./run_mcip.csh intel
/proj/ie/proj/CMAS/CMAQ/CMAQv5.3_branch_UNC3/mvapich2_2.3rc1_intel_17.2
Compiler is set to intel
/proj/ie/proj/CMAS/CMAQ/CMAQv5.3_branch_UNC3/mvapich2_2.3rc1_intel_17.2/PREP/mcip/scripts
Thu Apr 11 15:05:18 EDT 2019
rm: No match.
This program uses the EPA-AREAL/MCNC-EnvPgms/BAMS Models-3
I/O Applications Programming Interface, [I/O API] which is
built on top of the netCDF I/O library (Copyright 1993, 1996
University Corporation for Atmospheric Research/Unidata
Program) and the PVM parallel-programming library (from
Oak Ridge National Laboratory).
Copyright (C) 1992-2002 MCNC,
(C) 1992-2013 Carlie J. Coats, Jr.,
(C) 2003-2012 Baron Advanced Meteorological Systems, LLC, and
(C) 2014-2016 UNC Institute for the Environment.
Released under the GNU LGPL License, version 2.1. See URL
https://www.gnu.org/licenses/old-licenses/lgpl-2.1.html
for conditions of use.
ioapi-3.2: $Id: init3.F90 98 2018-04-05 14:35:07Z coats $
Version with PARMS3.EXT/PARAMETER::MXVARS3= 2048
netCDF version 4.1.3 of Jun 19 2018 15:50:22 $
EXECUTION_ID: mcip
==============================================================================
US EPA COMMUNITY MULTISCALE AIR QUALITY MODEL
METEOROLOGY-CHEMISTRY INTERFACE PROCESSOR
MCIP V4.5 FROZEN 03/16/2018
==============================================================================
- SUBROUTINE SETUP_WRFEM - READING WRF HEADER
WRF GRID DIMENSIONS (X,Y,Z) 103 83 35
------------------------------------------------------------------------------
USER OPTIONS AND GRID DEFINITIONS
Output start date = 2011-07-01-00:00:00.0000
Output end date = 2011-07-02-00:00:00.0000
Output interval = 60 minutes
LPV = 1: Will calculate and output 3D potential vorticity
LWOUT = 0: Will not output vertical velocity from the met model
LUVCOUT = 1: Will output u- and v-component winds on C grid
LSAT = 0: Will not use satellite adjustment of clouds for photolysis
LAI was found in the meteorology input file
FRACTIONAL LAND USE was found in the meteorology input file
FRACTIONAL LAND USE will be read from the WRF file
MONIN-OBUKHOV LENGTH was found in the meteorology input file
MONIN-OBUKHOV LENGTH will be recomputed by MCIP (Only for WRF P-X LSM)
AERODYNAMIC AND STOMATAL RESISTANCE was found in the meteorology input file
2-m TEMPERATURE was found in the meteorology input file
2-m MIXING RATIO was found in the meteorology input file
VEGETATION FRACTION was found in the meteorology input file
10-m WIND was found in the meteorology input file
CANOPY WETNESS was found in the meteorology input file
ROUGHNESS LENGTH was found in the meteorology input file
SOIL MOISTURE, TEMPERATURE, AND TYPE was found in the meteorology input file
SOIL MOISTURE, TEMPERATURE, AND TYPE will be in the output file
TURBULENT KINETIC ENERGY was NOT found in the meteorology input file
URBAN CANOPY MODEL (WRF ONLY) was NOT used in the meteorology model
3D RESOLVED CLOUD FRACTION was NOT found in the meteorology input file
HYBRID VERTICAL COORDINATE (WRF ONLY) was NOT used in the meteorology model
Met domain dimensions (col, row, lay): 103 83 35
MCIP X domain dimensions (col, row, lay): 102 82 35
Output domain dimensions (col, row, lay): 100 80 35
Output grid resolution: 12.0000000000000 km
Window domain origin on met domain (col,row): 1 , 1
Window domain far corner on met domain (col,row): 103 , 83
Cells and points for diagnostic prints
LPRT_COL LPRT_ROW 50 40
LPRT_XCOL LPRT_XROW 51 41
LPRT_METX LPRT_METY 51 41
IOAPI header variables:
GDTYP3D 2
GDNAM3D SE52BENCH
P_ALP3D 33.000
P_BET3D 45.000
P_GAM3D -97.000
XCENT3D -97.000
YCENT3D 40.000
XORIG3D 792000.000
YORIG3D -1080000.000
XCELL3D 12000.000
YCELL3D 12000.000
VGTYP3D 7
VGTOP3D 5000.000
VGLVS3D 1.00000 0.99750 0.99500 0.99000 0.98500 0.98000
0.97000 0.96000 0.95000 0.94000 0.93000 0.92000
0.91000 0.90000 0.88000 0.86000 0.84000 0.82000
0.80000 0.77000 0.74000 0.70000 0.65000 0.60000
0.55000 0.50000 0.45000 0.40000 0.35000 0.30000
0.25000 0.20000 0.15000 0.10000 0.05000 0.00000
------------------------------------------------------------------------------
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
~~~ Processing meteorology for time = 2011-07-01-00:00:00.0000
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
U -0.5964 -1.0403 -1.4599 -1.6978 -1.8199
-2.0807 -2.2600 -2.3977 -2.4544 -2.4656
-2.4317 -2.4188 -2.3964 -2.3666 -2.2259
-2.2015 -2.2499 -2.2323 -1.8278 -0.8676
1.0808 1.7251 0.3862 1.4014 3.9795
4.9639 4.0927 5.3105 7.8828 7.1929
6.8925 8.3963 10.0014 11.2059 -3.4874
V -1.1953 -1.9070 -2.3910 -2.5187 -2.4803
-2.4053 -2.3986 -2.3874 -2.4024 -2.4354
-2.4502 -2.4705 -2.4900 -2.5159 -2.4002
-2.4326 -2.5392 -2.4911 -2.2160 -2.5166
-4.5179 -4.7055 -6.6039 -8.7604 -10.0605
-8.8545 -7.6360 -6.5225 -5.1344 -8.6715
-7.9447 -7.0248 -4.6970 -6.3167 -7.6872
W -0.0013 -0.0068 -0.0088 -0.0094 -0.0093
-0.0093 -0.0090 -0.0089 -0.0090 -0.0093
-0.0090 -0.0091 -0.0093 -0.0090 -0.0089
-0.0086 -0.0091 -0.0087 -0.0071 -0.0038
-0.0047 0.0020 0.0146 0.0197 0.0084
-0.0096 -0.0158 -0.0003 0.0219 0.0285
0.0181 0.0030 -0.0022 0.0029 0.0054
-0.0011
PH 0.0000 14.6424 29.0047 57.2208 85.1313
112.8778 168.0043 222.7705 277.1376 331.0602
384.4995 437.4410 489.8475 541.6945 643.7048
743.5847 841.2031 936.4425 1029.3514 1164.3414
1293.2196 1450.1709 1647.4191 1869.0129 2107.1233
2377.2019 2688.4321 3040.7261 3434.5886 3844.6377
4247.0693 4712.2925 5303.2905 6377.2554 7210.1094
9052.6641
PHB 765.1553 962.6146 1160.4569 1557.3174 1955.7478
2355.7827 3160.6733 3972.0894 4790.1455 5614.9644
6446.6768 7285.3921 8131.2476 8984.3877 10713.0273
12472.4590 14263.9023 16088.6348 17948.0312 20805.4141
23749.4668 27820.1406 33164.7578 38828.4688 44855.9414
51302.3164 58236.8750 65748.4531 73953.8516 83011.2188
93142.9844 104677.5859 118132.9844 134406.7656 157094.2812
195157.4688
P 1326.0703 1312.6719 1300.0781 1289.3750 1280.6797
1265.7578 1246.5000 1226.9219 1207.9922 1188.6797
1168.3438 1148.8672 1129.0078 1099.9062 1060.5156
1021.5156 982.2031 943.2969 893.4297 836.0000
775.1797 715.8516 658.0469 602.6680 546.9102
491.5781 437.5586 384.9062 332.9883 281.2441
229.7676 178.5605 127.6582 76.9062 25.6348
PB 98966.5547 98731.3516 98378.5313 97908.1094 97437.6953
96732.0625 95791.2188 94850.3750 93909.5312 92968.6953
92027.8516 91087.0156 90146.1641 88734.9062 86853.2266
84971.5391 83089.8516 81208.1719 78856.0625 76033.5469
72740.6016 68506.8047 63802.6016 59098.3984 54394.1836
49689.9766 44985.7695 40281.5625 35577.3516 30873.1465
26168.9355 21464.7305 16760.5195 12056.3125 7352.1045
THETA 304.3011 305.0110 305.4361 305.6014 305.6483
305.6881 305.7391 305.7730 305.7908 305.8033
305.8008 305.7959 305.7897 305.7814 305.8174
305.8354 305.8592 305.9174 306.0104 306.1797
308.4571 311.1703 314.5753 316.9744 320.1039
323.7825 327.3818 330.8230 333.1943 334.8087
338.6066 343.4999 354.4267 377.0505 438.0953
T 304.5553 305.0494 305.1553 304.8987 304.5239
303.9268 303.1236 302.2969 301.4485 300.5883
299.7058 298.8154 297.9168 296.5580 294.7667
292.9295 291.0674 289.2069 286.8524 284.0107
282.4996 280.1269 277.4837 273.5389 269.7590
265.8854 261.2969 255.8285 248.6665 239.9319
231.4373 221.8293 213.2238 206.3851 208.0262
QVAPOR 0.0146 0.0123 0.0108 0.0101 0.0100
0.0099 0.0098 0.0098 0.0098 0.0098
0.0098 0.0098 0.0098 0.0098 0.0098
0.0098 0.0098 0.0098 0.0098 0.0097
0.0046 0.0019 0.0008 0.0009 0.0010
0.0008 0.0004 0.0003 0.0002 0.0002
0.0001 0.0001 0.0000 0.0000 0.0000
QCLOUD 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
QRAIN 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
QICE 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
QSNOW 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
QGRAUP 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
MU 1025.3440 Pa
MUB 94084.1641 Pa
T2 304.1827 K
Q2 0.0163 kg/kg
U10 -0.4512 m/s
V10 -0.9520 m/s
PSFC 100402.5469 Pa
MAPFAC_M 0.9992 dimensionless
MAPFAC_U 0.9991 dimensionless
MAPFAC_V 0.9992 dimensionless
HGT 77.9975 m
RAINC 0.0000 cm
RAINNC 0.0000 cm
SWDOWN 113.8828 W m^-2
GLW 399.8240 W m^-2
XLAT 33.4891 degrees
XLAT_U 33.4980 degrees
XLAT_V 33.4358 degrees
XLONG -81.9709 degrees
XLONG_U -82.0349 degrees
XLONG_V -81.9816 degrees
LU_INDEX 24 category
LANDMASK 1.0000 category
HFX -9.0076 W m^-2
LH 70.7717 W m^-2
UST 0.1068 m
MOL 30.9906 m
QFX 0.0000 kg m-2 s-1
PBLH 85.5292 m
RA 95.0254 s m^-1
RS 233.7984 s m^-1
LAI 2.9035 area/area
CANWAT 0.0000 kg m^-2
VEGF_PX 0.6907 area/area
ISLTYP 3 category
SMOIS 1 0.0399 m^3 m^-3
SMOIS 2 0.1896 m^3 m^-3
SMOIS 0.0399 0.1896
TSLB 1 303.5585 K
TSLB 2 299.1393 K
TSLB 303.5585 299.1393
TSK 303.5585 K
ALBEDO 0.1222 fraction
LANDUSEF 0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000
0.0000 0.0591 0.0000 0.0000 0.0000
0.0000 0.0000 0.1679 0.2715 0.1164
0.0577 0.0020 0.0818 0.0761 0.0103
0.0000 0.0095 0.0375 0.0000 0.0000
0.0000 0.0197 0.0192 0.0637 0.0074
**********************************************************************
*** SUBROUTINE: RDWRFEM
*** ERROR RETRIEVING VARIABLE FROM WRF FILE
*** VARIABLE = F
*** RCODE = NetCDF: Variable not found
**********************************************************************
*** ERROR ABORT in subroutine RDWRFEM
ABNORMAL TERMINATION IN RDWRFEM
Date and time 0:00:00 July 1, 2011 (2011182:000000)
Error running mcip
@lizadams
Copy link
Author

lizadams commented Apr 11, 2019

The variables that start with SNOW in the wrf input files are as follows:

grep SNOW  subset_wrfout_d01_2011-07-01_00:00:00.header
	float QSNOW(Time, bottom_top, south_north, west_east) ;
		QSNOW:FieldType = 104 ;
		QSNOW:MemoryOrder = "XYZ" ;
		QSNOW:description = "Snow mixing ratio" ;
		QSNOW:units = "kg kg-1" ;
		QSNOW:stagger = "" ;
		QSNOW:coordinates = "XLONG XLAT XTIME" ;
	float SNOWC(Time, south_north, west_east) ;
		SNOWC:FieldType = 104 ;
		SNOWC:MemoryOrder = "XY " ;
		SNOWC:description = "FLAG INDICATING SNOW COVERAGE (1 FOR SNOW COVER)" ;
		SNOWC:units = "" ;
		SNOWC:stagger = "" ;
		SNOWC:coordinates = "XLONG XLAT XTIME" ;

So, I don't have the SNOWH variable in the input file that comes with the BENCHMARK case in the following tar file: CMAQv5.2_Benchmark_SingleDay_Input_09_12_2017.tar.gz

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment