/[lmdze]/trunk/Sources/dyn3d/initdynav.f
ViewVC logotype

Diff of /trunk/Sources/dyn3d/initdynav.f

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

revision 52 by guez, Fri Sep 23 12:28:01 2011 UTC revision 56 by guez, Tue Jan 10 19:02:02 2012 UTC
# Line 2  module initdynav_m Line 2  module initdynav_m
2    
3    implicit none    implicit none
4    
5      integer histaveid
6    
7  contains  contains
8    
9    subroutine initdynav(day0, anne0, tstep, nq, fileid, t_ops, t_wrt)    subroutine initdynav(day0, anne0, tstep, nq, t_ops, t_wrt)
10    
11      ! From initdynav.F, version 1.1.1.1, 2004/05/19 12:53:05      ! From initdynav.F, version 1.1.1.1, 2004/05/19 12:53:05
12      ! L. Fairhead, LMD      ! L. Fairhead, LMD
# Line 23  contains Line 25  contains
25    
26      integer, intent(in):: day0, anne0 ! date de reference      integer, intent(in):: day0, anne0 ! date de reference
27      real, intent(in):: tstep ! frequence d'ecriture      real, intent(in):: tstep ! frequence d'ecriture
28        integer, intent(in):: nq ! nombre de traceurs
29      real, intent(in):: t_ops ! frequence de l'operation pour IOIPSL      real, intent(in):: t_ops ! frequence de l'operation pour IOIPSL
30      real, intent(in):: t_wrt ! frequence d'ecriture sur le fichier      real, intent(in):: t_wrt ! frequence d'ecriture sur le fichier
     integer, intent(out):: fileid ! ID du fichier netcdf cree  
     integer, intent(in):: nq ! nombre de traceurs  
31    
32      ! Variables locales      ! Variables locales
33      integer thoriid, zvertiid      integer thoriid, zvertiid
# Line 52  contains Line 53  contains
53      enddo      enddo
54    
55      call histbeg_totreg('dyn_hist_ave.nc', rlong(:, 1), rlat(1, :), 1, iip1, &      call histbeg_totreg('dyn_hist_ave.nc', rlong(:, 1), rlat(1, :), 1, iip1, &
56           1, jjp1, itau_dyn, zjulian, tstep, thoriid, fileid)           1, jjp1, itau_dyn, zjulian, tstep, thoriid, histaveid)
57      call histvert(fileid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &      call histvert(histaveid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &
58           zvertiid)           zvertiid)
59    
60      call histdef(fileid, 'u', 'vents u scalaires moyennes', &      call histdef(histaveid, 'u', 'vents u scalaires moyennes', &
61           'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
62           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
63      call histdef(fileid, 'v', 'vents v scalaires moyennes', &      call histdef(histaveid, 'v', 'vents v scalaires moyennes', &
64           'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
65           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
66      call histdef(fileid, 'temp', 'temperature moyennee', 'K', &      call histdef(histaveid, 'temp', 'temperature moyennee', 'K', &
67           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
68           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
69      call histdef(fileid, 'theta', 'temperature potentielle', 'K', &      call histdef(histaveid, 'theta', 'temperature potentielle', 'K', &
70           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
71           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
72      call histdef(fileid, 'phi', 'geopotentiel moyenne', '-', &      call histdef(histaveid, 'phi', 'geopotentiel moyenne', '-', &
73           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
74           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
75    
76      ! Traceurs      ! Traceurs
77      DO iq = 1, nq      DO iq = 1, nq
78         call histdef(fileid, ttext(iq), ttext(iq), '-', &         call histdef(histaveid, ttext(iq), ttext(iq), '-', &
79              iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &              iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
80              'ave(X)', t_ops, t_wrt)              'ave(X)', t_ops, t_wrt)
81      enddo      enddo
82    
83      call histdef(fileid, 'masse', 'masse', 'kg', &      call histdef(histaveid, 'masse', 'masse', 'kg', &
84           iip1, jjp1, thoriid, 1, 1, 1, -99, &           iip1, jjp1, thoriid, 1, 1, 1, -99, &
85           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
86      call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa', &      call histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', &
87           iip1, jjp1, thoriid, 1, 1, 1, -99, &           iip1, jjp1, thoriid, 1, 1, 1, -99, &
88           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
89      call histdef(fileid, 'phis', 'geopotentiel au sol', '-', &      call histdef(histaveid, 'phis', 'geopotentiel au sol', '-', &
90           iip1, jjp1, thoriid, 1, 1, 1, -99, &           iip1, jjp1, thoriid, 1, 1, 1, -99, &
91           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
92    
93      call histend(fileid)      call histend(histaveid)
94    
95    end subroutine initdynav    end subroutine initdynav
96    
Legend:
Removed from v.52  
changed lines
  Added in v.56
  ViewVC Help
Powered by ViewVC 1.1.21