Two-way coupled WRF-CMAQ process?

I am learning and running two-way coupled model, I can successfully run two-way coupled model with these setting (a. DO_SW_CAL = T, DO_LW_CAL = F; b. DO_SW_CAL = F, DO_LW_CAL = T; c.DO_SW_CAL = F, DO_LW_CAL = F), but can’t run with the setting (DO_SW_CAL = T, DO_LW_CAL = T) and have no errors. The attachment is configuration file. Can you tell me how to address this? Thanks in advance for your help!

neigae_20190620_sw_lw.csv (39.5 KB)

Hi gc13141112,

 Thanks for your interest in the WRF-CMAQ two-way coupled model. We did not fully tested the long wave option (LW_CAL) so it is not officially supported. Unfortunately the model crashed when LW_CAL and SW_CAL were both on. Were you using the test dataset? Anyway I will take a look of the code.

Cheers,
David

Hi David,
Thank you for your reply. In the twoway/phys folder, I can’t find twoway longwave code that looks like module_twoway_ra_rrtmg_lw.F. Also I can’t find any code related to the twoway CAM scheme. Please advise!

                                                                                                                Chao

Hi
I want to benchmark and run WRF-CMAQ TWO-WAY,but I don’t know how to do it? Could you guide me, please?
Thank you

Hi Ensiyeh,

Please download the required WRF-CMAQ model tarball which includes a step-by-step instruction to build the WRF-CMAQ two-way coupled model.

Hi David

Thank you for your response. I built WRF-CMAQ two-way model, but I don’t know how to benchmark it and how to run it? Could you guide me, please?
Thank you

Hi Ensiyeh,

I assume you are using the latest version WRF 4.1.1 and CMAQ 5.3 coupled model. You can download a one-day test dataset and a run script.

Hi David,

I built my model with WRF 3.8 and CMAQ 5.2. I downloaded WRFv38-CMAQv5.2-benchmark; But I don’t know how to benchmark and run this model. I know how to run CMAQ only. Could you guide me, please? Is there any user guidelines for running it?
Thank you

Hi Ensiyeh,
Please consider moving to the latest WRF-CMAQ release which uses WRFv4.1.1 and CMAQv5.3. A tarball with build and run scripts and a README file with instructions is available here (Google Drive): https://drive.google.com/drive/folders/1ru64xxZeRcE5buQ9fygV9jUHAQ9aDnL7

One day benchmark input and output datasets are available to download from here (Google Drive): https://drive.google.com/drive/folders/1apg7_LXU4Kpzx5EZgdc07Q3UTTpL-jah

The tarball and benchmark input data are also available from the EPA anonymous ftp server: ftp://newftp.epa.gov/exposure/CMAQ/V5_3/WRF-CMAQ_Coupled_Model/
The benchmark output data is not there yet but will be uploaded sometime next week.

Hi

Thank you. I download them. I can’t open 3rd link. Is that right?

Hi Ensiyeh,

If you have run CMAQ before, you can submit WRf-CMAQ job script the same way. The only thing you need to do is to modify the input, output, and source code paths in the run script.

The third link is correct and should be a publicly available link but perhaps it is being blocked by your system admin because it is anonymous ftp (not secure). The files on the ftp server are the same as what is posted on the Google Drive. If you were able to download them from the Google Drive folder, you do not need the third link. We provided this option for users that are unable to access Google Drive.

1 Like

Hi David,

Thank you. You mean I must just modify a scripts? Which script? I see just a scrip in twoway/scripts
Where I must put the inputs file?

Hi Ensiyeh,

When you unzip twoway.tar.gz, there is a run script in twoway/script. You can unzip the dataset and put it anywhere your want.

Hi David,

Thank you for your answer. I want to benchmark my model. I ran “twoway_model_run_script” and it gives me this:

mpirun -r ssh -np 8 wrf.exe
mpirun: Command not found.
0.0u 0.0s 0:00.00 0.0% 0+0k 0+0io 0pf+0w
date +Completed WRF at %m/%d/%y %H:%M:%S%n
Completed WRF at 09/04/19 18:57:52
mv: No match.
if ( F == T ) then
biogon: Undefined variable.

The log file is:
New Text Document.txt (19.4 KB)

There are another problems,too. Could you please check my logfile?
Could you guide me, please?
Thank you

Hi Ensiyeh,

According to the script, there is no biogon variable in the script. Please send me your script.

Hi David,

I don’t have access to the script that I modified, but I use this script:

New Text Document.txt (38.7 KB)

Hi Ensiyeh,

 Replace line 886:

    if ( $biogon ) then

with this line

if ( $CTM_BIOGEMIS == ‘Y’ ) then

Hi David,

Thank you for your answer. I did it but I couldn’t run it. Could you guide me, please?
The log file is:

set NPROCS = 8
set start_year = 2013
set end_year = 2013
set start_month = 8
set end_month = 8
set start_day = 1
set end_day = 1
set cycle = 0
set firstday = 20130728
set STTIME = 0
set NSTEPS = 240000
set TSTEP = 010000
set version = 3852
set MECH = cb6r3_ae6_aq
set resolution = 12km
setenv NZ 35
set ROOT_DIR = /home/ensiyeh/online
set WRF_DIR = /home/ensiyeh/online/WRFV38
set INPDIR = /home/ensiyeh/online/wrf_38_cmaq_52_input
set OMIpath = /home/ensiyeh/online/WRFV38/cmaq
set OUTPUT_ROOT = /home/ensiyeh/online
set NMLpath = /home/ensiyeh/online/WRFV38/cmaq
set NMLpath2 = /home/ensiyeh/online/WRFV38/cmaq
set NUM_LAND_USE_TYPE = 40
set cont_from_spinup_run = T
set radt = 20
setenv RUN_CMAQ_DRIVER T
setenv DO_SW_CAL T
setenv DO_LW_CAL F
setenv WRF_CMAQ_FREQ 5
setenv CREATE_PHYSICAL_FILE F
setenv FILE_TIME_STEP 10000
setenv WRF_COL_DIM 251
setenv WRF_ROW_DIM 185
setenv WRF_LAY_DIM 36
setenv CMAQ_COL_DIM 236
setenv CMAQ_ROW_DIM 170
setenv TWOWAY_DELTA_X 7
setenv TWOWAY_DELTA_Y 7
setenv WRF_LC_REF_LAT 40.0
setenv CTM_TURN_ON_PV N
setenv SD_TIME_SERIES F
setenv SD_SCOL 241
setenv SD_ECOL 248
setenv SD_SROW 160
setenv SD_EROW 169
setenv SD_CONC_SPCS NO2 NO O3 NO3 CO ASO4J ASO4I ANH4J ANH4I ANO3J
ANO3I AORGAJ AORGAI AORGPAJ AORGPAI AORGBJ AORGBI AECJ AECI A25J A25I
ACORS ASEAS ASOIL
setenv AVG_CONC_SPCS ALL
setenv CONC_SPCS
setenv CONC_BLEV_ELEV 1 35
setenv ACONC_BLEV_ELEV 1 1
setenv CTM_PROCAN N
setenv PA_BCOL_ECOL 30 100
setenv PA_BROW_EROW 10 50
setenv PA_BLEV_ELEV 1 5
setenv CTM_MAXSYNC 300
setenv CTM_MINSYNC 90
setenv SIGMA_SYNC_TOP 0.7
setenv CTM_ADV_CFL 0.95
setenv AVG_FILE_ENDTIME Y
setenv CTM_CKSUM Y
setenv CTM_WB_DUST N
setenv CTM_ERODE_AGLAND N
setenv CTM_WBDUST_BELD BELD3
setenv CTM_LTNG_NO N
setenv CTM_WVEL N
setenv KZMIN Y
setenv CTM_ILDEPV Y
setenv CTM_MOSAIC N
setenv CTM_ABFLUX N
setenv CTM_HGBIDI N
setenv CTM_SFC_HONO N
setenv CTM_PING N
setenv CTM_PT3DEMIS Y
setenv CTM_FST N
setenv CTM_BIOGEMIS N
setenv PX_VERSION Y
setenv BIOSW_YN N
setenv SUMMER_YN N
setenv AVG_FILE_ENDTIME N
setenv CTM_GRAV_SETL Y
setenv CORESHELL_OPTICS N
setenv MIE_OPTICS N
setenv LTNGPARAM N
setenv LTNG_ASSIM_FLAG N
setenv CTM_ZERO_PCSOA N
setenv IOAPI_LOG_WRITE F
setenv FL_ERR_STOP F
setenv IOAPI_OFFSET_64 Y
setenv CTM_EMISCHK N
setenv PROMPTFLAG F
setenv CLD_DIAG N
setenv CTM_PMDIAG Y
setenv CTM_APMDIAG Y
setenv CTM_AVISDIAG Y
setenv CTM_PHOTDIAG N
setenv CTM_SSEMDIAG N
setenv CTM_DUSTEM_DIAG N
setenv CTM_DEPV_FILE N
setenv B3GTS_DIAG N
setenv APT_DIAG N
setenv PT3DDIAG N
setenv PT3DFRAC N
setenv LTNGDIAG N
setenv VDIFF_DIAG_FILE N
setenv CTM_AOD N
setenv CTM_MGEMDIAG N
setenv REP_LAYER_MIN -1
setenv MPI_SM_POOL 16000
setenv MP_EAGER_LIMIT 65536
setenv MP_SINGLE_THREAD yes
setenv MP_STDOUTMODE ordered
setenv MP_LABELIO yes
setenv MP_SHARED_MEMORY yes
setenv MP_ADAPTER_USE shared
setenv MP_CPU_USE multiple
setenv MP_CSS_INTERRUPT yes
set OCEANfile = ssmask.4BEARPEX2.ncf
set OMIfile = OMI_1979_to_2015.dat
set OPTICSfile = PHOT_OPTICS.dat
set CSQYfile = CSQY_DATA_cb6r3_ae6_aq
set ICFILE = ICON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
set BCFILE = BCON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv RUNID v3852
set ICpath = /home/ensiyeh/online/wrf_38_cmaq_52_input/icbc
set BCpath = /home/ensiyeh/online/wrf_38_cmaq_52_input/icbc
set EMISpath = /home/ensiyeh/online/wrf_38_cmaq_52_input/emis
set IN_PTpath = /home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point
set IN_LTpath = /home/ensiyeh/online/wrf_38_cmaq_52_input/lightning
set METpath = /home/ensiyeh/online/wrf_38_cmaq_52_input/met_wrf
set LUpath = /home/ensiyeh/online/wrf_38_cmaq_52_input/land
set SZpath = /home/ensiyeh/online/wrf_38_cmaq_52_input/land
if ( ( T == T ) || ( F == T ) ) then
setenv CMAQ_WRF_FEEDBACK T
else
set rad_scheme = rrtmg_20
if ( T == T ) then
if ( F == T ) then
set feedback = sf
endif
else
set call_ratio = 5_1
set output_direct_name = output_12km_sf_rrtmg_20_5_1_v3852_2013
set EXEC = CCTM
set year = 2013
@ pre_year = 2013 - 1
@ cur_year = 2013
@ yr = 2013
@ start_hour = 0 / 10000
@ end_hour = ( 0 + 23 ) % 24
@ cur_year = 2013 - 1
while ( 2012 < 2013 )
@ cur_year++
if ( ( 2013 % 4 ) != 0 ) then
set leap_year_cur_year = 0
else if ( ( 2013 % 400 ) == 0 ) then
if ( 0 ) then
set numday_cur_year = ( 31 28 31 30 31 30 31 31 30 31 30 31 )
if ( ( 2012 % 4 ) != 0 ) then
if ( ( 2012 % 400 ) == 0 ) then
if ( ( 2012 % 100 ) == 0 ) then
set leap_year_pre_year = 1
endif
endif
if ( 2013 == 2013 ) then
@ cur_day_j = 0
@ temp_month = 1
while ( 1 < 8 )
@ cur_day_j = 0 + 31
@ temp_month++
end
while ( 2 < 8 )
@ cur_day_j = 31 + 28
@ temp_month++
end
while ( 3 < 8 )
@ cur_day_j = 59 + 31
@ temp_month++
end
while ( 4 < 8 )
@ cur_day_j = 90 + 30
@ temp_month++
end
while ( 5 < 8 )
@ cur_day_j = 120 + 31
@ temp_month++
end
while ( 6 < 8 )
@ cur_day_j = 151 + 30
@ temp_month++
end
while ( 7 < 8 )
@ cur_day_j = 181 + 31
@ temp_month++
end
while ( 8 < 8 )
@ cur_day_j = 212 + 1
@ pre_day_j = 213 - 1
else
if ( 213 == 1 ) then
if ( 2013 == 2013 ) then
@ cur_month = 8
else
if ( 2013 == 2013 ) then
@ cur_end_month = 8
else
@ cur_month = 8 - 1
while ( 7 < 8 )
@ cur_month++
@ cur_day = 0
@ cur_month_end_day = 31
if ( ( 2013 == 2013 ) && ( 8 == 8 ) ) then
@ cur_day = 1 - 1
endif
if ( ( 2013 == 2013 ) && ( 8 == 8 ) ) then
@ cur_month_end_day = 1
endif
while ( 0 < 1 )
@ cur_day++
set loc_pre_year = 2013
if ( 1 == 1 ) then
@ pre_month = 8 - 1
if ( 7 == 0 ) then
set yest = 31
set yesterday = printf "%d%2.2d%2.2d\n" $loc_pre_year $pre_month $yest
printf %d%2.2d%2.2d\n 2013 7 31
else
set today = printf "%d%2.2d%2.2d\n" $cur_year $cur_month $cur_day
printf %d%2.2d%2.2d\n 2013 8 1
set YYYYMMDD = 20130801
set G_DAY = 20130801
set today_j = printf "%d%3.3d\n" $cur_year $cur_day_j
printf %d%3.3d\n 2013 213
set yesterday_j = printf "%d%3.3d\n" $loc_pre_year $pre_day_j
printf %d%3.3d\n 2013 212
set wrf_tstep = 20
setenv CTM_APPL v3852_20130801
setenv CTM_APPL_yesterday v3852_20130731
if ( ! -d /home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013 ) then
set use_profile_init = 0
if ( 20130801 == 20130728 ) then
if ( T == T ) then
setenv FEEDBACK_RSTFLAG .false.
set cont_from_spinup_run = F
set use_profile_init = 1
else
setenv WRF_RSTFLAG .true.
setenv pxlsm_smois_init 0
endif
mkdir /home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/20130801
mkdir: cannot create directory
‘/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/20130801’:
File exists
set STDATE = 2013213
set APPL = 20130801
setenv GRIDDESC
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/GRIDDESC
setenv GRID_NAME WRF_CMAQ_2WAY
set DISP = delete
set OUTDIR = /home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013
set YEAR = 2013
set MONTH = printf "%2.2d\n" $cur_month
printf %2.2d\n 8
set day = 1
set EMISDATE = printf "%4d%2.2d%2.2d\n" ${YEAR} $cur_month $cur_day
printf %4d%2.2d%2.2d\n 2013 8 1
if ( 1 ) then
setenv INITIAL_RUN Y
else
set test = ls CTM_LOG_???.${APPL}
ls: No match.
if ( != ) then
if ( N == Y ) then
if ( N == Y ) then
if ( N == Y ) then
if ( N == Y ) then
setenv CSQY_DATA /home/ensiyeh/online/WRFV38/cmaq/CSQY_DATA_cb6r3_ae6_aq
if ( ! ( -e /home/ensiyeh/online/WRFV38/cmaq/CSQY_DATA_cb6r3_ae6_aq ) ) then
setenv OCEAN_1 /home/ensiyeh/online/wrf_38_cmaq_52_input/land/ssmask.4BEARPEX2.ncf
setenv EMISfile emis_mole_all_20130801_4BEARPEX2_cmaq_cb05e51_2013ej_v6_13i.ncf
setenv EMIS_1 /home/ensiyeh/online/wrf_38_cmaq_52_input/emis/emis_mole_all_20130801_4BEARPEX2_cmaq_cb05e51_2013ej_v6_13i.ncf
if ( Y == Y || Y == T ) then
setenv NPTGRPS 5
set iday = grep ^$STDATE $INPDIR/mergedates.2013.txt
grep ^2013213 /home/ensiyeh/online/wrf_38_cmaq_52_input/mergedates.2013.txt
setenv STK_GRPS_01
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/stack_groups_ptnonipm_12US1_2013ej_v6_13i.ncf
setenv STK_EMIS_01
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/inln_mole_ptnonipm_20130806_12US1_cmaq_cb05e51_2013ej_v6_13i.ncf
setenv STK_GRPS_02
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/stack_groups_othpt_12US1_2013ej_v6_13i.ncf
setenv STK_EMIS_02
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/inln_mole_othpt_20130806_12US1_cmaq_cb05e51_2013ej_v6_13i.ncf
setenv STK_GRPS_03
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/stack_groups_ptegu_12US1_2013ej_v6_13i.ncf
setenv STK_EMIS_03
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/inln_mole_ptegu_20130801_12US1_cmaq_cb05e51_2013ej_v6_13i.ncf
setenv STK_GRPS_04
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/stack_groups_pt_oilgas_12US1_2013ej_v6_13i.ncf
setenv STK_EMIS_04
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/inln_mole_pt_oilgas_20130806_12US1_cmaq_cb05e51_2013ej_v6_13i.ncf
setenv STK_GRPS_05
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/stack_groups_ptfire_20130801_12US1_2013ej_v6_13i.ncf
setenv STK_EMIS_05
/home/ensiyeh/online/wrf_38_cmaq_52_input/emis/inln_point/inln_mole_ptfire_20130801_12US1_cmaq_cb05e51_2013ej_v6_13i.ncf
unset echo
setenv LAYP_STDATE 2013213
setenv LAYP_STTIME 0
setenv LAYP_NSTEPS 240000
endif
setenv S_CGRID /home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_CGRID_v3852
20130801.nc
setenv CTM_CONC_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_CONC_v3852_20130801.nc
-v
setenv A_CONC_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_ACONC_v3852_20130801.nc
-v
setenv MEDIA_CONC
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_MEDIA_CONC_v3852_20130801.nc
-v
setenv CTM_DRY_DEP_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_DRYDEP_v3852_20130801.nc
-v
setenv CTM_DEPV_DIAG
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_DEPV_v3852_20130801.nc
-v
setenv CTM_PT3D_DIAG
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_PT3D_v3852_20130801.nc
-v
setenv B3GTS_S /home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_B3GTS_S_v3852_20130801.nc
-v
setenv SOILOUT /home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_SOILOUT_v3852_20130801.nc
setenv CTM_WET_DEP_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_WETDEP1_v3852_20130801.nc
-v
setenv CTM_WET_DEP_2
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_WETDEP2_v3852_20130801.nc
-v
setenv CTM_VIS_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_PMVIS_v3852_20130801.nc
-v
setenv CTM_AVIS_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_APMVIS_v3852_20130801.nc
-v
setenv CTM_PMDIAG_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_PMDIAG_v3852_20130801.nc
-v
setenv CTM_APMDIAG_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_APMDIAG_v3852_20130801.nc
-v
setenv CTM_RJ_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_PHOTDIAG1_v3852_20130801.nc
-v
setenv CTM_RJ_2
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_PHOTDIAG2_v3852_20130801.nc
-v
setenv CTM_SSEMIS_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_SSEMIS.v3852_20130801.nc
-v
setenv CTM_DUST_EMIS_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_DUSTEMIS.v3852_20130801.nc
-v
setenv CTM_IPR_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_PA_1_v3852_20130801.nc
-v
setenv CTM_IPR_2
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_PA_2_v3852_20130801.nc
-v
setenv CTM_IPR_3
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_PA_3_v3852_20130801.nc
-v
setenv CTM_IRR_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_IRR_1_v3852_20130801.nc
-v
setenv CTM_IRR_2
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_IRR_2_v3852_20130801.nc
-v
setenv CTM_IRR_3
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_IRR_3_v3852_20130801.nc
-v
setenv CTM_DRY_DEP_MOS
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_DDMOS_v3852_20130801.nc
-v
setenv CTM_DRY_DEP_FST
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_DDFST_v3852_20130801.nc
-v
setenv CTM_DEPV_MOS
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_DEPVFST_v3852_20130801.nc
-v
setenv CTM_DEPV_FST
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_DEPVMOS_v3852_20130801.nc
-v
setenv CTM_VDIFF_DIAG
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_VDIFF_DIAG_v3852_20130801.nc
-v
setenv CTM_VSED_DIAG
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_VSED_DIAG_v3852_20130801.nc
-v
setenv CTM_AOD_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_AOD_DIAG_v3852_20130801.nc
-v
setenv CTM_LTNGDIAG_1
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_LTNGHRLY_v3852_20130801.nc
-v
setenv CTM_LTNGDIAG_2
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/CCTM_LTNGCOL_v3852_20130801.nc
-v
setenv FLOOR_FILE
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/FLOOR_v3852_20130801.txt
setenv LAYER_FILE MET_CRO_3D
setenv GRID_DOT_2D BUFFERED
setenv GRID_CRO_2D BUFFERED
setenv MET_CRO_2D BUFFERED
setenv MET_DOT_3D BUFFERED
setenv MET_CRO_3D BUFFERED
@ n = 0
while ( 0 < 8 )
set name = printf "_%3.3d\n" $n
printf _%3.3d\n 0
setenv feed_back_000 BUFFERED
@ n++
end
while ( 1 < 8 )
set name = printf "_%3.3d\n" $n
printf _%3.3d\n 1
setenv feed_back_001 BUFFERED
@ n++
end
while ( 2 < 8 )
set name = printf "_%3.3d\n" $n
printf _%3.3d\n 2
setenv feed_back_002 BUFFERED
@ n++
end
while ( 3 < 8 )
set name = printf "_%3.3d\n" $n
printf _%3.3d\n 3
setenv feed_back_003 BUFFERED
@ n++
end
while ( 4 < 8 )
set name = printf "_%3.3d\n" $n
printf _%3.3d\n 4
setenv feed_back_004 BUFFERED
@ n++
end
while ( 5 < 8 )
set name = printf "_%3.3d\n" $n
printf _%3.3d\n 5
setenv feed_back_005 BUFFERED
@ n++
end
while ( 6 < 8 )
set name = printf "_%3.3d\n" $n
printf _%3.3d\n 6
setenv feed_back_006 BUFFERED
@ n++
end
while ( 7 < 8 )
set name = printf "_%3.3d\n" $n
printf _%3.3d\n 7
setenv feed_back_007 BUFFERED
@ n++
end
while ( 8 < 8 )
setenv PGRID_DOT_2D
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/GRID_DOT_2D20130801
-v
setenv PGRID_CRO_2D
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/GRID_CRO_2D20130801
-v
setenv PMET_CRO_2D
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/MET_CRO_2D20130801
-v
setenv PMET_DOT_3D
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/MET_DOT_2D20130801
-v
setenv PMET_CRO_3D
/home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013/MET_CRO_3D20130801
-v
if ( F == T ) then
setenv CTM_STDATE 2013213
setenv CTM_STTIME 0
setenv CTM_RUNLEN 240000
setenv CTM_TSTEP 010000
setenv CTM_PROGNAME CCTM
setenv INIT_GASC_1
/home/ensiyeh/online/wrf_38_cmaq_52_input/icbc/ICON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv INIT_AERO_1
/home/ensiyeh/online/wrf_38_cmaq_52_input/icbc/ICON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv INIT_NONR_1
/home/ensiyeh/online/wrf_38_cmaq_52_input/icbc/ICON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv INIT_TRAC_1
/home/ensiyeh/online/wrf_38_cmaq_52_input/icbc/ICON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv BNDY_GASC_1
/home/ensiyeh/online/wrf_38_cmaq_52_input/icbc/BCON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv BNDY_AERO_1
/home/ensiyeh/online/wrf_38_cmaq_52_input/icbc/BCON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv BNDY_NONR_1
/home/ensiyeh/online/wrf_38_cmaq_52_input/icbc/BCON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv BNDY_TRAC_1
/home/ensiyeh/online/wrf_38_cmaq_52_input/icbc/BCON.2013213.4BEARPEX2.2013ej_v6_13i.ncf
setenv OMI /home/ensiyeh/online/WRFV38/cmaq/OMI_1979_to_2015.dat
setenv OPTICS_DATA /home/ensiyeh/online/WRFV38/cmaq/PHOT_OPTICS.dat
cd /home/ensiyeh/online/output_12km_sf_rrtmg_20_5_1_v3852_2013
ln -s /home/ensiyeh/online/WRFV38/cmaq/GC_cb6r3_ae6_aq.nml gc_matrix_nml
ln: failed to create symbolic link ‘gc_matrix_nml’: File exists
ln -s /home/ensiyeh/online/WRFV38/cmaq/AE_cb6r3_ae6_aq.nml ae_matrix_nml
ln: failed to create symbolic link ‘ae_matrix_nml’: File exists
ln -s /home/ensiyeh/online/WRFV38/cmaq/NR_cb6r3_ae6_aq.nml nr_matrix_nml
ln: failed to create symbolic link ‘nr_matrix_nml’: File exists
ln -s /home/ensiyeh/online/WRFV38/cmaq/Species_Table_TR_0.nml tr_matrix_nml
ln: failed to create symbolic link ‘tr_matrix_nml’: File exists
if ( -f namelist.input ) rm -f namelist.input
rm -f namelist.input
cat
rm -f wrfbdy_d01 wrffdda_d01 wrfsfdda_d01 wrfinput_d01 wrflowinp_d01
ln -s /home/ensiyeh/online/wrf_38_cmaq_52_input/met_wrf/wrfbdy_d01 wrfbdy_d01
ln -s /home/ensiyeh/online/wrf_38_cmaq_52_input/met_wrf/wrffdda_d01 wrffdda_d01
ln -s /home/ensiyeh/online/wrf_38_cmaq_52_input/met_wrf/wrfsfdda_d01
wrfsfdda_d01
if ( .true. == .false. ) then
ln -s /home/ensiyeh/online/wrf_38_cmaq_52_input/met_wrf/wrfrst_d01_2013-08-01_00:00:00
wrfrst_d01_2013-08-01_00:00:00
ln: failed to create symbolic link ‘wrfrst_d01_2013-08-01_00:00:00’: File exists
endif
ln -s /home/ensiyeh/online/wrf_38_cmaq_52_input/met_wrf/wrflowinp_d01
wrflowinp_d01
if ( -f wrf.exe ) rm -f wrf.exe
rm -f wrf.exe
if ( -f ETAMPNEW_DATA ) rm -f ETAMPNEW_DATA
rm -f ETAMPNEW_DATA
if ( -f GENPARM.TBL ) rm -f GENPARM.TBL
rm -f GENPARM.TBL
if ( -f landFilenames ) rm -f landFilenames
rm -f landFilenames
if ( -f LANDUSE.TBL ) rm -f LANDUSE.TBL
rm -f LANDUSE.TBL
if ( -f RRTM_DATA ) rm -f RRTM_DATA
rm -f RRTM_DATA
if ( -f SOILPARM.TBL ) rm -f SOILPARM.TBL
rm -f SOILPARM.TBL
if ( -f tr49t67 ) rm -f tr49t67
rm -f tr49t67
if ( -f tr49t85 ) rm -f tr49t85
rm -f tr49t85
if ( -f tr67t85 ) rm -f tr67t85
rm -f tr67t85
if ( -f VEGPARM.TBL ) rm -f VEGPARM.TBL
rm -f VEGPARM.TBL
ln -s /home/ensiyeh/online/WRFV38/main/wrf.exe wrf.exe
ln -s /home/ensiyeh/online/WRFV38/test/em_real/ETAMPNEW_DATA ETAMPNEW_DATA
ln -s /home/ensiyeh/online/WRFV38/test/em_real/GENPARM.TBL GENPARM.TBL
ln -s /home/ensiyeh/online/WRFV38/test/em_real/landFilenames landFilenames
ln -s /home/ensiyeh/online/WRFV38/test/em_real/LANDUSE.TBL LANDUSE.TBL
ln -s /home/ensiyeh/online/WRFV38/test/em_real/RRTM_DATA RRTM_DATA
ln -s /home/ensiyeh/online/WRFV38/test/em_real/RRTMG_SW_DATA RRTMG_SW_DATA
ln: failed to create symbolic link ‘RRTMG_SW_DATA’: File exists
ln -s /home/ensiyeh/online/WRFV38/test/em_real/RRTMG_LW_DATA RRTMG_LW_DATA
ln: failed to create symbolic link ‘RRTMG_LW_DATA’: File exists
ln -s /home/ensiyeh/online/WRFV38/test/em_real/SOILPARM.TBL SOILPARM.TBL
ln -s /home/ensiyeh/online/WRFV38/test/em_real/tr49t67 tr49t67
ln -s /home/ensiyeh/online/WRFV38/test/em_real/tr49t85 tr49t85
ln -s /home/ensiyeh/online/WRFV38/test/em_real/tr67t85 tr67t85
ln -s /home/ensiyeh/online/WRFV38/test/em_real/VEGPARM.TBL VEGPARM.TBL
ln -s /home/ensiyeh/online/WRFV38/test/em_real/ozone_plev.formatted
ozone_plev.formatted
ln: failed to create symbolic link ‘ozone_plev.formatted’: File exists
ln -s /home/ensiyeh/online/WRFV38/test/em_real/ozone_lat.formatted
ozone_lat.formatted
ln: failed to create symbolic link ‘ozone_lat.formatted’: File exists
ln -s /home/ensiyeh/online/WRFV38/test/em_real/ozone.formatted ozone.formatted
ln: failed to create symbolic link ‘ozone.formatted’: File exists
date +Started wrf.exe at %m/%d/%y %H:%M:%S
Started wrf.exe at 09/05/19 17:21:34
set MPI = /home/ensiyeh/WRF-CMAQ/CMAQ-OFFLINE/CMAQ-5.2/lib/x86_64/intel/mpich/bin
set MPIRUN = /home/ensiyeh/WRF-CMAQ/CMAQ-OFFLINE/CMAQ-5.2/lib/x86_64/intel/mpich/bin/mpirun
mpirun -np 8 /home/ensiyeh/online/WRFV38/main/wrf.exe
mpirun: Command not found.
0.0u 0.0s 0:00.00 0.0% 0+0k 0+0io 0pf+0w
date +Completed WRF at %m/%d/%y %H:%M:%S%n
Completed WRF at 09/05/19 17:21:34

mv: No match.
if ( F == T ) then
if ( N == M-bM^@M^XYM-bM^@M^Y ) then
@ cur_day_j++
@ pre_day_j++
if ( 213 > 214 ) then
end
while ( 1 < 1 )
end
while ( 8 < 8 )
@ pre_year++
end
while ( 2013 < 2013 )

Thank you

You have an issue because the script cannot find the location or path of the mpirun command.
Try the following command

which mpirun

If you get the message “no mpirun in (…)”
then you need to set up your mpi environment
If you are running on a server that has Environment Modules use the following command to find out what modules are available:
Modules let you load predefined environments, and in this case they’re used to load different MPI and compiler environments.

module avail

If there is a module for openmpi and gcc such as the following module: openmpi_4.0.1/gcc_9.1.0 you can load that module using the following:

module load openmpi_4.0.1/gcc_9.1.0

If you do not have access to a server that allows you to load modules to configure your environment, then you will need to install openmpi or intel mpi and then once that is installed add the path to your mpirun command to your PATH.
You can try to follow these instructions: