/[lmdze]/trunk/libf/phylmd/ini_hist.f90
ViewVC logotype

Diff of /trunk/libf/phylmd/ini_hist.f90

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

revision 14 by guez, Tue Mar 4 14:00:42 2008 UTC revision 15 by guez, Fri Aug 1 15:24:12 2008 UTC
# Line 37  contains Line 37  contains
37      ENDDO      ENDDO
38    
39      CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)      CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)
40      CALL histbeg_totreg("histhf", iim, zx_lon(:, 1), (jjm + 1), &      CALL histbeg_totreg("histhf", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &
41           zx_lat(1, :), 1, iim, 1, (jjm + 1), itau_phy, zjulian, dtime,  &           (jjm + 1), itau_phy, zjulian, dtime, nhori, nid_hf)
          nhori, nid_hf)  
42    
43      CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb", &      CALL histvert(nid_hf, "presnivs", "Vertical levels", "mb", &
44           llm, presnivs/100., nvert)           llm, presnivs/100., nvert)
# Line 88  contains Line 87  contains
87      ENDDO      ENDDO
88    
89      CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)      CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)
90      CALL histbeg_totreg("histhf3d", iim, zx_lon(:, 1), (jjm + 1), &      CALL histbeg_totreg("histhf3d", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &
91           zx_lat(1, :), 1, iim, 1, (jjm + 1), itau_phy, zjulian, dtime,  &           (jjm + 1), itau_phy, zjulian, dtime, nhori, nid_hf3d)
          nhori, nid_hf3d)  
92    
93      CALL histvert(nid_hf3d, "presnivs", "Vertical levels", "mb", &      CALL histvert(nid_hf3d, "presnivs", "Vertical levels", "mb", &
94           llm, presnivs/100., nvert)           llm, presnivs/100., nvert)
# Line 98  contains Line 96  contains
96      ! Champs 3D:      ! Champs 3D:
97    
98      CALL histdef(nid_hf3d, "temp", "Air temperature", "K", &      CALL histdef(nid_hf3d, "temp", "Air temperature", "K", &
99           iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &           iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
100           "ave(X)", zstohf, zout)           "ave(X)", zstohf, zout)
101    
102      CALL histdef(nid_hf3d, "ovap", "Specific humidity", "kg/kg", &      CALL histdef(nid_hf3d, "ovap", "Specific humidity", "kg/kg", &
103           iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &           iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
104           "ave(X)", zstohf, zout)           "ave(X)", zstohf, zout)
105    
106      CALL histdef(nid_hf3d, "vitu", "Zonal wind", "m/s", &      CALL histdef(nid_hf3d, "vitu", "Zonal wind", "m/s", &
107           iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &           iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
108           "ave(X)", zstohf, zout)           "ave(X)", zstohf, zout)
109    
110      CALL histdef(nid_hf3d, "vitv", "Meridional wind", "m/s", &      CALL histdef(nid_hf3d, "vitv", "Meridional wind", "m/s", &
111           iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &           iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
112           "ave(X)", zstohf, zout)           "ave(X)", zstohf, zout)
113    
114      if (nbtr >= 3) then      if (nbtr >= 3) then
115         CALL histdef(nid_hf3d, "O3", "Ozone mass fraction", "?", iim, &         CALL histdef(nid_hf3d, "O3", "Ozone mass fraction", "?", iim, &
116              (jjm + 1), nhori, llm, 1, llm, nvert, 32, "ave(X)", zstohf, &              (jjm + 1), nhori, llm, 1, llm, nvert, "ave(X)", zstohf, &
117              zout)              zout)
118      end if      end if
119    
# Line 132  contains Line 130  contains
130      use dimens_m, only: iim, jjm, llm      use dimens_m, only: iim, jjm, llm
131      use dimphy, only: klon      use dimphy, only: klon
132      use temps, only: itau_phy, day_ref, annee_ref      use temps, only: itau_phy, day_ref, annee_ref
133      USE ioipsl, only: ymds2ju, histbeg_totreg, histvert, histend      USE ioipsl, only: ymds2ju, histbeg_totreg, histvert, histend, histdef
134      use phyetat0_m, only: rlon, rlat      use phyetat0_m, only: rlon, rlat
135        use clesphys, only: ecrit_day
136    
137      REAL, intent(in):: dtime ! pas temporel de la physique (s)      REAL, intent(in):: dtime ! pas temporel de la physique (s)
138      real, intent(in):: presnivs(:)      real, intent(in):: presnivs(:)
# Line 156  contains Line 155  contains
155            zx_lon(i, jjm + 1) = rlon(i+1)            zx_lon(i, jjm + 1) = rlon(i+1)
156         ENDDO         ENDDO
157         CALL gr_fi_ecrit(1, klon, iim, jjm + 1, rlat, zx_lat)         CALL gr_fi_ecrit(1, klon, iim, jjm + 1, rlat, zx_lat)
158         CALL histbeg_totreg("histday", iim, zx_lon(:, 1), jjm + 1, &         CALL histbeg_totreg("histday", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &
159              zx_lat(1, :), 1, iim, 1, jjm + 1, itau_phy, zjulian, dtime,  &              jjm + 1, itau_phy, zjulian, dtime, nhori, nid_day)
             nhori, nid_day)  
160         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", &         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", &
161              llm, presnivs/100., nvert)              llm, presnivs/100., nvert)
162           call histdef(nid_day, "Sigma_O3_Royer", &
163                "column-density of ozone, in a cell, from Royer", "DU", &
164                pxsize=iim, pysize=jjm+1, phoriid=nhori, pzsize=llm, par_oriz=1, &
165                par_szz=llm, pzid=nvert, popp="ave(X)", pfreq_opp=dtime, &
166                pfreq_wrt=real(ecrit_day))
167         CALL histend(nid_day)         CALL histend(nid_day)
168      ENDIF                     ! fin de test sur ok_journe      ENDIF
169    
170    end subroutine ini_histday    end subroutine ini_histday
171    
# Line 206  contains Line 208  contains
208            zx_lon(i, (jjm + 1)) = rlon(i+1)            zx_lon(i, (jjm + 1)) = rlon(i+1)
209         ENDDO         ENDDO
210         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)         CALL gr_fi_ecrit(1, klon, iim, (jjm + 1), rlat, zx_lat)
211         CALL histbeg_totreg("histins", iim, zx_lon(:, 1), (jjm + 1), &         CALL histbeg_totreg("histins", zx_lon(:, 1), zx_lat(1, :), 1, iim, 1, &
212              zx_lat(1, :), 1, iim, 1, (jjm + 1), itau_phy, zjulian, dtime, &              jjm + 1, itau_phy, zjulian, dtime, nhori, nid_ins)
             nhori, nid_ins)  
213         write(*, *)'Inst ', itau_phy, zjulian         write(*, *)'Inst ', itau_phy, zjulian
214         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &         CALL histvert(nid_ins, "presnivs", "Vertical levels", "mb", &
215              llm, presnivs/100., nvert)              llm, presnivs/100., nvert)
216    
217         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &         CALL histdef(nid_ins, "phis", "Surface geop. height", "-", &
218              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
219              "once", zsto, zout)              "once", zsto, zout)
220    
221         CALL histdef(nid_ins, "aire", "Grid area", "-", &         CALL histdef(nid_ins, "aire", "Grid area", "-", &
222              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
223              "once", zsto, zout)              "once", zsto, zout)
224    
225         ! Champs 2D:         ! Champs 2D:
226    
227         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K", &         CALL histdef(nid_ins, "tsol", "Surface Temperature", "K", &
228              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
229              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
230    
231         CALL histdef(nid_ins, "t2m", "Temperature 2m", "K", &         CALL histdef(nid_ins, "t2m", "Temperature 2m", "K", &
232              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
233              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
234    
235         CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg", &         CALL histdef(nid_ins, "q2m", "Specific humidity 2m", "Kg/Kg", &
236              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
237              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
238    
239         CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s", &         CALL histdef(nid_ins, "u10m", "Vent zonal 10m", "m/s", &
240              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
241              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
242    
243         CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s", &         CALL histdef(nid_ins, "v10m", "Vent meridien 10m", "m/s", &
244              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
245              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
246    
247         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", &         CALL histdef(nid_ins, "psol", "Surface Pressure", "Pa", &
248              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
249              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
250    
251         CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day", &         CALL histdef(nid_ins, "plul", "Large-scale Precip.", "mm/day", &
252              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
253              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
254    
255         CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day", &         CALL histdef(nid_ins, "pluc", "Convective Precip.", "mm/day", &
256              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
257              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
258    
259         CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-", &         CALL histdef(nid_ins, "cdrm", "Momentum drag coef.", "-", &
260              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
261              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
262    
263         CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-", &         CALL histdef(nid_ins, "cdrh", "Heat drag coef.", "-", &
264              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
265              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
266    
267         CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",  &         CALL histdef(nid_ins, "precip", "Precipitation Totale liq+sol",  &
268              "kg/(s*m2)", &              "kg/(s*m2)", &
269              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
270              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
271    
272         CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)", &         CALL histdef(nid_ins, "snow", "Snow fall", "kg/(s*m2)", &
273              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
274              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
275    
276         !        CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",         !        CALL histdef(nid_ins, "snow_mass", "Snow Mass", "kg/m2",
277         !    .                iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,         !    .                iim, (jjm + 1), nhori, 1, 1, 1, -99,
278         !    .                "inst(X)", zsto, zout)         !    .                "inst(X)", zsto, zout)
279    
280         CALL histdef(nid_ins, "topl", "OLR", "W/m2", &         CALL histdef(nid_ins, "topl", "OLR", "W/m2", &
281              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
282              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
283    
284         CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)", &         CALL histdef(nid_ins, "evap", "Evaporation", "kg/(s*m2)", &
285              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
286              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
287    
288         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", &         CALL histdef(nid_ins, "sols", "Solar rad. at surf.", "W/m2", &
289              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
290              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
291    
292         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", &         CALL histdef(nid_ins, "soll", "IR rad. at surface", "W/m2", &
293              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
294              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
295    
296         CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",  &         CALL histdef(nid_ins, "solldown", "Down. IR rad. at surface",  &
297              "W/m2", iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              "W/m2", iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
298              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
299    
300         CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2", &         CALL histdef(nid_ins, "bils", "Surf. total heat flux", "W/m2", &
301              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
302              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
303    
304         CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2", &         CALL histdef(nid_ins, "sens", "Sensible heat flux", "W/m2", &
305              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
306              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
307    
308         CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2", &         CALL histdef(nid_ins, "fder", "Heat flux derivation", "W/m2", &
309              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
310              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
311    
312         CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s", &         CALL histdef(nid_ins, "dtsvdfo", "Boundary-layer dTs(o)", "K/s", &
313              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
314              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
315    
316         CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s", &         CALL histdef(nid_ins, "dtsvdft", "Boundary-layer dTs(t)", "K/s", &
317              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
318              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
319    
320         CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s", &         CALL histdef(nid_ins, "dtsvdfg", "Boundary-layer dTs(g)", "K/s", &
321              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
322              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
323    
324         CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", &         CALL histdef(nid_ins, "dtsvdfi", "Boundary-layer dTs(g)", "K/s", &
325              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
326              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
327    
328         DO nsrf = 1, nbsrf         DO nsrf = 1, nbsrf
329    
330            call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &            call histdef(nid_ins, "pourc_"//clnsurf(nsrf),  &
331                 "% "//clnsurf(nsrf), "%",   &                 "% "//clnsurf(nsrf), "%",   &
332                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
333                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
334    
335            call histdef(nid_ins, "fract_"//clnsurf(nsrf),  &            call histdef(nid_ins, "fract_"//clnsurf(nsrf),  &
336                 "Fraction "//clnsurf(nsrf), "1",   &                 "Fraction "//clnsurf(nsrf), "1",   &
337                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
338                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
339    
340            call histdef(nid_ins, "sens_"//clnsurf(nsrf),  &            call histdef(nid_ins, "sens_"//clnsurf(nsrf),  &
341                 "Sensible heat flux "//clnsurf(nsrf), "W/m2",   &                 "Sensible heat flux "//clnsurf(nsrf), "W/m2",   &
342                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
343                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
344    
345            call histdef(nid_ins, "tsol_"//clnsurf(nsrf),  &            call histdef(nid_ins, "tsol_"//clnsurf(nsrf),  &
346                 "Surface Temperature"//clnsurf(nsrf), "W/m2",   &                 "Surface Temperature"//clnsurf(nsrf), "W/m2",   &
347                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
348                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
349    
350            call histdef(nid_ins, "lat_"//clnsurf(nsrf),  &            call histdef(nid_ins, "lat_"//clnsurf(nsrf),  &
351                 "Latent heat flux "//clnsurf(nsrf), "W/m2",   &                 "Latent heat flux "//clnsurf(nsrf), "W/m2",   &
352                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
353                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
354    
355            call histdef(nid_ins, "taux_"//clnsurf(nsrf),  &            call histdef(nid_ins, "taux_"//clnsurf(nsrf),  &
356                 "Zonal wind stress"//clnsurf(nsrf), "Pa", &                 "Zonal wind stress"//clnsurf(nsrf), "Pa", &
357                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
358                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
359    
360            call histdef(nid_ins, "tauy_"//clnsurf(nsrf),  &            call histdef(nid_ins, "tauy_"//clnsurf(nsrf),  &
361                 "Meridional xind stress "//clnsurf(nsrf), "Pa",   &                 "Meridional xind stress "//clnsurf(nsrf), "Pa",   &
362                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
363                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
364    
365            call histdef(nid_ins, "albe_"//clnsurf(nsrf),  &            call histdef(nid_ins, "albe_"//clnsurf(nsrf),  &
366                 "Albedo "//clnsurf(nsrf), "-",   &                 "Albedo "//clnsurf(nsrf), "-",   &
367                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
368                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
369    
370            call histdef(nid_ins, "rugs_"//clnsurf(nsrf),  &            call histdef(nid_ins, "rugs_"//clnsurf(nsrf),  &
371                 "rugosite "//clnsurf(nsrf), "-",   &                 "rugosite "//clnsurf(nsrf), "-",   &
372                 iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &                 iim, (jjm + 1), nhori, 1, 1, 1, -99, &
373                 "inst(X)", zsto, zout)                 "inst(X)", zsto, zout)
374            !XXX            !XXX
375         END DO         END DO
376         CALL histdef(nid_ins, "rugs", "rugosity", "-", &         CALL histdef(nid_ins, "rugs", "rugosity", "-", &
377              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
378              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
379    
380         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &         CALL histdef(nid_ins, "albs", "Surface albedo", "-", &
381              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
382              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
383         CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-", &         CALL histdef(nid_ins, "albslw", "Surface albedo LW", "-", &
384              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32,  &              iim, (jjm + 1), nhori, 1, 1, 1, -99,  &
385              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
386    
387         !IM cf. AM 081204 BEG         !IM cf. AM 081204 BEG
388         ! HBTM2         ! HBTM2
389         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &         CALL histdef(nid_ins, "s_pblh", "Boundary Layer Height", "m", &
390              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
391              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
392    
393         CALL histdef(nid_ins, "s_pblt", "T at Boundary Layer Height",  &         CALL histdef(nid_ins, "s_pblt", "T at Boundary Layer Height",  &
394              "K", &              "K", &
395              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
396              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
397    
398         CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &         CALL histdef(nid_ins, "s_lcl", "Condensation level", "m", &
399              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
400              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
401    
402         CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for ABL", "J/m2", &         CALL histdef(nid_ins, "s_capCL", "Conv avlbl pot ener for ABL", "J/m2", &
403              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
404              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
405    
406         CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &         CALL histdef(nid_ins, "s_oliqCL", "Liq Water in BL", "kg/m2", &
407              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
408              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
409    
410         CALL histdef(nid_ins, "s_cteiCL", "Instability criteria (ABL)", "K", &         CALL histdef(nid_ins, "s_cteiCL", "Instability criteria (ABL)", "K", &
411              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
412              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
413    
414         CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &         CALL histdef(nid_ins, "s_therm", "Exces du thermique", "K", &
415              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
416              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
417    
418         CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2", &         CALL histdef(nid_ins, "s_trmb1", "deep_cape(HBTM2)", "J/m2", &
419              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
420              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
421    
422         CALL histdef(nid_ins, "s_trmb2", "inhibition (HBTM2)", "J/m2", &         CALL histdef(nid_ins, "s_trmb2", "inhibition (HBTM2)", "J/m2", &
423              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
424              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
425    
426         CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m", &         CALL histdef(nid_ins, "s_trmb3", "Point Omega (HBTM2)", "m", &
427              iim, (jjm + 1), nhori, 1, 1, 1, -99, 32, &              iim, (jjm + 1), nhori, 1, 1, 1, -99, &
428              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
429    
430         !IM cf. AM 081204 END         !IM cf. AM 081204 END
# Line 431  contains Line 432  contains
432         ! Champs 3D:         ! Champs 3D:
433    
434         CALL histdef(nid_ins, "temp", "Temperature", "K", &         CALL histdef(nid_ins, "temp", "Temperature", "K", &
435              iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
436              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
437    
438         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", &         CALL histdef(nid_ins, "vitu", "Zonal wind", "m/s", &
439              iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
440              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
441    
442         CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", &         CALL histdef(nid_ins, "vitv", "Merid wind", "m/s", &
443              iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
444              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
445    
446         CALL histdef(nid_ins, "geop", "Geopotential height", "m", &         CALL histdef(nid_ins, "geop", "Geopotential height", "m", &
447              iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
448              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
449    
450         CALL histdef(nid_ins, "pres", "Air pressure", "Pa", &         CALL histdef(nid_ins, "pres", "Air pressure", "Pa", &
451              iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
452              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
453    
454         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", &         CALL histdef(nid_ins, "dtvdf", "Boundary-layer dT", "K/s", &
455              iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
456              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
457    
458         CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", &         CALL histdef(nid_ins, "dqvdf", "Boundary-layer dQ", "Kg/Kg/s", &
459              iim, (jjm + 1), nhori, llm, 1, llm, nvert, 32, &              iim, (jjm + 1), nhori, llm, 1, llm, nvert, &
460              "inst(X)", zsto, zout)              "inst(X)", zsto, zout)
461    
462         CALL histend(nid_ins)         CALL histend(nid_ins)
Legend:
Removed from v.14  
changed lines
  Added in v.15
  ViewVC Help
Powered by ViewVC 1.1.21