diff --git a/examples/PMillennium/PMillennium-1536/getIC.sh b/examples/PMillennium/PMillennium-1536/getIC.sh
new file mode 100755
index 0000000000000000000000000000000000000000..0ba78a0fd11025b12a7a9a7be9de47c2b5ad2898
--- /dev/null
+++ b/examples/PMillennium/PMillennium-1536/getIC.sh
@@ -0,0 +1,2 @@
+#!/bin/bash
+wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/PMillennium/PMill-1536.hdf5
diff --git a/examples/PMillennium/PMillennium-1536/p-mill-1536.yml b/examples/PMillennium/PMillennium-1536/p-mill-1536.yml
index 254f2a4397e7b2a6c6ef22dd50f4544d550e992a..ea44572f90eb69330c973946cb5533a4c58b2c82 100644
--- a/examples/PMillennium/PMillennium-1536/p-mill-1536.yml
+++ b/examples/PMillennium/PMillennium-1536/p-mill-1536.yml
@@ -45,7 +45,7 @@ Gravity:
   
 # Parameters related to the initial conditions
 InitialConditions:
-  file_name:                   ics.hdf5
+  file_name:                   PMill-1536.hdf5
   periodic:                    1
   cleanup_h_factors:           1    
   cleanup_velocity_factors:    1  
diff --git a/examples/PMillennium/PMillennium-384/getIC.sh b/examples/PMillennium/PMillennium-384/getIC.sh
new file mode 100755
index 0000000000000000000000000000000000000000..eef95aac218a55e3e432fd1afb59ee7823c930fd
--- /dev/null
+++ b/examples/PMillennium/PMillennium-384/getIC.sh
@@ -0,0 +1,2 @@
+#!/bin/bash
+wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/PMillennium/PMill-384.hdf5
diff --git a/examples/PMillennium/PMillennium-384/p-mill-384.yml b/examples/PMillennium/PMillennium-384/p-mill-384.yml
index 4aede77c0c1a8a6818c95c318364150ede919a01..0040f580c8e67182949875d7d85cee2f75851654 100644
--- a/examples/PMillennium/PMillennium-384/p-mill-384.yml
+++ b/examples/PMillennium/PMillennium-384/p-mill-384.yml
@@ -41,11 +41,11 @@ Gravity:
   theta:                  0.5          
   comoving_softening:     0.08333  # 83.333 kpc = 1/25 mean inter-particle separation
   max_physical_softening: 0.08333  # 83.333 kpc = 1/25 mean inter-particle separation
-  mesh_side_length:       128
+  mesh_side_length:       64
   
 # Parameters related to the initial conditions
 InitialConditions:
-  file_name:                   ics.hdf5
+  file_name:                   PMill-384.hdf5
   periodic:                    1
   cleanup_h_factors:           1    
   cleanup_velocity_factors:    1  
diff --git a/examples/PMillennium/PMillennium-768/getIC.sh b/examples/PMillennium/PMillennium-768/getIC.sh
new file mode 100755
index 0000000000000000000000000000000000000000..6c020a5bb4269d4fd4284c617cfd94ab19e87326
--- /dev/null
+++ b/examples/PMillennium/PMillennium-768/getIC.sh
@@ -0,0 +1,2 @@
+#!/bin/bash
+wget http://virgodb.cosma.dur.ac.uk/swift-webstorage/ICs/PMillennium/PMill-768.hdf5
diff --git a/examples/PMillennium/PMillennium-768/p-mill-768.yml b/examples/PMillennium/PMillennium-768/p-mill-768.yml
index a70c9c70831af9c237a466165b25b6300df69336..5ba97af72513dca5fdccd216f0ede78d0e279b0a 100644
--- a/examples/PMillennium/PMillennium-768/p-mill-768.yml
+++ b/examples/PMillennium/PMillennium-768/p-mill-768.yml
@@ -41,11 +41,11 @@ Gravity:
   theta:                  0.5          
   comoving_softening:     0.041666  # 41.6666 kpc = 1/25 mean inter-particle separation
   max_physical_softening: 0.041666  # 41.6666 kpc = 1/25 mean inter-particle separation
-  mesh_side_length:       256
+  mesh_side_length:       128
   
 # Parameters related to the initial conditions
 InitialConditions:
-  file_name:                   ics.hdf5
+  file_name:                   PMill-768.hdf5
   periodic:                    1
   cleanup_h_factors:           1    
   cleanup_velocity_factors:    1  
diff --git a/examples/PMillennium/README b/examples/PMillennium/README
new file mode 100644
index 0000000000000000000000000000000000000000..e4755bc214d438a282521b77a34ad853d0f36871
--- /dev/null
+++ b/examples/PMillennium/README
@@ -0,0 +1,9 @@
+Initial conditions for the P-Millennium simulation (Baugh et al. 2018)
+at various resolution. These are all DMONLY runs in a 800^3 Mpc^3
+volume with the cosmology given by Planck-13:
+
+ Om = 0.307, Ol = 0.693, h = 0.6777, Ob = 0.0455
+
+The ICs exist at different resolution. The Millennium simulation
+(Springel 2005) has a resolution in between the 1536^3 and 3072^3
+examples given here.