Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlsplt_loc.f90

Timestamp:

Aug 2, 2024, 9:58:25 PM (3 months ago)

Author:

abarral

Message:

Put dimensions.h and paramet.h into modules

File:

• LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlsplt_loc.f90 (modified) (18 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3dmem/vlsplt_loc.f90

@@ -4 +4 @@
 
   ! 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 parallel_lmdz
@@ -16 +16 @@
           min_qParent, min_qMass, min_ratio ! MVals et CRisi
   USE lmdz_iniprint, ONLY: lunout, prt_level
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
-  !
-  !
+
+
+
+
   !   Arguments:
   !   ----------
@@ -29 +31 @@
   REAL :: w(ijb_u:ije_u, llm)
   INTEGER :: iq ! CRisi
-  !
+
   !  Local
   !   ---------
-  !
+
   INTEGER :: ij, l, j, i, iju, ijq, indu(ijnb_u), niju
   INTEGER :: n0, iadvplus(ijb_u:ije_u, llm), nl(llm)
-  !
+
   REAL :: new_m, zu_m, zdum(ijb_u:ije_u, llm)
   REAL :: sigu(ijb_u:ije_u), dxq(ijb_u:ije_u, llm), dxqu(ijb_u:ije_u)
@@ -367 +369 @@
 
 SUBROUTINE vly_loc(q, pente_max, masse, masse_adv_v, iq)
-  !
+
   ! Auteurs:   P.Le Van, F.Hourdin, F.Forget
-  !
+
   !    ********************************************************************
   ! Shema  d'advection " pseudo amont " .
@@ -375 +377 @@
   ! 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 parallel_lmdz
@@ -385 +387 @@
   USE lmdz_comgeom
 
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
-  !
-  !
+
+
+
+
   !   Arguments:
   !   ----------
@@ -397 +401 @@
   REAL :: q(ijb_u:ije_u, llm, nqtot), dq(ijb_u:ije_u, llm)
   INTEGER :: iq ! CRisi
-  !
+
   !  Local
   !   ---------
-  !
+
   INTEGER :: i, ij, l
-  !
+
   REAL :: airej2, airejjm, airescb(iim), airesch(iim)
   REAL :: dyq(ijb_u:ije_u, llm), dyqv(ijb_v:ije_v), zdvm(ijb_u:ije_u, llm)
@@ -456 +460 @@
   ENDIF
 
-  !
+
   ! PRINT*,'CALCUL EN LATITUDE'
 
   !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
   DO l = 1, llm
-    !
+
     !   --------------------------------
     !  CALCUL EN LATITUDE
@@ -596 +600 @@
     ! 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.
@@ -605 +609 @@
     !    &   dyqv(ismin(iim,dyqv(ip1jm-iip1+1),1)+ip1jm-iip1+1).le.0.)
     !    &   apps=0.
-    !
+
     !   limitation des pentes aux poles
     ! DO ij=1,iip1
@@ -611 +615 @@
     !    dyq(ip1jm+ij)=apps*dyq(ip1jm+ij)
     ! ENDDO
-    !
+
     !   test
     !  DO ij=1,iip1
@@ -620 +624 @@
     !     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.)
@@ -632 +636 @@
     !    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.)
@@ -797 +801 @@
 
 SUBROUTINE vlz_loc(q, pente_max, masse, w, ijb_x, ije_x, iq)
-  !
+
   ! Auteurs:   P.Le Van, F.Hourdin, F.Forget
-  !
+
   !    ********************************************************************
   ! Shema  d'advection " pseudo amont " .
@@ -805 +809 @@
   !    q,pbaru,pbarv,w sont des arguments d'entree  pour le s-pg ....
   ! dq            sont des arguments de sortie pour le s-pg ....
-  !
-  !
+
+
   !   --------------------------------------------------------------------
   USE parallel_lmdz
@@ -815 +819 @@
 
 
+USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
   !
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
-  !
-  !
+
+
+
+
   !   Arguments:
   !   ----------
@@ -827 +833 @@
   REAL :: w(ijb_u:ije_u, llm + 1, nqtot)
   INTEGER :: iq
-  !
+
   !  Local
   !   ---------
-  !
+
   INTEGER :: i, ij, l, j, ii
 
@@ -837 +843 @@
   INTEGER, SAVE :: countcfl
   !$OMP THREADPRIVATE(countcfl)
-  !
+
   REAL :: newmasse