/[lmdze]/trunk/dyn3d/bilan_dyn.f
ViewVC logotype

Diff of /trunk/dyn3d/bilan_dyn.f

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

revision 15 by guez, Fri Aug 1 15:24:12 2008 UTC revision 39 by guez, Tue Jan 25 15:11:05 2011 UTC
# Line 9  c   Prevoir en champ nq+1 le diagnostiqu Line 9  c   Prevoir en champ nq+1 le diagnostiqu
9  c   en faisant Qzon=Cv T + L * ...  c   en faisant Qzon=Cv T + L * ...
10  c             vQ..A=Cp T + L * ...  c             vQ..A=Cp T + L * ...
11    
12        USE IOIPSL        USE histcom
13          use calendar
14          use histwrite_m
15        use dimens_m        use dimens_m
16        use paramet_m        use paramet_m
17        use comconst        use comconst
# Line 19  c             vQ..A=Cp T + L * ... Line 20  c             vQ..A=Cp T + L * ...
20        use temps        use temps
21        use iniprint        use iniprint
22        use inigrads_m, only: inigrads        use inigrads_m, only: inigrads
23          use nr_util, only: pi
24    
25        IMPLICIT NONE        IMPLICIT NONE
26    
# Line 141  cym      character*10 zunites(ntr,nQ) Line 143  cym      character*10 zunites(ntr,nQ)
143  c   Initialisation du fichier contenant les moyennes zonales.  c   Initialisation du fichier contenant les moyennes zonales.
144  c   ---------------------------------------------------------  c   ---------------------------------------------------------
145    
       character*10 infile  
   
146        integer fileid        integer fileid
147        integer thoriid, zvertiid        integer thoriid, zvertiid
148        save fileid        save fileid
# Line 151  c   ------------------------------------ Line 151  c   ------------------------------------
151    
152  C   Variables locales  C   Variables locales
153  C  C
       integer tau0  
154        real zjulian        real zjulian
155        character*3 str        character*3 str
156        character*10 ctrac        character*10 ctrac
# Line 187  c   ncum est la frequence de stokage en Line 186  c   ncum est la frequence de stokage en
186    
187          if (i_sortie.eq.1) then          if (i_sortie.eq.1) then
188           file='dynzon'           file='dynzon'
189           call inigrads(ifile           call inigrads(ifile ,(/0./),180./pi,0.,0.,rlatv,-90.,90.,
190       s  ,(/0./),180./pi,0.,0.,rlatv,-90.,90.,180./pi       $        180./pi ,presnivs,1. ,dt_cum,file,'dyn_zon ')
      s  ,presnivs,1.  
      s  ,dt_cum,file,'dyn_zon ')  
191          endif          endif
192    
193          nom(itemp)='T'          nom(itemp)='T'
# Line 213  c   ncum est la frequence de stokage en Line 210  c   ncum est la frequence de stokage en
210  c   Initialisation du fichier contenant les moyennes zonales.  c   Initialisation du fichier contenant les moyennes zonales.
211  c   ---------------------------------------------------------  c   ---------------------------------------------------------
212    
       infile='dynzon'  
   
213        zan = annee_ref        zan = annee_ref
214        dayref = day_ref        dayref = day_ref
215        CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)        CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
       tau0 = itau_dyn  
216                
217        rlong=0.        rlong=0.
218        rlatg=rlatv*180./pi        rlatg=rlatv*180./pi
219                
220        call histbeg_totreg(infile, rlong(:1), rlatg,        call histbeg_totreg('dynzon', rlong(:1), rlatg,
221       .             1, 1, 1, jjm,       .             1, 1, 1, jjm,
222       .             tau0, zjulian, dt_cum, thoriid, fileid)       .             itau_dyn, zjulian, dt_cum, thoriid, fileid)
223    
224  C  C
225  C  Appel a histvert pour la grille verticale  C  Appel a histvert pour la grille verticale
Legend:
Removed from v.15  
changed lines
  Added in v.39
  ViewVC Help
Powered by ViewVC 1.1.21