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) |