Run the coupled case

Case initialization

To test the coupled case, enter runCase.init folder

[ruisun@acc00]~/scripps_kaust_model/coupler/L3.C1.coupled_RS2012_ring$ cd runCase.init/
[ruisun@acc00]~/.../coupler/L3.C1.coupled_RS2012_ring/runCase.init$ ls
Allclean  Allrun  namelist.input  README

Initialize the test case

[ruisun@acc00]~/.../coupler/L3.C1.coupled_RS2012_ring/runCase.init$ ./Allrun

Run the coupled case

Note

I assume the user is familiar with both MITgcm and WRF. Before running the coupled simulations. Please make sure that the WRF and MITgcm setups are the same as the coupled model, including data,``data.exf``,``data.cal`` for MITgcm and namelist.input for WRF.

Enter the runCase folder and run the test case:

[ruisun@acc00]~/.../coupler/L3.C1.coupled_RS2012_ring/runCase.init$ cd ../runCase/
[ruisun@acc00]~/.../coupler/L3.C1.coupled_RS2012_ring/runCase$ ./Allrun

If the coupled code is compiled successfully, the coupled model will run and generate coupled simulation results.

Something more about the coupled model

In the Allrun file, MATLAB is required for pre-processing (the python scripts are available in the L3.C3 case for SHAHEEN at KAUST):

\3 ln -sf ../caseInput/* .
\4
\5 cp ../../../WRFV3911_AO/main/wrf.exe .
\6 cp ../runCase.init/wrfbdy_d01 .
\7 cp ../runCase.init/wrfinput_d01 .
\8 cp ../runCase.init/wrflowinp_d01 .
\9 matlab -nodisplay < updateLowinp.m
10
11 cp namelist.input.set namelist.input
12 mpirun -np 4 wrf.exe
13 matlab -nodisplay < updateHFlux.m
14
15 cp namelist.input.run namelist.input
16 cp ../coupledCode/esmf_application .
17 echo "running coupled MITgcm--WRF simulation.."
18 mpirun -np 4 esmf_application &> log.esmf

The namelist.rc file controls the coupled run:

Debuglevel: 0
## mode 1 = sequential mode
## mode 2 = concurrent mode
coupleMode: 1
cpuOCN:     4
cpuATM:     4

StartYear:  2012
StartMonth:   06
StartDay:     01
StartHour:    00
StartMinute:  00
StartSecond:  00

EndYear:    2012
EndMonth:     06
EndDay:       01
EndHour:      06
EndMinute:    00
EndSecond:    00

## EsmStepSeconds is the coupling time step
EsmStepSeconds: 60
## ATMStepSeconds is the time step for WRF
ATMStepSeconds: 60
## OCNStepSeconds is the time step for MITgcm
OCNStepSeconds: 60
## The wave model is not activated
WAVStepSeconds: 60

In namelist.input file, I added the following lines to control the input and output from ESMF:

## auxinput5_interval_s sets the coupling interval for ESMF inputs
auxinput5_inname                    = 'wrfin_esmf',
auxinput5_interval_s                = 60,
auxinput5_end_d                     = 60,
io_form_auxinput5                   = 7,
## auxhist5_interval_s sets the coupling interval for ESMF outputs
auxhist5_outname                    = 'wrfout_esmf',
auxhist5_interval_s                 = 60,
auxhist5_end_d                      = 60,
io_form_auxhist5                    = 7,