--- trunk/Sources/phylmd/stdlevvar.f	2016/03/22 16:31:39	188
+++ trunk/Sources/phylmd/stdlevvar.f	2017/10/16 12:35:41	225
@@ -9,58 +9,40 @@
 
     ! From LMDZ4/libf/phylmd/stdlevvar.F90, version 1.3 2005/05/25 13:10:09
 
-    use coefcdrag_m, only: coefcdrag
-    USE suphec_m, ONLY: rg, rkappa
-    use screenp_m, only: screenp
-
     ! Objet : calcul de la température et de l'humidité relative à 2 m
     ! et du module du vent à 10 m à partir des relations de
     ! Dyer-Businger et des équations de Louis.
 
     ! Reference: Hess, Colman and McAvaney (1995) 
 
-    ! Author: I. Musat, 01.07.2002
+    ! Author: I. Musat, July 1st, 2002
+
+    use coefcdrag_m, only: coefcdrag
+    USE suphec_m, ONLY: rg, rkappa
+    use screenc_m, only: screenc
+    use screenp_m, only: screenp
 
     INTEGER, intent(in):: klon
     ! dimension de la grille physique (= nb_pts_latitude X nb_pts_longitude)
 
-    INTEGER, intent(in):: knon
-    ! knon----input-I- nombre de points pour un type de surface
-    INTEGER, intent(in):: nsrf
-    ! nsrf----input-I- indice pour le type de surface; voir indicesol.inc
-    LOGICAL, intent(in):: zxli
-    ! zxli----input-L- TRUE si calcul des cdrags selon Laurent Li
-    REAL, dimension(klon), intent(in):: u1
-    ! u1------input-R- vent zonal au 1er niveau du modele
-    REAL, dimension(klon), intent(in):: v1
-    ! v1------input-R- vent meridien au 1er niveau du modele
-    REAL, dimension(klon), intent(in):: t1
-    ! t1------input-R- temperature de l'air au 1er niveau du modele
-    REAL, dimension(klon), intent(in):: q1
-    ! q1------input-R- humidite relative au 1er niveau du modele
-    REAL, dimension(klon), intent(in):: z1
-    ! z1------input-R- geopotentiel au 1er niveau du modele
-    REAL, dimension(klon), intent(in):: ts1
-    ! ts1-----input-R- temperature de l'air a la surface
-    REAL, dimension(klon), intent(in):: qsurf
-     ! qsurf---input-R- humidite relative a la surface
-    REAL, dimension(klon), intent(in):: rugos
-    ! rugos---input-R- rugosite
-   REAL, dimension(klon), intent(in):: psol
-    ! psol----input-R- pression au sol
-   REAL, dimension(klon), intent(in):: pat1
-    ! pat1----input-R- pression au 1er niveau du modele
-
-    REAL, dimension(klon), intent(out):: t_2m
-    ! t_2m---output-R- temperature de l'air a 2m
-    REAL, dimension(klon), intent(out):: q_2m
-    ! q_2m---output-R- humidite relative a 2m
-    REAL, dimension(klon), intent(out):: t_10m
-    ! t_10m--output-R- temperature de l'air a 10m
-    REAL, dimension(klon), intent(out):: q_10m
-    ! q_10m--output-R- humidite specifique a 10m
-    REAL, dimension(klon), intent(out):: u_10m
-    ! u_10m--output-R- vitesse du vent a 10m
+    INTEGER, intent(in):: knon ! nombre de points pour un type de surface
+    INTEGER, intent(in):: nsrf ! indice pour le type de surface
+    LOGICAL, intent(in):: zxli ! calcul des cdrags selon Laurent Li
+    REAL, intent(in):: u1(:) ! (knon) vent zonal au 1er niveau du modele
+    REAL, intent(in):: v1(:) ! (knon) vent meridien au 1er niveau du modele
+    REAL, intent(in):: t1 (klon) ! temperature de l'air au 1er niveau du modele
+    REAL, intent(in):: q1(klon) ! humidite relative au 1er niveau du modele
+    REAL, intent(in):: z1 (klon) ! geopotentiel au 1er niveau du modele
+    REAL, intent(in):: ts1(klon) ! temperature de l'air a la surface
+    REAL, intent(in):: qsurf(klon) ! humidite relative a la surface
+    REAL, intent(in):: rugos(klon) ! rugosite
+    REAL, intent(in):: psol(klon) ! pression au sol
+    REAL, intent(in):: pat1(klon) ! pression au 1er niveau du modele
+    REAL, intent(out):: t_2m(klon) ! temperature de l'air a 2m
+    REAL, intent(out):: q_2m(klon) ! humidite relative a 2m
+    REAL, intent(out):: t_10m(klon) ! temperature de l'air a 10m
+    REAL, intent(out):: q_10m(klon) ! humidite specifique a 10m
+    REAL, intent(out):: u_10m(klon) ! vitesse du vent a 10m
     REAL, intent(out):: ustar(klon) ! u*
 
     ! Local:
@@ -175,10 +157,8 @@
     ! First aproximation of variables at zref  
 
     zref = 10.0
-    CALL screenp(klon, knon, speed, tpot, q1, &
-         ts1, qsurf, rugos, lmon, &
-         ustar, testar, qstar, zref, &
-         delu, delte, delq)
+    CALL screenp(klon, knon, speed, tpot, q1, ts1, qsurf, rugos, lmon, ustar, &
+         testar, qstar, zref, delu, delte, delq)
 
     DO i = 1, knon
        u_zref(i) = delu(i)
@@ -193,11 +173,9 @@
 
     DO n = 1, niter
        okri=.TRUE.
-       CALL screenc(klon, knon, nsrf, zxli, &
-            u_zref, temp, q_zref, zref, &
-            ts1, qsurf, rugos, psol, &
-            ustar, testar, qstar, okri, ri1, &
-            pref, delu, delte, delq)
+       CALL screenc(klon, knon, nsrf, zxli, u_zref, temp, q_zref, zref, ts1, &
+            qsurf, rugos, psol, ustar, testar, qstar, okri, ri1, pref, delu, &
+            delte, delq)
 
        DO i = 1, knon
           u_zref(i) = delu(i)
@@ -209,9 +187,7 @@
 
     DO i = 1, knon
        u_zref_c(i) = u_zref(i)
-
        u_10m(i) = u_zref_p(i) * ok_pred(i) + u_zref_c(i) * ok_corr(i)
-
        q_zref_c(i) = q_zref(i)
        temp_c(i) = temp(i)
        t_10m(i) = temp_p(i) * ok_pred(i) + temp_c(i) * ok_corr(i)