--- trunk/libf/dyn3d/interpost.f	2008/02/27 13:16:39	3
+++ trunk/dyn3d/interpost.f90	2014/03/05 14:38:41	81
@@ -1,45 +1,47 @@
-!
-! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/interpost.F,v 1.1.1.1 2004/05/19 12:53:06 lmdzadmin Exp $
-!
-        subroutine interpost(q,qppm)
-
-       use dimens_m
-      use paramet_m
-      use comconst
-      use comvert
-      use comgeom
-       implicit none
 
+! $Header: /home/cvsroot/LMDZ4/libf/dyn3d/interpost.F,v 1.1.1.1 2004/05/19
+! 12:53:06 lmdzadmin Exp $
 
+SUBROUTINE interpost(q, qppm)
 
-c Arguments   
-      real   q(iip1,jjp1,llm)
-      real   qppm(iim,jjp1,llm)
-c Local
-      integer l,i,j
-  
-c RE-INVERSION DES NIVEAUX
-c le programme ppm3d travaille avec une 3ème coordonnée inversée par rapport
-c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
-c On passe donc des niveaux de Lin à ceux du LMDZ
-           
-        do l=1,llm
-          do j=1,jjp1
-             do i=1,iim
-                 q(i,j,l)=qppm(i,j,llm-l+1)
-             enddo
-          enddo
-         enddo
-            
-c BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D
-
-         do l=1,llm
-           do j=1,jjp1
-            q(iip1,j,l)=q(1,j,l)
-           enddo
-         enddo
-  
-      
-       return
+  USE dimens_m
+  USE paramet_m
+  USE comconst
+  USE disvert_m
+  USE comgeom
+  IMPLICIT NONE
 
-       end
+
+
+  ! Arguments
+  REAL q(iip1, jjp1, llm)
+  REAL qppm(iim, jjp1, llm)
+  ! Local
+  INTEGER l, i, j
+
+  ! RE-INVERSION DES NIVEAUX
+  ! le programme ppm3d travaille avec une 3ème coordonnée inversée par
+  ! rapport
+  ! de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
+  ! On passe donc des niveaux de Lin à ceux du LMDZ
+
+  DO l = 1, llm
+    DO j = 1, jjp1
+      DO i = 1, iim
+        q(i, j, l) = qppm(i, j, llm-l+1)
+      END DO
+    END DO
+  END DO
+
+  ! BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D
+
+  DO l = 1, llm
+    DO j = 1, jjp1
+      q(iip1, j, l) = q(1, j, l)
+    END DO
+  END DO
+
+
+  RETURN
+
+END SUBROUTINE interpost