/[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 212 by guez, Thu Jan 12 12:31:31 2017 UTC revision 253 by guez, Tue Jan 23 15:49:10 2018 UTC
# Line 8  contains Line 8  contains
8    
9    subroutine inithist(tstep, nq, t_ops, t_wrt)    subroutine inithist(tstep, nq, 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
13    
14      ! Routine d'initialisation des écritures des fichiers histoires      ! Routine d'initialisation des écritures des fichiers histoires au
15      ! LMDZ au format IOIPSL.      ! format IOIPSL.
16    
17      USE dimens_m, ONLY: jjm, llm      USE dimens_m, ONLY: jjm, llm
18      USE disvert_m, ONLY: presnivs      USE disvert_m, ONLY: presnivs
19      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv      use dynetat0_m, only: day_ref, annee_ref, rlatu, rlatv, rlonu, rlonv
20      USE histbeg_totreg_m, ONLY : histbeg_totreg      use histbeg_totreg_m, only: histbeg_totreg
21      USE histdef_m, ONLY : histdef      USE histdef_m, ONLY: histdef
22      USE histend_m, ONLY : histend      USE histend_m, ONLY: histend
23      USE histvert_m, ONLY : histvert      USE histvert_m, ONLY: histvert
24      USE iniadvtrac_m, ONLY: ttext      USE iniadvtrac_m, ONLY: ttext
25      USE nr_util, ONLY: pi      USE nr_util, ONLY: pi
26      USE paramet_m, ONLY: iip1, jjp1      USE paramet_m, ONLY: iip1, jjp1
27      USE temps, ONLY: itau_dyn      USE temps, ONLY: itau_dyn
28      USE ymds2ju_m, ONLY: ymds2ju      use ymds2ju_m, ONLY: ymds2ju
29    
30      real, intent(in):: tstep ! durée du pas de temps en secondes      real, intent(in):: tstep ! durée du pas de temps en secondes
31      integer, intent(in):: nq ! nombre de traceurs      integer, intent(in):: nq ! nombre de traceurs
32      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
33      real, intent(in):: t_wrt ! fréquence d'écriture sur le fichier      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  
39    
40      !-----------------------------------------------------------------------      !-----------------------------------------------------------------------
41    
42      CALL ymds2ju(annee_ref, 1, day_ref, 0.0, zjulian)      print *, "Call sequence information: inithist"
43        CALL ymds2ju(annee_ref, 1, day_ref, 0., julian)
44    
45      do jj = 1, jjp1      call histbeg_totreg("dyn_histu.nc", rlonu * 180. / pi, rlatu * 180. / pi, &
46         do ii = 1, iip1           1, iip1, 1, jjp1, itau_dyn, julian, tstep, uhoriid, histuid)
           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), 1, iip1, &      call histbeg_totreg('dyn_histv.nc', rlonv * 180. / pi, rlatv * 180. / pi, &
53           1, jjm, itau_dyn, zjulian, tstep, vhoriid, histvid)           1, iip1, 1, jjm, itau_dyn, julian, tstep, vhoriid, histvid)
   
     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 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  
54    
55      call histvert(histid, 'presnivs', 'Niveaux pression','mb', presnivs/100., &      call histbeg_totreg("dyn_hist.nc", rlonv * 180. / pi, rlatu * 180. / pi, &
56           zvertiid,'down')           1, iip1, 1, jjp1, itau_dyn, julian, tstep, thoriid, histid)
     call histvert(histvid, 'presnivs', 'Niveaux pression','mb', &  
          presnivs/100., zvertiid,'down')  
     call histvert(histuid, 'presnivs', 'Niveaux pression','mb', &  
          presnivs/100., zvertiid,'down')  
57    
58      ! Appels a histdef pour la definition des variables a sauvegarder      call histvert(histid, 'presnivs', 'Niveaux pression', 'mb', presnivs/100., &
59             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, &      call histdef(histuid, 'u', 'vent u', 'm/s', iip1, jjp1, uhoriid, llm, 1, &
66           llm, zvertiid, 'inst(X)', t_ops, t_wrt)           llm, zvertiid, 'inst(X)', t_ops, t_wrt)
# Line 98  contains Line 72  contains
72           llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)           llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
73    
74      ! Traceurs      ! Traceurs
75      DO iq=1, nq      DO iq = 1, nq
76         call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &         call histdef(histid, ttext(iq), ttext(iq), '-', iip1, jjp1, thoriid, &
77              llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)              llm, 1, llm, zvertiid, 'inst(X)', t_ops, t_wrt)
78      enddo      enddo
Legend:
Removed from v.212  
changed lines
  Added in v.253
  ViewVC Help
Powered by ViewVC 1.1.21