Ignore:
Timestamp:
Aug 2, 2024, 9:58:25 PM (3 months ago)
Author:
abarral
Message:

Put dimensions.h and paramet.h into modules

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/branches/Amaury_dev/libf/dyn3dmem/inithist_loc.F90

    r5158 r5159  
    1717  USE lmdz_comgeom
    1818
     19  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
     20  USE lmdz_paramet
    1921  IMPLICIT NONE
    2022
    21   !
     23
    2224  !   Routine d'initialisation des ecritures des fichiers histoires LMDZ
    2325  !   au format IOIPSL
    24   !
     26
    2527  !   Appels succesifs des routines: histbeg
    2628  !                              histhori
     
    2830  !                              histdef
    2931  !                              histend
    30   !
     32
    3133  !   Entree:
    32   !
     34
    3335  !  day0,anne0: date de reference
    3436  !  tstep: duree du pas de temps en seconde
     
    3638  !  t_wrt: frequence d'ecriture sur le fichier
    3739  !  nq: nombre de traceurs
    38   !
    39   !
     40
     41
    4042  !   L. Fairhead, LMD, 03/99
    41   !
     43
    4244  ! =====================================================================
    43   !
     45
    4446  !   Declarations
    45   INCLUDE "dimensions.h"
    46   INCLUDE "paramet.h"
     47
     48
    4749
    4850  !   Arguments
    49   !
     51
    5052  INTEGER :: day0, anne0
    5153  REAL :: tstep, t_ops, t_wrt
     
    5355  ! This routine needs IOIPSL
    5456  !   Variables locales
    55   !
     57
    5658  INTEGER :: tau0
    5759  REAL :: zjulian
     
    8082  IF (adjust) return
    8183
    82   !
     84
    8385  !  Initialisations
    84   !
     86
    8587  pi = 4. * atan (1.)
    86   !
     88
    8789  !  Appel a histbeg: creation du fichier netcdf et initialisations diverses
    8890  !
     
    200202          llm, presnivs / 100., zvertiidu, 'down')
    201203
    202   !
     204
    203205  ! -------------------------------------------------------------
    204206  !  Appels a histdef pour la definition des variables a sauvegarder
    205207  ! -------------------------------------------------------------
    206   !
     208
    207209  !  Vents U
    208   !
     210
    209211  jjn = jj_nb
    210212  CALL histdef(histuid, 'u', 'vent u', &
     
    212214          32, 'inst(X)', t_ops, t_wrt)
    213215
    214   !
     216
    215217  !  Vents V
    216   !
     218
    217219  IF (pole_sud) jjn = jj_nb - 1
    218220  CALL histdef(histvid, 'v', 'vent v', &
     
    220222          32, 'inst(X)', t_ops, t_wrt)
    221223
    222   !
     224
    223225  !  Temperature
    224   !
     226
    225227  jjn = jj_nb
    226228  CALL histdef(histid, 'temp', 'temperature', 'K', &
    227229          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
    228230          32, 'inst(X)', t_ops, t_wrt)
    229   !
     231
    230232  !  Temperature potentielle
    231   !
     233
    232234  CALL histdef(histid, 'theta', 'temperature potentielle', 'K', &
    233235          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
     
    235237
    236238
    237   !
     239
    238240  !  Geopotentiel
    239   !
     241
    240242  CALL histdef(histid, 'phi', 'geopotentiel', '-', &
    241243          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
    242244          32, 'inst(X)', t_ops, t_wrt)
    243   !
     245
    244246  !  Traceurs
    245   !
     247
    246248  !    DO iq=1,nqtot
    247249  !      CALL histdef(histid, tracers(iq)%name,
     
    250252  ! .             32, 'inst(X)', t_ops, t_wrt)
    251253  !    enddo
    252   !
     254
    253255  !  Masse
    254   !
     256
    255257  CALL histdef(histid, 'masse', 'masse', 'kg', &
    256258          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
    257259          32, 'inst(X)', t_ops, t_wrt)
    258   !
     260
    259261  !  Pression au sol
    260   !
     262
    261263  CALL histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', &
    262264          iip1, jjn, thoriid, 1, 1, 1, -99, &
    263265          32, 'inst(X)', t_ops, t_wrt)
    264   !
     266
    265267  !  Geopotentiel au sol
    266   !
     268
    267269  !  CALL histdef(histid, 'phis', 'geopotentiel au sol', '-',
    268270  ! .             iip1, jjn, thoriid, 1, 1, 1, -99,
    269271  ! .             32, 'inst(X)', t_ops, t_wrt)
    270   !
     272
    271273  !  Fin
    272   !
     274
    273275  CALL histend(histid)
    274276  CALL histend(histuid)
Note: See TracChangeset for help on using the changeset viewer.