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

Diff of /trunk/dyn3d/inithist.f90

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

trunk/libf/bibio/inithist.f90 revision 56 by guez, Tue Jan 10 19:02:02 2012 UTC trunk/dyn3d/inithist.f revision 261 by guez, Wed Mar 7 13:33:15 2018 UTC
# Line 2  module inithist_m Line 2  module inithist_m
2    
3    implicit none    implicit none
4    
5  contains    integer histid, histvid, histuid
   
   subroutine inithist(day0, anne0, tstep, nq, t_ops, t_wrt)  
6    
7      ! From inithist.F, version 1.1.1.1 2004/05/19 12:53:05  contains
8    
9      ! Routine d'initialisation des écritures des fichiers histoires LMDZ    subroutine inithist(tstep, nq, t_ops, t_wrt)
     ! au format IOIPSL  
     ! Appels successifs des routines : histbeg, histhori, histver,  
     ! histdef, histend  
   
     ! Entrées :  
     ! day0, anne0: date de référence  
     ! tstep : durée du pas de temps en secondes  
     ! t_ops : fréquence de l'opération pour IOIPSL  
     ! t_wrt : fréquence d'écriture sur le fichier  
     ! nq : nombre de traceurs  
10    
11        ! 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
13    
14      USE calendar      ! Routine d'initialisation des écritures des fichiers histoires au
15      use com_io_dyn, only: histid, histvid, histuid      ! format IOIPSL.
16      use histcom  
17      use dimens_m      USE dimens_m, ONLY: jjm, llm
18      use paramet_m      USE disvert_m, ONLY: presnivs
19      use comconst      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv
20      use comvert      use histbeg_totreg_m, only: histbeg_totreg
21      use logic      USE histdef_m, ONLY: histdef
22      use comgeom      USE histend_m, ONLY: histend
23      use serre      USE histvert_m, ONLY: histvert
24      use temps      USE iniadvtrac_m, ONLY: ttext
25      use ener      USE nr_util, ONLY: pi
26      use iniadvtrac_m      USE paramet_m, ONLY: iip1, jjp1
27      use nr_util, only: pi      USE temps, ONLY: itau_dyn
28        use ymds2ju_m, ONLY: ymds2ju
29      ! Arguments  
30      integer day0, anne0      real, intent(in):: tstep ! durée du pas de temps en secondes
31      real, intent(in):: tstep, t_ops, t_wrt      integer, intent(in):: nq ! nombre de traceurs
32      integer nq      real, intent(in):: t_ops ! fréquence de l'opération pour IOIPSL
33        real, intent(in):: t_wrt ! fréquence d'écriture sur le fichier
34    
35      ! Variables locales      ! Local:
36      real zjulian      real julian
37      integer iq      integer iq
     real rlong(iip1, jjp1), rlat(iip1, jjp1)  
38      integer uhoriid, vhoriid, thoriid, zvertiid      integer uhoriid, vhoriid, thoriid, zvertiid
     integer ii, jj  
     integer zan, dayref  
39    
40      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
41    
42      ! Appel a histbeg: creation du fichier netcdf et initialisations diverses      print *, "Call sequence information: inithist"
43        CALL ymds2ju(annee_ref, 1, day_ref, 0., julian)
44    
45      zan = anne0      call histbeg_totreg("dyn_histu.nc", rlonu * 180. / pi, rlatu * 180. / pi, &
46      dayref = day0           1, iip1, 1, jjp1, itau_dyn, julian, tstep, uhoriid, histuid)
     CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)  
   
     do jj = 1, jjp1  
        do ii = 1, iip1  
           rlong(ii, jj) = rlonu(ii) * 180. / pi  
           rlat(ii, jj) = rlatu(jj) * 180. / pi  
        enddo  
     enddo  
   
     call histbeg_totreg("dyn_histu.nc", rlong(:,1), rlat(1,:), 1, iip1, 1, &  
          jjp1, itau_dyn, zjulian, tstep, uhoriid, histuid)  
   
     do jj = 1, jjm  
        do ii = 1, iip1  
           rlong(ii, jj) = rlonv(ii) * 180. / pi  
           rlat(ii, jj) = rlatv(jj) * 180. / pi  
        enddo  
     enddo  
47    
48      ! 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,
49      ! 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
50      ! un meme fichier)      ! un meme fichier)
51    
52      call histbeg_totreg('dyn_histv.nc', rlong(:, 1), rlat(1, :jjm), &      call histbeg_totreg('dyn_histv.nc', rlonv * 180. / pi, rlatv * 180. / pi, &
53           1, iip1, 1, jjm, &           1, iip1, 1, jjm, itau_dyn, julian, tstep, vhoriid, histvid)
54           itau_dyn, zjulian, tstep, vhoriid, histvid)  
55      !      call histbeg_totreg("dyn_hist.nc", rlonv * 180. / pi, rlatu * 180. / pi, &
56      ! Appel a histhori pour rajouter les autres grilles horizontales           1, iip1, 1, jjp1, itau_dyn, julian, tstep, thoriid, histid)
57      !  
58      do jj = 1, jjp1      call histvert(histid, 'presnivs', 'Niveaux pression', 'mb', presnivs/100., &
59         do ii = 1, iip1           zvertiid, 'down')
60            rlong(ii, jj) = rlonv(ii) * 180. / pi      call histvert(histvid, 'presnivs', 'Niveaux pression', 'mb', &
61            rlat(ii, jj) = rlatu(jj) * 180. / pi           presnivs/100., zvertiid, 'down')
62         enddo      call histvert(histuid, 'presnivs', 'Niveaux pression', 'mb', &
63      enddo           presnivs/100., zvertiid, 'down')
64    
65        call histdef(histuid, 'u', 'vent u', 'm/s', iip1, jjp1, uhoriid, llm, 1, &
66             llm, zvertiid, 'inst(X)', t_ops, t_wrt)
67        call histdef(histvid, 'v', 'vent v', 'm/s', iip1, jjm, vhoriid, llm, 1, &
68             llm, zvertiid, 'inst(X)', t_ops, t_wrt)
69        call histdef(histid, 'temp', 'temperature', 'K', iip1, jjp1, &
70             thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
71        call histdef(histid, 'theta', 'temperature potentielle', 'K', iip1, jjp1, &
72             thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
73        call histdef(histid, 'phi', 'geopotentiel', '-', iip1, jjp1, thoriid, &
74             llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
75    
     call histbeg_totreg("dyn_hist.nc", rlong(:, 1), rlat(1, :), 1, iip1, 1, &  
          jjp1, itau_dyn, zjulian, tstep, thoriid, histid)  
     !  
     ! Appel a histvert pour la grille verticale  
     !  
       call histvert(histid, 'presnivs', 'Niveaux pression','mb', llm, &  
            presnivs/100., zvertiid,'down')  
       call histvert(histvid, 'presnivs', 'Niveaux pression','mb', llm, &  
            presnivs/100., zvertiid,'down')  
       call histvert(histuid, 'presnivs', 'Niveaux pression','mb', llm, &  
            presnivs/100., zvertiid,'down')  
     !  
     ! Appels a histdef pour la definition des variables a sauvegarder  
     !  
     ! Vents U  
     !  
     call histdef(histuid, 'u', 'vent u', 'm/s', &  
          iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     ! Vents V  
     !  
     call histdef(histvid, 'v', 'vent v', 'm/s', &  
          iip1, jjm, vhoriid, llm, 1, llm, zvertiid, &  
          'inst(X)', t_ops, t_wrt)  
   
     !  
     ! Temperature potentielle  
     !  
     call histdef(histid, 'teta', 'temperature potentielle', '-', &  
          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     ! Geopotentiel  
     !  
     call histdef(histid, 'phi', 'geopotentiel', '-', &  
          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &  
          'inst(X)', t_ops, t_wrt)  
     !  
76      ! Traceurs      ! Traceurs
77      !      DO iq = 1, nq
78      DO iq=1, nq         call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &
79         call histdef(histid, ttext(iq), ttext(iq), '-', &              llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
             iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &  
             'inst(X)', t_ops, t_wrt)  
80      enddo      enddo
81      !  
82      ! Masse      call histdef(histid, 'masse', 'masse', 'kg', iip1, jjp1, thoriid, llm, 1, &
83      !           llm, zvertiid, 'inst(X)', t_ops, t_wrt)
84      call histdef(histid, 'masse', 'masse', 'kg', &      call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', iip1, jjp1, &
85           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
86           'inst(X)', t_ops, t_wrt)  
     !  
     ! Pression au sol  
     !  
     call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', &  
          iip1, jjp1, thoriid, 1, 1, 1, -99, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     ! Geopotentiel au sol  
     !  
     call histdef(histid, 'phis', 'geopotentiel au sol', '-', &  
          iip1, jjp1, thoriid, 1, 1, 1, -99, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     ! Fin  
     !  
87      call histend(histid)      call histend(histid)
88      call histend(histuid)      call histend(histuid)
89      call histend(histvid)      call histend(histvid)
Legend:
Removed from v.56  
changed lines
  Added in v.261
  ViewVC Help
Powered by ViewVC 1.1.21