--- trunk/phylmd/indicesol.f90	2013/11/15 18:45:49	76
+++ trunk/phylmd/indicesol.f	2018/08/02 17:23:07	301
@@ -3,11 +3,19 @@
   implicit none
 
   INTEGER, parameter:: nbsrf=4
-  INTEGER, parameter:: is_oce=3 ! ocean
-  INTEGER, parameter:: is_sic = 4 ! sea ice
+
   INTEGER, parameter:: is_ter = 1 ! land
+  ! Note: the first surface type must be land. pbl_surface must call
+  ! the land interface first because the ocean needs the run-off.
+
   INTEGER, parameter:: is_lic = 2 ! land ice
-  REAL, parameter:: epsfra = 1E-05
+
+  INTEGER, parameter:: is_oce = 3 ! ocean
+  ! Ocean must be before sea ice.
+  
+  INTEGER, parameter:: is_sic = 4 ! sea ice
+
+  REAL, parameter:: epsfra = 1E-5
   CHARACTER(len=*), parameter:: clnsurf(nbsrf) = (/'ter', 'lic', 'oce', 'sic'/)
 
 end module indicesol