/[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 40 by guez, Tue Feb 22 13:49:36 2011 UTC revision 62 by guez, Thu Jul 26 14:37:37 2012 UTC
# Line 1  Line 1 
1  module initdynav_m  module initdynav_m
2    
   ! This module is clean: no C preprocessor directive, no include line  
   
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 14  contains Line 14  contains
14      ! au format IOIPSL. Initialisation du fichier histoire moyenne.      ! au format IOIPSL. Initialisation du fichier histoire moyenne.
15    
16      use calendar, ONLY: ymds2ju      use calendar, ONLY: ymds2ju
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 iniadvtrac_m, ONLY : ttext      USE histdef_m, ONLY: histdef
22      USE nr_util, ONLY : pi      USE histend_m, ONLY: histend
23      USE paramet_m, ONLY : iip1, jjp1      USE histvert_m, ONLY: histvert
24      USE temps, ONLY : itau_dyn      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    
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 horiid, zvertiid
37      real zjulian      real julian
38      integer iq      integer iq
     real rlong(iip1, jjp1), rlat(iip1, jjp1)  
39      integer ii, jj      integer ii, jj
     integer zan, dayref  
40    
41      !----------------------------------------------------      !----------------------------------------------------
42    
43      print *, "Call sequence information: initdynav"      print *, "Call sequence information: initdynav"
44    
45      zan = anne0      CALL ymds2ju(anne0, 1, day0, 0., julian)
46      dayref = day0      call histbeg_totreg('dyn_hist_ave.nc', rlonv * 180. / pi, &
47      CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)           rlatu * 180. / pi, 1, iip1, 1, jjp1, itau_dyn, julian, tstep, &
48             horiid, histaveid)
49      do jj = 1, jjp1      call histvert(histaveid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &
        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_ave.nc', rlong(:, 1), rlat(1, :), 1, iip1, &  
          1, jjp1, itau_dyn, zjulian, tstep, thoriid, fileid)  
     call histvert(fileid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &  
50           zvertiid)           zvertiid)
51    
52      call histdef(fileid, 'u', 'vents u scalaires moyennes', &      call histdef(histaveid, 'u', 'vents u scalaires moyennes', 'm/s', iip1, &
53           'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           jjp1, horiid, llm, 1, llm, zvertiid, 'ave(X)', t_ops, t_wrt)
54           'ave(X)', t_ops, t_wrt)      call histdef(histaveid, 'v', 'vents v scalaires moyennes', 'm/s', iip1, &
55      call histdef(fileid, 'v', 'vents v scalaires moyennes', &           jjp1, horiid, llm, 1, llm, zvertiid, 'ave(X)', t_ops, t_wrt)
56           'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &      call histdef(histaveid, 'temp', 'temperature moyennee', 'K', iip1, jjp1, &
57           'ave(X)', t_ops, t_wrt)           horiid, llm, 1, llm, zvertiid, 'ave(X)', t_ops, t_wrt)
58      call histdef(fileid, 'temp', 'temperature moyennee', 'K', &      call histdef(histaveid, 'theta', 'temperature potentielle', 'K', iip1, &
59           iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &           jjp1, horiid, llm, 1, llm, zvertiid, 'ave(X)', t_ops, t_wrt)
60           'ave(X)', t_ops, t_wrt)      call histdef(histaveid, 'phi', 'geopotentiel moyenne', '-', iip1, jjp1, &
61      call histdef(fileid, 'theta', 'temperature potentielle', 'K', &           horiid, llm, 1, llm, zvertiid, 'ave(X)', t_ops, t_wrt)
          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &  
          'ave(X)', t_ops, t_wrt)  
     call histdef(fileid, 'phi', 'geopotentiel moyenne', '-', &  
          iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &  
          'ave(X)', t_ops, t_wrt)  
62    
63      ! Traceurs      ! Traceurs
64      DO iq = 1, nq      DO iq = 1, nq
65         call histdef(fileid, ttext(iq), ttext(iq), '-', &         call histdef(histaveid, ttext(iq), ttext(iq), '-', iip1, jjp1, &
66              iip1, jjp1, thoriid, llm, 1, llm, zvertiid, &              horiid, llm, 1, llm, zvertiid, 'ave(X)', t_ops, t_wrt)
             'ave(X)', t_ops, t_wrt)  
67      enddo      enddo
68    
69      call histdef(fileid, 'masse', 'masse', 'kg', &      call histdef(histaveid, 'masse', 'masse', 'kg', iip1, jjp1, horiid, 1, &
70           iip1, jjp1, thoriid, 1, 1, 1, -99, &           1, 1, -99, 'ave(X)', t_ops, t_wrt)
71           'ave(X)', t_ops, t_wrt)      call histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', iip1, &
72      call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa', &           jjp1, horiid, 1, 1, 1, -99, 'ave(X)', t_ops, t_wrt)
73           iip1, jjp1, thoriid, 1, 1, 1, -99, &      call histdef(histaveid, 'phis', 'geopotentiel au sol', '-', iip1, jjp1, &
74           'ave(X)', t_ops, t_wrt)           horiid, 1, 1, 1, -99, 'ave(X)', t_ops, t_wrt)
     call histdef(fileid, 'phis', 'geopotentiel au sol', '-', &  
          iip1, jjp1, thoriid, 1, 1, 1, -99, &  
          'ave(X)', t_ops, t_wrt)  
75    
76      call histend(fileid)      call histend(histaveid)
77    
78    end subroutine initdynav    end subroutine initdynav
79    
Legend:
Removed from v.40  
changed lines
  Added in v.62
  ViewVC Help
Powered by ViewVC 1.1.21