/[lmdze]/trunk/phylmd/ini_histday.f
ViewVC logotype

Diff of /trunk/phylmd/ini_histday.f

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

trunk/libf/phylmd/ini_histday.f90 revision 61 by guez, Fri Apr 20 14:58:43 2012 UTC trunk/phylmd/ini_histday.f90 revision 76 by guez, Fri Nov 15 18:45:49 2013 UTC
# Line 18  contains Line 18  contains
18      use phyetat0_m, only: rlon, rlat      use phyetat0_m, only: rlon, rlat
19      use clesphys, only: ecrit_day      use clesphys, only: ecrit_day
20      use grid_change, only: gr_phy_write_2d      use grid_change, only: gr_phy_write_2d
21      use comvert, only: presnivs      use disvert_m, only: presnivs
22    
23      REAL, intent(in):: dtime ! pas temporel de la physique (s)      REAL, intent(in):: dtime ! pas temporel de la physique (s)
24      logical, intent(in):: ok_journe      logical, intent(in):: ok_journe
# Line 38  contains Line 38  contains
38         CALL histbeg_totreg("histday", rlon(2: iim+1), zx_lat(1, :), 1, iim, &         CALL histbeg_totreg("histday", rlon(2: iim+1), zx_lat(1, :), 1, iim, &
39              1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_day)              1, jjm + 1, itau_phy, zjulian, dtime, nhori, nid_day)
40         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", &         CALL histvert(nid_day, "presnivs", "Vertical levels", "mb", &
41              llm, presnivs/100., nvert)              presnivs/100., nvert)
42         if (nq <= 4) then         if (nq <= 4) then
43            call histdef(nid_day, "Sigma_O3_Royer", &            call histdef(nid_day, "Sigma_O3_Royer", &
44                 "column-density of ozone, in a cell, from Royer", "DU", &                 "column-density of ozone, in a cell, from Royer", "DU", &
45                 xsize=iim, ysize=jjm+1, horiid=nhori, pzsize=llm, par_oriz=1, &                 xsize=iim, ysize=jjm+1, horiid=nhori, pzsize=llm, oriz=1, &
46                 par_szz=llm, pzid=nvert, popp="ave(X)", pfreq_opp=dtime, &                 szz=llm, zid=nvert, opp="ave(X)", pfreq_opp=dtime, &
47                 pfreq_wrt=real(ecrit_day))                 pfreq_wrt=real(ecrit_day))
48         end if         end if
49         CALL histend(nid_day)         CALL histend(nid_day)
Legend:
Removed from v.61  
changed lines
  Added in v.76
  ViewVC Help
Powered by ViewVC 1.1.21