Changeset 5159 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initfluxsto.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/dyn3d_common/initfluxsto.f90 (modified) (10 diffs)

LMDZ6/branches/Amaury_dev/libf/dyn3d_common/initfluxsto.f90

@@ -13 +13 @@
   USE lmdz_comgeom
 
+  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
+  USE lmdz_paramet
   IMPLICIT NONE
 
-  !
+
   !   Routine d'initialisation des ecritures des fichiers histoires LMDZ
   !   au format IOIPSL
-  !
+
   !   Appels succesifs des routines: histbeg
   !                              histhori
@@ -24 +26 @@
   !                              histdef
   !                              histend
-  !
+
   !   Entree:
-  !
+
   !  infile: nom du fichier histoire a creer
   !  day0,anne0: date de reference
@@ -32 +34 @@
   !  t_ops: frequence de l'operation pour IOIPSL
   !  t_wrt: frequence d'ecriture sur le fichier
-  !
+
   !   Sortie:
   !  fileid: ID du fichier netcdf cree
   !  filevid:ID du fichier netcdf pour la grille v
-  !
+
   !   L. Fairhead, LMD, 03/99
-  !
+
   ! =====================================================================
-  !
+
   !   Declarations
-  INCLUDE "dimensions.h"
-  INCLUDE "paramet.h"
+
+
 
   !   Arguments
-  !
+
   CHARACTER(LEN = *) :: infile
   REAL :: tstep, t_ops, t_wrt
@@ -53 +55 @@
   ! This routine needs IOIPSL to work
   !   Variables locales
-  !
+
   REAL :: nivd(1)
   INTEGER :: tau0
@@ -65 +67 @@
   INTEGER :: zan, idayref
   LOGICAL :: ok_sync
-  !
+
   !  Initialisations
-  !
+
   pi = 4. * atan (1.)
   str = 'q  '
   ctrac = 'traceur   '
   ok_sync = .TRUE.
-  !
+
   !  Appel a histbeg: creation du fichier netcdf et initialisations diverses
   !
@@ -91 +93 @@
           1, iip1, 1, jjp1, &
           tau0, zjulian, tstep, uhoriid, fileid)
-  !
+
   !  Creation du fichier histoire pour la grille en V (oblige pour l'instant,
   !  IOIPSL ne permet pas de grilles avec des nombres de point differents dans
@@ -112 +114 @@
           tau0, zjulian, tstep, dhoriid, filedid)
 
-  !
+
   !  Appel a histhori pour rajouter les autres grilles horizontales
-  !
+
   DO jj = 1, jjp1
     DO ii = 1, iip1
@@ -125 +127 @@
           'Grille points scalaires', thoriid)
 
-  !
+
   !  Appel a histvert pour la grille verticale
-  !
+
   CALL histvert(fileid, 'sig_s', 'Niveaux sigma', &
           'sigma_level', &
@@ -141 +143 @@
           1, nivd, dvertiid)
 
-  !
+
   !  Appels a histdef pour la definition des variables a sauvegarder
 
@@ -165 +167 @@
 
 
-  !
+
   ! Masse
-  !
+
   CALL histdef(fileid, 'masse', 'Masse', 'kg', &
           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
           32, 'inst(X)', t_ops, t_wrt)
-  !
+
   !  Pbaru
-  !
+
   CALL histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s', &
           iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, &
           32, 'inst(X)', t_ops, t_wrt)
 
-  !
+
   !  Pbarv
-  !
+
   CALL histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s', &
           iip1, jjm, vhoriid, llm, 1, llm, zvertiid, &
           32, 'inst(X)', t_ops, t_wrt)
+
+  !  w
+
+  CALL histdef(fileid, 'w', 'flx de masse vert', 'kg m/s', &
+          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
+          32, 'inst(X)', t_ops, t_wrt)
+
+
+  !  Temperature potentielle
+
+  CALL histdef(fileid, 'teta', 'temperature potentielle', '-', &
+          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
+          32, 'inst(X)', t_ops, t_wrt)
   !
-  !  w
-  !
-  CALL histdef(fileid, 'w', 'flx de masse vert', 'kg m/s', &
-          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
-          32, 'inst(X)', t_ops, t_wrt)
-
-  !
-  !  Temperature potentielle
-  !
-  CALL histdef(fileid, 'teta', 'temperature potentielle', '-', &
-          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
-          32, 'inst(X)', t_ops, t_wrt)
-  !
-
-  !
+
+
   ! Geopotentiel
-  !
+
   CALL histdef(fileid, 'phi', 'geopotentiel instantane', '-', &
           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
           32, 'inst(X)', t_ops, t_wrt)
-  !
+
   !  Fin
-  !
+
   CALL histend(fileid)
   CALL histend(filevid)