Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d/vlsplt.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/vlsplt.F90 (modified) (18 diffs)

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

@@ -1 @@
-!
+
 ! $Id$
 !
@@ -6 +6 @@
   USE infotrac, ONLY: nqtot, tracers
   USE lmdz_ssum_scopy, ONLY: scopy
-  !
+
   ! Auteurs:   P.Le Van, F.Hourdin, F.Forget
-  !
+
   !    ********************************************************************
   ! Shema  d'advection " pseudo amont " .
   !    ********************************************************************
   ! q,pbaru,pbarv,w sont des arguments d'entree  pour le s-pg ....
-  !
+
   !   pente_max facteur de limitation des pentes: 2 en general
   !                                           0 pour un schema amont
   !   pbaru,pbarv,w flux de masse en u ,v ,w
   !   pdt pas de temps
-  !
+
   !   --------------------------------------------------------------------
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
-
-  !
+
+
+
+
   !   Arguments:
   !   ----------
@@ -33 +35 @@
   REAL :: w(ip1jmp1, llm), pdt
   INTEGER :: iq ! CRisi
-  !
+
   !  Local
   !   ---------
-  !
+
   INTEGER :: ij, l
-  !
+
   REAL :: zm(ip1jmp1, llm, nqtot)
   REAL :: mu(ip1jmp1, llm)
@@ -110 +112 @@
 
   ! Auteurs:   P.Le Van, F.Hourdin, F.Forget
-  !
+
   !    ********************************************************************
   ! Shema  d'advection " pseudo amont " .
   !    ********************************************************************
   ! nq,iq,q,pbaru,pbarv,w sont des arguments d'entree  pour le s-pg ....
-  !
-  !
+
+
   !   --------------------------------------------------------------------
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
-  !
-  !
+
+
+
+
   !   Arguments:
   !   ----------
@@ -130 +134 @@
   REAL :: q(ip1jmp1, llm, nqtot)
   INTEGER :: iq ! CRisi
-  !
+
   !  Local
   !   ---------
-  !
+
   INTEGER :: ij, l, j, i, iju, ijq, indu(ip1jmp1), niju
   INTEGER :: n0, iadvplus(ip1jmp1, llm), nl(llm)
-  !
+
   REAL :: new_m, zu_m, zdum(ip1jmp1, llm)
   REAL :: dxq(ip1jmp1, llm), dxqu(ip1jmp1)
@@ -424 +428 @@
   USE lmdz_comgeom
 
- !
+
   ! Auteurs:   P.Le Van, F.Hourdin, F.Forget
-  !
+
   !    ********************************************************************
   ! Shema  d'advection " pseudo amont " .
@@ -432 +436 @@
   ! q,masse_adv_v,w sont des arguments d'entree  pour le s-pg ....
   ! dq            sont des arguments de sortie pour le s-pg ....
-  !
-  !
+
+
   !   --------------------------------------------------------------------
   USE comconst_mod, ONLY: pi
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
-  !
-  !
+
+
+
+
   !   Arguments:
   !   ----------
@@ -448 +454 @@
   REAL :: q(ip1jmp1, llm, nqtot)
   INTEGER :: iq ! CRisi
-  !
+
   !  Local
   !   ---------
-  !
+
   INTEGER :: i, ij, l
-  !
+
   REAL :: airej2, airejjm, airescb(iim), airesch(iim)
   REAL :: dyq(ip1jmp1, llm), dyqv(ip1jm)
@@ -493 +499 @@
   ENDIF
 
-  !
+
   !PRINT*,'CALCUL EN LATITUDE'
 
   DO l = 1, llm
-    !
+
     !   --------------------------------
     !  CALCUL EN LATITUDE
@@ -590 +596 @@
     ! appn=min(pente_max/appn,1.)
     ! apps=min(pente_max/apps,1.)
-    !
-    !
+
+
     !   cas ou on a un extremum au pole
-    !
+
     ! IF(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.)
     !    &   appn=0.
@@ -599 +605 @@
     !    &   dyqv(ismin(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1).le.0.)
     !    &   apps=0.
-    !
+
     !   limitation des pentes aux poles
     ! DO ij=1,iip1
@@ -605 +611 @@
     !    dyq(ip1jm+ij)=apps*dyq(ip1jm+ij)
     ! ENDDO
-    !
+
     !   test
     !  DO ij=1,iip1
@@ -614 +620 @@
     !     dyq(ij)=dyq(ij)*cos(rlatu((ij-1)/iip1+1))
     !  ENDDO
-    !
+
     ! changement 10 07 96
     ! IF(dyqv(ismin(iim,dyqv,1))*dyqv(ismax(iim,dyqv,1)).le.0.)
@@ -626 +632 @@
     !    ENDDO
     ! ENDIF
-    !
+
     ! IF(dyqv(ismax(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1)*
     !    & dyqv(ismin(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1).le.0.)
@@ -745 +751 @@
   USE infotrac, ONLY: nqtot, tracers, & ! CRisi
           min_qParent, min_qMass, min_ratio ! MVals et CRisi
-  !
+
   ! Auteurs:   P.Le Van, F.Hourdin, F.Forget
-  !
+
   !    ********************************************************************
   ! Shema  d'advection " pseudo amont " .
@@ -754 +760 @@
   ! dq            sont des arguments de sortie pour le s-pg ....
   !   --------------------------------------------------------------------
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
-  !
-  !
+
+
+
+
   !   Arguments:
   !   ----------
@@ -766 +774 @@
   REAL :: w(ip1jmp1, llm + 1)
   INTEGER :: iq
-  !
+
   !  Local
   !   ---------
-  !
+
   INTEGER :: ij, l
-  !
+
   REAL :: wq(ip1jmp1, llm + 1), newmasse
 
@@ -886 +894 @@
 
 SUBROUTINE minmaxq(zq, qmin, qmax, comment)
-
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
+  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
 
   CHARACTER(LEN = 20) :: comment