/[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 61 by guez, Fri Apr 20 14:58:43 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 15  contains Line 17  contains
17      USE comvert, ONLY : nivsigs      USE comvert, ONLY : nivsigs
18      USE comgeom, ONLY : rlatu, rlonv      USE comgeom, ONLY : rlatu, rlonv
19      USE dimens_m, ONLY : llm      USE dimens_m, ONLY : llm
20      USE histcom, ONLY: histbeg_totreg, histdef, histend, histvert      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      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
# Line 23  contains Line 28  contains
28    
29      integer, intent(in):: day0, anne0 ! date de reference      integer, intent(in):: day0, anne0 ! date de reference
30      real, intent(in):: tstep ! frequence d'ecriture      real, intent(in):: tstep ! frequence d'ecriture
31        integer, intent(in):: nq ! nombre de traceurs
32      real, intent(in):: t_ops ! frequence de l'operation pour IOIPSL      real, intent(in):: t_ops ! frequence de l'operation pour IOIPSL
33      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  
34    
35      ! Variables locales      ! Variables locales
36      integer thoriid, zvertiid      integer thoriid, zvertiid
# Line 52  contains Line 56  contains
56      enddo      enddo
57    
58      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, &
59           1, jjp1, itau_dyn, zjulian, tstep, thoriid, fileid)           1, jjp1, itau_dyn, zjulian, tstep, thoriid, histaveid)
60      call histvert(fileid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &      call histvert(histaveid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &
61           zvertiid)           zvertiid)
62    
63      call histdef(fileid, 'u', 'vents u scalaires moyennes', &      call histdef(histaveid, 'u', 'vents u 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, 'v', 'vents v scalaires moyennes', &      call histdef(histaveid, 'v', 'vents v scalaires moyennes', &
67           'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
68           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
69      call histdef(fileid, 'temp', 'temperature moyennee', 'K', &      call histdef(histaveid, 'temp', 'temperature moyennee', '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, 'theta', 'temperature potentielle', 'K', &      call histdef(histaveid, 'theta', 'temperature potentielle', 'K', &
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      call histdef(fileid, 'phi', 'geopotentiel moyenne', '-', &      call histdef(histaveid, 'phi', 'geopotentiel moyenne', '-', &
76           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
77           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
78    
79      ! Traceurs      ! Traceurs
80      DO iq = 1, nq      DO iq = 1, nq
81         call histdef(fileid, ttext(iq), ttext(iq), '-', &         call histdef(histaveid, ttext(iq), ttext(iq), '-', &
82              iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &              iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &
83              'ave(X)', t_ops, t_wrt)              'ave(X)', t_ops, t_wrt)
84      enddo      enddo
85    
86      call histdef(fileid, 'masse', 'masse', 'kg', &      call histdef(histaveid, 'masse', 'masse', 'kg', &
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, 'ps', 'pression naturelle au sol', 'Pa', &      call histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', &
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      call histdef(fileid, 'phis', 'geopotentiel au sol', '-', &      call histdef(histaveid, 'phis', 'geopotentiel au sol', '-', &
93           iip1, jjp1, thoriid, 1, 1, 1, -99, &           iip1, jjp1, thoriid, 1, 1, 1, -99, &
94           'ave(X)', t_ops, t_wrt)           'ave(X)', t_ops, t_wrt)
95    
96      call histend(fileid)      call histend(histaveid)
97    
98    end subroutine initdynav    end subroutine initdynav
99    
Legend:
Removed from v.52  
changed lines
  Added in v.61
  ViewVC Help
Powered by ViewVC 1.1.21