diff --git a/srcF/sfx_api_inmcm.f90 b/srcF/sfx_api_inmcm.f90
index c04567e1af0112818cb722ddb0f29d341ba19faa..c1354a6ea39f2d3630ebeca755081349f895b99d 100644
--- a/srcF/sfx_api_inmcm.f90
+++ b/srcF/sfx_api_inmcm.f90
@@ -42,7 +42,7 @@ contains
         meteo%depth = depth_inm
         meteo%lai = lai_inm
         meteo%surface_type = IVEG_sfx
-
+        write(*,*) 'surface_type, IVEG_sfx', meteo%surface_type, IVEG_sfx
     end subroutine inmcm_to_sfx_in_cell
     ! --------------------------------------------------------------------------------
     
diff --git a/srcF/sfx_sheba_noniterative.f90 b/srcF/sfx_sheba_noniterative.f90
index 35e3b7ca50b6bac1b5d51bb2343d01e1f291f9fb..c7912efee0a4025ac5fbc703185c693e93c08abb 100644
--- a/srcF/sfx_sheba_noniterative.f90
+++ b/srcF/sfx_sheba_noniterative.f90
@@ -154,7 +154,7 @@ contains
             h = meteo%h(i), &
             U = meteo%U(i), dT = meteo%dT(i), Tsemi = meteo%Tsemi(i), dQ = meteo%dQ(i), &
             z0_m = meteo%z0_m(i), depth=meteo%depth(i), lai=meteo%lai(i), surface_type=meteo%surface_type(i))
-
+            write(*,*) 'get_flux', meteo%surface_type(i)
             call get_surface_fluxes(sfx_cell, meteo_cell, numerics)
 
             call push_sfx_data(sfx, sfx_cell, i)
@@ -188,6 +188,7 @@ contains
         real :: z0_m    !< surface aerodynamic roughness (should be < 0 for water bodies surface)
         real :: depth   
         real :: lai
+        integer :: surface_type
         ! ----------------------------------------------------------------------------
 
         ! --- local variables
@@ -221,7 +222,7 @@ contains
 
         real Cm, Ct             !< transfer coeff. for (momentum) & (heat) [n/d]
 
-        integer surface_type    !< surface type = (ocean || land)
+        
 
         real fval               !< just a shortcut for partial calculations