Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j.F90

Timestamp:

Aug 2, 2024, 9:58:25 PM (7 weeks ago)

Author:

abarral

Message:

Put dimensions.h and paramet.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j.F90 (modified) (8 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d/tetaleveli1j.F90

@@ -6 +6 @@
 
   ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
 
-  INCLUDE "dimensions.h"
+
 
   !================================================================
-  !
+
   ! Interpoler des champs 3-D u, v et g du modele a un niveau de
   ! pression donnee (pres)
-  !
+
   ! INPUT:  ilon ----- nombre de points
   !     ilev ----- nombre de couches
@@ -22 +24 @@
   !     Qgcm ----- champ GCM
   !     Qpres ---- champ interpolle au niveau pres
-  !
+
   !================================================================
-  !
+
   !   arguments :
   !   -----------
@@ -43 +45 @@
   ! REAL ptop, pbot, aist(klon), aisb(klon)
   !
-  INCLUDE "paramet.h"
-  !
+
+
   INTEGER :: lt(ip1jm), lb(ip1jm)
   REAL :: ptop, pbot, aist(ip1jm), aisb(ip1jm)
@@ -51 +53 @@
 
   INTEGER :: i, k
-  !
+
   ! PRINT*,'tetalevel pres=',pres
   !=====================================================================
@@ -60 +62 @@
 
     ! Chercher les 2 couches les plus proches du niveau a obtenir
-    !
+
     ! Eventuellement, faire l'extrapolation a partir des deux couches
     ! les plus basses ou les deux couches les plus hautes:
@@ -87 +89 @@
       END DO
     END DO
-    !
+
     ! Interpolation lineaire:
-    !
+
     DO i = 1, ilon
       ! interpolation en logarithme de pression:
-      !
+
       ! ...   Modif . P. Le Van    ( 20/01/98) ....
       !   Modif Frédéric Hourdin (3/01/02)
@@ -98 +100 @@
       IF(pgcm(i, lb(i))==0.OR. &
               pgcm(i, lt(i))==0.) THEN
-        !
+
         PRINT*, 'i,lb,lt,2pgcm,pres', i, lb(i), &
                 lt(i), pgcm(i, lb(i)), pgcm(i, lt(i)), pres
-        !
+
       ENDIF
-      !
+
       aist(i) = LOG(pgcm(i, lb(i)) / pres) &
               / LOG(pgcm(i, lb(i)) / pgcm(i, lt(i)))
@@ -121 +123 @@
     !IM  $   Qgcm(i,lt(i)),aist(i),Qpres(i)
   enddo
-  !
+
   ! Je mets les vents a zero quand je rencontre une montagne
   DO i = 1, ilon