/[lmdze]/trunk/libf/phylmd/initphysto.f
ViewVC logotype

Diff of /trunk/libf/phylmd/initphysto.f

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

revision 14 by guez, Mon Jul 21 16:05:07 2008 UTC revision 15 by guez, Fri Aug 1 15:24:12 2008 UTC
# Line 97  C Line 97  C
97           CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat)           CALL gr_fi_ecrit(1,klon,iim,jjm+1,rlat,zx_lat)
98    
99    
100        call histbeg_totreg(infile, iim, zx_lon(:,1), jjm+1, zx_lat(1,:),        call histbeg_totreg(infile, zx_lon(:,1), zx_lat(1,:),
101       .             1, iim, 1, jjm+1,       .             1, iim, 1, jjm+1,
102       .             tau0, zjulian, tstep, nhoriid, fileid)       .             tau0, zjulian, tstep, nhoriid, fileid)
103    
# Line 118  C Line 118  C
118          write(*,*) 'apres histvert ds initphysto'          write(*,*) 'apres histvert ds initphysto'
119    
120         CALL histdef(fileid, "phis", "Surface geop. height", "-",         CALL histdef(fileid, "phis", "Surface geop. height", "-",
121       .                iim,jjm+1,nhoriid, 1,1,1, -99, 32,       .                iim,jjm+1,nhoriid, 1,1,1, -99,
122       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
123  c  c
124          write(*,*) 'apres phis ds initphysto'          write(*,*) 'apres phis ds initphysto'
125    
126           CALL histdef(fileid, "aire", "Grid area", "-",           CALL histdef(fileid, "aire", "Grid area", "-",
127       .                iim,jjm+1,nhoriid, 1,1,1, -99, 32,       .                iim,jjm+1,nhoriid, 1,1,1, -99,
128       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
129           write(*,*) 'apres aire ds initphysto'           write(*,*) 'apres aire ds initphysto'
130    
131           CALL histdef(fileid, "dtime", "tps phys ", "s",           CALL histdef(fileid, "dtime", "tps phys ", "s",
132       .                1,1,nhoriid, 1,1,1, -99, 32,       .                1,1,nhoriid, 1,1,1, -99,
133       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
134                    
135           CALL histdef(fileid, "istphy", "tps stock", "s",           CALL histdef(fileid, "istphy", "tps stock", "s",
136       .                1,1,nhoriid, 1,1,1, -99, 32,       .                1,1,nhoriid, 1,1,1, -99,
137       .                "once", t_ops, t_wrt)       .                "once", t_ops, t_wrt)
138    
139  C T  C T
140  C  C
141        call histdef(fileid, 't', 'Temperature', 'K',        call histdef(fileid, 't', 'Temperature', 'K',
142       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
143       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
144          write(*,*) 'apres t ds initphysto'          write(*,*) 'apres t ds initphysto'
145  C mfu  C mfu
146  C  C
147        call histdef(fileid, 'mfu', 'flx m. pan. mt', 'kg m/s',        call histdef(fileid, 'mfu', 'flx m. pan. mt', 'kg m/s',
148       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
149       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
150          write(*,*) 'apres mfu ds initphysto'          write(*,*) 'apres mfu ds initphysto'
151  C  C
152  C mfd  C mfd
153  C  C
154        call histdef(fileid, 'mfd', 'flx m. pan. des', 'kg m/s',        call histdef(fileid, 'mfd', 'flx m. pan. des', 'kg m/s',
155       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
156       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
157    
158  C  C
159  C en_u  C en_u
160  C  C
161        call histdef(fileid, 'en_u', 'flx ent pan mt', 'kg m/s',        call histdef(fileid, 'en_u', 'flx ent pan mt', 'kg m/s',
162       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
163       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
164                 write(*,*) 'apres en_u ds initphysto'                 write(*,*) 'apres en_u ds initphysto'
165  C  C
166  C de_u  C de_u
167  C  C
168        call histdef(fileid, 'de_u', 'flx det pan mt', 'kg m/s',        call histdef(fileid, 'de_u', 'flx det pan mt', 'kg m/s',
169       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
170       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
171    
172  C  C
173  C en_d  C en_d
174  C  C
175        call histdef(fileid, 'en_d', 'flx ent pan dt', 'kg m/s',        call histdef(fileid, 'en_d', 'flx ent pan dt', 'kg m/s',
176       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
177       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
178  C  C
179    
180  C  C
# Line 182  C de_d Line 182  C de_d
182  C  C
183        call histdef(fileid, 'de_d', 'flx det pan dt', 'kg m/s',        call histdef(fileid, 'de_d', 'flx det pan dt', 'kg m/s',
184       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
185       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
186    
187  c coefh frac_impa,frac_nucl  c coefh frac_impa,frac_nucl
188                    
189          call histdef(fileid, "coefh", " ", " ",          call histdef(fileid, "coefh", " ", " ",
190       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
191       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
192    
193  c abderrahmane le 16 09 02  c abderrahmane le 16 09 02
194          call histdef(fileid, "fm_th", " ", " ",          call histdef(fileid, "fm_th", " ", " ",
195       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
196       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
197    
198          call histdef(fileid, "en_th", " ", " ",          call histdef(fileid, "en_th", " ", " ",
199       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
200       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
201  c fin aj  c fin aj
202                    
203          write(*,*) 'apres coefh ds initphysto'            write(*,*) 'apres coefh ds initphysto'  
204    
205          call histdef(fileid, 'frac_impa', ' ', ' ',          call histdef(fileid, 'frac_impa', ' ', ' ',
206       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
207       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
208                    
209          call histdef(fileid, 'frac_nucl', ' ', ' ',          call histdef(fileid, 'frac_nucl', ' ', ' ',
210       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,       .             iim, jjm+1, nhoriid, llm, 1, llm, zvertiid,
211       .             32, 'inst(X)', t_ops, t_wrt)       .              'inst(X)', t_ops, t_wrt)
212    
213  c  c
214  c pyu1  c pyu1
215  c  c
216        CALL histdef(fileid, "pyu1", " ", " ",        CALL histdef(fileid, "pyu1", " ", " ",
217       .                iim,jjm+1,nhoriid, 1,1,1, -99, 32,       .                iim,jjm+1,nhoriid, 1,1,1, -99,
218       .                "inst(X)", t_ops, t_wrt)       .                "inst(X)", t_ops, t_wrt)
219    
220  c  c
221  c pyv1  c pyv1
222  c  c
223          CALL histdef(fileid, "pyv1", " ", " ",          CALL histdef(fileid, "pyv1", " ", " ",
224       .                iim,jjm+1,nhoriid, 1,1,1, -99, 32,       .                iim,jjm+1,nhoriid, 1,1,1, -99,
225       .                "inst(X)", t_ops, t_wrt)       .                "inst(X)", t_ops, t_wrt)
226                    
227          write(*,*) 'apres pyv1 ds initphysto'          write(*,*) 'apres pyv1 ds initphysto'
# Line 229  c Line 229  c
229  c ftsol1  c ftsol1
230  c  c
231          call histdef(fileid, "ftsol1", " ", " ",          call histdef(fileid, "ftsol1", " ", " ",
232       .             iim, jjm+1, nhoriid, 1, 1,1, -99,32,       .             iim, jjm+1, nhoriid, 1, 1,1, -99,
233       .             "inst(X)", t_ops, t_wrt)       .             "inst(X)", t_ops, t_wrt)
234    
235  c  c
236  c ftsol2  c ftsol2
237  c  c
238          call histdef(fileid, "ftsol2", " ", " ",          call histdef(fileid, "ftsol2", " ", " ",
239       .             iim, jjm+1, nhoriid, 1, 1,1, -99,32,       .             iim, jjm+1, nhoriid, 1, 1,1, -99,
240       .             "inst(X)", t_ops, t_wrt)       .             "inst(X)", t_ops, t_wrt)
241    
242  c  c
# Line 244  c ftsol3 Line 244  c ftsol3
244  c  c
245          call histdef(fileid, "ftsol3", " ", " ",          call histdef(fileid, "ftsol3", " ", " ",
246       .             iim, jjm+1, nhoriid, 1, 1,1, -99,       .             iim, jjm+1, nhoriid, 1, 1,1, -99,
247       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
248    
249  c  c
250  c ftsol4  c ftsol4
251  c  c
252          call histdef(fileid, "ftsol4", " ", " ",          call histdef(fileid, "ftsol4", " ", " ",
253       .             iim, jjm+1, nhoriid, 1, 1,1, -99,       .             iim, jjm+1, nhoriid, 1, 1,1, -99,
254       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
255                    
256  c  c
257  c rain  c rain
258  c  c
259          call histdef(fileid, "rain", " ", " ",          call histdef(fileid, "rain", " ", " ",
260       .             iim, jjm+1, nhoriid, 1, 1,1, -99,       .             iim, jjm+1, nhoriid, 1, 1,1, -99,
261       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
262    
263  c  c
264  c psrf1  c psrf1
265  c  c
266          call histdef(fileid, "psrf1", " ", " ",          call histdef(fileid, "psrf1", " ", " ",
267       .             iim, jjm+1, nhoriid, 1, 1, 1, -99,       .             iim, jjm+1, nhoriid, 1, 1, 1, -99,
268       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
269                    
270  c  c
271  c psrf2  c psrf2
272  c  c
273          call histdef(fileid, "psrf2", " ", " ",          call histdef(fileid, "psrf2", " ", " ",
274       .             iim, jjm+1, nhoriid, 1, 1, 1, -99,       .             iim, jjm+1, nhoriid, 1, 1, 1, -99,
275       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
276    
277  c  c
278  c psrf3  c psrf3
279  c  c
280          call histdef(fileid, "psrf3", " ", " ",          call histdef(fileid, "psrf3", " ", " ",
281       .             iim, jjm+1, nhoriid, 1, 1, 1, -99,       .             iim, jjm+1, nhoriid, 1, 1, 1, -99,
282       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
283    
284  c  c
285  c psrf4  c psrf4
286  c  c
287          call histdef(fileid, "psrf4", " ", " ",          call histdef(fileid, "psrf4", " ", " ",
288       .             iim, jjm+1, nhoriid, 1, 1, 1, -99,       .             iim, jjm+1, nhoriid, 1, 1, 1, -99,
289       .             32, "inst(X)", t_ops, t_wrt)       .              "inst(X)", t_ops, t_wrt)
290                    
291          write(*,*) 'avant histend ds initphysto'                  write(*,*) 'avant histend ds initphysto'        
292    
Legend:
Removed from v.14  
changed lines
  Added in v.15
  ViewVC Help
Powered by ViewVC 1.1.21