Aquaplanet-slab-ocean.md

@@ -120,10 +120,14 @@ Not yet there, but getting close...... Namelist parameters and input forcing fil
 **Modifications to namelist settings** 
 1. Under **runctl**, please add \
    `echo   "    l_orbvsop87 = .false."                           >>  namelist.echam` \
-   `echo   "    lmlo = .true."                                   >>  namelist.echam` 
+   `echo   "    lmlo = .true."                                   >>  namelist.echam`
+   `echo   "    lmlo_ice = .false."                                   >>  namelist.echam` 
+
+ 
 2. Under **radctl**, please add \
    `echo   "    ldiur = false"                                   >>  namelist.echam` \
    `echo   "    cecc = 0.0"                                      >>  namelist.echam` \
+   `echo   "    cobld = 0.0"                                     >>  namelist.echam` \
    Please also comment out as follows:
    ```bash
    #if [[ ${read_ghg} == false ]]; then
@@ -170,7 +174,7 @@ Not yet there, but getting close...... Namelist parameters and input forcing fil
 
 **Modifications to linking input files**
 
-While many input forcing files need to change, this is done by changing the link/pointing. However, to run with slab ocean, you'll need an extra file, the `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 `T63_heatflux_test.nc`. There are also two other files that are not necessarily readily available, `T63GR15_surf_aqua_irad_noseason_final.nc` and `T63_ozone_historical_1850.nc` \
+While many input forcing files need to change, this is done by changing the link/pointing. However, to run with slab ocean, you'll need an extra file, the `heatflux` or `qflx` file. Assume for the moment you have this file `T63_heatflux_test.nc`. There are also two other files that are not necessarily readily available, `T63GR15_surf_aqua_irad_noseason_final.nc` and `T63_ozone_historical_1850.nc` \
 Below, unless otherwise stated, `#` means to comment these lines out and be sure to add the new pointers. 
 
 1. 3d initial file of atmosphere (spectral, no dependence on lsmask) 
@@ -249,6 +253,12 @@ Below, unless otherwise stated, `#` means to comment these lines out and be sure
      get_file ${srfmod} input ${ludir}/cover_fract_${res_srf}_${ntiles}tiles_2015.nc  cover_fract.${startyear}.nc
    ```
 
+A little more on preparations of these input files, particularly for aquaplanet slab ocean.
+
+- The 'qlfx' file is normally obtained by first running a slab ocean aquaplanet with zero flux correction in order to get the climatological SST. The SST is then used to run a pure aquaplanet run with prescribed SST so that we can then compute the surface flux. This surface flux then requires post-processing to get a climatological (12 months) zonally invariant and hemispherically symmetric flux file, that is then used as an input to the slab ocean part.
+- Ozone file needs a vertical profile to ensure radiative properties of the stratosphere and troposphere, but for an aquaplanet slab ocean, we would need zonally invariant and hemispherically symmetric ozone concentrations. 
+- Solar irradiance must have no seasons. Here, we ensure eccentricity and obliquity is also set to zero. However, all these do not account for varying zenith angle of the sun. 
+
 # Submit job
 ```bash
 cd experiments/aqua0101/scripts