/[lmdze]/trunk/dyn3d/inithist.f
ViewVC logotype

Diff of /trunk/dyn3d/inithist.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 254 by guez, Mon Feb 5 10:39:38 2018 UTC revision 302 by guez, Thu Sep 6 13:19:51 2018 UTC
# Line 6  module inithist_m Line 6  module inithist_m
6    
7  contains  contains
8    
9    subroutine inithist(tstep, nq, t_ops, t_wrt)    subroutine inithist(t_ops, t_wrt)
10    
11      ! From inithist.F, version 1.1.1.1, 2004/05/19 12:53:05      ! From inithist.F, version 1.1.1.1, 2004/05/19 12:53:05
12      ! L. Fairhead, LMD, 03/99      ! L. Fairhead, LMD, 03/99
# Line 14  contains Line 14  contains
14      ! Routine d'initialisation des écritures des fichiers histoires au      ! Routine d'initialisation des écritures des fichiers histoires au
15      ! format IOIPSL.      ! format IOIPSL.
16    
17      USE dimens_m, ONLY: jjm, llm      use comconst, only: dtvr
18        USE dimensions, ONLY: jjm, llm, nqmx
19      USE disvert_m, ONLY: presnivs      USE disvert_m, ONLY: presnivs
20      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv
21      use histbeg_totreg_m, only: histbeg_totreg      use histbeg_totreg_m, only: histbeg_totreg
# Line 27  contains Line 28  contains
28      USE temps, ONLY: itau_dyn      USE temps, ONLY: itau_dyn
29      use ymds2ju_m, ONLY: ymds2ju      use ymds2ju_m, ONLY: ymds2ju
30    
     real, intent(in):: tstep ! durée du pas de temps en secondes  
     integer, intent(in):: nq ! nombre de traceurs  
31      real, intent(in):: t_ops ! fréquence de l'opération pour IOIPSL      real, intent(in):: t_ops ! fréquence de l'opération pour IOIPSL
32      real, intent(in):: t_wrt ! fréquence d'écriture sur le fichier      real, intent(in):: t_wrt ! fréquence d'écriture sur le fichier
33    
# Line 43  contains Line 42  contains
42      CALL ymds2ju(annee_ref, 1, day_ref, 0., julian)      CALL ymds2ju(annee_ref, 1, day_ref, 0., julian)
43    
44      call histbeg_totreg("dyn_histu.nc", rlonu * 180. / pi, rlatu * 180. / pi, &      call histbeg_totreg("dyn_histu.nc", rlonu * 180. / pi, rlatu * 180. / pi, &
45           1, iip1, 1, jjp1, itau_dyn, julian, tstep, uhoriid, histuid)           1, iip1, 1, jjp1, itau_dyn, julian, dtvr, uhoriid, histuid)
46    
47      ! Creation du fichier histoire pour la grille en V (oblige pour l'instant,      ! Creation du fichier histoire pour la grille en V (oblige pour l'instant,
48      ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans      ! IOIPSL ne permet pas de grilles avec des nombres de point differents dans
49      ! un meme fichier)      ! un meme fichier)
50    
51      call histbeg_totreg('dyn_histv.nc', rlonv * 180. / pi, rlatv * 180. / pi, &      call histbeg_totreg('dyn_histv.nc', rlonv * 180. / pi, rlatv * 180. / pi, &
52           1, iip1, 1, jjm, itau_dyn, julian, tstep, vhoriid, histvid)           1, iip1, 1, jjm, itau_dyn, julian, dtvr, vhoriid, histvid)
53    
54      call histbeg_totreg("dyn_hist.nc", rlonv * 180. / pi, rlatu * 180. / pi, &      call histbeg_totreg("dyn_hist.nc", rlonv * 180. / pi, rlatu * 180. / pi, &
55           1, iip1, 1, jjp1, itau_dyn, julian, tstep, thoriid, histid)           1, iip1, 1, jjp1, itau_dyn, julian, dtvr, thoriid, histid)
56    
57      call histvert(histid, 'presnivs', 'Niveaux pression', 'mb', presnivs/100., &      call histvert(histid, 'presnivs', 'Niveaux pression', 'mb', presnivs/100., &
58           zvertiid, 'down')           zvertiid, 'down')
# Line 66  contains Line 65  contains
65           llm, zvertiid, 'inst(X)', t_ops, t_wrt)           llm, zvertiid, 'inst(X)', t_ops, t_wrt)
66      call histdef(histvid, 'v', 'vent v', 'm/s', iip1, jjm, vhoriid, llm, 1, &      call histdef(histvid, 'v', 'vent v', 'm/s', iip1, jjm, vhoriid, llm, 1, &
67           llm, zvertiid, 'inst(X)', t_ops, t_wrt)           llm, zvertiid, 'inst(X)', t_ops, t_wrt)
68      call histdef(histid, 'teta', 'temperature potentielle', '-', iip1, jjp1, &      call histdef(histid, 'temp', 'temperature', 'K', iip1, jjp1, &
69           thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)           thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
70      call histdef(histid, 'phi', 'geopotentiel', '-', iip1, jjp1, thoriid, &      call histdef(histid, 'theta', 'temperature potentielle', 'K', iip1, jjp1, &
71             thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
72        call histdef(histid, 'phi', 'geopotential', 'm2 s-2', iip1, jjp1, thoriid, &
73           llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)           llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
74    
75      ! Traceurs      ! Traceurs
76      DO iq = 1, nq      DO iq = 1, nqmx
77         call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &         call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &
78              llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)              llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
79      enddo      enddo
# Line 81  contains Line 82  contains
82           llm, zvertiid, 'inst(X)', t_ops, t_wrt)           llm, zvertiid, 'inst(X)', t_ops, t_wrt)
83      call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', iip1, jjp1, &      call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', iip1, jjp1, &
84           thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)           thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
     call histdef(histid, 'phis', 'geopotentiel au sol', '-', iip1, jjp1, &  
          thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)  
85    
86      call histend(histid)      call histend(histid)
87      call histend(histuid)      call histend(histuid)
Legend:
Removed from v.254  
changed lines
  Added in v.302
  ViewVC Help
Powered by ViewVC 1.1.21