Changeset 5368 for LMDZ6/trunk/libf/phylmd/dyn1d/m_1dconv_mod_h.f90

Timestamp:

Dec 3, 2024, 12:40:26 PM (37 hours ago)

Author:

abarral

Message:

(WIP) Turn implicit into explicit declarations
Turn 1dconv.h into a module to that end

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/m_1dconv_mod_h.f90 (moved) (moved from LMDZ6/trunk/libf/phylmd/dyn1d/1Dconv.h) (11 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/m_1dconv_mod_h.f90

@@ -1 @@
-!
-! $Id$
-!
+MODULE m_1dconv_mod_h
+  IMPLICIT NONE; PRIVATE
+  PUBLIC get_uvd, copie, get_uvd2
+
+  REAL play(100)  !pression en Pa au milieu de chaque couche GCM
+  INTEGER JM(100) !pression en Pa au milieu de chaque couche GCM
+  REAL coef1(100) !coefficient d interpolation
+  REAL coef2(100) !coefficient d interpolation
+  INTEGER klev
+
+  INTEGER nblvlm !nombre de niveau de pression du mesoNH
+  REAL playm(100)  !pression en Pa au milieu de chaque couche Meso-NH
+  REAL hplaym(100) !pression en hPa milieux des couches Meso-NH
+
+CONTAINS
+
         subroutine get_uvd(itap,dtime,file_forctl,file_fordat,                  &
      &       ht,hq,hw,hu,hv,hthturb,hqturb,                                     &
@@ -13 +26 @@
 ! pouvoir calculer la convergence et le cisaillement dans la physiq
 !cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-
-
-      INTEGER klev
-      REAL play(100)  !pression en Pa au milieu de chaque couche GCM
-      INTEGER JM(100) !pression en Pa au milieu de chaque couche GCM
-      REAL coef1(100) !coefficient d interpolation
-      REAL coef2(100) !coefficient d interpolation
-
-      INTEGER nblvlm !nombre de niveau de pression du mesoNH
-      REAL playm(100)  !pression en Pa au milieu de chaque couche Meso-NH
-      REAL hplaym(100) !pression en hPa milieux des couches Meso-NH
-
       integer i,j,k,ll,in
 
       CHARACTER*80 file_forctl,file_fordat
-
-      COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev
-      COMMON/com2_phys_gcss/playm,hplaym,nblvlm
 
 !======================================================================
@@ -66 +63 @@
         integer itap
         real dtime
-        real ht(100)
-        real hq(100)
-        real hu(100)
-        real hv(100)
-        real hw(100)
-        real hthturb(100)
-        real hqturb(100)
+        real ht(:)
+        real hq(:)
+        real hu(:)
+        real hv(:)
+        real hw(:)
+        real hthturb(:)
+        real hqturb(:)
         real Ts, Ts_subr
         logical imp_fcg
@@ -141 +138 @@
 ! variables destinees a la lecture du pas de temps du fichier de donnees
 !---------------------------------------------------------------------
-       character*80 aaa,atemps,spaces,apasmax
+       character*80 aaa,atemps,apasmax
        integer nch,imn,ipa
-!---------------------------------------------------------------------
-!  procedures appelees
-        external rdgrads    !lire en iterant dans forcing.dat
 !---------------------------------------------------------------------
                print*,'le pas itap est:',itap
@@ -496 +490 @@
 
 !cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-! cette routine remplit les COMMON com1_phys_gcss et com2_phys_gcss.h
+! cette routine remplit les variables du module
 !ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-      INTEGER klev !nombre de niveau de pression du GCM
-      REAL play(100)  !pression en Pa au milieu de chaque couche GCM
-      INTEGER JM(100)
-      REAL coef1(100)   !coefficient d interpolation
-      REAL coef2(100)   !coefficient d interpolation
-
-      INTEGER nblvlm !nombre de niveau de pression du mesoNH
-      REAL playm(100) !pression en Pa au milieu de chaque couche Meso-NH
-      REAL hplaym(100)!pression en hecto-Pa des milieux de couche Meso-NH
-
-      COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev
-      COMMON/com2_phys_gcss/playm,hplaym,nblvlm
 
       integer k,klevgcm
@@ -578 +559 @@
 !ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
 !
-      INTEGER nblvlm !nombre de niveau de pression du mesoNH
-      REAL playm(100)  !pression en Pa milieu de chaque couche Meso-NH
-      REAL hplaym(100) !pression en hPa des milieux de couche Meso-NH
-      COMMON/com2_phys_gcss/playm,hplaym,nblvlm
 
       INTEGER i,lu,mlz,mlzh
@@ -588 +565 @@
 
       character*4 a
-      character*80 aaa,anblvl,spaces
+      character*80 aaa,anblvl
       integer nch
 
@@ -718 +695 @@
 ! meme.
 !ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc
-
-      INTEGER klev    !nombre de niveau de pression du GCM
-      REAL play(100)  !pression en Pa au milieu de chaque couche GCM
-      INTEGER JM(100)
-      REAL coef1(100) !coefficient d interpolation
-      REAL coef2(100) !coefficient d interpolation
-
-      INTEGER nblvlm !nombre de niveau de pression du mesoNH
-      REAL playm(100) !pression en Pa milieu de chaque couche Meso-NH
-      REAL hplaym(100)!pression en hPa des milieux de couche Meso-NH
-
-      COMMON/com1_phys_gcss/play,coef1,coef2,JM,klev
-      COMMON/com2_phys_gcss/playm,hplaym,nblvlm
 
       REAL psol
@@ -823 +787 @@
       END
 
-      CHARACTER*(*) FUNCTION SPACES(STR,NSPACE)
+      CHARACTER*(80) FUNCTION SPACES(STR,NSPACE)
 !
 ! CERN PROGLIB# M433    SPACES          .VERSION KERNFOR  4.14  860211
@@ -832 +796 @@
 !-    replaced by NSPACE blanks inside the string STR
 !
-      CHARACTER*(*), INTENT(OUT) :: str
-      INTEGER, INTENT(INOUT) :: nspace
+      CHARACTER*(80), INTENT(OUT) :: str
+      INTEGER :: nspace
       INTEGER :: iblank, inonbl, ispace, lenspa, i, lens
 !
@@ -881 +845 @@
 !
   999 END
+
+END MODULE m_1dconv_mod_h