--- trunk/phylmd/physiq.f 2018/09/11 11:08:38 305 +++ trunk/phylmd/physiq.f 2018/09/11 12:23:47 306 @@ -1004,21 +1004,6 @@ CALL histwrite_phy("dtsvdfg", d_ts(:, is_lic)) CALL histwrite_phy("dtsvdfi", d_ts(:, is_sic)) CALL histwrite_phy("zxfqcalving", sum(fqcalving * pctsrf, dim = 2)) - - DO nsrf = 1, nbsrf - CALL histwrite_phy("pourc_"//clnsurf(nsrf), pctsrf(:, nsrf) * 100.) - CALL histwrite_phy("fract_"//clnsurf(nsrf), pctsrf(:, nsrf)) - CALL histwrite_phy("sens_"//clnsurf(nsrf), flux_t(:, nsrf)) - CALL histwrite_phy("lat_"//clnsurf(nsrf), fluxlat(:, nsrf)) - CALL histwrite_phy("tsol_"//clnsurf(nsrf), ftsol(:, nsrf)) - CALL histwrite_phy("taux_"//clnsurf(nsrf), flux_u(:, nsrf)) - CALL histwrite_phy("tauy_"//clnsurf(nsrf), flux_v(:, nsrf)) - CALL histwrite_phy("rugs_"//clnsurf(nsrf), frugs(:, nsrf)) - CALL histwrite_phy("albe_"//clnsurf(nsrf), falbe(:, nsrf)) - CALL histwrite_phy("u10m_"//clnsurf(nsrf), u10m_srf(:, nsrf)) - CALL histwrite_phy("v10m_"//clnsurf(nsrf), v10m_srf(:, nsrf)) - END DO - CALL histwrite_phy("albs", albsol) CALL histwrite_phy("tro3", wo * dobson_u * 1e3 / zmasse / rmo3 * md) CALL histwrite_phy("rugs", zxrugs) @@ -1029,12 +1014,6 @@ CALL histwrite_phy("s_oliqCL", s_oliqCL) CALL histwrite_phy("s_cteiCL", s_cteiCL) CALL histwrite_phy("s_therm", s_therm) - - if (conv_emanuel) then - CALL histwrite_phy("ptop", ema_pct) - CALL histwrite_phy("dnwd0", - mp) - end if - CALL histwrite_phy("temp", t_seri) CALL histwrite_phy("vitu", u_seri) CALL histwrite_phy("vitv", v_seri) @@ -1048,6 +1027,26 @@ CALL histwrite_phy("dtlw0", - cool0 / 86400.) CALL histwrite_phy("msnow", sum(fsnow * pctsrf, dim = 2)) call histwrite_phy("qsurf", sum(fqsurf * pctsrf, dim = 2)) + call histwrite_phy("flat", zxfluxlat) + + DO nsrf = 1, nbsrf + CALL histwrite_phy("pourc_"//clnsurf(nsrf), pctsrf(:, nsrf) * 100.) + CALL histwrite_phy("fract_"//clnsurf(nsrf), pctsrf(:, nsrf)) + CALL histwrite_phy("sens_"//clnsurf(nsrf), flux_t(:, nsrf)) + CALL histwrite_phy("lat_"//clnsurf(nsrf), fluxlat(:, nsrf)) + CALL histwrite_phy("tsol_"//clnsurf(nsrf), ftsol(:, nsrf)) + CALL histwrite_phy("taux_"//clnsurf(nsrf), flux_u(:, nsrf)) + CALL histwrite_phy("tauy_"//clnsurf(nsrf), flux_v(:, nsrf)) + CALL histwrite_phy("rugs_"//clnsurf(nsrf), frugs(:, nsrf)) + CALL histwrite_phy("albe_"//clnsurf(nsrf), falbe(:, nsrf)) + CALL histwrite_phy("u10m_"//clnsurf(nsrf), u10m_srf(:, nsrf)) + CALL histwrite_phy("v10m_"//clnsurf(nsrf), v10m_srf(:, nsrf)) + END DO + + if (conv_emanuel) then + CALL histwrite_phy("ptop", ema_pct) + CALL histwrite_phy("dnwd0", - mp) + end if if (ok_instan) call histsync(nid_ins)