Aquaplanet-slab-ocean.md

@@ -114,6 +114,41 @@ cd util/running/tools
 Create_TASKS.frm -F setup_mpiesm-as_aqua --id aqua0101 mpiesm-as
 ```
 
+### Edits to run script
+Not yet there, but getting close...... Namelist parameters and input forcing files must be specified in accordance to aquaplanet slab ocean requirements. 
+Namelist settings you'll need to add:
+1. Under **runctl** \
+   `echo   "    l_orbvsop87 = .false."                           >>  namelist.echam` \
+   `echo   "    lmlo = .true."                                   >>  namelist.echam` 
+2. Under **radctl** \
+   `echo   "    ldiur = false"                                   >>  namelist.echam` \
+   `echo   "    cecc = 0.0"                                      >>  namelist.echam` \
+   Please also comment out as follows:
+   ```bash
+   #if [[ ${read_ghg} == false ]]; then
+     echo "    ighg = 0"                                        >>  namelist.echam
+     echo "    icfc = 2"                                        >>  namelist.echam
+     echo "    ich4 = 3"                                        >>  namelist.echam
+   #else
+   #  echo "    ighg = 1"                                        >>  namelist.echam
+   #  echo "    icfc = 4"                                        >>  namelist.echam
+   #  echo "    ich4 = 4"                                        >>  namelist.echam
+   #fi
+
+   #if [ "${piControl}" = "true" ]; then
+   #  echo "    co2vmr = 284.3169860840e-06"                     >>  namelist.echam
+     echo "    co2vmr = 280.0e-06"                              >>  namelist.echam
+     echo "    ch4vmr = 808.2490234375e-09"                     >>  namelist.echam
+     echo "    n2ovmr = 273.0210571289e-09"                     >>  namelist.echam
+   #  echo "    yr_perp = 1850"                                  >>  namelist.echam
+   #fi
+   ```
+
+
+In fact, to run with slab ocean, you'll need a `heatflux` or `qflx` file. This is normally obtained by running a pure aquaplanet run first, and post-processing the output to get a climatological (maybe also zonally symmetric) flux file, that is then used as an input to the slab ocean part. Assume for the moment you have this file and its called `T63_heatflux_test.nc`
+
+
+