--- trunk/phylmd/Interface_surf/calcul_fluxs.f	2018/09/18 15:14:40	308
+++ trunk/phylmd/Interface_surf/calcul_fluxs.f	2018/12/03 17:52:21	311
@@ -4,9 +4,9 @@
 
 contains
 
-  SUBROUTINE calcul_fluxs(tsurf, p1lay, cal, beta, cdragh, ps, qsurf, &
-       radsol, t1lay, q1lay, u1lay, v1lay, tAcoef, qAcoef, tBcoef, &
-       qBcoef, tsurf_new, evap, fluxlat, flux_t, dflux_s, dflux_l, dif_grnd)
+  SUBROUTINE calcul_fluxs(tsurf, p1lay, cal, beta, cdragh, ps, qsurf, radsol, &
+       t1lay, q1lay, u1lay, v1lay, tAcoef, qAcoef, tBcoef, qBcoef, tsurf_new, &
+       evap, fluxlat, flux_t, dflux_s, dflux_l, dif_grnd)
 
     ! Cette routine calcule les flux en h et q à l'interface et une
     ! température de surface.
@@ -32,11 +32,11 @@
     real, intent(IN):: cal(:) ! (knon) capacité calorifique du sol
     real, intent(IN):: beta(:) ! (knon) évaporation réelle
     real, intent(IN):: cdragh(:) ! (knon) coefficient d'échange
-    real, intent(IN):: ps(:) ! (knon) pression au sol
+    real, intent(IN):: ps(:) ! (knon) pression au sol, en Pa
     real, intent(OUT):: qsurf(:) ! (knon) humidité de l'air au-dessus du sol
 
     real, intent(IN):: radsol(:) ! (knon)
-    ! rayonnement net au sol (longwave + shortwave)
+    ! net downward radiative (longwave + shortwave) flux at the surface
 
     real, intent(IN):: dif_grnd ! coefficient de diffusion vers le sol profond
     real, intent(IN):: t1lay(:), q1lay(:), u1lay(:), v1lay(:) ! (knon)
@@ -51,7 +51,7 @@
     real, intent(OUT):: evap(:) ! (knon)
 
     real, intent(OUT):: fluxlat(:), flux_t(:) ! (knon)
-    ! flux de chaleurs latente et sensible
+    ! flux de chaleurs latente et sensible, en W m-2
 
     real, intent(OUT):: dflux_s(:), dflux_l(:) ! (knon)
     ! dérivées des flux de chaleurs sensible et latente par rapport à