/[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 18 by guez, Thu Aug 7 12:29:13 2008 UTC trunk/Sources/dyn3d/inithist.f revision 253 by guez, Tue Jan 23 15:49:10 2018 UTC
# Line 1  Line 1 
1  module inithist_m  module inithist_m
2    
   ! This module is clean: no C preprocessor directive, no include line  
   
3    implicit none    implicit none
4    
5  contains    integer histid, histvid, histuid
6    
7    subroutine inithist(day0, anne0, tstep, nq, fileid, filevid, infile, t_ops, &  contains
        t_wrt)  
8    
9      ! From inithist.F,v 1.1.1.1 2004/05/19 12:53:05    subroutine inithist(tstep, nq, t_ops, t_wrt)
10    
11      !   Routine d'initialisation des ecritures des fichiers histoires LMDZ      ! From inithist.F, version 1.1.1.1, 2004/05/19 12:53:05
12      !   au format IOIPSL      ! L. Fairhead, LMD, 03/99
13    
14      !   Appels succesifs des routines: histbeg      ! Routine d'initialisation des écritures des fichiers histoires au
15      !                                  histhori      ! format IOIPSL.
     !                                  histver  
     !                                  histdef  
     !                                  histend  
   
     !   Entree:  
     !      infile: nom du fichier histoire a creer  
     !      day0,anne0: date de reference  
     !      tstep: duree du pas de temps en seconde  
     !      t_ops: frequence de l'operation pour IOIPSL  
     !      t_wrt: frequence d'ecriture sur le fichier  
     !      nq: nombre de traceurs  
   
     !   Sortie:  
     !      fileid: ID du fichier netcdf cree  
     !      filevid:ID du fichier netcdf pour la grille v  
   
     !   L. Fairhead, LMD, 03/99  
   
     USE IOIPSL  
     use dimens_m  
     use paramet_m  
     use comconst  
     use comvert  
     use logic  
     use comgeom  
     use serre  
     use temps  
     use ener  
     use iniadvtrac_m  
   
     !   Arguments  
     character(len=*) infile  
     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        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      integer tau0      ! 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      !  Initialisations      print *, "Call sequence information: inithist"
43        CALL ymds2ju(annee_ref, 1, day_ref, 0., julian)
44    
45      pi = 4. * atan (1.)      call histbeg_totreg("dyn_histu.nc", rlonu * 180. / pi, rlatu * 180. / pi, &
46             1, iip1, 1, jjp1, itau_dyn, julian, tstep, uhoriid, histuid)
47    
48      !  Appel a histbeg: creation du fichier netcdf et initialisations diverses      ! 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
50      zan = anne0      ! un meme fichier)
51      dayref = day0  
52      CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)      call histbeg_totreg('dyn_histv.nc', rlonv * 180. / pi, rlatv * 180. / pi, &
53      tau0 = itau_dyn           1, iip1, 1, jjm, itau_dyn, julian, tstep, vhoriid, histvid)
54    
55      do jj = 1, jjp1      call histbeg_totreg("dyn_hist.nc", rlonv * 180. / pi, rlatu * 180. / pi, &
56         do ii = 1, iip1           1, iip1, 1, jjp1, itau_dyn, julian, tstep, thoriid, histid)
57            rlong(ii,jj) = rlonu(ii) * 180. / pi  
58            rlat(ii,jj) = rlatu(jj) * 180. / pi      call histvert(histid, 'presnivs', 'Niveaux pression', 'mb', presnivs/100., &
59         enddo           zvertiid, 'down')
60        call histvert(histvid, 'presnivs', 'Niveaux pression', 'mb', &
61             presnivs/100., zvertiid, 'down')
62        call histvert(histuid, 'presnivs', 'Niveaux pression', 'mb', &
63             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, 'teta', 'temperature potentielle', '-', iip1, jjp1, &
70             thoriid, llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
71        call histdef(histid, 'phi', 'geopotentiel', '-', iip1, jjp1, thoriid, &
72             llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
73    
74        ! Traceurs
75        DO iq = 1, nq
76           call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &
77                llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
78      enddo      enddo
79    
80      call histbeg_totreg(infile, rlong(:,1), rlat(1,:), &      call histdef(histid, 'masse', 'masse', 'kg', iip1, jjp1, thoriid, llm, 1, &
81           1, iip1, 1, jjp1, &           llm, zvertiid, 'inst(X)', t_ops, t_wrt)
82           tau0, zjulian, tstep, uhoriid, fileid)      call histdef(histid, 'ps', 'pression naturelle au sol', 'Pa', iip1, jjp1, &
83      !           thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
84      !  Creation du fichier histoire pour la grille en V (oblige pour l'instant,      call histdef(histid, 'phis', 'geopotentiel au sol', '-', iip1, jjp1, &
85      !  IOIPSL ne permet pas de grilles avec des nombres de point differents dans           thoriid, 1, 1, 1, -99, 'inst(X)', t_ops, t_wrt)
86      !  un meme fichier)  
87        call histend(histid)
88      do jj = 1, jjm      call histend(histuid)
89         do ii = 1, iip1      call histend(histvid)
           rlong(ii,jj) = rlonv(ii) * 180. / pi  
           rlat(ii,jj) = rlatv(jj) * 180. / pi  
        enddo  
     enddo  
   
     call histbeg_totreg('dyn_histv.nc', rlong(:,1), rlat(1,:jjm), &  
          1, iip1, 1, jjm, &  
          tau0, zjulian, tstep, vhoriid, filevid)  
     !  
     !  Appel a histhori pour rajouter les autres grilles horizontales  
     !  
     do jj = 1, jjp1  
        do ii = 1, iip1  
           rlong(ii,jj) = rlonv(ii) * 180. / pi  
           rlat(ii,jj) = rlatu(jj) * 180. / pi  
        enddo  
     enddo  
   
     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)  
     !  
     !  Traceurs  
     !  
     DO iq=1,nq  
        call histdef(fileid, ttext(iq),  ttext(iq), '-', &  
             iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &  
             'inst(X)', t_ops, t_wrt)  
     enddo  
     !  
     !  Masse  
     !  
     call histdef(fileid, 'masse', 'masse', 'kg', &  
          iip1, jjp1, thoriid, 1, 1, 1, -99, &  
          'inst(X)', t_ops, t_wrt)  
     !  
     !  Pression au sol  
     !  
     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.18  
changed lines
  Added in v.253
  ViewVC Help
Powered by ViewVC 1.1.21