diff --git a/srcF/sfx_run.f90 b/srcF/sfx_run.f90
index 6304f7d572f07efdadf26899ffdd9cbd896ce52a..784ecb839319cd737488052a3642beab58eadcf0 100644
--- a/srcF/sfx_run.f90
+++ b/srcF/sfx_run.f90
@@ -174,6 +174,8 @@ contains
             call get_surface_fluxes_vec_sheba_noit(sfx, meteo, numerics_sheba_noit, num)
         else if (model == model_most_snow) then
             call get_surface_fluxes_vec_most_snow(sfx, meteo, numerics_most_snow, num)
+        else if (model == model_sheba_coare) then
+            call get_surface_fluxes_vec_sheba_coare(sfx, meteo, numerics_sheba_coare, num)
         end if
     
     
@@ -260,7 +262,7 @@ contains
                 write(*, *) ' --help'
                 write(*, *) '    print usage options'
                 write(*, *) ' --model [key]'
-                write(*, *) '    key = esm (default) || log || most || sheba || sheba_noit || most_snow'
+                write(*, *) '    key = esm (default) || log || most || sheba || sheba_noit || most_snow || sheba_coare'
                 write(*, *) ' --dataset [key]'
                 write(*, *) '    key = mosaic (default) || irgason || sheba'
                 write(*, *) '        = lake || papa || toga || user [filename] [param]'