/[lmdze]/trunk/libf/bibio/initfluxsto.f90
ViewVC logotype

Diff of /trunk/libf/bibio/initfluxsto.f90

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

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 27 by guez, Thu Mar 25 14:29:07 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 tstep, t_ops, t_wrt
52        integer fileid, filevid,filedid        integer fileid, filevid,filedid
# Line 58  C Line 54  C
54        real nivd(1)        real nivd(1)
55    
56  C   Variables locales  C   Variables locales
 C  
       integer tau0  
57        real zjulian        real zjulian
58        character*3 str        character*3 str
59        character*10 ctrac        character*10 ctrac
60        integer iq        integer iq
61        real rlong(iip1,jjp1), rlat(iip1,jjp1),rl(1,1)        real rlong(iip1,jjp1), rlat(iip1,jjp1)
62        integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid        integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid
63        integer ii,jj        integer ii,jj
64        integer zan, idayref        integer zan, idayref
# Line 83  C Line 77  C
77        zan = annee_ref        zan = annee_ref
78        idayref = day_ref        idayref = day_ref
79        CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)        CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
       tau0 = itau_dyn  
80                    
81          do jj = 1, jjp1          do jj = 1, jjp1
82          do ii = 1, iip1          do ii = 1, iip1
# Line 92  C Line 85  C
85          enddo          enddo
86        enddo        enddo
87    
88        call histbeg_totreg(infile, iip1, rlong(:,1), jjp1, rlat(1,:),        call histbeg_totreg('fluxstoke', rlong(:,1), rlat(1,:),
89       .             1, iip1, 1, jjp1,       .             1, iip1, 1, jjp1,
90       .             tau0, zjulian, tstep, uhoriid, fileid)       .             itau_dyn, zjulian, tstep, uhoriid, fileid)
91  C  C
92  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,
93  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 101  C  un meme fichier)
101          enddo          enddo
102        enddo        enddo
103    
104        call histbeg_totreg('fluxstokev.nc', iip1, rlong(:,1), jjm,        call histbeg_totreg('fluxstokev.nc', rlong(:,1),
105       .             rlat(1,:),1, iip1, 1, jjm,       .             rlat(1,:jjm),1, iip1, 1, jjm,
106       .             tau0, zjulian, tstep, vhoriid, filevid)       .             itau_dyn, zjulian, tstep, vhoriid, filevid)
107                    
108          rl(1,1) = 1.            call histbeg_totreg('defstoke.nc', (/1./), (/1./),
       call histbeg_regular('defstoke.nc', 1, rl, 1, rl,  
109       .             1, 1, 1, 1,       .             1, 1, 1, 1,
110       .             tau0, zjulian, tstep, dhoriid, filedid)       .             itau_dyn, zjulian, tstep, dhoriid, filedid)
111    
112  C  C
113  C  Appel a histhori pour rajouter les autres grilles horizontales  C  Appel a histhori pour rajouter les autres grilles horizontales
# Line 127  C Line 119  C
119          enddo          enddo
120        enddo        enddo
121    
122        call histhori(fileid, iip1, rlong, jjp1, rlat, 'scalar',        call histhori_regular(fileid, iip1, rlong, jjp1, rlat, 'scalar',
123       .              'Grille points scalaires', thoriid)       .              'Grille points scalaires', thoriid)
124                    
125  C  C
# Line 150  C Line 142  C
142  C  Appels a histdef pour la definition des variables a sauvegarder  C  Appels a histdef pour la definition des variables a sauvegarder
143                    
144          CALL histdef(fileid, "phis", "Surface geop. height", "-",          CALL histdef(fileid, "phis", "Surface geop. height", "-",
145       .                iip1,jjp1,thoriid, 1,1,1, -99, 32,       .                iip1,jjp1,thoriid, 1,1,1, -99,
146       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
147    
148           CALL histdef(fileid, "aire", "Grid area", "-",           CALL histdef(fileid, "aire", "Grid area", "-",
149       .                iip1,jjp1,thoriid, 1,1,1, -99, 32,       .                iip1,jjp1,thoriid, 1,1,1, -99,
150       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
151                    
152          CALL histdef(filedid, "dtvr", "tps dyn", "s",          CALL histdef(filedid, "dtvr", "tps dyn", "s",
153       .                1,1,dhoriid, 1,1,1, -99, 32,       .                1,1,dhoriid, 1,1,1, -99,
154       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
155                    
156           CALL histdef(filedid, "istdyn", "tps stock", "s",           CALL histdef(filedid, "istdyn", "tps stock", "s",
157       .                1,1,dhoriid, 1,1,1, -99, 32,       .                1,1,dhoriid, 1,1,1, -99,
158       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
159                    
160           CALL histdef(filedid, "istphy", "tps stock phy", "s",           CALL histdef(filedid, "istphy", "tps stock phy", "s",
161       .                1,1,dhoriid, 1,1,1, -99, 32,       .                1,1,dhoriid, 1,1,1, -99,
162       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
163    
164    
# Line 175  C Masse Line 167  C Masse
167  C  C
168        call histdef(fileid, 'masse', 'Masse', 'kg',        call histdef(fileid, 'masse', 'Masse', 'kg',
169       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
170       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
171  C  C
172  C  Pbaru  C  Pbaru
173  C  C
174        call histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s',        call histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s',
175       .             iip1, jjp1, uhoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, uhoriid, llm, 1, llm, zvertiid,
176       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
177    
178  C  C
179  C  Pbarv  C  Pbarv
180  C  C
181        call histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s',        call histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s',
182       .             iip1, jjm, vhoriid, llm, 1, llm, zvertiid,       .             iip1, jjm, vhoriid, llm, 1, llm, zvertiid,
183       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
184  C  C
185  C  w  C  w
186  C  C
187        call histdef(fileid, 'w', 'flx de masse vert', 'kg m/s',        call histdef(fileid, 'w', 'flx de masse vert', 'kg m/s',
188       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
189       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
190    
191  C  C
192  C  Temperature potentielle  C  Temperature potentielle
193  C  C
194        call histdef(fileid, 'teta', 'temperature potentielle', '-',        call histdef(fileid, 'teta', 'temperature potentielle', '-',
195       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
196       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
197  C  C
198    
199  C  C
# Line 209  C Geopotentiel Line 201  C Geopotentiel
201  C  C
202        call histdef(fileid, 'phi', 'geopotentiel instantane', '-',        call histdef(fileid, 'phi', 'geopotentiel instantane', '-',
203       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
204       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
205  C  C
206  C  Fin  C  Fin
207  C  C
Legend:
Removed from v.3  
changed lines
  Added in v.27
  ViewVC Help
Powered by ViewVC 1.1.21