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

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

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

revision 14 by guez, Mon Mar 31 13:58:05 2008 UTC revision 15 by guez, Fri Aug 1 15:24:12 2008 UTC
# Line 62  C Line 62  C
62        character*3 str        character*3 str
63        character*10 ctrac        character*10 ctrac
64        integer iq        integer iq
65        real rlong(iip1,jjp1), rlat(iip1,jjp1),rl(1,1)        real rlong(iip1,jjp1), rlat(iip1,jjp1)
66        integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid        integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid
67        integer ii,jj        integer ii,jj
68        integer zan, idayref        integer zan, idayref
# Line 90  C Line 90  C
90          enddo          enddo
91        enddo        enddo
92    
93        call histbeg_totreg(infile, iip1, rlong(:,1), jjp1, rlat(1,:),        call histbeg_totreg(infile, rlong(:,1), rlat(1,:),
94       .             1, iip1, 1, jjp1,       .             1, iip1, 1, jjp1,
95       .             tau0, zjulian, tstep, uhoriid, fileid)       .             tau0, zjulian, tstep, uhoriid, fileid)
96  C  C
# Line 106  C  un meme fichier) Line 106  C  un meme fichier)
106          enddo          enddo
107        enddo        enddo
108    
109        call histbeg_totreg('fluxstokev.nc', iip1, rlong(:,1), jjm,        call histbeg_totreg('fluxstokev.nc', rlong(:,1),
110       .             rlat(1,:),1, iip1, 1, jjm,       .             rlat(1,:jjm),1, iip1, 1, jjm,
111       .             tau0, zjulian, tstep, vhoriid, filevid)       .             tau0, zjulian, tstep, vhoriid, filevid)
112                    
113          rl(1,1) = 1.            call histbeg_totreg('defstoke.nc', (/1./), (/1./),
       call histbeg_regular('defstoke.nc', 1, rl, 1, rl,  
114       .             1, 1, 1, 1,       .             1, 1, 1, 1,
115       .             tau0, zjulian, tstep, dhoriid, filedid)       .             tau0, zjulian, tstep, dhoriid, filedid)
116    
# Line 125  C Line 124  C
124          enddo          enddo
125        enddo        enddo
126    
127        call histhori(fileid, iip1, rlong, jjp1, rlat, 'scalar',        call histhori_regular(fileid, iip1, rlong, jjp1, rlat, 'scalar',
128       .              'Grille points scalaires', thoriid)       .              'Grille points scalaires', thoriid)
129                    
130  C  C
# Line 148  C Line 147  C
147  C  Appels a histdef pour la definition des variables a sauvegarder  C  Appels a histdef pour la definition des variables a sauvegarder
148                    
149          CALL histdef(fileid, "phis", "Surface geop. height", "-",          CALL histdef(fileid, "phis", "Surface geop. height", "-",
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(fileid, "aire", "Grid area", "-",           CALL histdef(fileid, "aire", "Grid area", "-",
154       .                iip1,jjp1,thoriid, 1,1,1, -99, 32,       .                iip1,jjp1,thoriid, 1,1,1, -99,
155       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
156                    
157          CALL histdef(filedid, "dtvr", "tps dyn", "s",          CALL histdef(filedid, "dtvr", "tps dyn", "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, "istdyn", "tps stock", "s",           CALL histdef(filedid, "istdyn", "tps stock", "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           CALL histdef(filedid, "istphy", "tps stock phy", "s",           CALL histdef(filedid, "istphy", "tps stock phy", "s",
166       .                1,1,dhoriid, 1,1,1, -99, 32,       .                1,1,dhoriid, 1,1,1, -99,
167       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
168    
169    
# Line 173  C Masse Line 172  C Masse
172  C  C
173        call histdef(fileid, 'masse', 'Masse', 'kg',        call histdef(fileid, 'masse', 'Masse', 'kg',
174       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
175       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
176  C  C
177  C  Pbaru  C  Pbaru
178  C  C
179        call histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s',        call histdef(fileid, 'pbaru', 'flx de masse zonal', 'kg m/s',
180       .             iip1, jjp1, uhoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, uhoriid, llm, 1, llm, zvertiid,
181       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
182    
183  C  C
184  C  Pbarv  C  Pbarv
185  C  C
186        call histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s',        call histdef(filevid, 'pbarv', 'flx de masse mer', 'kg m/s',
187       .             iip1, jjm, vhoriid, llm, 1, llm, zvertiid,       .             iip1, jjm, vhoriid, llm, 1, llm, zvertiid,
188       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
189  C  C
190  C  w  C  w
191  C  C
192        call histdef(fileid, 'w', 'flx de masse vert', 'kg m/s',        call histdef(fileid, 'w', 'flx de masse vert', 'kg m/s',
193       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
194       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
195    
196  C  C
197  C  Temperature potentielle  C  Temperature potentielle
198  C  C
199        call histdef(fileid, 'teta', 'temperature potentielle', '-',        call histdef(fileid, 'teta', 'temperature potentielle', '-',
200       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
201       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
202  C  C
203    
204  C  C
# Line 207  C Geopotentiel Line 206  C Geopotentiel
206  C  C
207        call histdef(fileid, 'phi', 'geopotentiel instantane', '-',        call histdef(fileid, 'phi', 'geopotentiel instantane', '-',
208       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,       .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
209       .             32, 'inst(X)', t_ops, t_wrt)       .             'inst(X)', t_ops, t_wrt)
210  C  C
211  C  Fin  C  Fin
212  C  C
Legend:
Removed from v.14  
changed lines
  Added in v.15
  ViewVC Help
Powered by ViewVC 1.1.21