--- trunk/dyn3d/groupeun.f	2013/11/15 18:45:49	76
+++ trunk/dyn3d/groupeun.f90	2014/03/05 14:38:41	81
@@ -1,59 +1,60 @@
-!
-! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/groupeun.F,v 1.1.1.1 2004/05/19 12:53:07 lmdzadmin Exp $
-!
-      subroutine groupeun(jjmax,llmax,q)
-      use dimens_m
-      use paramet_m
-      use comconst
-      use comgeom
-      implicit none
-
-
-      integer jjmax,llmax
-      real q(iip1,jjmax,llmax)
-
-      integer ngroup
-      parameter (ngroup=3)
-
-      real airen,airecn,qn
-      real aires,airecs,qs
-
-      integer i,j,l,ig,j1,j2,i0,jd
-
-Champs 3D
-      jd=jjp1-jjmax
-      do l=1,llm
-      j1=1+jd
-      j2=2
-      do ig=1,ngroup
-         do j=j1-jd,j2-jd
-            do i0=1,iim,2**(ngroup-ig+1)
-               airen=0.
-               airecn=0.
-               qn=0.
-               aires=0.
-               airecs=0.
-               qs=0.
-               do i=i0,i0+2**(ngroup-ig+1)-1
-                  airen=airen+aire_2d(i,j)
-                  aires=aires+aire_2d(i,jjp1-j+1)
-                  qn=qn+q(i,j,l)
-                  qs=qs+q(i,jjp1-j+1-jd,l)
-               enddo
-               airecn=0.
-               airecs=0.
-               do i=i0,i0+2**(ngroup-ig+1)-1
-                  q(i,j,l)=qn*aire_2d(i,j)/airen
-                  q(i,jjp1-j+1-jd,l)=qs*aire_2d(i,jjp1-j+1)/aires
-               enddo
-            enddo
-            q(iip1,j,l)=q(1,j,l)
-            q(iip1,jjp1-j+1-jd,l)=q(1,jjp1-j+1-jd,l)
-         enddo
-         j1=j2+1
-         j2=j2+2**ig
-      enddo
-      enddo
 
-      return
-      end
+! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/groupeun.F,v 1.1.1.1 2004/05/19
+! 12:53:07 lmdzadmin Exp $
+
+SUBROUTINE groupeun(jjmax, llmax, q)
+  USE dimens_m
+  USE paramet_m
+  USE comconst
+  USE comgeom
+  IMPLICIT NONE
+
+
+  INTEGER jjmax, llmax
+  REAL q(iip1, jjmax, llmax)
+
+  INTEGER ngroup
+  PARAMETER (ngroup=3)
+
+  REAL airen, airecn, qn
+  REAL aires, airecs, qs
+
+  INTEGER i, j, l, ig, j1, j2, i0, jd
+
+  ! hamps 3D
+  jd = jjp1 - jjmax
+  DO l = 1, llm
+    j1 = 1 + jd
+    j2 = 2
+    DO ig = 1, ngroup
+      DO j = j1 - jd, j2 - jd
+        DO i0 = 1, iim, 2**(ngroup-ig+1)
+          airen = 0.
+          airecn = 0.
+          qn = 0.
+          aires = 0.
+          airecs = 0.
+          qs = 0.
+          DO i = i0, i0 + 2**(ngroup-ig+1) - 1
+            airen = airen + aire_2d(i, j)
+            aires = aires + aire_2d(i, jjp1-j+1)
+            qn = qn + q(i, j, l)
+            qs = qs + q(i, jjp1-j+1-jd, l)
+          END DO
+          airecn = 0.
+          airecs = 0.
+          DO i = i0, i0 + 2**(ngroup-ig+1) - 1
+            q(i, j, l) = qn*aire_2d(i, j)/airen
+            q(i, jjp1-j+1-jd, l) = qs*aire_2d(i, jjp1-j+1)/aires
+          END DO
+        END DO
+        q(iip1, j, l) = q(1, j, l)
+        q(iip1, jjp1-j+1-jd, l) = q(1, jjp1-j+1-jd, l)
+      END DO
+      j1 = j2 + 1
+      j2 = j2 + 2**ig
+    END DO
+  END DO
+
+  RETURN
+END SUBROUTINE groupeun