/[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 61 by guez, Fri Apr 20 14:58:43 2012 UTC revision 62 by guez, Thu Jul 26 14:37:37 2012 UTC
# 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 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    
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
# Line 33  contains Line 33  contains
33      real, intent(in):: t_wrt ! frequence d'ecriture sur le fichier      real, intent(in):: t_wrt ! frequence d'ecriture sur le fichier
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)
     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_ave.nc', rlong(:, 1), rlat(1, :), 1, iip1, &  
          1, jjp1, itau_dyn, zjulian, tstep, thoriid, histaveid)  
49      call histvert(histaveid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &      call histvert(histaveid, 'sigss', 'Niveaux sigma', 'Pa', llm, nivsigs, &
50           zvertiid)           zvertiid)
51    
52      call histdef(histaveid, '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(histaveid, '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(histaveid, '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(histaveid, '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(histaveid, '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(histaveid, 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(histaveid, '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(histaveid, '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(histaveid, 'phis', 'geopotentiel au sol', '-', &  
          iip1, jjp1, thoriid, 1, 1, 1, -99, &  
          'ave(X)', t_ops, t_wrt)  
75    
76      call histend(histaveid)      call histend(histaveid)
77    
Legend:
Removed from v.61  
changed lines
  Added in v.62
  ViewVC Help
Powered by ViewVC 1.1.21