--- trunk/Sources/phylmd/screenc.f	2017/11/15 13:56:45	246
+++ trunk/phylmd/screenc.f	2018/07/11 17:06:09	275
@@ -19,7 +19,7 @@
 
     ! I. Musat, 01.07.2002
 
-    use coefcdrag_m, only: coefcdrag
+    use cdrag_m, only: cdrag
     use SUPHEC_M, only: RG
 
     INTEGER, intent(in):: klon
@@ -46,8 +46,7 @@
     REAL, dimension(klon), intent(in):: qstar
     ! qstar---input-R- facteur d'echelle pour l'humidite relative
 
-    REAL, dimension(klon), intent(out):: pref
-    ! pref----input-R- pression au niveau de reference
+    REAL, intent(out):: pref(:) ! (knon) pression au niveau de reference
     REAL, dimension(klon), intent(out):: delu
     ! delu----input-R- anomalie du vent par rapport au 1er niveau
     REAL, dimension(klon), intent(out):: delte
@@ -58,7 +57,7 @@
 
     ! Local:
     INTEGER i
-    REAL, dimension(klon):: cdram, cdrah, cdran, zri1, gref
+    REAL, dimension(knon):: cdram, cdrah, gref
 
     !------------------------------------------------------------------------- 
 
@@ -68,9 +67,8 @@
 
     ! Richardson at reference level
 
-    CALL coefcdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), &
-         gref(:knon), psol(:knon), ts, qsurf, rugos, cdram, cdrah, cdran, &
-         zri1, pref)
+    CALL cdrag(nsrf, speed(:knon), temp(:knon), q_zref(:knon), gref(:knon), &
+         psol(:knon), ts(:knon), qsurf(:knon), rugos(:knon), cdram, cdrah, pref)
 
     DO i = 1, knon
        delu(i) = ustar(i) / sqrt(cdram(i))