--- trunk/libf/dyn3d/tau2alpha.f90	2011/01/25 15:11:05	39
+++ trunk/dyn3d/tau2alpha.f	2014/07/15 13:43:24	102
@@ -1,23 +1,29 @@
 module tau2alpha_m
 
+    USE paramet_m, ONLY : iip1, jjp1
+    USE dimens_m, ONLY : jjm
+
     IMPLICIT NONE
 
-    REAL lat_min_guide, lat_max_guide
+    private iip1, jjp1, jjm
+
+    REAL dxdys(iip1, jjp1), dxdyu(iip1, jjp1), dxdyv(iip1, jjm)
 
 contains
 
   SUBROUTINE tau2alpha(type, pim, pjm, factt, taumin, taumax, alpha)
 
-    USE dimens_m, ONLY : iim, jjm
-    USE paramet_m, ONLY : iip1, jjp1
-    USE nr_util, ONLY : pi
     USE comgeom, ONLY : cu_2d, cv_2d, rlatu, rlatv
+    use conf_guide_m, only: lat_min_guide, lat_max_guide
+    use dump2d_m, only: dump2d
+    USE dimens_m, ONLY : iim
+    USE nr_util, ONLY : pi
     USE serre, ONLY : clat, clon, grossismx, grossismy
 
     !   arguments :
     INTEGER type
     INTEGER pim, pjm
-    REAL factt, taumin, taumax
+    REAL, intent(in):: factt, taumin, taumax
     REAL dxdy_, alpha(pim, pjm)
     REAL dxdy_min, dxdy_max
 
@@ -31,10 +37,9 @@
     DATA first/ .TRUE./
 
     REAL zdx(iip1, jjp1), zdy(iip1, jjp1)
-
     REAL zlat
-    REAL dxdys(iip1, jjp1), dxdyu(iip1, jjp1), dxdyv(iip1, jjm)
-    COMMON /comdxdy/dxdys, dxdyu, dxdyv
+
+    !------------------------------------------------------------
 
     IF (first) THEN
        DO j = 2, jjm
@@ -67,7 +72,7 @@
        END DO
        DO j = 1, jjm
           DO i = 1, iip1
-             dxdyv(i, j) = 0.5*(dxdys(i, j)+dxdys(i+1, j))
+             dxdyv(i, j) = 0.5*(dxdys(i, j)+dxdys(i, j + 1))
           END DO
        END DO