Ignore:
Timestamp:
Aug 2, 2024, 9:58:25 PM (7 weeks 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/initdynav_loc.f90

    r5158 r5159  
    1818  USE lmdz_comgeom
    1919
     20  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
     21  USE lmdz_paramet
    2022  IMPLICIT NONE
    2123
    22   !
     24
    2325  !   Routine d'initialisation des ecritures des fichiers histoires LMDZ
    2426  !   au format IOIPSL. Initialisation du fichier histoire moyenne.
    25   !
     27
    2628  !   Appels succesifs des routines: histbeg
    2729  !                              histhori
     
    2931  !                              histdef
    3032  !                              histend
    31   !
     33
    3234  !   Entree:
    33   !
     35
    3436  !  day0,anne0: date de reference
    3537  !  tstep : frequence d'ecriture
    3638  !  t_ops: frequence de l'operation pour IOIPSL
    3739  !  t_wrt: frequence d'ecriture sur le fichier
    38   !
     40
    3941  !   Sortie:
    4042  !  fileid: ID du fichier netcdf cree
    41   !
     43
    4244  !   L. Fairhead, LMD, 03/99
    43   !
     45
    4446  ! =====================================================================
    45   !
     47
    4648  !   Declarations
    47   INCLUDE "dimensions.h"
    48   INCLUDE "paramet.h"
     49
     50
    4951
    5052  !   Arguments
    51   !
     53
    5254  INTEGER(kind = 4) day0, anne0
    5355  REAL :: tstep, t_ops, t_wrt
     
    5557  ! This routine needs IOIPSL
    5658  !   Variables locales
    57   !
     59
    5860  INTEGER :: tau0
    5961  REAL :: zjulian
     
    8284  IF (adjust) return
    8385
    84   !
     86
    8587  !  Initialisations
    86   !
     88
    8789  pi = 4. * atan (1.)
    88   !
     90
    8991  !  Appel a histbeg: creation du fichier netcdf et initialisations diverses
    9092  !
     
    192194
    193195
    194   !
     196
    195197  !  Appel a histvert pour la grille verticale
    196   !
     198
    197199  CALL histvert(histaveid, 'presnivs', 'Niveaux Pression&
    198200          &     approximatifs', 'mb', llm, presnivs / 100., zvertiid, 'down')
     
    202204          &     approximatifs', 'mb', llm, presnivs / 100., zvertiidu, 'down')
    203205
    204   !
     206
    205207  !  Appels a histdef pour la definition des variables a sauvegarder
    206   !
     208
    207209  !  Vents U
    208   !
     210
    209211  jjn = jj_nb
    210212  CALL histdef(histuaveid, 'u', 'vent u moyen ', &
     
    212214          32, 'ave(X)', t_ops, t_wrt)
    213215
    214   !
     216
    215217  !  Vents V
    216   !
     218
    217219  IF (pole_sud) jjn = jj_nb - 1
    218220  CALL histdef(histvaveid, 'v', 'vent v moyen', &
     
    220222          32, 'ave(X)', t_ops, t_wrt)
    221223
    222   !
     224
    223225  !  Temperature
    224   !
     226
    225227  jjn = jj_nb
    226228  CALL histdef(histaveid, 'temp', 'temperature moyenne', 'K', &
    227229          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
    228230          32, 'ave(X)', t_ops, t_wrt)
    229   !
     231
    230232  !  Temperature potentielle
    231   !
     233
    232234  CALL histdef(histaveid, 'theta', 'temperature potentielle', 'K', &
    233235          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
     
    235237
    236238
    237   !
     239
    238240  !  Geopotentiel
    239   !
     241
    240242  CALL histdef(histaveid, 'phi', 'geopotentiel moyen', '-', &
    241243          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
    242244          32, 'ave(X)', t_ops, t_wrt)
    243   !
     245
    244246  !  Traceurs
    245   !
     247
    246248  !    DO iq=1,nqtot
    247249  !      CALL histdef(histaveid, tracers(iq)%name,
     
    250252  ! .             32, 'ave(X)', t_ops, t_wrt)
    251253  !    enddo
    252   !
     254
    253255  !  Masse
    254   !
     256
    255257  CALL histdef(histaveid, 'masse', 'masse moyenne', 'kg', &
    256258          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
    257259          32, 'ave(X)', t_ops, t_wrt)
    258   !
     260
    259261  !  Pression au sol
    260   !
     262
    261263  CALL histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', &
    262264          iip1, jjn, thoriid, 1, 1, 1, -99, &
    263265          32, 'ave(X)', t_ops, t_wrt)
    264   !
     266
    265267  !  Geopotentiel au sol
    266   !
     268
    267269  !  CALL histdef(histaveid, 'phis', 'geopotentiel au sol', '-',
    268270  ! .             iip1, jjn, thoriid, 1, 1, 1, -99,
    269271  ! .             32, 'ave(X)', t_ops, t_wrt)
    270   !
     272
    271273  !  Fin
    272   !
     274
    273275  CALL histend(histaveid)
    274276  CALL histend(histuaveid)
Note: See TracChangeset for help on using the changeset viewer.