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

Diff of /trunk/dyn3d/inithist.f

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

trunk/libf/bibio/inithist.f90 revision 52 by guez, Fri Sep 23 12:28:01 2011 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, fileid, filevid, 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  
   
     ! Sorties :  
     ! fileid : ID du fichier Netcdf créé  
     ! filevid : ID du fichier Netcdf pour la grille v  
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 histcom      ! format IOIPSL.
     use dimens_m  
     use paramet_m  
     use comconst  
     use comvert  
     use logic  
     use comgeom  
     use serre  
     use temps  
     use ener  
     use iniadvtrac_m  
     use nr_util, only: pi  
   
     ! Arguments  
     integer day0, anne0  
     real, intent(in):: tstep, t_ops, t_wrt  
     integer fileid, filevid  
     integer nq  
16    
17      ! Variables locales      USE dimens_m, ONLY: jjm, llm
18      real zjulian      USE disvert_m, ONLY: presnivs
19        use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv
20        use histbeg_totreg_m, only: histbeg_totreg
21        USE histdef_m, ONLY: histdef
22        USE histend_m, ONLY: histend
23        USE histvert_m, ONLY: histvert
24        USE iniadvtrac_m, ONLY: ttext
25        USE nr_util, ONLY: pi
26        USE paramet_m, ONLY: iip1, jjp1
27        USE temps, ONLY: itau_dyn
28        use ymds2ju_m, ONLY: ymds2ju
29    
30        real, intent(in):: tstep ! durée du pas de temps en secondes
31        integer, intent(in):: nq ! nombre de traceurs
32        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        ! Local:
36        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  
47    
     call histbeg_totreg("dyn_hist.nc", rlong(:, 1), rlat(1, :), &  
          1, iip1, 1, jjp1, &  
          itau_dyn, zjulian, tstep, uhoriid, fileid)  
     !  
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      do jj = 1, jjm      call histbeg_totreg('dyn_histv.nc', rlonv * 180. / pi, rlatv * 180. / pi, &
53         do ii = 1, iip1           1, iip1, 1, jjm, itau_dyn, julian, tstep, vhoriid, histvid)
           rlong(ii, jj) = rlonv(ii) * 180. / pi  
           rlat(ii, jj) = rlatv(jj) * 180. / pi  
        enddo  
     enddo  
54    
55      call histbeg_totreg('dyn_histv.nc', rlong(:, 1), rlat(1, :jjm), &      call histbeg_totreg("dyn_hist.nc", rlonv * 180. / pi, rlatu * 180. / pi, &
56           1, iip1, 1, jjm, &           1, iip1, 1, jjp1, itau_dyn, julian, tstep, thoriid, histid)
57           itau_dyn, zjulian, tstep, vhoriid, filevid)  
58      !      call histvert(histid, 'presnivs', 'Niveaux pression', 'mb', presnivs/100., &
59      ! Appel a histhori pour rajouter les autres grilles horizontales           zvertiid, 'down')
60      !      call histvert(histvid, 'presnivs', 'Niveaux pression', 'mb', &
61      do jj = 1, jjp1           presnivs/100., zvertiid, 'down')
62         do ii = 1, iip1      call histvert(histuid, 'presnivs', 'Niveaux pression', 'mb', &
63            rlong(ii, jj) = rlonv(ii) * 180. / pi           presnivs/100., zvertiid, 'down')
64            rlat(ii, jj) = rlatu(jj) * 180. / pi  
65         enddo      call histdef(histuid, 'u', 'vent u', 'm/s', iip1, jjp1, uhoriid, llm, 1, &
66      enddo           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 histhori_regular(fileid, iip1, rlong, jjp1, rlat, 'scalar', &  
          'Grille points scalaires', thoriid)  
     !  
     ! Appel a histvert pour la grille verticale  
     !  
     call histvert(fileid, 'sig_s', 'Niveaux sigma', '-', &  
          llm, nivsigs, zvertiid)  
     ! Pour le fichier V  
     call histvert(filevid, 'sig_s', 'Niveaux sigma', '-', &  
          llm, nivsigs, zvertiid)  
     !  
     ! Appels a histdef pour la definition des variables a sauvegarder  
     !  
     ! Vents U  
     !  
     call histdef(fileid, 'ucov', 'vents u covariants', 'm/s', &  
          iip1, jjp1, uhoriid, llm, 1, llm, zvertiid, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     ! Vents V  
     !  
     call histdef(filevid, 'vcov', 'vents v covariants', 'm/s', &  
          iip1, jjm, vhoriid, llm, 1, llm, zvertiid, &  
          'inst(X)', t_ops, t_wrt)  
   
     !  
     ! Temperature potentielle  
     !  
     call histdef(fileid, 'teta', 'temperature potentielle', '-', &  
          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     ! Geopotentiel  
     !  
     call histdef(fileid, 'phi', 'geopotentiel instantane', '-', &  
          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(fileid, 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(fileid, 'masse', 'masse', 'kg', &      call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', iip1, jjp1, &
85           iip1, jjp1, thoriid, 1, 1, 1, -99, &           thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
86           'inst(X)', t_ops, t_wrt)  
87      !      call histend(histid)
88      ! Pression au sol      call histend(histuid)
89      !      call histend(histvid)
     call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa', &  
          iip1, jjp1, thoriid, 1, 1, 1, -99, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     ! Pression au sol  
     !  
     call histdef(fileid, 'phis', 'geopotentiel au sol', '-', &  
          iip1, jjp1, thoriid, 1, 1, 1, -99, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     ! Fin  
     !  
     call histend(fileid)  
     call histend(filevid)  
90    
91    end subroutine inithist    end subroutine inithist
92    
Legend:
Removed from v.52  
changed lines
  Added in v.261
  ViewVC Help
Powered by ViewVC 1.1.21