/[lmdze]/trunk/bibio/initfluxsto.f
ViewVC logotype

Diff of /trunk/bibio/initfluxsto.f

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 28 by guez, Fri Mar 26 18:33:04 2010 UTC
# Line 2  Line 2 
2  ! $Header: /home/cvsroot/LMDZ4/libf/bibio/initfluxsto.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $  ! $Header: /home/cvsroot/LMDZ4/libf/bibio/initfluxsto.F,v 1.1.1.1 2004/05/19 12:53:05 lmdzadmin Exp $
3  !  !
4        subroutine initfluxsto        subroutine initfluxsto
5       .  (infile,tstep,t_ops,t_wrt,nq,       .  (tstep,t_ops,t_wrt,nq,
6       .                    fileid,filevid,filedid)       .                    fileid,filevid,filedid)
7    
8         USE IOIPSL         USE IOIPSL
        USE histcom  
   
9    
10  C  C
11  C   Routine d'initialisation des ecritures des fichiers histoires LMDZ  C   Routine d'initialisation des ecritures des fichiers histoires LMDZ
# Line 21  C                                  histe Line 19  C                                  histe
19  C  C
20  C   Entree:  C   Entree:
21  C  C
 C      infile: nom du fichier histoire a creer  
22  C      day0,anne0: date de reference  C      day0,anne0: date de reference
23  C      tstep: duree du pas de temps en seconde  C      tstep: duree du pas de temps en seconde
24  C      t_ops: frequence de l'operation pour IOIPSL  C      t_ops: frequence de l'operation pour IOIPSL
# Line 50  C   Declarations Line 47  C   Declarations
47    
48  C   Arguments  C   Arguments
49  C  C
       character*(*) infile  
50        integer*4 itau        integer*4 itau
51        real tstep, t_ops, t_wrt        real, intent(in):: tstep
52          real t_ops, t_wrt
53        integer fileid, filevid,filedid        integer fileid, filevid,filedid
54        integer nq,ndex(1)        integer nq,ndex(1)
55        real nivd(1)        real nivd(1)
56    
57  C   Variables locales  C   Variables locales
 C  
       integer tau0  
58        real zjulian        real zjulian
59        character*3 str        character*3 str
60        character*10 ctrac        character*10 ctrac
61        integer iq        integer iq
62        real rlong(iip1,jjp1), rlat(iip1,jjp1),rl(1,1)        real rlong(iip1,jjp1), rlat(iip1,jjp1)
63        integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid        integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid
64        integer ii,jj        integer ii,jj
65        integer zan, idayref        integer zan, idayref
# Line 83  C Line 78  C
78        zan = annee_ref        zan = annee_ref
79        idayref = day_ref        idayref = day_ref
80        CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)        CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
       tau0 = itau_dyn  
81                    
82          do jj = 1, jjp1          do jj = 1, jjp1
83          do ii = 1, iip1          do ii = 1, iip1
# Line 92  C Line 86  C
86          enddo          enddo
87        enddo        enddo
88    
89        call histbeg_totreg(infile, iip1, rlong(:,1), jjp1, rlat(1,:),        call histbeg_totreg('fluxstoke', rlong(:,1), rlat(1,:),
90       .             1, iip1, 1, jjp1,       .             1, iip1, 1, jjp1,
91       .             tau0, zjulian, tstep, uhoriid, fileid)       .             itau_dyn, zjulian, tstep, uhoriid, fileid)
92  C  C
93  C  Creation du fichier histoire pour la grille en V (oblige pour l'instant,  C  Creation du fichier histoire pour la grille en V (oblige pour l'instant,
94  C  IOIPSL ne permet pas de grilles avec des nombres de point differents dans  C  IOIPSL ne permet pas de grilles avec des nombres de point differents dans
# Line 108  C  un meme fichier) Line 102  C  un meme fichier)
102          enddo          enddo
103        enddo        enddo
104    
105        call histbeg_totreg('fluxstokev.nc', iip1, rlong(:,1), jjm,        call histbeg_totreg('fluxstokev.nc', rlong(:,1),
106       .             rlat(1,:),1, iip1, 1, jjm,       .             rlat(1,:jjm),1, iip1, 1, jjm,
107       .             tau0, zjulian, tstep, vhoriid, filevid)       .             itau_dyn, zjulian, tstep, vhoriid, filevid)
108                    
109          rl(1,1) = 1.            call histbeg_totreg('defstoke.nc', (/1./), (/1./),
       call histbeg_regular('defstoke.nc', 1, rl, 1, rl,  
110       .             1, 1, 1, 1,       .             1, 1, 1, 1,
111       .             tau0, zjulian, tstep, dhoriid, filedid)       .             itau_dyn, zjulian, tstep, dhoriid, filedid)
112    
113  C  C
114  C  Appel a histhori pour rajouter les autres grilles horizontales  C  Appel a histhori pour rajouter les autres grilles horizontales
# Line 127  C Line 120  C
120          enddo          enddo
121        enddo        enddo
122    
123        call histhori(fileid, iip1, rlong, jjp1, rlat, 'scalar',        call histhori_regular(fileid, iip1, rlong, jjp1, rlat, 'scalar',
124       .              'Grille points scalaires', thoriid)       .              'Grille points scalaires', thoriid)
125                    
126  C  C
# Line 150  C Line 143  C
143  C  Appels a histdef pour la definition des variables a sauvegarder  C  Appels a histdef pour la definition des variables a sauvegarder
144                    
145          CALL histdef(fileid, "phis", "Surface geop. height", "-",          CALL histdef(fileid, "phis", "Surface geop. height", "-",
146       .                iip1,jjp1,thoriid, 1,1,1, -99, 32,       .                iip1,jjp1,thoriid, 1,1,1, -99,
147       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
148    
149           CALL histdef(fileid, "aire", "Grid area", "-",           CALL histdef(fileid, "aire", "Grid area", "-",
150       .                iip1,jjp1,thoriid, 1,1,1, -99, 32,       .                iip1,jjp1,thoriid, 1,1,1, -99,
151       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
152                    
153          CALL histdef(filedid, "dtvr", "tps dyn", "s",          CALL histdef(filedid, "dtvr", "tps dyn", "s",
154       .                1,1,dhoriid, 1,1,1, -99, 32,       .                1,1,dhoriid, 1,1,1, -99,
155       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
156                    
157           CALL histdef(filedid, "istdyn", "tps stock", "s",           CALL histdef(filedid, "istdyn", "tps stock", "s",
158       .                1,1,dhoriid, 1,1,1, -99, 32,       .                1,1,dhoriid, 1,1,1, -99,
159       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
160                    
161           CALL histdef(filedid, "istphy", "tps stock phy", "s",           CALL histdef(filedid, "istphy", "tps stock phy", "s",
162       .                1,1,dhoriid, 1,1,1, -99, 32,       .                1,1,dhoriid, 1,1,1, -99,
163       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
164    
165    
# Line 175  C Masse Line 168  C Masse
168  C  C
169        call histdef(fileid, 'masse', 'Masse', 'kg',        call histdef(fileid, 'masse', 'Masse', 'kg',
170       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
171       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
172  C  C
173  C  Pbaru  C  Pbaru
174  C  C
175        call histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s',        call histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s',
176       .             iip1, jjp1, uhoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, uhoriid, llm, 1, llm, zvertiid,
177       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
178    
179  C  C
180  C  Pbarv  C  Pbarv
181  C  C
182        call histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s',        call histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s',
183       .             iip1, jjm, vhoriid, llm, 1, llm, zvertiid,       .             iip1, jjm, vhoriid, llm, 1, llm, zvertiid,
184       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
185  C  C
186  C  w  C  w
187  C  C
188        call histdef(fileid, 'w', 'flx de masse vert', 'kg m/s',        call histdef(fileid, 'w', 'flx de masse vert', 'kg m/s',
189       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
190       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
191    
192  C  C
193  C  Temperature potentielle  C  Temperature potentielle
194  C  C
195        call histdef(fileid, 'teta', 'temperature potentielle', '-',        call histdef(fileid, 'teta', 'temperature potentielle', '-',
196       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
197       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
198  C  C
199    
200  C  C
# Line 209  C Geopotentiel Line 202  C Geopotentiel
202  C  C
203        call histdef(fileid, 'phi', 'geopotentiel instantane', '-',        call histdef(fileid, 'phi', 'geopotentiel instantane', '-',
204       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
205       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
206  C  C
207  C  Fin  C  Fin
208  C  C
Legend:
Removed from v.3  
changed lines
  Added in v.28
  ViewVC Help
Powered by ViewVC 1.1.21