GCC Code Coverage Report


Directory: ./
File: phys/phys_output_write_mod.f90
Date: 2022-01-11 19:19:34
Exec Total Coverage
Lines: 942 1198 78.6%
Branches: 1101 1446 76.1%

Line Branch Exec Source
1 !
2 ! $Id: phys_output_write_mod.F90 3956 2021-07-06 07:16:14Z jyg $
3 !
4 MODULE phys_output_write_mod
5
6 USE phytrac_mod, ONLY : d_tr_cl, d_tr_th, d_tr_cv, d_tr_lessi_impa, &
7 d_tr_lessi_nucl, d_tr_insc, d_tr_bcscav, d_tr_evapls, d_tr_ls, &
8 d_tr_trsp, d_tr_sscav, d_tr_sat, d_tr_uscav, flux_tr_dry
9
10 ! Author: Abderrahmane IDELKADI (original include file)
11 ! Author: Laurent FAIRHEAD (transformation to module/subroutine)
12 ! Author: Ulysse GERARD (effective implementation)
13
14 CONTAINS
15
16 ! ug Routine pour définir (lors du premier passage) ET sortir les variables
17 23568 SUBROUTINE phys_output_write(itap, pdtphys, paprs, pphis, &
18 481 pplay, lmax_th, aerosol_couple, &
19 ok_ade, ok_aie, ok_volcan, ivap, iliq, isol, ok_sync, &
20 ptconv, read_climoz, clevSTD, ptconvth, &
21 481 d_u, d_t, qx, d_qx, zmasse, flag_aerosol, flag_aerosol_strat, ok_cdnc)
22
23 ! This subroutine does the actual writing of diagnostics that were
24 ! defined and initialised in phys_output_mod.F90
25
26 USE dimphy, ONLY: klon, klev, klevp1
27 USE infotrac_phy, ONLY: nbtr
28 USE mod_phys_lmdz_para, ONLY: is_north_pole_phy,is_south_pole_phy
29 USE mod_grid_phy_lmdz, ONLY : nbp_lon, nbp_lat
30 USE time_phylmdz_mod, ONLY: day_step_phy, start_time, itau_phy
31 USE vertical_layers_mod, ONLY : ap, bp, aps, bps
32 USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, &
33 o_longitude, o_latitude, &
34 o_Ahyb, o_Bhyb,o_Ahyb_bounds, o_Bhyb_bounds, &
35 o_Ahyb_mid, o_Bhyb_mid,o_Ahyb_mid_bounds, o_Bhyb_mid_bounds, &
36 is_ave, is_sic, o_contfracATM, o_contfracOR, &
37 o_aireTER, o_flat, o_slp, o_ptstar, o_pt0, o_tsol, &
38 o_t2m, o_t2m_min, o_t2m_max, &
39 o_t2m_min_mon, o_t2m_max_mon, &
40 o_nt2mout, o_nt2moutfg, &
41 o_nq2mout, o_nq2moutfg, &
42 o_nu2mout, o_nu2moutfg, &
43 o_q2m, o_ustar, o_u10m, o_v10m, &
44 o_wind10m, o_wind10max, o_wind100m, o_gusts, o_sicf, &
45 o_loadfactor_wind_onshore, o_loadfactor_wind_offshore, &
46 o_psol, o_mass, o_qsurf, o_qsol, &
47 o_precip, o_rain_fall, o_rain_con, o_ndayrain, o_plul, o_pluc, o_plun, &
48 o_snow, o_msnow, o_fsnow, o_evap, o_ep,o_epmax_diag, & ! epmax_cape
49 o_tops, o_tops0, o_topl, o_topl0, &
50 o_SWupTOA, o_SWupTOAclr, o_SWupTOAcleanclr, o_SWdnTOA, o_fdiffSWdnSFC, &
51 o_SWdnTOAclr, o_nettop, o_SWup200, &
52 o_SWup200clr, o_SWdn200, o_SWdn200clr, &
53 o_LWup200, o_LWup200clr, o_LWdn200, &
54 o_LWdn200clr, o_sols, o_sols0, &
55 o_soll, o_radsol, o_soll0, o_SWupSFC, &
56 o_SWupSFCclr, o_SWupSFCcleanclr, o_SWdnSFC, o_SWdnSFCclr, o_SWdnSFCcleanclr, &
57 o_LWupSFC, o_LWdnSFC, o_LWupSFCclr, &
58 o_LWdnSFCclr, o_LWupTOAcleanclr, o_LWdnSFCcleanclr, o_bils, o_bils_diss, &
59 o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &
60 o_bils_latent, o_bils_enthalp, o_sens, &
61 o_fder, o_ffonte, o_fqcalving, o_fqfonte, o_mrroli, o_runofflic, &
62 o_taux, o_tauy, o_snowsrf, o_qsnow, &
63 o_snowhgt, o_toice, o_sissnow, o_runoff, &
64 o_albslw3, o_pourc_srf, o_fract_srf, &
65 o_taux_srf, o_tauy_srf, o_tsol_srf, &
66 o_evappot_srf, o_ustar_srf, o_u10m_srf, &
67 o_v10m_srf, o_t2m_srf, o_evap_srf, &
68 o_sens_srf, o_lat_srf, o_flw_srf, &
69 o_fsw_srf, o_wbils_srf, o_wbilo_srf, &
70 o_wevap_srf, o_wrain_srf, o_wsnow_srf, &
71 o_tke_srf, o_tke_max_srf,o_dltpbltke_srf, o_wstar, &
72 o_l_mixmin,o_l_mix, &
73 o_cdrm, o_cdrh, o_cldl, o_cldm, o_cldh, &
74 o_cldt, o_JrNt, o_cldljn, o_cldmjn, &
75 o_cldhjn, o_cldtjn, o_cldq, o_lwp, o_iwp, &
76 o_ue, o_ve, o_uq, o_vq, o_cape, o_pbase, &
77 o_uwat, o_vwat, &
78 o_ptop, o_fbase, o_plcl, o_plfc, &
79 o_wbeff, o_convoccur, o_cape_max, o_upwd, o_ep,o_epmax_diag, &
80 o_Mipsh, o_Ma, &
81 o_dnwd, o_dnwd0, o_ftime_deepcv, o_ftime_con, o_mc, &
82 o_prw, o_prlw, o_prsw, o_s_pblh, o_s_pblt, o_s_lcl, &
83 o_s_therm, o_uSTDlevs, o_vSTDlevs, &
84 o_wSTDlevs, o_zSTDlevs, o_qSTDlevs, &
85 o_tSTDlevs, epsfra, o_t_oce_sic, &
86 o_ale_bl, o_alp_bl, o_ale_wk, o_alp_wk, &
87 o_dtvdf_x , o_dtvdf_w , o_dqvdf_x , o_dqvdf_w , &
88 o_sens_x , o_sens_w , o_flat_x , o_flat_w , &
89 o_delta_tsurf, o_delta_tsurf_srf, &
90 o_cdragh_x , o_cdragh_w , o_cdragm_x , o_cdragm_w , &
91 o_kh , o_kh_x , o_kh_w , &
92 o_ale, o_alp, o_cin, o_WAPE, o_wake_h, o_cv_gen, o_wake_dens, &
93 o_wake_s, o_wake_deltat, o_wake_deltaq, &
94 o_wake_omg, o_dtwak, o_dqwak, o_dqwak2d, o_Vprecip, &
95 o_qtaa, o_Clwaa, &
96 o_ftd, o_fqd, o_wdtrainA, o_wdtrainS, o_wdtrainM, &
97 o_n2, o_s2, o_proba_notrig, &
98 o_random_notrig, o_ale_bl_stat, &
99 o_ale_bl_trig, o_alp_bl_det, &
100 o_alp_bl_fluct_m, o_alp_bl_fluct_tke, &
101 o_alp_bl_conv, o_alp_bl_stat, &
102 o_slab_qflux, o_tslab, o_slab_bils, &
103 o_slab_bilg, o_slab_sic, o_slab_tice, &
104 o_slab_hdiff, o_slab_ekman, o_slab_gm, &
105 o_weakinv, o_dthmin, o_cldtau, &
106 o_cldemi, o_pr_con_l, o_pr_con_i, &
107 o_pr_lsc_l, o_pr_lsc_i, o_re, o_fl, &
108 o_rh2m, &
109 o_qsat2m, o_tpot, o_tpote, o_SWnetOR, &
110 o_LWdownOR, o_snowl, &
111 o_solldown, o_dtsvdfo, o_dtsvdft, &
112 o_dtsvdfg, o_dtsvdfi, o_z0m, o_z0h, o_od443aer, o_od550aer, &
113 o_dryod550aer, o_od865aer, o_abs550aer, o_od550lt1aer, &
114 o_sconcso4, o_sconcno3, o_sconcoa, o_sconcbc, &
115 o_sconcss, o_sconcdust, o_concso4, o_concno3, &
116 o_concoa, o_concbc, o_concss, o_concdust, &
117 o_loadso4, o_loadoa, o_loadbc, o_loadss, &
118 o_loaddust, o_loadno3, o_tausumaero, &
119 o_drytausumaero, o_tausumaero_lw, &
120 o_topswad, o_topswad0, o_solswad, o_solswad0, &
121 o_toplwad, o_toplwad0, o_sollwad, o_sollwad0, &
122 o_swtoaas_nat, o_swsrfas_nat, &
123 o_swtoacs_nat, o_swtoaas_ant, &
124 o_swsrfas_ant, o_swtoacs_ant, &
125 o_swsrfcs_ant, o_swtoacf_nat, &
126 o_swsrfcf_nat, o_swtoacf_ant, &
127 o_swsrfcs_nat, o_swsrfcf_ant, &
128 o_swtoacf_zero, o_swsrfcf_zero, &
129 o_topswai, o_solswai, o_toplwai, o_sollwai, o_scdnc, &
130 o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&
131 o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
132 o_lcc3dstra, o_icc3dcon, o_icc3dstra, &
133 o_cldicemxrat, o_cldwatmxrat, o_reffclwtop, o_ec550aer, &
134 o_lwcon, o_iwcon, o_temp, o_theta, &
135 o_ovapinit, o_ovap, o_oliq, o_ocond, o_geop, &
136 o_vitu, o_vitv, o_vitw, o_pres, o_paprs, &
137 o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, &
138 o_rnebls, o_rneblsvol, o_rhum, o_rhl, o_rhi, o_ozone, o_ozone_light, &
139 o_duphy, o_dtphy, o_dqphy, o_dqphy2d, o_dqlphy, o_dqlphy2d, &
140 o_dqsphy, o_dqsphy2d, o_albe_srf, o_z0m_srf, o_z0h_srf, &
141 o_ages_srf, o_snow_srf, o_alb1, o_alb2, o_tke, o_tke_dissip, &
142 o_tke_max, o_kz, o_kz_max, o_clwcon, &
143 o_dtdyn, o_dqdyn, o_dqdyn2d, o_dqldyn, o_dqldyn2d, &
144 o_dqsdyn, o_dqsdyn2d, o_dudyn, o_dvdyn, &
145 o_dtcon, o_tntc, o_ducon, o_dvcon, &
146 o_dqcon, o_dqcon2d, o_tnhusc, o_tnhusc, o_dtlsc, &
147 o_dtlschr, o_dqlsc, o_dqlsc2d, o_beta_prec, &
148 o_dtlscth, o_dtlscst, o_dqlscth, o_dqlscth2d, &
149 o_dqlscst, o_dqlscst2d, o_plulth, o_plulst, &
150 o_ptconvth, o_lmaxth, o_dtvdf, &
151 o_dtdis, o_dqvdf, o_dqvdf2d, o_dteva, o_dqeva, o_dqeva2d, &
152 o_ptconv, o_ratqs, o_dtthe, &
153 o_duthe, o_dvthe, o_ftime_th, &
154 o_f_th, o_e_th, o_w_th, o_q_th, &
155 o_a_th, o_cloudth_sth, o_cloudth_senv, &
156 o_cloudth_sigmath, o_cloudth_sigmaenv, &
157 o_d_th, o_f0_th, o_zmax_th, &
158 o_dqthe, o_dqthe2d, o_dtajs, o_dqajs, o_dqajs2d, o_dtswr, &
159 o_dtsw0, o_dtlwr, o_dtlw0, o_dtec, &
160 o_duvdf, o_dvvdf, o_duoro, o_dvoro, &
161 o_dtoro, o_dulif, o_dvlif, o_dtlif, &
162 o_du_gwd_hines, o_dv_gwd_hines, o_dthin, o_dqch4, o_rsu, &
163 o_du_gwd_front, o_dv_gwd_front, &
164 o_east_gwstress, o_west_gwstress, &
165 o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, o_rsucsaf, o_rsdcsaf, &
166 o_rlucs, o_rldcs, o_tnt, o_tntr, &
167 o_tntscpbl, o_tnhus, o_tnhusscpbl, &
168 o_evu, o_h2o, o_mcd, o_dmc, o_ref_liq, &
169 o_ref_ice, o_rsut4co2, o_rlut4co2, &
170 o_rsutcs4co2, o_rlutcs4co2, o_rsu4co2, &
171 o_rlu4co2, o_rsucs4co2, o_rlucs4co2, &
172 o_rsd4co2, o_rld4co2, o_rsdcs4co2, &
173 o_rldcs4co2, o_tnondef, o_ta, o_zg, &
174 o_hus, o_hur, o_ua, o_va, o_wap, &
175 o_psbg, o_tro3, o_tro3_daylight, &
176 o_uxv, o_vxq, o_vxT, o_wxq, o_vxphi, &
177 o_wxT, o_uxu, o_vxv, o_TxT, o_trac, &
178 o_dtr_vdf, o_dtr_the, o_dtr_con, &
179 o_dtr_lessi_impa, o_dtr_lessi_nucl, &
180 o_dtr_insc, o_dtr_bcscav, o_dtr_evapls, &
181 o_dtr_ls, o_dtr_trsp, o_dtr_sscav, o_dtr_dry, &
182 o_dtr_sat, o_dtr_uscav, o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
183 o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &
184 o_ustr_gwd_front,o_vstr_gwd_front, &
185 o_sens_prec_liq_oce, o_sens_prec_liq_sic, &
186 o_sens_prec_sol_oce, o_sens_prec_sol_sic, &
187 o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
188 o_lat_prec_sol_oce, o_lat_prec_sol_sic, &
189 o_sza, &
190 ! Marine
191 o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &
192 o_map_deltaz_hc, o_map_pcld_hc, o_map_tcld_hc, &
193 o_map_emis_hist, o_map_iwp_hist, o_map_deltaz_hist, &
194 o_map_rad_hist, &
195 o_map_emis_Cb, o_map_pcld_Cb, o_map_tcld_Cb, &
196 o_map_emis_ThCi, o_map_pcld_ThCi, o_map_tcld_ThCi, &
197 o_map_emis_Anv, o_map_pcld_Anv, o_map_tcld_Anv, &
198 o_map_ntot, o_map_hc,o_map_hist,o_map_Cb,o_map_ThCi,o_map_Anv, &
199 o_alt_tropo, &
200 ! Tropopause
201 o_p_tropopause, o_z_tropopause, o_t_tropopause, &
202 o_col_O3_strato, o_col_O3_tropo, &
203 !--interactive CO2
204 o_flx_co2_ocean, o_flx_co2_ocean_cor, &
205 o_flx_co2_land, o_flx_co2_land_cor, &
206 o_flx_co2_ff, o_flx_co2_bb, &
207 o_delta_sst, o_delta_sal, o_ds_ns, o_dt_ns, o_dter, o_dser, o_tkt, &
208 o_tks, o_taur, o_sss
209
210
211 USE phys_output_ctrlout_mod, ONLY: o_heat_volc, o_cool_volc !NL
212 USE phys_state_var_mod, ONLY: heat_volc, cool_volc !NL
213
214 USE phys_state_var_mod, ONLY: pctsrf, rain_fall, snow_fall, &
215 qsol, z0m, z0h, fevap, agesno, &
216 nday_rain, rain_con, snow_con, &
217 topsw, toplw, toplw0, swup, swdn, solswfdiff, &
218 topsw0, swupc0, swdnc0, swup0, swdn0, SWup200, SWup200clr, &
219 SWdn200, SWdn200clr, LWup200, LWup200clr, &
220 LWdn200, LWdn200clr, solsw, solsw0, sollw, &
221 radsol, swradcorr, sollw0, sollwdown, sollw, gustiness, &
222 sollwdownclr, lwdnc0, lwdn0, ftsol, ustar, u10m, &
223 v10m, pbl_tke, wake_delta_pbl_TKE, &
224 delta_tsurf, &
225 wstar, cape, ema_pcb, ema_pct, &
226 ema_cbmf, Mipsh, Ma, fm_therm, ale_bl, alp_bl, ale, &
227 alp, cin, wake_pe, wake_dens, cv_gen, wake_s, wake_deltat, &
228 wake_deltaq, ftd, fqd, ale_bl_trig, albsol1, &
229 ale_wake, ale_bl_stat, &
230 rnebcon, wo, falb1, albsol2, coefh, clwcon0, &
231 ratqs, entr_therm, zqasc, detr_therm, f0, &
232 lwup, lwdn, lwupc0, lwup0, coefm, &
233 swupp, lwupp, swupc0p, swup0p, lwupc0p, lwup0p, swdnp, lwdnp, &
234 swdnc0p, swdn0p, lwdnc0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
235 vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
236 vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
237 T2sumSTD, nlevSTD, du_gwd_rando, du_gwd_front, &
238 ulevSTD, vlevSTD, wlevSTD, philevSTD, qlevSTD, tlevSTD, &
239 rhlevSTD, O3STD, O3daySTD, uvSTD, vqSTD, vTSTD, wqSTD, vphiSTD, &
240 wTSTD, u2STD, v2STD, T2STD, missing_val_nf90, delta_sal, ds_ns, &
241 dt_ns, delta_sst
242
243 USE phys_local_var_mod, ONLY: zxfluxlat, slp, ptstar, pt0, zxtsol, zt2m, &
244 zn2mout, t2m_min_mon, t2m_max_mon, evap, &
245 l_mixmin,l_mix, tke_dissip, &
246 zu10m, zv10m, zq2m, zustar, zxqsurf, &
247 rain_lsc, rain_num, snow_lsc, bils, sens, fder, &
248 zxffonte, zxfqcalving, zxfqfonte, zxrunofflic, fluxu, &
249 fluxv, zxsnow, qsnow, snowhgt, to_ice, &
250 sissnow, runoff, albsol3_lic, evap_pot, &
251 t2m, fluxt, fluxlat, fsollw, fsolsw, &
252 wfbils, wfbilo, wfevap, wfrain, wfsnow, &
253 cdragm, cdragh, cldl, cldm, &
254 cldh, cldt, JrNt, & ! only output names: cldljn,cldmjn,cldhjn,cldtjn
255 cldq, flwp, fiwp, ue, ve, uq, vq, &
256 uwat, vwat, &
257 plcl, plfc, wbeff, convoccur, upwd, dnwd, dnwd0, prw, prlw, prsw, &
258 s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
259 vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &
260 twriteSTD, alp_wake, &
261 !! dtvdf_x ,dtvdf_w ,dqvdf_x ,dqvdf_w , &
262 d_t_vdf_x ,d_t_vdf_w ,d_q_vdf_x ,d_q_vdf_w , &
263 sens_x ,sens_w ,zxfluxlat_x,zxfluxlat_w, &
264 cdragh_x ,cdragh_w ,cdragm_x ,cdragm_w , &
265 kh ,kh_x ,kh_w , &
266 wake_h, &
267 wake_omg, d_t_wake, d_q_wake, Vprecip, qtaa, Clw, &
268 wdtrainA, wdtrainS, wdtrainM, n2, s2, proba_notrig, &
269 random_notrig, &
270 alp_bl_det, alp_bl_fluct_m, alp_bl_conv, &
271 alp_bl_stat, alp_bl_fluct_tke, slab_wfbils, &
272 weak_inversion, dthmin, cldtau, cldemi, &
273 pmflxr, pmflxs, prfl, psfl, re, fl, rh2m, &
274 qsat2m, tpote, tpot, d_ts, od443aer, od550aer, dryod550aer, &
275 od865aer, abs550aer, od550lt1aer, sconcso4, sconcno3, &
276 sconcoa, sconcbc, sconcss, sconcdust, concso4, concno3, &
277 concoa, concbc, concss, concdust, loadso4, &
278 loadoa, loadbc, loadss, loaddust, loadno3, tausum_aero, drytausum_aero, &
279 topswad_aero, topswad0_aero, solswad_aero, &
280 solswad0_aero, topsw_aero, solsw_aero, &
281 topsw0_aero, solsw0_aero, topswcf_aero, &
282 solswcf_aero, topswai_aero, solswai_aero, &
283 toplwad_aero, toplwad0_aero, sollwad_aero, &
284 sollwad0_aero, toplwai_aero, sollwai_aero, &
285 scdnc, cldncl, reffclws, reffclwc, cldnvi, stratomask,&
286 lcc, lcc3d, lcc3dcon, lcc3dstra, &
287 icc3dcon, icc3dstra, zfice, reffclwtop, &
288 ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
289 ql_seri, qs_seri, tr_seri, &
290 zphi, u_seri, v_seri, omega, cldfra, &
291 rneb, rnebjn, rneblsvol, zx_rh, zx_rhl, zx_rhi, d_t_dyn, &
292 d_q_dyn, d_ql_dyn, d_qs_dyn, &
293 d_q_dyn2d, d_ql_dyn2d, d_qs_dyn2d, &
294 d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &
295 d_u_ajs, d_v_ajs, &
296 d_u_con, d_v_con, d_q_con, d_q_ajs, d_t_lsc, &
297 d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0, &
298 d_t_eva, d_q_lsc, beta_prec, d_t_lscth, &
299 d_t_lscst, d_q_lscth, d_q_lscst, plul_th, &
300 plul_st, d_t_vdf, d_t_diss, d_q_vdf, d_q_eva, &
301 zw2, fraca, zmax_th, d_q_ajsb, d_t_ec, d_u_vdf, &
302 d_v_vdf, d_u_oro, d_v_oro, d_t_oro, d_u_lif, &
303 d_v_lif, d_t_lif, du_gwd_hines, dv_gwd_hines, d_t_hin, &
304 dv_gwd_rando, dv_gwd_front, &
305 east_gwstress, west_gwstress, &
306 d_q_ch4, pmfd, pmfu, ref_liq, ref_ice, rhwriteSTD, &
307 ep, epmax_diag, & ! epmax_cape
308 p_tropopause, t_tropopause, z_tropopause
309
310
311 USE carbon_cycle_mod, ONLY: fco2_ff, fco2_bb, fco2_land, fco2_ocean
312 USE carbon_cycle_mod, ONLY: fco2_ocean_cor, fco2_land_cor
313
314 USE phys_output_var_mod, ONLY: vars_defined, snow_o, zfra_o, bils_diss, &
315 bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
316 itau_con, nfiles, clef_files, nid_files, dryaod_diag, &
317 zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
318 zustr_gwd_front, zvstr_gwd_front, sza_o, &
319 sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
320 cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
321 ! Marine
322 map_prop_hc, map_prop_hist, &
323 map_emis_hc,map_iwp_hc,map_deltaz_hc,&
324 map_pcld_hc,map_tcld_hc,&
325 map_emis_hist,map_iwp_hist,map_deltaz_hist,&
326 map_rad_hist,&
327 map_ntot,map_hc,map_hist,&
328 map_Cb,map_ThCi,map_Anv,&
329 map_emis_Cb,map_pcld_Cb,map_tcld_Cb,&
330 map_emis_ThCi,map_pcld_ThCi,map_tcld_ThCi,&
331 map_emis_Anv,map_pcld_Anv,map_tcld_Anv, &
332 alt_tropo, &
333 !Ionela
334 ok_4xCO2atm, dter, dser, tkt, tks, taur, sss
335
336 USE ocean_slab_mod, ONLY: nslay, tslab, slab_bilg, tice, seaice, &
337 slab_ekman,slab_hdiff,slab_gm,dt_ekman, dt_hdiff, dt_gm, dt_qflux
338 USE pbl_surface_mod, ONLY: snow
339 USE indice_sol_mod, ONLY: nbsrf
340 USE infotrac_phy, ONLY: nqtot, nqo, type_trac, tname, niadv
341 USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg
342 USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, landice_opt
343 USE aero_mod, ONLY: naero_tot, id_STRAT_phy
344 USE ioipsl, ONLY: histend, histsync
345 USE iophy, ONLY: set_itau_iophy, histwrite_phy
346 USE netcdf, ONLY: nf90_fill_real
347 USE print_control_mod, ONLY: prt_level,lunout
348
349
350 USE phys_cal_mod, ONLY : mth_len
351
352 USE YOESW, ONLY : RSUN
353 USE tracinca_mod, ONLY: config_inca
354 use config_ocean_skin_m, only: activate_ocean_skin
355
356 USE vertical_layers_mod, ONLY: presnivs
357
358 IMPLICIT NONE
359
360 INCLUDE "clesphys.h"
361 INCLUDE "thermcell.h"
362 INCLUDE "compbl.h"
363 INCLUDE "YOMCST.h"
364
365 ! Input
366 INTEGER :: itap, ivap, iliq, isol, read_climoz
367 INTEGER, DIMENSION(klon) :: lmax_th
368 LOGICAL :: aerosol_couple, ok_sync
369 LOGICAL :: ok_ade, ok_aie, ok_volcan
370 LOGICAL, DIMENSION(klon, klev) :: ptconv, ptconvth
371 REAL :: pdtphys
372 CHARACTER (LEN=4), DIMENSION(nlevSTD) :: clevSTD
373 962 REAL, DIMENSION(klon,nlevSTD) :: zx_tmp_fi3d_STD
374 REAL, DIMENSION(klon) :: pphis
375 REAL, DIMENSION(klon, klev) :: pplay, d_u, d_t
376 REAL, DIMENSION(klon, klev+1) :: paprs
377 REAL, DIMENSION(klon,klev,nqtot) :: qx, d_qx
378 REAL, DIMENSION(klon, klev) :: zmasse
379 INTEGER :: flag_aerosol_strat
380 INTEGER :: flag_aerosol
381 LOGICAL :: ok_cdnc
382 REAL, DIMENSION(3) :: freq_moyNMC
383
384 ! Local
385 INTEGER :: itau_w
386 INTEGER :: i, iinit, iinitend=1, iff, iq, iiq, nsrf, k, ll, naero
387 962 REAL, DIMENSION (klon) :: zx_tmp_fi2d, zpt_conv2d, wind100m
388 962 REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv
389 962 REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1
390 REAL, DIMENSION (klon,NSW) :: zx_tmp_fi3dsp
391 CHARACTER (LEN=4) :: bb2
392 962 INTEGER, DIMENSION(nbp_lon*nbp_lat) :: ndex2d
393 962 INTEGER, DIMENSION(nbp_lon*nbp_lat*klev) :: ndex3d
394 REAL, PARAMETER :: dobson_u = 2.1415e-05 ! Dobson unit, in kg m-2
395 ! REAL, PARAMETER :: missing_val=nf90_fill_real
396 962 REAL, DIMENSION(klev+1,2) :: Ahyb_bounds, Bhyb_bounds
397 962 REAL, DIMENSION(klev,2) :: Ahyb_mid_bounds, Bhyb_mid_bounds
398 INTEGER :: ilev
399 INTEGER, SAVE :: kmax_100m
400 !$OMP THREADPRIVATE(kmax_100m)
401 REAL :: x
402 REAL :: missing_val
403 REAL, PARAMETER :: un_jour=86400.
404 INTEGER ISW
405 CHARACTER*1 ch1
406 CHARACTER*20 varname
407
408 962 REAL,DIMENSION(klon,klev) :: z, dz
409 962 REAL,DIMENSION(klon) :: zrho, zt
410
411 INTEGER :: nqup
412
413 ! On calcul le nouveau tau:
414 481 itau_w = itau_phy + itap
415 ! On le donne à iophy pour que les histwrite y aient accès:
416 481 CALL set_itau_iophy(itau_w)
417
418 ! IF (.NOT.vars_defined) THEN
419 481 iinitend = 1
420 ! ELSE
421 ! iinitend = 1
422 ! ENDIF
423
424
425 missing_val=missing_val_nf90
426
427
2/2
✓ Branch 0 taken 1 times.
✓ Branch 1 taken 480 times.
481 IF (.NOT.vars_defined) THEN
428 1 kmax_100m=1
429
2/2
✓ Branch 0 taken 38 times.
✓ Branch 1 taken 1 times.
39 DO k=1, klev-1
430
2/2
✓ Branch 0 taken 3 times.
✓ Branch 1 taken 35 times.
39 IF (presnivs(k).GT.0.97*101325.) kmax_100m = k !--finding out max level for 100 m with a good margin
431 ENDDO
432 ENDIF
433
434 481 Ahyb_bounds(1,1) = 0.
435 481 Ahyb_bounds(1,2) = aps(1)
436 481 Bhyb_bounds(1,1) = 1.
437 481 Bhyb_bounds(1,2) = bps(1)
438
2/2
✓ Branch 0 taken 18278 times.
✓ Branch 1 taken 481 times.
18759 DO ilev=2,klev
439 18278 Ahyb_bounds(ilev,1) = aps(ilev-1)
440 18278 Ahyb_bounds(ilev,2) = aps(ilev)
441 18278 Bhyb_bounds(ilev,1) = bps(ilev-1)
442 18759 Bhyb_bounds(ilev,2) = bps(ilev)
443 ENDDO
444 481 Ahyb_bounds(klev+1,1) = aps(klev)
445 481 Ahyb_bounds(klev+1,2) = 0.
446 481 Bhyb_bounds(klev+1,1) = bps(klev)
447 481 Bhyb_bounds(klev+1,2) = 0.
448
449
2/2
✓ Branch 0 taken 18759 times.
✓ Branch 1 taken 481 times.
19240 DO ilev=1, klev
450 18759 Ahyb_mid_bounds(ilev,1) = ap(ilev)
451 18759 Ahyb_mid_bounds(ilev,2) = ap(ilev+1)
452 18759 Bhyb_mid_bounds(ilev,1) = bp(ilev)
453 19240 Bhyb_mid_bounds(ilev,2) = bp(ilev+1)
454 END DO
455
456 ! ug la boucle qui suit ne sert qu'une fois, pour l'initialisation, sinon il n'y a toujours qu'un seul passage:
457
2/2
✓ Branch 0 taken 481 times.
✓ Branch 1 taken 481 times.
962 DO iinit=1, iinitend
458 ! print *,'IFF iinit=', iinit, iinitend
459 ! On procède à l'écriture ou à la définition des nombreuses variables:
460 !!! Champs 1D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
461 481 CALL histwrite_phy(o_phis, pphis)
462
463
2/2
✓ Branch 0 taken 478114 times.
✓ Branch 1 taken 481 times.
478595 zx_tmp_fi2d = cell_area
464
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (is_north_pole_phy) then
465 481 zx_tmp_fi2d(1) = cell_area(1)/nbp_lon
466 ENDIF
467
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (is_south_pole_phy) then
468 481 zx_tmp_fi2d(klon) = cell_area(klon)/nbp_lon
469 ENDIf
470 481 CALL histwrite_phy(o_aire, zx_tmp_fi2d)
471
472
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
473
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
474 477600 zx_tmp_fi2d(i)=pctsrf(i,is_ter)+pctsrf(i,is_lic)
475 ENDDO
476 ENDIF
477
478 481 CALL histwrite_phy(o_contfracATM, zx_tmp_fi2d)
479 481 CALL histwrite_phy(o_contfracOR, pctsrf(:,is_ter))
480 !
481
482 !!! Champs 2D !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
483 ! Simulateur AIRS
484
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (ok_airs) then
485 CALL histwrite_phy(o_alt_tropo,alt_tropo)
486
487 CALL histwrite_phy(o_map_prop_hc,map_prop_hc)
488 CALL histwrite_phy(o_map_prop_hist,map_prop_hist)
489
490 CALL histwrite_phy(o_map_emis_hc,map_emis_hc)
491 CALL histwrite_phy(o_map_iwp_hc,map_iwp_hc)
492 CALL histwrite_phy(o_map_deltaz_hc,map_deltaz_hc)
493 CALL histwrite_phy(o_map_pcld_hc,map_pcld_hc)
494 CALL histwrite_phy(o_map_tcld_hc,map_tcld_hc)
495
496 CALL histwrite_phy(o_map_emis_hist,map_emis_hist)
497 CALL histwrite_phy(o_map_iwp_hist,map_iwp_hist)
498 CALL histwrite_phy(o_map_deltaz_hist,map_deltaz_hist)
499
500 CALL histwrite_phy(o_map_ntot,map_ntot)
501 CALL histwrite_phy(o_map_hc,map_hc)
502 CALL histwrite_phy(o_map_hist,map_hist)
503
504 CALL histwrite_phy(o_map_Cb,map_Cb)
505 CALL histwrite_phy(o_map_ThCi,map_ThCi)
506 CALL histwrite_phy(o_map_Anv,map_Anv)
507
508 CALL histwrite_phy(o_map_emis_Cb,map_emis_Cb)
509 CALL histwrite_phy(o_map_pcld_Cb,map_pcld_Cb)
510 CALL histwrite_phy(o_map_tcld_Cb,map_tcld_Cb)
511
512 CALL histwrite_phy(o_map_emis_ThCi,map_emis_ThCi)
513 CALL histwrite_phy(o_map_pcld_ThCi,map_pcld_ThCi)
514 CALL histwrite_phy(o_map_tcld_ThCi,map_tcld_ThCi)
515
516 CALL histwrite_phy(o_map_emis_Anv,map_emis_Anv)
517 CALL histwrite_phy(o_map_pcld_Anv,map_pcld_Anv)
518 CALL histwrite_phy(o_map_tcld_Anv,map_tcld_Anv)
519 ENDIF
520
521 481 CALL histwrite_phy(o_sza, sza_o)
522 481 CALL histwrite_phy(o_flat, zxfluxlat)
523 481 CALL histwrite_phy(o_ptstar, ptstar)
524 481 CALL histwrite_phy(o_pt0, pt0)
525 481 CALL histwrite_phy(o_slp, slp)
526 481 CALL histwrite_phy(o_tsol, zxtsol)
527 481 CALL histwrite_phy(o_t2m, zt2m)
528 481 CALL histwrite_phy(o_t2m_min, zt2m)
529 481 CALL histwrite_phy(o_t2m_max, zt2m)
530 481 CALL histwrite_phy(o_t2m_max_mon, t2m_max_mon)
531 481 CALL histwrite_phy(o_t2m_min_mon, t2m_min_mon)
532
533
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
534
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
535 477600 zx_tmp_fi2d(i)=real(zn2mout(i,1))
536 ENDDO
537 ENDIF
538 481 CALL histwrite_phy(o_nt2mout, zx_tmp_fi2d)
539
540
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
541
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
542 477600 zx_tmp_fi2d(i)=real(zn2mout(i,2))
543 ENDDO
544 ENDIF
545 481 CALL histwrite_phy(o_nt2moutfg, zx_tmp_fi2d)
546
547
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
548
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
549 477600 zx_tmp_fi2d(i)=real(zn2mout(i,3))
550 ENDDO
551 ENDIF
552 481 CALL histwrite_phy(o_nq2mout, zx_tmp_fi2d)
553
554
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
555
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
556 477600 zx_tmp_fi2d(i)=real(zn2mout(i,4))
557 ENDDO
558 ENDIF
559 481 CALL histwrite_phy(o_nq2moutfg, zx_tmp_fi2d)
560
561
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
562
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
563 477600 zx_tmp_fi2d(i)=real(zn2mout(i,5))
564 ENDDO
565 ENDIF
566 481 CALL histwrite_phy(o_nu2mout, zx_tmp_fi2d)
567
568
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
569
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
570 477600 zx_tmp_fi2d(i)=real(zn2mout(i,6))
571 ENDDO
572 ENDIF
573 481 CALL histwrite_phy(o_nu2moutfg, zx_tmp_fi2d)
574
575
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
576
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
577 477600 zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
578 ENDDO
579 ENDIF
580 481 CALL histwrite_phy(o_wind10m, zx_tmp_fi2d)
581
582
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
583
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
584 477600 zx_tmp_fi2d(i)=SQRT(zu10m(i)*zu10m(i)+zv10m(i)*zv10m(i))
585 ENDDO
586 ENDIF
587 481 CALL histwrite_phy(o_wind10max, zx_tmp_fi2d)
588
589 481 CALL histwrite_phy(o_gusts, gustiness)
590
591
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
592
2/2
✓ Branch 0 taken 1440 times.
✓ Branch 1 taken 480 times.
1920 DO k = 1, kmax_100m !--we could stop much lower
593
2/2
✓ Branch 0 taken 1431360 times.
✓ Branch 1 taken 1440 times.
1432800 zrho(:) = pplay(:,k)/t_seri(:,k)/RD ! air density in kg/m3
594
2/2
✓ Branch 0 taken 1431360 times.
✓ Branch 1 taken 1440 times.
1432800 dz(:,k) = (paprs(:,k)-paprs(:,k+1))/zrho(:)/RG ! layer thickness in m
595
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 960 times.
1920 IF (k==1) THEN
596
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 z(:,1) = (paprs(:,1)-pplay(:,1))/zrho(:)/RG ! altitude middle of first layer in m
597
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 477120 times.
477600 zt(:) = dz(:,1) ! altitude top of first layer in m
598 ELSE
599
2/2
✓ Branch 0 taken 954240 times.
✓ Branch 1 taken 960 times.
955200 z(:,k) = zt(:) + (paprs(:,k)-pplay(:,k))/zrho(:)/RG ! altitude middle of layer k in m
600
2/2
✓ Branch 0 taken 960 times.
✓ Branch 1 taken 954240 times.
955200 zt(:) = zt(:) + dz(:,k) ! altitude top of layer k in m
601 ENDIF
602 ENDDO
603
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 wind100m(:)=missing_val
604
2/2
✓ Branch 0 taken 960 times.
✓ Branch 1 taken 480 times.
1440 DO k=1, kmax_100m-1 !--we could stop much lower
605
2/2
✓ Branch 0 taken 954240 times.
✓ Branch 1 taken 960 times.
955680 DO i=1,klon
606
4/4
✓ Branch 0 taken 591328 times.
✓ Branch 1 taken 362912 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 114208 times.
955200 IF (z(i,k).LT.100..AND.z(i,k+1).GE.100.) THEN
607 wind100m(i)=SQRT( (u_seri(i,k)+(100.-z(i,k))/(z(i,k+1)-z(i,k))*(u_seri(i,k+1)-u_seri(i,k)))**2.0 + &
608 477120 (v_seri(i,k)+(100.-z(i,k))/(z(i,k+1)-z(i,k))*(v_seri(i,k+1)-v_seri(i,k)))**2.0 )
609 ENDIF
610 ENDDO
611 ENDDO
612 ENDIF
613 481 CALL histwrite_phy(o_wind100m, wind100m)
614
615
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
616 !--polynomial fit for 14,Vestas,1074,V136/3450 kW windmill - Olivier
617
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1,klon
618
3/4
✓ Branch 0 taken 186720 times.
✓ Branch 1 taken 290400 times.
✓ Branch 2 taken 186720 times.
✗ Branch 3 not taken.
477600 IF (pctsrf(i,is_ter).GT.0.05 .AND. wind100m(i).NE.missing_val) THEN
619 x=wind100m(i)
620
3/4
✓ Branch 0 taken 128617 times.
✓ Branch 1 taken 58103 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 128617 times.
186720 IF (x.LE.3.0 .OR. x.GE.22.5) THEN
621 58103 zx_tmp_fi2d(i)=0.0
622
2/2
✓ Branch 0 taken 1876 times.
✓ Branch 1 taken 126741 times.
128617 ELSE IF (x.GE.10.0) THEN
623 1876 zx_tmp_fi2d(i)=1.0
624 ELSE
625 zx_tmp_fi2d(i)= 10.73 + x*(-14.69 + x*(8.339 + x*(-2.59 + x*(0.4893 + x*(-0.05898 + x*(0.004627 + &
626 126741 x*(-0.0002352 + x*(7.478e-06 + x*(-1.351e-07 + x*(1.059e-09))))))))))
627 126741 zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
628 ENDIF
629 ELSE
630 290400 zx_tmp_fi2d(i)=missing_val
631 ENDIF
632 ENDDO
633 ENDIF
634 481 CALL histwrite_phy(o_loadfactor_wind_onshore, zx_tmp_fi2d)
635
636
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
637 !--polynomial fit for 14,Vestas,867,V164/8000 kW - Olivier
638
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1,klon
639
3/4
✓ Branch 0 taken 325441 times.
✓ Branch 1 taken 151679 times.
✓ Branch 2 taken 325441 times.
✗ Branch 3 not taken.
477600 IF (pctsrf(i,is_oce).GT.0.05 .AND. wind100m(i).NE.missing_val) THEN
640 x=wind100m(i)
641
3/4
✓ Branch 0 taken 255879 times.
✓ Branch 1 taken 69562 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 255879 times.
325441 IF (x.LE.3.0 .OR. x.GE.25.5) THEN
642 69562 zx_tmp_fi2d(i)=0.0
643
2/2
✓ Branch 0 taken 3378 times.
✓ Branch 1 taken 252501 times.
255879 ELSE IF (x.GE.12.5) THEN
644 3378 zx_tmp_fi2d(i)=1.0
645 ELSE
646 zx_tmp_fi2d(i)= 20.59 + x*(-22.39 + x*(10.25 + x*(-2.601 + x*(0.4065 + x*(-0.04099 + x*(0.002716 + &
647 252501 x*(-0.0001175 + x*(3.195e-06 + x*(-4.959e-08 + x*(3.352e-10))))))))))
648 252501 zx_tmp_fi2d(i)=MIN(MAX(zx_tmp_fi2d(i),0.0),1.0)
649 ENDIF
650 ELSE
651 151679 zx_tmp_fi2d(i)=missing_val
652 ENDIF
653 ENDDO
654 ENDIF
655 481 CALL histwrite_phy(o_loadfactor_wind_offshore, zx_tmp_fi2d)
656
657
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
658
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i = 1, klon
659 477600 zx_tmp_fi2d(i) = pctsrf(i,is_sic)
660 ENDDO
661 ENDIF
662 481 CALL histwrite_phy(o_sicf, zx_tmp_fi2d)
663 481 CALL histwrite_phy(o_q2m, zq2m)
664
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d = zustar
665 481 CALL histwrite_phy(o_ustar, zx_tmp_fi2d)
666 481 CALL histwrite_phy(o_u10m, zu10m)
667 481 CALL histwrite_phy(o_v10m, zv10m)
668
669
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
670
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i = 1, klon
671 477600 zx_tmp_fi2d(i) = paprs(i,1)
672 ENDDO
673 ENDIF
674 481 CALL histwrite_phy(o_psol, zx_tmp_fi2d)
675 481 CALL histwrite_phy(o_mass, zmasse)
676 481 CALL histwrite_phy(o_qsurf, zxqsurf)
677
678
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (.NOT. ok_veget) THEN
679 481 CALL histwrite_phy(o_qsol, qsol)
680 ENDIF
681
682
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
683
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i = 1, klon
684 477600 zx_tmp_fi2d(i) = rain_fall(i) + snow_fall(i)
685 ENDDO
686 ENDIF
687
688 481 CALL histwrite_phy(o_precip, zx_tmp_fi2d)
689 481 CALL histwrite_phy(o_rain_fall, rain_fall)
690 481 CALL histwrite_phy(o_ndayrain, nday_rain)
691
692 ! epmax_cape:
693 ! CALL histwrite_phy(o_epmax_diag, epmax_diag)
694 481 CALL histwrite_phy(o_ep, ep)
695
696
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
697
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i = 1, klon
698 477600 zx_tmp_fi2d(i) = rain_lsc(i) + snow_lsc(i)
699 ENDDO
700 ENDIF
701 481 CALL histwrite_phy(o_plul, zx_tmp_fi2d)
702 481 CALL histwrite_phy(o_plun, rain_num)
703
704
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
705
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i = 1, klon
706 477600 zx_tmp_fi2d(i) = rain_con(i) + snow_con(i)
707 ENDDO
708 ENDIF
709 481 CALL histwrite_phy(o_pluc, zx_tmp_fi2d)
710 481 CALL histwrite_phy(o_rain_con, rain_con)
711 481 CALL histwrite_phy(o_snow, snow_fall)
712 481 CALL histwrite_phy(o_msnow, zxsnow)
713 481 CALL histwrite_phy(o_fsnow, zfra_o)
714 481 CALL histwrite_phy(o_evap, evap)
715
716
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
717
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = topsw*swradcorr
718 ENDIF
719 481 CALL histwrite_phy(o_tops, zx_tmp_fi2d)
720
721
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
722
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = topsw0*swradcorr
723 ENDIF
724 481 CALL histwrite_phy(o_tops0, zx_tmp_fi2d)
725
726 481 CALL histwrite_phy(o_topl, toplw)
727 481 CALL histwrite_phy(o_topl0, toplw0)
728
729
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
730
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swup(:,klevp1)*swradcorr(:)
731 ENDIF
732 481 CALL histwrite_phy(o_SWupTOA, zx_tmp_fi2d)
733
734
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
735
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swup0(:,klevp1)*swradcorr(:)
736 ENDIF
737 481 CALL histwrite_phy(o_SWupTOAclr, zx_tmp_fi2d)
738
739
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
740
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swupc0(:,klevp1)*swradcorr(:)
741 ENDIF
742 481 CALL histwrite_phy(o_SWupTOAcleanclr, zx_tmp_fi2d)
743
744
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
745
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swdn(:,klevp1)*swradcorr(:)
746 ENDIF
747 481 CALL histwrite_phy(o_SWdnTOA, zx_tmp_fi2d)
748
749
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
750
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swdn0(:,klevp1)*swradcorr(:)
751 ENDIF
752 481 CALL histwrite_phy(o_SWdnTOAclr, zx_tmp_fi2d)
753
754
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
755
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = topsw(:)*swradcorr(:)-toplw(:)
756 ENDIF
757 481 CALL histwrite_phy(o_nettop, zx_tmp_fi2d)
758
759
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
760
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = SWup200*swradcorr
761 ENDIF
762 481 CALL histwrite_phy(o_SWup200, zx_tmp_fi2d)
763
764
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
765
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = SWup200clr*swradcorr
766 ENDIF
767 481 CALL histwrite_phy(o_SWup200clr, zx_tmp_fi2d)
768
769
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
770
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = SWdn200*swradcorr
771 ENDIF
772 481 CALL histwrite_phy(o_SWdn200, zx_tmp_fi2d)
773
774
775
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
776
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = SWdn200clr*swradcorr
777 ENDIF
778 481 CALL histwrite_phy(o_SWdn200clr, zx_tmp_fi2d)
779
780 481 CALL histwrite_phy(o_LWup200, LWup200)
781 481 CALL histwrite_phy(o_LWup200clr, LWup200clr)
782 481 CALL histwrite_phy(o_LWdn200, LWdn200)
783 481 CALL histwrite_phy(o_LWdn200clr, LWdn200clr)
784
785
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
786
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = solsw*swradcorr
787 ENDIF
788 481 CALL histwrite_phy(o_sols, zx_tmp_fi2d)
789
790
791
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
792
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = solsw0*swradcorr
793 ENDIF
794 481 CALL histwrite_phy(o_sols0, zx_tmp_fi2d)
795 481 CALL histwrite_phy(o_soll, sollw)
796 481 CALL histwrite_phy(o_soll0, sollw0)
797 481 CALL histwrite_phy(o_radsol, radsol)
798
799
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
800
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swup(:,1)*swradcorr(:)
801 ENDIF
802 481 CALL histwrite_phy(o_SWupSFC, zx_tmp_fi2d)
803
804
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
805
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swup0(:,1)*swradcorr(:)
806 ENDIF
807 481 CALL histwrite_phy(o_SWupSFCclr, zx_tmp_fi2d)
808
809
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
810
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swupc0(:,1)*swradcorr(:)
811 ENDIF
812 481 CALL histwrite_phy(o_SWupSFCcleanclr, zx_tmp_fi2d)
813
814
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
815
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swdn(:,1)*swradcorr(:)
816 ENDIF
817 481 CALL histwrite_phy(o_SWdnSFC, zx_tmp_fi2d)
818
819
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
820
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swdn0(:,1)*swradcorr(:)
821 ENDIF
822 481 CALL histwrite_phy(o_SWdnSFCclr, zx_tmp_fi2d)
823
824
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
825
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = swdnc0(:,1)*swradcorr(:)
826 ENDIF
827 481 CALL histwrite_phy(o_SWdnSFCcleanclr, zx_tmp_fi2d)
828
829 481 CALL histwrite_phy(o_fdiffSWdnSFC, solswfdiff)
830
831
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
832
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:)=sollwdown(:)-sollw(:)
833 ENDIF
834 481 CALL histwrite_phy(o_LWupSFC, zx_tmp_fi2d)
835 481 CALL histwrite_phy(o_LWdnSFC, sollwdown)
836
837
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
838
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 sollwdownclr(1:klon) = -1.*lwdn0(1:klon,1)
839
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(1:klon)=sollwdownclr(1:klon)-sollw0(1:klon)
840 ENDIF
841 481 CALL histwrite_phy(o_LWupSFCclr, zx_tmp_fi2d)
842 481 CALL histwrite_phy(o_LWdnSFCclr, sollwdownclr)
843
844
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
845
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = lwupc0(:,klevp1)
846 ENDIF
847 481 CALL histwrite_phy(o_LWupTOAcleanclr, zx_tmp_fi2d)
848
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
849
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = -1.*lwdnc0(:,1)
850 ENDIF
851 481 CALL histwrite_phy(o_LWdnSFCcleanclr, zx_tmp_fi2d)
852
853 481 CALL histwrite_phy(o_bils, bils)
854 481 CALL histwrite_phy(o_bils_diss, bils_diss)
855 481 CALL histwrite_phy(o_bils_ec, bils_ec)
856 481 CALL histwrite_phy(o_bils_ech, bils_ech)
857 481 CALL histwrite_phy(o_bils_tke, bils_tke)
858 481 CALL histwrite_phy(o_bils_kinetic, bils_kinetic)
859 481 CALL histwrite_phy(o_bils_latent, bils_latent)
860 481 CALL histwrite_phy(o_bils_enthalp, bils_enthalp)
861
862
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
863
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(1:klon)=-1*sens(1:klon)
864 ENDIF
865 481 CALL histwrite_phy(o_sens, zx_tmp_fi2d)
866 481 CALL histwrite_phy(o_fder, fder)
867 481 CALL histwrite_phy(o_ffonte, zxffonte)
868 481 CALL histwrite_phy(o_fqcalving, zxfqcalving)
869 481 CALL histwrite_phy(o_fqfonte, zxfqfonte)
870
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
871
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(1:klon)=(zxfqfonte(1:klon)+rain_fall(1:klon))*pctsrf(1:klon,is_lic)
872 ENDIF
873 481 CALL histwrite_phy(o_mrroli, zx_tmp_fi2d)
874 481 CALL histwrite_phy(o_runofflic, zxrunofflic)
875
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
876
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d=0.
877
2/2
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 480 times.
2400 DO nsrf=1,nbsrf
878
2/2
✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
1910880 zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxu(:,1,nsrf)
879 ENDDO
880 ENDIF
881 481 CALL histwrite_phy(o_taux, zx_tmp_fi2d)
882
883
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
884
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d=0.
885
2/2
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 480 times.
2400 DO nsrf=1,nbsrf
886
2/2
✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
1910880 zx_tmp_fi2d(:)=zx_tmp_fi2d(:)+pctsrf(:,nsrf)*fluxv(:,1,nsrf)
887 ENDDO
888 ENDIF
889 481 CALL histwrite_phy(o_tauy, zx_tmp_fi2d)
890
891 ! Etienne: test sorties pour compil sur JZ
892 ! IF (landice_opt .GE. 1) THEN
893 ! CALL histwrite_phy(o_snowsrf, snow_o)
894 ! CALL histwrite_phy(o_qsnow, qsnow)
895 ! CALL histwrite_phy(o_snowhgt,snowhgt)
896 ! CALL histwrite_phy(o_toice,to_ice)
897 ! CALL histwrite_phy(o_sissnow,sissnow)
898 ! CALL histwrite_phy(o_runoff,runoff)
899 ! CALL histwrite_phy(o_albslw3,albsol3_lic)
900 ! ENDIF
901
902
2/2
✓ Branch 0 taken 1924 times.
✓ Branch 1 taken 481 times.
2405 DO nsrf = 1, nbsrf
903
904
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)*100.
905 1924 CALL histwrite_phy(o_pourc_srf(nsrf), zx_tmp_fi2d)
906
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = pctsrf( 1 : klon, nsrf)
907 1924 CALL histwrite_phy(o_fract_srf(nsrf), zx_tmp_fi2d)
908
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxu( 1 : klon, 1, nsrf)
909 1924 CALL histwrite_phy(o_taux_srf(nsrf), zx_tmp_fi2d)
910
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxv( 1 : klon, 1, nsrf)
911 1924 CALL histwrite_phy(o_tauy_srf(nsrf), zx_tmp_fi2d)
912
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = ftsol( 1 : klon, nsrf)
913 1924 CALL histwrite_phy(o_tsol_srf(nsrf), zx_tmp_fi2d)
914
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = evap_pot( 1 : klon, nsrf)
915 1924 CALL histwrite_phy(o_evappot_srf(nsrf), zx_tmp_fi2d)
916
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = ustar(1 : klon, nsrf)
917 1924 CALL histwrite_phy(o_ustar_srf(nsrf), zx_tmp_fi2d)
918
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = u10m(1 : klon, nsrf)
919 1924 CALL histwrite_phy(o_u10m_srf(nsrf), zx_tmp_fi2d)
920
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = v10m(1 : klon, nsrf)
921 1924 CALL histwrite_phy(o_v10m_srf(nsrf), zx_tmp_fi2d)
922
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = t2m(1 : klon, nsrf)
923 1924 CALL histwrite_phy(o_t2m_srf(nsrf), zx_tmp_fi2d)
924
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = fevap(1 : klon, nsrf)
925 1924 CALL histwrite_phy(o_evap_srf(nsrf), zx_tmp_fi2d)
926
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxt( 1 : klon, 1, nsrf)
927 1924 CALL histwrite_phy(o_sens_srf(nsrf), zx_tmp_fi2d)
928
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = fluxlat( 1 : klon, nsrf)
929 1924 CALL histwrite_phy(o_lat_srf(nsrf), zx_tmp_fi2d)
930
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsollw( 1 : klon, nsrf)
931 1924 CALL histwrite_phy(o_flw_srf(nsrf), zx_tmp_fi2d)
932
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, nsrf)
933 1924 CALL histwrite_phy(o_fsw_srf(nsrf), zx_tmp_fi2d)
934
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfbils( 1 : klon, nsrf)
935 1924 CALL histwrite_phy(o_wbils_srf(nsrf), zx_tmp_fi2d)
936
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfbilo( 1 : klon, nsrf)
937 1924 CALL histwrite_phy(o_wbilo_srf(nsrf), zx_tmp_fi2d)
938
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfevap( 1 : klon, nsrf)
939 1924 CALL histwrite_phy(o_wevap_srf(nsrf), zx_tmp_fi2d)
940
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfrain( 1 : klon, nsrf)
941 1924 CALL histwrite_phy(o_wrain_srf(nsrf), zx_tmp_fi2d)
942
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = wfsnow( 1 : klon, nsrf)
943 1924 CALL histwrite_phy(o_wsnow_srf(nsrf), zx_tmp_fi2d)
944
945
1/2
✓ Branch 0 taken 1924 times.
✗ Branch 1 not taken.
1924 IF (iflag_pbl > 1) THEN
946 1924 CALL histwrite_phy(o_tke_srf(nsrf), pbl_tke(:,1:klev,nsrf))
947 !CALL histwrite_phy(o_l_mix(nsrf), l_mix(:,1:klev,nsrf))
948 1924 CALL histwrite_phy(o_l_mixmin(nsrf), l_mixmin(:,1:klev,nsrf))
949 1924 CALL histwrite_phy(o_tke_max_srf(nsrf), pbl_tke(:,1:klev,nsrf))
950
951
952 ENDIF
953 !jyg<
954
3/6
✓ Branch 0 taken 1924 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 1924 times.
✗ Branch 3 not taken.
✓ Branch 4 taken 1924 times.
✗ Branch 5 not taken.
2405 IF (iflag_pbl > 1 .AND. iflag_wake>=1 .AND. iflag_pbl_split >=1) THEN
955 1924 CALL histwrite_phy(o_dltpbltke_srf(nsrf), wake_delta_pbl_TKE(:,1:klev,nsrf))
956 ENDIF
957 !>jyg
958 ! IF (iflag_pbl > 1 .AND. ifl_pbltree >=1 ) THEN
959 ! CALL histwrite_phy(o_treedrg_srf(nsrf), treedrg(:,1:klev,nsrf))
960 ! ENDIF
961
962 ENDDO
963
964
965
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_pbl > 1) THEN
966
4/4
✓ Branch 0 taken 18759 times.
✓ Branch 1 taken 481 times.
✓ Branch 2 taken 18646446 times.
✓ Branch 3 taken 18759 times.
18665686 zx_tmp_fi3d=0.
967
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
968
2/2
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 480 times.
2400 DO nsrf=1,nbsrf
969
2/2
✓ Branch 0 taken 74880 times.
✓ Branch 1 taken 1920 times.
77280 DO k=1,klev
970 zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
971
2/2
✓ Branch 0 taken 74430720 times.
✓ Branch 1 taken 74880 times.
74507520 +pctsrf(:,nsrf)*tke_dissip(:,k,nsrf)
972 ENDDO
973 ENDDO
974 ENDIF
975
976 481 CALL histwrite_phy(o_tke_dissip, zx_tmp_fi3d)
977 ENDIF
978
979
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_liq_o(1 : klon, 1)
980 481 CALL histwrite_phy(o_sens_prec_liq_oce, zx_tmp_fi2d)
981
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_liq_o(1 : klon, 2)
982 481 CALL histwrite_phy(o_sens_prec_liq_sic, zx_tmp_fi2d)
983
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_sol_o(1 : klon, 1)
984 481 CALL histwrite_phy(o_sens_prec_sol_oce, zx_tmp_fi2d)
985
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_sol_o(1 : klon, 2)
986 481 CALL histwrite_phy(o_sens_prec_sol_sic, zx_tmp_fi2d)
987
988
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_liq_o(1 : klon, 1)
989 481 CALL histwrite_phy(o_lat_prec_liq_oce, zx_tmp_fi2d)
990
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_liq_o(1 : klon, 2)
991 481 CALL histwrite_phy(o_lat_prec_liq_sic, zx_tmp_fi2d)
992
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_sol_o(1 : klon, 1)
993 481 CALL histwrite_phy(o_lat_prec_sol_oce, zx_tmp_fi2d)
994
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_sol_o(1 : klon, 2)
995 481 CALL histwrite_phy(o_lat_prec_sol_sic, zx_tmp_fi2d)
996
997
2/2
✓ Branch 0 taken 2405 times.
✓ Branch 1 taken 481 times.
2886 DO nsrf=1,nbsrf+1
998 2886 CALL histwrite_phy(o_wstar(nsrf), wstar(1 : klon, nsrf))
999 ENDDO
1000
1001 481 CALL histwrite_phy(o_cdrm, cdragm)
1002 481 CALL histwrite_phy(o_cdrh, cdragh)
1003 481 CALL histwrite_phy(o_cldl, cldl)
1004 481 CALL histwrite_phy(o_cldm, cldm)
1005 481 CALL histwrite_phy(o_cldh, cldh)
1006 481 CALL histwrite_phy(o_cldt, cldt)
1007 481 CALL histwrite_phy(o_JrNt, JrNt)
1008
1009
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d=cldl*JrNt
1010 481 CALL histwrite_phy(o_cldljn, zx_tmp_fi2d)
1011
1012
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d=cldm*JrNt
1013 481 CALL histwrite_phy(o_cldmjn, zx_tmp_fi2d)
1014
1015
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d=cldh*JrNt
1016 481 CALL histwrite_phy(o_cldhjn, zx_tmp_fi2d)
1017
1018
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d=cldt*JrNt
1019 481 CALL histwrite_phy(o_cldtjn, zx_tmp_fi2d)
1020
1021 481 CALL histwrite_phy(o_cldq, cldq)
1022
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1:klon) = flwp(1:klon)
1023 481 CALL histwrite_phy(o_lwp, zx_tmp_fi2d)
1024
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1:klon) = fiwp(1:klon)
1025 481 CALL histwrite_phy(o_iwp, zx_tmp_fi2d)
1026 481 CALL histwrite_phy(o_ue, ue)
1027 481 CALL histwrite_phy(o_ve, ve)
1028 481 CALL histwrite_phy(o_uq, uq)
1029 481 CALL histwrite_phy(o_vq, vq)
1030 481 CALL histwrite_phy(o_uwat, uwat)
1031 481 CALL histwrite_phy(o_vwat, vwat)
1032
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_con.GE.3) THEN ! sb
1033 481 CALL histwrite_phy(o_cape, cape)
1034 481 CALL histwrite_phy(o_pbase, ema_pcb)
1035 481 CALL histwrite_phy(o_ptop, ema_pct)
1036 481 CALL histwrite_phy(o_fbase, ema_cbmf)
1037
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_con /= 30) THEN
1038 481 CALL histwrite_phy(o_plcl, plcl)
1039 481 CALL histwrite_phy(o_plfc, plfc)
1040 481 CALL histwrite_phy(o_wbeff, wbeff)
1041 481 CALL histwrite_phy(o_convoccur, convoccur)
1042 ENDIF
1043
1044 481 CALL histwrite_phy(o_cape_max, cape)
1045
1046 481 CALL histwrite_phy(o_upwd, upwd)
1047 481 CALL histwrite_phy(o_Ma, Ma)
1048 481 CALL histwrite_phy(o_dnwd, dnwd)
1049 481 CALL histwrite_phy(o_dnwd0, dnwd0)
1050 !! The part relative to the frequency of occurence of convection
1051 !! is now grouped with the part relative to thermals and shallow
1052 !! convection (output of the 3 fields: ftime_deepcv, ftime_th and
1053 !! ftime_con).
1054
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1055
1/2
✓ Branch 0 taken 480 times.
✗ Branch 1 not taken.
480 IF (iflag_thermals>=1)THEN
1056
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 zx_tmp_fi3d=-dnwd+dnwd0+upwd+fm_therm(:,1:klev)
1057 ELSE
1058 zx_tmp_fi3d=-dnwd+dnwd0+upwd
1059 ENDIF
1060 ENDIF
1061 481 CALL histwrite_phy(o_mc, zx_tmp_fi3d)
1062 ENDIF !iflag_con .GE. 3
1063 481 CALL histwrite_phy(o_prw, prw)
1064 481 CALL histwrite_phy(o_prlw, prlw)
1065 481 CALL histwrite_phy(o_prsw, prsw)
1066 481 CALL histwrite_phy(o_s_pblh, s_pblh)
1067 481 CALL histwrite_phy(o_s_pblt, s_pblt)
1068 481 CALL histwrite_phy(o_s_lcl, s_lcl)
1069 481 CALL histwrite_phy(o_s_therm, s_therm)
1070 !IM : Les champs suivants (s_capCL, s_oliqCL, s_cteiCL, s_trmb1, s_trmb2, s_trmb3) ne sont pas definis dans HBTM.F
1071 ! IF (o_s_capCL%flag(iff)<=lev_files(iff)) THEN
1072 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1073 ! $o_s_capCL%name,itau_w,s_capCL)
1074 ! ENDIF
1075 ! IF (o_s_oliqCL%flag(iff)<=lev_files(iff)) THEN
1076 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1077 ! $o_s_oliqCL%name,itau_w,s_oliqCL)
1078 ! ENDIF
1079 ! IF (o_s_cteiCL%flag(iff)<=lev_files(iff)) THEN
1080 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1081 ! $o_s_cteiCL%name,itau_w,s_cteiCL)
1082 ! ENDIF
1083 ! IF (o_s_trmb1%flag(iff)<=lev_files(iff)) THEN
1084 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1085 ! $o_s_trmb1%name,itau_w,s_trmb1)
1086 ! ENDIF
1087 ! IF (o_s_trmb2%flag(iff)<=lev_files(iff)) THEN
1088 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1089 ! $o_s_trmb2%name,itau_w,s_trmb2)
1090 ! ENDIF
1091 ! IF (o_s_trmb3%flag(iff)<=lev_files(iff)) THEN
1092 ! CALL histwrite_phy(nid_files(iff),clef_stations(iff),
1093 ! $o_s_trmb3%name,itau_w,s_trmb3)
1094 ! ENDIF
1095
1096
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (.NOT.ok_all_xml) THEN
1097 ! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX:
1098 ! Champs interpolles sur des niveaux de pression
1099
2/2
✓ Branch 0 taken 4810 times.
✓ Branch 1 taken 481 times.
5291 DO iff=1, nfiles
1100 ll=0
1101
2/2
✓ Branch 0 taken 81770 times.
✓ Branch 1 taken 4810 times.
87061 DO k=1, nlevSTD
1102 81770 bb2=clevSTD(k)
1103 IF (bb2.EQ."850".OR.bb2.EQ."700".OR. &
1104 bb2.EQ."500".OR.bb2.EQ."200".OR. &
1105 bb2.EQ."100".OR. &
1106
14/14
✓ Branch 0 taken 76960 times.
✓ Branch 1 taken 4810 times.
✓ Branch 2 taken 72150 times.
✓ Branch 3 taken 4810 times.
✓ Branch 4 taken 67340 times.
✓ Branch 5 taken 4810 times.
✓ Branch 6 taken 62530 times.
✓ Branch 7 taken 4810 times.
✓ Branch 8 taken 57720 times.
✓ Branch 9 taken 4810 times.
✓ Branch 10 taken 52910 times.
✓ Branch 11 taken 4810 times.
✓ Branch 12 taken 4810 times.
✓ Branch 13 taken 48100 times.
86580 bb2.EQ."50".OR.bb2.EQ."10") THEN
1107
1108 ! a refaire correctement !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1109 33670 ll=ll+1
1110 33670 CALL histwrite_phy(o_uSTDlevs(ll),uwriteSTD(:,k,iff), iff)
1111 33670 CALL histwrite_phy(o_vSTDlevs(ll),vwriteSTD(:,k,iff), iff)
1112 33670 CALL histwrite_phy(o_wSTDlevs(ll),wwriteSTD(:,k,iff), iff)
1113 33670 CALL histwrite_phy(o_zSTDlevs(ll),phiwriteSTD(:,k,iff), iff)
1114 33670 CALL histwrite_phy(o_qSTDlevs(ll),qwriteSTD(:,k,iff), iff)
1115 33670 CALL histwrite_phy(o_tSTDlevs(ll),twriteSTD(:,k,iff), iff)
1116
1117 ENDIF !(bb2.EQ."850".OR.bb2.EQ."700".OR.
1118 ENDDO
1119 ENDDO
1120 ENDIF
1121
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1122
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
1123
4/4
✓ Branch 0 taken 114143 times.
✓ Branch 1 taken 362977 times.
✓ Branch 2 taken 36383 times.
✓ Branch 3 taken 77760 times.
477120 IF (pctsrf(i,is_oce).GT.epsfra.OR. &
1124 480 pctsrf(i,is_sic).GT.epsfra) THEN
1125 zx_tmp_fi2d(i) = (ftsol(i, is_oce) * pctsrf(i,is_oce)+ &
1126 ftsol(i, is_sic) * pctsrf(i,is_sic))/ &
1127 399360 (pctsrf(i,is_oce)+pctsrf(i,is_sic))
1128 ELSE
1129 77760 zx_tmp_fi2d(i) = 273.15
1130 ENDIF
1131 ENDDO
1132 ENDIF
1133 481 CALL histwrite_phy(o_t_oce_sic, zx_tmp_fi2d)
1134
1135 ! Couplage convection-couche limite
1136
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_con.GE.3) THEN
1137
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_coupl>=1) THEN
1138 481 CALL histwrite_phy(o_ale_bl, ale_bl)
1139 481 CALL histwrite_phy(o_alp_bl, alp_bl)
1140 ENDIF !iflag_coupl>=1
1141 ENDIF !(iflag_con.GE.3)
1142 ! Wakes
1143
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_con.EQ.3) THEN
1144 481 CALL histwrite_phy(o_Mipsh, Mipsh)
1145
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_wake>=1) THEN
1146 481 CALL histwrite_phy(o_ale_wk, ale_wake)
1147 481 CALL histwrite_phy(o_alp_wk, alp_wake)
1148
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_pbl_split>=1) THEN
1149 !! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_x(1:klon,1:klev)/pdtphys
1150 !! CALL histwrite_phy(o_dtvdf_x ,zx_tmp_fi3d)
1151 !! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_w(1:klon,1:klev)/pdtphys
1152 !! CALL histwrite_phy(o_dtvdf_w ,zx_tmp_fi3d)
1153 !! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_x(1:klon,1:klev)/pdtphys
1154 !! CALL histwrite_phy(o_dqvdf_x ,zx_tmp_fi3d)
1155 !! IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_w(1:klon,1:klev)/pdtphys
1156 !
1157
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf_x(1:klon,1:klev)/pdtphys
1158 481 CALL histwrite_phy(o_dtvdf_x ,zx_tmp_fi3d)
1159
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf_w(1:klon,1:klev)/pdtphys
1160 481 CALL histwrite_phy(o_dtvdf_w ,zx_tmp_fi3d)
1161
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf_x(1:klon,1:klev)/pdtphys
1162 481 CALL histwrite_phy(o_dqvdf_x ,zx_tmp_fi3d)
1163
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf_w(1:klon,1:klev)/pdtphys
1164 !
1165 481 CALL histwrite_phy(o_dqvdf_w ,zx_tmp_fi3d)
1166
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1:klon)=-1*sens_x(1:klon)
1167 481 CALL histwrite_phy(o_sens_x ,zx_tmp_fi2d)
1168
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1:klon)=-1*sens_w(1:klon)
1169 481 CALL histwrite_phy(o_sens_w ,zx_tmp_fi2d)
1170 481 CALL histwrite_phy(o_flat_x ,zxfluxlat_x)
1171 481 CALL histwrite_phy(o_flat_w ,zxfluxlat_w)
1172
2/2
✓ Branch 0 taken 478114 times.
✓ Branch 1 taken 481 times.
478595 zx_tmp_fi2d=0.
1173
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1174
2/2
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 480 times.
2400 DO nsrf=1,nbsrf
1175 zx_tmp_fi2d(:)=zx_tmp_fi2d(:) &
1176
2/2
✓ Branch 0 taken 1908480 times.
✓ Branch 1 taken 1920 times.
1910880 +pctsrf(:,nsrf)*delta_tsurf(:,nsrf)
1177 ENDDO
1178 ENDIF
1179 481 CALL histwrite_phy(o_delta_tsurf,zx_tmp_fi2d)
1180 481 CALL histwrite_phy(o_cdragh_x ,cdragh_x )
1181 481 CALL histwrite_phy(o_cdragh_w ,cdragh_w )
1182 481 CALL histwrite_phy(o_cdragm_x ,cdragm_x )
1183 481 CALL histwrite_phy(o_cdragm_w ,cdragm_w )
1184 481 CALL histwrite_phy(o_kh ,kh )
1185 481 CALL histwrite_phy(o_kh_x ,kh_x )
1186 481 CALL histwrite_phy(o_kh_w ,kh_w )
1187 ENDIF ! (iflag_pbl_split>=1)
1188 481 CALL histwrite_phy(o_ale, ale)
1189 481 CALL histwrite_phy(o_alp, alp)
1190 481 CALL histwrite_phy(o_cin, cin)
1191 481 CALL histwrite_phy(o_WAPE, wake_pe)
1192 481 CALL histwrite_phy(o_cv_gen, cv_gen)
1193 481 CALL histwrite_phy(o_wake_h, wake_h)
1194 481 CALL histwrite_phy(o_wake_dens, wake_dens)
1195 481 CALL histwrite_phy(o_wake_s, wake_s)
1196 481 CALL histwrite_phy(o_wake_deltat, wake_deltat)
1197 481 CALL histwrite_phy(o_wake_deltaq, wake_deltaq)
1198 481 CALL histwrite_phy(o_wake_omg, wake_omg)
1199
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_wake(1:klon,1:klev) &
1200
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 /pdtphys
1201 481 CALL histwrite_phy(o_dtwak, zx_tmp_fi3d)
1202
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_wake(1:klon,1:klev)/pdtphys
1203 481 CALL histwrite_phy(o_dqwak, zx_tmp_fi3d)
1204
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1205 481 CALL histwrite_phy(o_dqwak2d, zx_tmp_fi2d)
1206 ENDIF ! iflag_wake>=1
1207 481 CALL histwrite_phy(o_ftd, ftd)
1208 481 CALL histwrite_phy(o_fqd, fqd)
1209 ENDIF !(iflag_con.EQ.3)
1210
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_con.EQ.3.OR.iflag_con.EQ.30) THEN
1211 ! sortie RomP convection descente insaturee iflag_con=30
1212 ! etendue a iflag_con=3 (jyg)
1213 481 CALL histwrite_phy(o_Vprecip, Vprecip)
1214 481 CALL histwrite_phy(o_qtaa, qtaa)
1215 481 CALL histwrite_phy(o_clwaa, clw)
1216 481 CALL histwrite_phy(o_wdtrainA, wdtrainA)
1217 481 CALL histwrite_phy(o_wdtrainS, wdtrainS)
1218 481 CALL histwrite_phy(o_wdtrainM, wdtrainM)
1219 ENDIF !(iflag_con.EQ.3.or.iflag_con.EQ.30)
1220 !!! nrlmd le 10/04/2012
1221
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_trig_bl>=1) THEN
1222 481 CALL histwrite_phy(o_n2, n2)
1223 481 CALL histwrite_phy(o_s2, s2)
1224 481 CALL histwrite_phy(o_proba_notrig, proba_notrig)
1225 481 CALL histwrite_phy(o_random_notrig, random_notrig)
1226 481 CALL histwrite_phy(o_ale_bl_stat, ale_bl_stat)
1227 481 CALL histwrite_phy(o_ale_bl_trig, ale_bl_trig)
1228 ENDIF !(iflag_trig_bl>=1)
1229
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_clos_bl>=1) THEN
1230 481 CALL histwrite_phy(o_alp_bl_det, alp_bl_det)
1231 481 CALL histwrite_phy(o_alp_bl_fluct_m, alp_bl_fluct_m)
1232 CALL histwrite_phy(o_alp_bl_fluct_tke, &
1233 481 alp_bl_fluct_tke)
1234 481 CALL histwrite_phy(o_alp_bl_conv, alp_bl_conv)
1235 481 CALL histwrite_phy(o_alp_bl_stat, alp_bl_stat)
1236 ENDIF !(iflag_clos_bl>=1)
1237 !!! fin nrlmd le 10/04/2012
1238 ! Output of slab ocean variables
1239
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (type_ocean=='slab ') THEN
1240 CALL histwrite_phy(o_slab_bils, slab_wfbils)
1241 IF (nslay.EQ.1) THEN
1242 IF (vars_defined) zx_tmp_fi2d(:)=tslab(:,1)
1243 CALL histwrite_phy(o_tslab, zx_tmp_fi2d)
1244 IF (vars_defined) zx_tmp_fi2d(:)=dt_qflux(:,1)
1245 CALL histwrite_phy(o_slab_qflux, zx_tmp_fi2d)
1246 ELSE
1247 CALL histwrite_phy(o_tslab, tslab(:,1:nslay))
1248 CALL histwrite_phy(o_slab_qflux, dt_qflux(:,1:nslay))
1249 ENDIF
1250 IF (version_ocean=='sicINT') THEN
1251 CALL histwrite_phy(o_slab_bilg, slab_bilg)
1252 CALL histwrite_phy(o_slab_tice, tice)
1253 CALL histwrite_phy(o_slab_sic, seaice)
1254 ENDIF
1255 IF (slab_gm) THEN
1256 CALL histwrite_phy(o_slab_gm, dt_gm(:,1:nslay))
1257 ENDIF
1258 IF (slab_hdiff) THEN
1259 IF (nslay.EQ.1) THEN
1260 IF (vars_defined) zx_tmp_fi2d(:)=dt_hdiff(:,1)
1261 CALL histwrite_phy(o_slab_hdiff, zx_tmp_fi2d)
1262 ELSE
1263 CALL histwrite_phy(o_slab_hdiff, dt_hdiff(:,1:nslay))
1264 ENDIF
1265 ENDIF
1266 IF (slab_ekman.GT.0) THEN
1267 IF (nslay.EQ.1) THEN
1268 IF (vars_defined) zx_tmp_fi2d(:)=dt_ekman(:,1)
1269 CALL histwrite_phy(o_slab_ekman, zx_tmp_fi2d)
1270 ELSE
1271 CALL histwrite_phy(o_slab_ekman, dt_ekman(:,1:nslay))
1272 ENDIF
1273 ENDIF
1274 ENDIF !type_ocean == force/slab
1275 481 CALL histwrite_phy(o_weakinv, weak_inversion)
1276 481 CALL histwrite_phy(o_dthmin, dthmin)
1277 481 CALL histwrite_phy(o_cldtau, cldtau)
1278 481 CALL histwrite_phy(o_cldemi, cldemi)
1279 481 CALL histwrite_phy(o_pr_con_l, pmflxr(:,1:klev))
1280 481 CALL histwrite_phy(o_pr_con_i, pmflxs(:,1:klev))
1281 481 CALL histwrite_phy(o_pr_lsc_l, prfl(:,1:klev))
1282 481 CALL histwrite_phy(o_pr_lsc_i, psfl(:,1:klev))
1283 481 CALL histwrite_phy(o_re, re)
1284 481 CALL histwrite_phy(o_fl, fl)
1285
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1286
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1, klon
1287
2/2
✓ Branch 0 taken 159860 times.
✓ Branch 1 taken 317260 times.
477600 IF (zt2m(i).LE.273.15) then
1288 159860 zx_tmp_fi2d(i)=MAX(0.,rh2m(i)*100.)
1289 ELSE
1290 317260 zx_tmp_fi2d(i)=MAX(0.,MIN(100.,rh2m(i)*100.))
1291 ENDIF
1292 ENDDO
1293 ENDIF
1294 481 CALL histwrite_phy(o_rh2m, zx_tmp_fi2d)
1295
1296 ! IF (vars_defined) THEN
1297 ! DO i=1, klon
1298 ! zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
1299 ! ENDDO
1300 ! ENDIF
1301 ! CALL histwrite_phy(o_rh2m_min, zx_tmp_fi2d)
1302
1303 ! IF (vars_defined) THEN
1304 ! DO i=1, klon
1305 ! zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
1306 ! ENDDO
1307 ! ENDIF
1308 ! CALL histwrite_phy(o_rh2m_max, zx_tmp_fi2d)
1309
1310 481 CALL histwrite_phy(o_qsat2m, qsat2m)
1311 481 CALL histwrite_phy(o_tpot, tpot)
1312 481 CALL histwrite_phy(o_tpote, tpote)
1313
4/4
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 477120 times.
✓ Branch 3 taken 480 times.
477601 IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
1314 481 CALL histwrite_phy(o_SWnetOR, zx_tmp_fi2d)
1315 481 CALL histwrite_phy(o_LWdownOR, sollwdown)
1316 481 CALL histwrite_phy(o_snowl, snow_lsc)
1317 481 CALL histwrite_phy(o_solldown, sollwdown)
1318 481 CALL histwrite_phy(o_dtsvdfo, d_ts(:,is_oce))
1319 481 CALL histwrite_phy(o_dtsvdft, d_ts(:,is_ter))
1320 481 CALL histwrite_phy(o_dtsvdfg, d_ts(:,is_lic))
1321 481 CALL histwrite_phy(o_dtsvdfi, d_ts(:,is_sic))
1322 481 CALL histwrite_phy(o_z0m, z0m(:,nbsrf+1))
1323 481 CALL histwrite_phy(o_z0h, z0h(:,nbsrf+1))
1324
1325 ! od550 per species
1326 !--OLIVIER
1327 !This is warranted by treating INCA aerosols as offline aerosols
1328
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (flag_aerosol.GT.0) THEN
1329 CALL histwrite_phy(o_od443aer, od443aer)
1330 CALL histwrite_phy(o_od550aer, od550aer)
1331 CALL histwrite_phy(o_od865aer, od865aer)
1332 CALL histwrite_phy(o_abs550aer, abs550aer)
1333 CALL histwrite_phy(o_od550lt1aer, od550lt1aer)
1334 CALL histwrite_phy(o_sconcso4, sconcso4)
1335 CALL histwrite_phy(o_sconcno3, sconcno3)
1336 CALL histwrite_phy(o_sconcoa, sconcoa)
1337 CALL histwrite_phy(o_sconcbc, sconcbc)
1338 CALL histwrite_phy(o_sconcss, sconcss)
1339 CALL histwrite_phy(o_sconcdust, sconcdust)
1340 CALL histwrite_phy(o_concso4, concso4)
1341 CALL histwrite_phy(o_concno3, concno3)
1342 CALL histwrite_phy(o_concoa, concoa)
1343 CALL histwrite_phy(o_concbc, concbc)
1344 CALL histwrite_phy(o_concss, concss)
1345 CALL histwrite_phy(o_concdust, concdust)
1346 CALL histwrite_phy(o_loadso4, loadso4)
1347 CALL histwrite_phy(o_loadoa, loadoa)
1348 CALL histwrite_phy(o_loadbc, loadbc)
1349 CALL histwrite_phy(o_loadss, loadss)
1350 CALL histwrite_phy(o_loaddust, loaddust)
1351 CALL histwrite_phy(o_loadno3, loadno3)
1352 CALL histwrite_phy(o_dryod550aer, dryod550aer)
1353 DO naero = 1, naero_tot-1
1354 CALL histwrite_phy(o_drytausumaero(naero),drytausum_aero(:,naero))
1355 END DO
1356 ENDIF
1357 !--STRAT AER
1358
2/4
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
✗ Branch 2 not taken.
✓ Branch 3 taken 481 times.
481 IF (flag_aerosol.GT.0.OR.flag_aerosol_strat.GT.0) THEN
1359 DO naero = 1, naero_tot
1360 481 CALL histwrite_phy(o_tausumaero(naero),tausum_aero(:,2,naero))
1361 END DO
1362 ENDIF
1363
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (flag_aerosol_strat.GT.0) THEN
1364 CALL histwrite_phy(o_tausumaero_lw,tausum_aero(:,6,id_STRAT_phy))
1365 ENDIF
1366
1367 481 CALL histwrite_phy(o_p_tropopause, p_tropopause)
1368 481 CALL histwrite_phy(o_t_tropopause, t_tropopause)
1369 481 CALL histwrite_phy(o_z_tropopause, z_tropopause)
1370
1371 ! ThL -- In the following, we assume read_climoz == 1
1372
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1373
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = 0.0 ! Computation for strato, added ThL
1374
2/2
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
19200 DO k=1, klev
1375
2/2
✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
18626880 zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * stratomask(:,k) * 1.e3
1376 END DO
1377 ENDIF
1378 481 CALL histwrite_phy(o_col_O3_strato, zx_tmp_fi2d) ! Added ThL
1379
1380
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1381
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d = 0.0 ! Computation for tropo, added ThL
1382
2/2
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
19200 DO k=1, klev
1383
2/2
✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
18626880 zx_tmp_fi2d(:) = zx_tmp_fi2d(:) + wo(:,k,1) * (1.0-stratomask(:,k)) * 1.e3
1384 END DO
1385 ENDIF
1386 481 CALL histwrite_phy(o_col_O3_tropo, zx_tmp_fi2d) ! Added ThL
1387 ! end add ThL
1388
1389 !NL
1390
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
481 IF (ok_volcan .AND. ok_ade) THEN
1391 DO k=1, klev
1392 IF (vars_defined) zx_tmp_fi3d(:,k)=heat_volc(:,k)*swradcorr(:)
1393 ENDDO
1394 CALL histwrite_phy(o_heat_volc, zx_tmp_fi3d)
1395 DO k=1, klev
1396 IF (vars_defined) zx_tmp_fi3d(:,k)=cool_volc(:,k)
1397 ENDDO
1398 CALL histwrite_phy(o_cool_volc, zx_tmp_fi3d)
1399 ENDIF
1400
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (ok_ade) THEN
1401 IF (vars_defined) zx_tmp_fi2d(:)=topswad_aero*swradcorr
1402 CALL histwrite_phy(o_topswad, zx_tmp_fi2d)
1403
1404 IF (vars_defined) zx_tmp_fi2d(:)=topswad0_aero*swradcorr
1405 CALL histwrite_phy(o_topswad0, zx_tmp_fi2d)
1406
1407 IF (vars_defined) zx_tmp_fi2d(:)=solswad_aero*swradcorr
1408 CALL histwrite_phy(o_solswad, zx_tmp_fi2d)
1409
1410 IF (vars_defined) zx_tmp_fi2d(:)=solswad0_aero*swradcorr
1411 CALL histwrite_phy(o_solswad0, zx_tmp_fi2d)
1412
1413 CALL histwrite_phy(o_toplwad, toplwad_aero)
1414 CALL histwrite_phy(o_toplwad0, toplwad0_aero)
1415 CALL histwrite_phy(o_sollwad, sollwad_aero)
1416 CALL histwrite_phy(o_sollwad0, sollwad0_aero)
1417 !====MS forcing diagnostics
1418 !ym warning : topsw_aero, solsw_aero, topsw0_aero, solsw0_aero are not defined by model
1419 !ym => init to 0 in radlwsw_m.F90 ztopsw_aero, zsolsw_aero, ztopsw0_aero, zsolsw0_aero
1420
1421 IF (vars_defined) zx_tmp_fi2d(:)=topsw_aero(:,1)*swradcorr(:)
1422 CALL histwrite_phy(o_swtoaas_nat,zx_tmp_fi2d)
1423 IF (vars_defined) zx_tmp_fi2d(:)=solsw_aero(:,1)*swradcorr(:)
1424 CALL histwrite_phy(o_swsrfas_nat,zx_tmp_fi2d)
1425 IF (vars_defined) zx_tmp_fi2d(:)=topsw0_aero(:,1)*swradcorr(:)
1426 CALL histwrite_phy(o_swtoacs_nat,zx_tmp_fi2d)
1427 IF (vars_defined) zx_tmp_fi2d(:)=solsw0_aero(:,1)*swradcorr(:)
1428 CALL histwrite_phy(o_swsrfcs_nat,zx_tmp_fi2d)
1429 !ant
1430 IF (vars_defined) zx_tmp_fi2d(:)=topsw_aero(:,2)*swradcorr(:)
1431 CALL histwrite_phy(o_swtoaas_ant,zx_tmp_fi2d)
1432 IF (vars_defined) zx_tmp_fi2d(:)=solsw_aero(:,2)*swradcorr(:)
1433 CALL histwrite_phy(o_swsrfas_ant,zx_tmp_fi2d)
1434 IF (vars_defined) zx_tmp_fi2d(:)=topsw0_aero(:,2)*swradcorr(:)
1435 CALL histwrite_phy(o_swtoacs_ant,zx_tmp_fi2d)
1436 IF (vars_defined) zx_tmp_fi2d(:)=solsw0_aero(:,2)*swradcorr(:)
1437 CALL histwrite_phy(o_swsrfcs_ant,zx_tmp_fi2d)
1438 !cf
1439 IF (.not. aerosol_couple) THEN
1440 IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,1)*swradcorr(:)
1441 CALL histwrite_phy(o_swtoacf_nat,zx_tmp_fi2d)
1442 IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,1)*swradcorr(:)
1443 CALL histwrite_phy(o_swsrfcf_nat,zx_tmp_fi2d)
1444 IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,2)*swradcorr(:)
1445 CALL histwrite_phy(o_swtoacf_ant,zx_tmp_fi2d)
1446 IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,2)*swradcorr(:)
1447 CALL histwrite_phy(o_swsrfcf_ant,zx_tmp_fi2d)
1448 IF (vars_defined) zx_tmp_fi2d(:)=topswcf_aero(:,3)*swradcorr(:)
1449 CALL histwrite_phy(o_swtoacf_zero,zx_tmp_fi2d)
1450 IF (vars_defined) zx_tmp_fi2d(:)=solswcf_aero(:,3)*swradcorr(:)
1451 CALL histwrite_phy(o_swsrfcf_zero,zx_tmp_fi2d)
1452 ENDIF
1453 !====MS forcing diagnostics
1454 ENDIF
1455
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (ok_aie) THEN
1456 IF (vars_defined) zx_tmp_fi2d(:)= topswai_aero*swradcorr
1457 CALL histwrite_phy(o_topswai, zx_tmp_fi2d)
1458
1459 IF (vars_defined) zx_tmp_fi2d(:)=toplwai_aero*swradcorr
1460 CALL histwrite_phy(o_toplwai, zx_tmp_fi2d)
1461
1462 IF (vars_defined) zx_tmp_fi2d(:)=solswai_aero*swradcorr
1463 CALL histwrite_phy(o_solswai, zx_tmp_fi2d)
1464
1465 IF (vars_defined) zx_tmp_fi2d(:)=sollwai_aero*swradcorr
1466 CALL histwrite_phy(o_sollwai, zx_tmp_fi2d)
1467 ENDIF
1468
1/4
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
✗ Branch 2 not taken.
✗ Branch 3 not taken.
481 IF (flag_aerosol.GT.0.AND.ok_cdnc) THEN
1469 CALL histwrite_phy(o_scdnc, scdnc)
1470 CALL histwrite_phy(o_cldncl, cldncl)
1471 CALL histwrite_phy(o_reffclws, reffclws)
1472 CALL histwrite_phy(o_reffclwc, reffclwc)
1473 CALL histwrite_phy(o_cldnvi, cldnvi)
1474 CALL histwrite_phy(o_lcc, lcc)
1475 CALL histwrite_phy(o_lcc3d, lcc3d)
1476 CALL histwrite_phy(o_lcc3dcon, lcc3dcon)
1477 CALL histwrite_phy(o_lcc3dstra, lcc3dstra)
1478 CALL histwrite_phy(o_icc3dcon, icc3dcon)
1479 CALL histwrite_phy(o_icc3dstra, icc3dstra)
1480 CALL histwrite_phy(o_cldicemxrat, zfice)
1481 IF (vars_defined) zx_tmp_fi3d(:,:)=1-zfice(:,:)
1482 CALL histwrite_phy(o_cldwatmxrat, zx_tmp_fi3d)
1483 CALL histwrite_phy(o_reffclwtop, reffclwtop)
1484 ENDIF
1485 ! Champs 3D:
1486
2/4
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
✗ Branch 2 not taken.
✓ Branch 3 taken 481 times.
481 IF (ok_ade .OR. ok_aie) then
1487 CALL histwrite_phy(o_ec550aer, ec550aer)
1488 ENDIF
1489 481 CALL histwrite_phy(o_lwcon, flwc)
1490 481 CALL histwrite_phy(o_iwcon, fiwc)
1491 481 CALL histwrite_phy(o_temp, t_seri)
1492 481 CALL histwrite_phy(o_theta, theta)
1493 481 CALL histwrite_phy(o_ovapinit, qx(:,:,ivap))
1494 481 CALL histwrite_phy(o_ovap, q_seri)
1495 481 CALL histwrite_phy(o_oliq, ql_seri)
1496
1497
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d = ql_seri+qs_seri
1498 481 CALL histwrite_phy(o_ocond, zx_tmp_fi3d)
1499
1500 481 CALL histwrite_phy(o_geop, zphi)
1501 481 CALL histwrite_phy(o_vitu, u_seri)
1502 481 CALL histwrite_phy(o_vitv, v_seri)
1503 481 CALL histwrite_phy(o_vitw, omega)
1504 481 CALL histwrite_phy(o_pres, pplay)
1505 481 CALL histwrite_phy(o_paprs, paprs(:,1:klev))
1506
1507
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d = zphi/RG
1508 481 CALL histwrite_phy(o_zfull,zx_tmp_fi3d)
1509
1510
1511
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (flag_aerosol_strat.EQ.2) THEN
1512 CALL histwrite_phy(o_stratomask, stratomask)
1513 ENDIF
1514
1515
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1516
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi3d(:,1)= pphis(:)/RG
1517
2/2
✓ Branch 0 taken 18240 times.
✓ Branch 1 taken 480 times.
18720 DO k = 2, klev
1518
2/2
✓ Branch 0 taken 18130560 times.
✓ Branch 1 taken 18240 times.
18149280 DO i = 1, klon
1519 zx_tmp_fi3d(i,k) = zphi(i,k-1)/RG + &
1520 (zphi(i,k)-zphi(i,k-1))/RG * &
1521 18148800 (paprs(i,k)-pplay(i,k-1))/(pplay(i,k)-pplay(i,k-1))
1522 ENDDO
1523 ENDDO
1524 ENDIF
1525 481 CALL histwrite_phy(o_zhalf, zx_tmp_fi3d)
1526 481 CALL histwrite_phy(o_rneb, cldfra)
1527 481 CALL histwrite_phy(o_rnebcon, rnebcon)
1528 481 CALL histwrite_phy(o_rnebls, rneb)
1529 481 CALL histwrite_phy(o_rneblsvol, rneblsvol)
1530
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1531
2/2
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
19200 DO k=1, klev
1532
2/2
✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
18626880 DO i=1, klon
1533 18626400 zx_tmp_fi3d(i,k)=cldfra(i,k)*JrNt(i)
1534 ENDDO
1535 ENDDO
1536 ENDIF
1537 481 CALL histwrite_phy(o_rnebjn, zx_tmp_fi3d)
1538 481 CALL histwrite_phy(o_rhum, zx_rh)
1539
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_ice_thermo .GT. 0) THEN
1540
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d = zx_rhl * 100.
1541 481 CALL histwrite_phy(o_rhl, zx_tmp_fi3d)
1542
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d = zx_rhi * 100.
1543 481 CALL histwrite_phy(o_rhi, zx_tmp_fi3d)
1544 ENDIF
1545
1546
1547
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d = wo(:, :, 1) * dobson_u * 1e3 / zmasse / rmo3 * rmd
1548 481 CALL histwrite_phy(o_ozone, zx_tmp_fi3d)
1549
1550
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (read_climoz == 2) THEN
1551 IF (vars_defined) zx_tmp_fi3d = wo(:, :, 2) * dobson_u * 1e3 / zmasse / rmo3 * rmd
1552 CALL histwrite_phy(o_ozone_light, zx_tmp_fi3d)
1553 ENDIF
1554
1555 481 CALL histwrite_phy(o_duphy, d_u)
1556
1557 481 CALL histwrite_phy(o_dtphy, d_t)
1558
1559 481 CALL histwrite_phy(o_dqphy, d_qx(:,:,ivap))
1560
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,ivap),zmasse,zx_tmp_fi2d)
1561 481 CALL histwrite_phy(o_dqphy2d, zx_tmp_fi2d)
1562
1563 481 CALL histwrite_phy(o_dqlphy, d_qx(:,:,iliq))
1564
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,iliq),zmasse,zx_tmp_fi2d)
1565 481 CALL histwrite_phy(o_dqlphy2d, zx_tmp_fi2d)
1566
1567
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (nqo.EQ.3) THEN
1568 481 CALL histwrite_phy(o_dqsphy, d_qx(:,:,isol))
1569
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,d_qx(:,:,isol),zmasse,zx_tmp_fi2d)
1570 481 CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d)
1571 ELSE
1572 zx_tmp_fi3d=0.0
1573 CALL histwrite_phy(o_dqsphy, zx_tmp_fi3d)
1574 zx_tmp_fi2d=0.0
1575 CALL histwrite_phy(o_dqsphy2d, zx_tmp_fi2d)
1576 ENDIF
1577
1578
2/2
✓ Branch 0 taken 1924 times.
✓ Branch 1 taken 481 times.
2405 DO nsrf=1, nbsrf
1579
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
1580 1924 CALL histwrite_phy(o_albe_srf(nsrf), zx_tmp_fi2d)
1581
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = z0m( 1 : klon, nsrf)
1582 1924 CALL histwrite_phy(o_z0m_srf(nsrf), zx_tmp_fi2d)
1583
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = z0h( 1 : klon, nsrf)
1584 1924 CALL histwrite_phy(o_z0h_srf(nsrf), zx_tmp_fi2d)
1585
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = agesno( 1 : klon, nsrf)
1586 1924 CALL histwrite_phy(o_ages_srf(nsrf), zx_tmp_fi2d)
1587
4/4
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 4 times.
✓ Branch 2 taken 1908480 times.
✓ Branch 3 taken 1920 times.
1910404 IF (vars_defined) zx_tmp_fi2d(1 : klon) = snow( 1 : klon, nsrf)
1588 2405 CALL histwrite_phy(o_snow_srf(nsrf), zx_tmp_fi2d)
1589 ENDDO !nsrf=1, nbsrf
1590 481 CALL histwrite_phy(o_alb1, albsol1)
1591 481 CALL histwrite_phy(o_alb2, albsol2)
1592 !FH Sorties pour la couche limite
1593
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_pbl>1) THEN
1594
4/4
✓ Branch 0 taken 18759 times.
✓ Branch 1 taken 481 times.
✓ Branch 2 taken 18646446 times.
✓ Branch 3 taken 18759 times.
18665686 zx_tmp_fi3d=0.
1595
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1596
2/2
✓ Branch 0 taken 1920 times.
✓ Branch 1 taken 480 times.
2400 DO nsrf=1,nbsrf
1597
2/2
✓ Branch 0 taken 74880 times.
✓ Branch 1 taken 1920 times.
77280 DO k=1,klev
1598 zx_tmp_fi3d(:,k)=zx_tmp_fi3d(:,k) &
1599
2/2
✓ Branch 0 taken 74430720 times.
✓ Branch 1 taken 74880 times.
74507520 +pctsrf(:,nsrf)*pbl_tke(:,k,nsrf)
1600 ENDDO
1601 ENDDO
1602 ENDIF
1603 481 CALL histwrite_phy(o_tke, zx_tmp_fi3d)
1604 481 CALL histwrite_phy(o_tke_max, zx_tmp_fi3d)
1605
1606 ENDIF
1607
1608 481 CALL histwrite_phy(o_kz, coefh(:,:,is_ave))
1609
1610 481 CALL histwrite_phy(o_kz_max, coefh(:,:,is_ave))
1611
1612 481 CALL histwrite_phy(o_clwcon, clwcon0)
1613 481 CALL histwrite_phy(o_dtdyn, d_t_dyn)
1614
1615 481 CALL histwrite_phy(o_dqdyn, d_q_dyn)
1616
1617 481 CALL histwrite_phy(o_dqdyn2d,d_q_dyn2d)
1618
1619 481 CALL histwrite_phy(o_dqldyn, d_ql_dyn)
1620
1621 481 CALL histwrite_phy(o_dqldyn2d, d_ql_dyn2d)
1622
1623 481 CALL histwrite_phy(o_dqsdyn, d_qs_dyn)
1624
1625 481 CALL histwrite_phy(o_dqsdyn2d, d_qs_dyn2d)
1626
1627 481 CALL histwrite_phy(o_dudyn, d_u_dyn)
1628 481 CALL histwrite_phy(o_dvdyn, d_v_dyn)
1629
1630
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1631
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys
1632 ENDIF
1633 481 CALL histwrite_phy(o_dtcon, zx_tmp_fi3d)
1634
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (iflag_thermals.EQ.0) THEN
1635 IF (vars_defined) THEN
1636 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
1637 d_t_ajsb(1:klon,1:klev)/pdtphys
1638 ENDIF
1639 CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
1640
2/4
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 481 times.
✗ Branch 3 not taken.
481 ELSE IF(iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
1641
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1642 zx_tmp_fi3d(1:klon,1:klev)=d_t_con(1:klon,1:klev)/pdtphys + &
1643 d_t_ajs(1:klon,1:klev)/pdtphys + &
1644
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_t_wake(1:klon,1:klev)/pdtphys
1645 ENDIF
1646 481 CALL histwrite_phy(o_tntc, zx_tmp_fi3d)
1647 ENDIF
1648
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_con(1:klon,1:klev)/pdtphys
1649 481 CALL histwrite_phy(o_ducon, zx_tmp_fi3d)
1650
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_con(1:klon,1:klev)/pdtphys
1651 481 CALL histwrite_phy(o_dvcon, zx_tmp_fi3d)
1652
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1653 481 CALL histwrite_phy(o_dqcon, zx_tmp_fi3d)
1654
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1655 481 CALL histwrite_phy(o_dqcon2d, zx_tmp_fi2d)
1656
1657
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (iflag_thermals.EQ.0) THEN
1658 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys
1659 CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
1660
2/4
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 481 times.
✗ Branch 3 not taken.
481 ELSE IF (iflag_thermals.GE.1.AND.iflag_wake.EQ.1) THEN
1661
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1662 zx_tmp_fi3d(1:klon,1:klev)=d_q_con(1:klon,1:klev)/pdtphys + &
1663 d_q_ajs(1:klon,1:klev)/pdtphys + &
1664
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_q_wake(1:klon,1:klev)/pdtphys
1665 ENDIF
1666 481 CALL histwrite_phy(o_tnhusc, zx_tmp_fi3d)
1667 ENDIF
1668
1669
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lsc(1:klon,1:klev)/pdtphys
1670 481 CALL histwrite_phy(o_dtlsc, zx_tmp_fi3d)
1671
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) zx_tmp_fi3d(1:klon, 1:klev)=(d_t_lsc(1:klon,1:klev)+ &
1672
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_t_eva(1:klon,1:klev))/pdtphys
1673 481 CALL histwrite_phy(o_dtlschr, zx_tmp_fi3d)
1674
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys
1675 481 CALL histwrite_phy(o_dqlsc, zx_tmp_fi3d)
1676
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1677 481 CALL histwrite_phy(o_dqlsc2d, zx_tmp_fi2d)
1678
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=beta_prec(1:klon,1:klev)
1679 481 CALL histwrite_phy(o_beta_prec, zx_tmp_fi3d)
1680 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1681 ! Sorties specifiques a la separation thermiques/non thermiques
1682
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_thermals>=1) THEN
1683
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscth(1:klon,1:klev)/pdtphys
1684 481 CALL histwrite_phy(o_dtlscth, zx_tmp_fi3d)
1685
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lscst(1:klon,1:klev)/pdtphys
1686 481 CALL histwrite_phy(o_dtlscst, zx_tmp_fi3d)
1687
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscth(1:klon,1:klev)/pdtphys
1688 481 CALL histwrite_phy(o_dqlscth, zx_tmp_fi3d)
1689
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1690 481 CALL histwrite_phy(o_dqlscth2d, zx_tmp_fi2d)
1691
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_lscst(1:klon,1:klev)/pdtphys
1692 481 CALL histwrite_phy(o_dqlscst, zx_tmp_fi3d)
1693
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1694 481 CALL histwrite_phy(o_dqlscst2d, zx_tmp_fi2d)
1695 481 CALL histwrite_phy(o_plulth, plul_th)
1696 481 CALL histwrite_phy(o_plulst, plul_st)
1697
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1698
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 DO i=1,klon
1699
2/2
✓ Branch 0 taken 474257280 times.
✓ Branch 1 taken 477120 times.
474734880 zx_tmp_fi2d(1:klon)=lmax_th(:)
1700 ENDDO
1701 ENDIF
1702 481 CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d)
1703
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1704
2/2
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
19200 DO k=1,klev
1705
2/2
✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
18626880 DO i=1,klon
1706
2/2
✓ Branch 0 taken 1305081 times.
✓ Branch 1 taken 17302599 times.
18626400 IF (ptconvth(i,k)) THEN
1707 1305081 zx_tmp_fi3d(i,k)=1.
1708 ELSE
1709 17302599 zx_tmp_fi3d(i,k)=0.
1710 ENDIF
1711 ENDDO
1712 ENDDO
1713 ENDIF
1714 481 CALL histwrite_phy(o_ptconvth, zx_tmp_fi3d)
1715 ENDIF ! iflag_thermals>=1
1716 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1717
4/4
✓ Branch 0 taken 18759 times.
✓ Branch 1 taken 481 times.
✓ Branch 2 taken 18646446 times.
✓ Branch 3 taken 18759 times.
18665686 zpt_conv = 0.
1718
6/6
✓ Branch 0 taken 18759 times.
✓ Branch 1 taken 481 times.
✓ Branch 2 taken 18646446 times.
✓ Branch 3 taken 18759 times.
✓ Branch 4 taken 1038140 times.
✓ Branch 5 taken 17608306 times.
18665686 WHERE (ptconv) zpt_conv = 1.
1719 481 CALL histwrite_phy(o_ptconv, zpt_conv)
1720 !! IF (vars_defined) zx_tmp_fi2d=float(itau_con)/float(itap)
1721 !! CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
1722
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1723
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zpt_conv2d(:) = 0.
1724
2/2
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
19200 DO k=1,klev
1725
4/4
✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
✓ Branch 2 taken 1038140 times.
✓ Branch 3 taken 17569540 times.
18626880 WHERE (ptconv(:,k)) zpt_conv2d(:) = 1.
1726 ENDDO
1727 ENDIF
1728 481 CALL histwrite_phy(o_ftime_deepcv, zpt_conv2d)
1729
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1730
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = 0.
1731
2/2
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
19200 DO k=1,klev
1732
4/4
✓ Branch 0 taken 18607680 times.
✓ Branch 1 taken 18720 times.
✓ Branch 2 taken 1305081 times.
✓ Branch 3 taken 17302599 times.
18626880 WHERE (ptconvth(:,k)) zx_tmp_fi2d(:) = 1.
1733 ENDDO
1734 ENDIF
1735 481 CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
1736
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1737
2/2
✓ Branch 0 taken 477120 times.
✓ Branch 1 taken 480 times.
477600 zx_tmp_fi2d(:) = max(zx_tmp_fi2d(:),zpt_conv2d(:))
1738 ENDIF
1739 481 CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
1740 !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
1741
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
1742 481 CALL histwrite_phy(o_dtvdf, zx_tmp_fi3d)
1743
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_diss(1:klon,1:klev)/pdtphys
1744 481 CALL histwrite_phy(o_dtdis, zx_tmp_fi3d)
1745
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf(1:klon,1:klev)/pdtphys
1746 481 CALL histwrite_phy(o_dqvdf, zx_tmp_fi3d)
1747
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1748 481 CALL histwrite_phy(o_dqvdf2d, zx_tmp_fi2d)
1749
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_eva(1:klon,1:klev)/pdtphys
1750 481 CALL histwrite_phy(o_dteva, zx_tmp_fi3d)
1751
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_eva(1:klon,1:klev)/pdtphys
1752 481 CALL histwrite_phy(o_dqeva, zx_tmp_fi3d)
1753
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1754 481 CALL histwrite_phy(o_dqeva2d, zx_tmp_fi2d)
1755 481 CALL histwrite_phy(o_ratqs, ratqs)
1756
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1757 zx_tmp_fi3d(1:klon,1:klev)=d_t_ajs(1:klon,1:klev)/pdtphys - &
1758
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_t_ajsb(1:klon,1:klev)/pdtphys
1759 ENDIF
1760 481 CALL histwrite_phy(o_dtthe, zx_tmp_fi3d)
1761
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1762
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys
1763 ENDIF
1764 481 CALL histwrite_phy(o_duthe, zx_tmp_fi3d)
1765
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1766
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 zx_tmp_fi3d(1:klon,1:klev)=d_v_ajs(1:klon,1:klev)/pdtphys
1767 ENDIF
1768 481 CALL histwrite_phy(o_dvthe, zx_tmp_fi3d)
1769
1770
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_thermals>=1) THEN
1771 ! Pour l instant 0 a y reflichir pour les thermiques
1772 ! regroupe avec ftime_deepcv et ftime_con
1773 !!zx_tmp_fi2d=0.
1774 !!CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
1775 481 CALL histwrite_phy(o_f_th, fm_therm)
1776 481 CALL histwrite_phy(o_e_th, entr_therm)
1777 481 CALL histwrite_phy(o_w_th, zw2)
1778 481 CALL histwrite_phy(o_q_th, zqasc)
1779 481 CALL histwrite_phy(o_a_th, fraca)
1780 481 CALL histwrite_phy(o_cloudth_sth, cloudth_sth)
1781 481 CALL histwrite_phy(o_cloudth_senv, cloudth_senv)
1782 481 CALL histwrite_phy(o_cloudth_sigmath, cloudth_sigmath)
1783 481 CALL histwrite_phy(o_cloudth_sigmaenv, cloudth_sigmaenv)
1784 481 CALL histwrite_phy(o_d_th, detr_therm)
1785 481 CALL histwrite_phy(o_f0_th, f0)
1786 481 CALL histwrite_phy(o_zmax_th, zmax_th)
1787
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1788 zx_tmp_fi3d(1:klon,1:klev)=d_q_ajs(1:klon,1:klev)/pdtphys - &
1789
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_q_ajsb(1:klon,1:klev)/pdtphys
1790 ENDIF
1791 481 CALL histwrite_phy(o_dqthe, zx_tmp_fi3d)
1792
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1793 481 CALL histwrite_phy(o_dqthe2d, zx_tmp_fi2d)
1794 ENDIF !iflag_thermals
1795
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ajsb(1:klon,1:klev)/pdtphys
1796 481 CALL histwrite_phy(o_dtajs, zx_tmp_fi3d)
1797
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
1798 481 CALL histwrite_phy(o_dqajs, zx_tmp_fi3d)
1799
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
1800 481 CALL histwrite_phy(o_dqajs2d, zx_tmp_fi2d)
1801
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys
1802 481 CALL histwrite_phy(o_dtswr, zx_tmp_fi3d)
1803
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_sw0(1:klon,1:klev)/pdtphys
1804 481 CALL histwrite_phy(o_dtsw0, zx_tmp_fi3d)
1805
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lwr(1:klon,1:klev)/pdtphys
1806 481 CALL histwrite_phy(o_dtlwr, zx_tmp_fi3d)
1807
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lw0(1:klon,1:klev)/pdtphys
1808 481 CALL histwrite_phy(o_dtlw0, zx_tmp_fi3d)
1809
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)/pdtphys
1810 481 CALL histwrite_phy(o_dtec, zx_tmp_fi3d)
1811
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_vdf(1:klon,1:klev)/pdtphys
1812 481 CALL histwrite_phy(o_duvdf, zx_tmp_fi3d)
1813
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_vdf(1:klon,1:klev)/pdtphys
1814 481 CALL histwrite_phy(o_dvvdf, zx_tmp_fi3d)
1815
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (ok_orodr) THEN
1816
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_oro(1:klon,1:klev)/pdtphys
1817 481 CALL histwrite_phy(o_duoro, zx_tmp_fi3d)
1818
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_oro(1:klon,1:klev)/pdtphys
1819 481 CALL histwrite_phy(o_dvoro, zx_tmp_fi3d)
1820
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_oro(1:klon,1:klev)/pdtphys
1821 481 CALL histwrite_phy(o_dtoro, zx_tmp_fi3d)
1822 ENDIF
1823
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (ok_orolf) THEN
1824
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
1825 481 CALL histwrite_phy(o_dulif, zx_tmp_fi3d)
1826
1827
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
1828 481 CALL histwrite_phy(o_dvlif, zx_tmp_fi3d)
1829
1830
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys
1831 481 CALL histwrite_phy(o_dtlif, zx_tmp_fi3d)
1832 ENDIF
1833
1834
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (ok_hines) THEN
1835 IF (vars_defined) zx_tmp_fi3d=du_gwd_hines/pdtphys
1836 CALL histwrite_phy(o_du_gwd_hines, zx_tmp_fi3d)
1837
1838 IF (vars_defined) zx_tmp_fi3d= dv_gwd_hines/pdtphys
1839 CALL histwrite_phy(o_dv_gwd_hines, zx_tmp_fi3d)
1840
1841 IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_hin(1:klon,1:klev)/pdtphys
1842 CALL histwrite_phy(o_dthin, zx_tmp_fi3d)
1843 CALL histwrite_phy(o_ustr_gwd_hines, zustr_gwd_hines)
1844 CALL histwrite_phy(o_vstr_gwd_hines, zvstr_gwd_hines)
1845 ENDIF
1846
1847
2/4
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
✓ Branch 2 taken 481 times.
✗ Branch 3 not taken.
481 IF (.not. ok_hines .and. ok_gwd_rando) THEN
1848
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d=du_gwd_front / pdtphys
1849 481 CALL histwrite_phy(o_du_gwd_front, zx_tmp_fi3d)
1850
1851
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d=dv_gwd_front / pdtphys
1852 481 CALL histwrite_phy(o_dv_gwd_front, zx_tmp_fi3d)
1853
1854 481 CALL histwrite_phy(o_ustr_gwd_front, zustr_gwd_front)
1855 481 CALL histwrite_phy(o_vstr_gwd_front, zvstr_gwd_front)
1856 ENDIF
1857
1858
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (ok_gwd_rando) THEN
1859
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d=du_gwd_rando / pdtphys
1860 481 CALL histwrite_phy(o_du_gwd_rando, zx_tmp_fi3d)
1861
1862
6/6
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
✓ Branch 2 taken 18720 times.
✓ Branch 3 taken 480 times.
✓ Branch 4 taken 18607680 times.
✓ Branch 5 taken 18720 times.
18626881 IF (vars_defined) zx_tmp_fi3d=dv_gwd_rando / pdtphys
1863 481 CALL histwrite_phy(o_dv_gwd_rando, zx_tmp_fi3d)
1864 481 CALL histwrite_phy(o_ustr_gwd_rando, zustr_gwd_rando)
1865 481 CALL histwrite_phy(o_vstr_gwd_rando, zvstr_gwd_rando)
1866 481 CALL histwrite_phy(o_east_gwstress, east_gwstress )
1867 481 CALL histwrite_phy(o_west_gwstress, west_gwstress )
1868 ENDIF
1869
1870
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (ok_qch4) THEN
1871 IF (vars_defined) zx_tmp_fi3d=d_q_ch4 / pdtphys
1872 CALL histwrite_phy(o_dqch4, zx_tmp_fi3d)
1873 ENDIF
1874
1875
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1876
2/2
✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 480 times.
19680 DO k=1, klevp1
1877
2/2
✓ Branch 0 taken 19084800 times.
✓ Branch 1 taken 19200 times.
19104480 zx_tmp_fi3d1(:,k)=swup(:,k)*swradcorr(:)
1878 ENDDO
1879 ENDIF
1880
1881 481 CALL histwrite_phy(o_rsu, zx_tmp_fi3d1)
1882
1883
1884
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1885
2/2
✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 480 times.
19680 DO k=1, klevp1
1886
2/2
✓ Branch 0 taken 19084800 times.
✓ Branch 1 taken 19200 times.
19104480 zx_tmp_fi3d1(:,k)=swdn(:,k)*swradcorr(:)
1887 ENDDO
1888 ENDIF
1889
1890 481 CALL histwrite_phy(o_rsd, zx_tmp_fi3d1)
1891
1892
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1893
2/2
✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 480 times.
19680 DO k=1, klevp1
1894
2/2
✓ Branch 0 taken 19084800 times.
✓ Branch 1 taken 19200 times.
19104480 zx_tmp_fi3d1(:,k)=swup0(:,k)*swradcorr(:)
1895 ENDDO
1896 ENDIF
1897
1898 481 CALL histwrite_phy(o_rsucs, zx_tmp_fi3d1)
1899
1900
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1901
2/2
✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 480 times.
19680 DO k=1, klevp1
1902
2/2
✓ Branch 0 taken 19084800 times.
✓ Branch 1 taken 19200 times.
19104480 zx_tmp_fi3d1(:,k)=swupc0(:,k)*swradcorr(:)
1903 ENDDO
1904 ENDIF
1905 481 CALL histwrite_phy(o_rsucsaf, zx_tmp_fi3d1)
1906
1907
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1908
2/2
✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 480 times.
19680 DO k=1, klevp1
1909
2/2
✓ Branch 0 taken 19084800 times.
✓ Branch 1 taken 19200 times.
19104480 zx_tmp_fi3d1(:,k)=swdn0(:,k)*swradcorr(:)
1910 ENDDO
1911 ENDIF
1912 481 CALL histwrite_phy(o_rsdcs, zx_tmp_fi3d1)
1913
1914
1915
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1916
2/2
✓ Branch 0 taken 19200 times.
✓ Branch 1 taken 480 times.
19680 DO k=1, klevp1
1917
2/2
✓ Branch 0 taken 19084800 times.
✓ Branch 1 taken 19200 times.
19104480 zx_tmp_fi3d1(:,k)=swdnc0(:,k)*swradcorr(:)
1918 ENDDO
1919 ENDIF
1920 481 CALL histwrite_phy(o_rsdcsaf, zx_tmp_fi3d1)
1921
1922 481 CALL histwrite_phy(o_rlu, lwup)
1923 481 CALL histwrite_phy(o_rld, lwdn)
1924 481 CALL histwrite_phy(o_rlucs, lwup0)
1925 481 CALL histwrite_phy(o_rldcs, lwdn0)
1926
1927
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1928 zx_tmp_fi3d(1:klon,1:klev)=d_t(1:klon,1:klev)+ &
1929
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_t_dyn(1:klon,1:klev)
1930 ENDIF
1931 481 CALL histwrite_phy(o_tnt, zx_tmp_fi3d)
1932
1933
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1934 zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys + &
1935
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_t_lwr(1:klon,1:klev)/pdtphys
1936 ENDIF
1937 481 CALL histwrite_phy(o_tntr, zx_tmp_fi3d)
1938
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1939 zx_tmp_fi3d(1:klon,1:klev)= (d_t_lsc(1:klon,1:klev)+ &
1940 d_t_eva(1:klon,1:klev)+ &
1941
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_t_vdf(1:klon,1:klev))/pdtphys
1942 ENDIF
1943 481 CALL histwrite_phy(o_tntscpbl, zx_tmp_fi3d)
1944
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1945 zx_tmp_fi3d(1:klon,1:klev)=d_qx(1:klon,1:klev,ivap)+ &
1946
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_q_dyn(1:klon,1:klev)
1947 ENDIF
1948 481 CALL histwrite_phy(o_tnhus, zx_tmp_fi3d)
1949
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1950 zx_tmp_fi3d(1:klon,1:klev)=d_q_lsc(1:klon,1:klev)/pdtphys+ &
1951
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 d_q_eva(1:klon,1:klev)/pdtphys
1952 ENDIF
1953 481 CALL histwrite_phy(o_tnhusscpbl, zx_tmp_fi3d)
1954 481 CALL histwrite_phy(o_evu, coefm(:,:,is_ave))
1955
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1956 zx_tmp_fi3d(1:klon,1:klev)=q_seri(1:klon,1:klev)+ &
1957
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 ql_seri(1:klon,1:klev)
1958 ENDIF
1959 481 CALL histwrite_phy(o_h2o, zx_tmp_fi3d)
1960
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (iflag_con >= 3) THEN
1961
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1962 zx_tmp_fi3d(1:klon,1:klev)=-1 * (dnwd(1:klon,1:klev)+ &
1963
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 dnwd0(1:klon,1:klev))
1964 ENDIF
1965 481 CALL histwrite_phy(o_mcd, zx_tmp_fi3d)
1966
2/2
✓ Branch 0 taken 480 times.
✓ Branch 1 taken 1 times.
481 IF (vars_defined) THEN
1967 zx_tmp_fi3d(1:klon,1:klev)=upwd(1:klon,1:klev) + &
1968
4/4
✓ Branch 0 taken 18720 times.
✓ Branch 1 taken 480 times.
✓ Branch 2 taken 18607680 times.
✓ Branch 3 taken 18720 times.
18626880 dnwd(1:klon,1:klev)+ dnwd0(1:klon,1:klev)
1969 ENDIF
1970 481 CALL histwrite_phy(o_dmc, zx_tmp_fi3d)
1971 ELSE IF (iflag_con == 2) THEN
1972 CALL histwrite_phy(o_mcd, pmfd)
1973 IF (vars_defined) zx_tmp_fi3d = pmfu + pmfd
1974 CALL histwrite_phy(o_dmc, zx_tmp_fi3d)
1975 ENDIF
1976 481 CALL histwrite_phy(o_ref_liq, ref_liq)
1977 481 CALL histwrite_phy(o_ref_ice, ref_ice)
1978 !
1979
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 481 times.
481 IF (ok_4xCO2atm) THEN
1980 IF (vars_defined) zx_tmp_fi2d(:) = swupp(:,klevp1)*swradcorr(:)
1981 CALL histwrite_phy(o_rsut4co2, zx_tmp_fi2d)
1982 IF (vars_defined) zx_tmp_fi2d(:) = lwupp(:,klevp1)
1983 CALL histwrite_phy(o_rlut4co2, zx_tmp_fi2d)
1984 IF (vars_defined) zx_tmp_fi2d(:) = swup0p(:,klevp1)*swradcorr(:)
1985 CALL histwrite_phy(o_rsutcs4co2, zx_tmp_fi2d)
1986 IF (vars_defined) zx_tmp_fi2d(:) = lwup0p(:,klevp1)
1987 CALL histwrite_phy(o_rlutcs4co2, zx_tmp_fi2d)
1988 IF (vars_defined) THEN
1989 DO k=1, klevp1
1990 zx_tmp_fi3d1(:,k)=swupp(:,k)*swradcorr(:)
1991 ENDDO
1992 ENDIF
1993 CALL histwrite_phy(o_rsu4co2, zx_tmp_fi3d1)
1994 IF (vars_defined) THEN
1995 DO k=1, klevp1
1996 zx_tmp_fi3d1(:,k)=swup0p(:,k)*swradcorr(:)
1997 ENDDO
1998 ENDIF
1999 CALL histwrite_phy(o_rsucs4co2, zx_tmp_fi3d1)
2000 IF (vars_defined) THEN
2001 DO k=1, klevp1
2002 zx_tmp_fi3d1(:,k)=swdnp(:,k)*swradcorr(:)
2003 ENDDO
2004 ENDIF
2005 CALL histwrite_phy(o_rsd4co2, zx_tmp_fi3d1)
2006 IF (vars_defined) THEN
2007 DO k=1, klevp1
2008 zx_tmp_fi3d1(:,k)=swdn0p(:,k)*swradcorr(:)
2009 ENDDO
2010 ENDIF
2011 CALL histwrite_phy(o_rsdcs4co2, zx_tmp_fi3d1)
2012 CALL histwrite_phy(o_rlu4co2, lwupp)
2013 CALL histwrite_phy(o_rlucs4co2, lwup0p)
2014 CALL histwrite_phy(o_rld4co2, lwdnp)
2015 CALL histwrite_phy(o_rldcs4co2, lwdn0p)
2016 ENDIF !ok_4xCO2atm
2017 !!!!!!!!!!!! Sorties niveaux de pression NMC !!!!!!!!!!!!!!!!!!!!
2018
1/2
✓ Branch 0 taken 481 times.
✗ Branch 1 not taken.
481 IF (.NOT.ok_all_xml) THEN
2019 ! ATTENTION, LES ANCIENS HISTWRITE ONT ETES CONSERVES EN ATTENDANT MIEUX:
2020 ! Champs interpolles sur des niveaux de pression
2021
2/2
✓ Branch 0 taken 1443 times.
✓ Branch 1 taken 481 times.
1924 DO iff=7, nfiles-1 !--OB: here we deal with files 7,8,9
2022
2023 1443 CALL histwrite_phy(o_tnondef,tnondef(:,:,iff-6),iff)
2024 1443 CALL histwrite_phy(o_ta,twriteSTD(:,:,iff-6),iff)
2025 1443 CALL histwrite_phy(o_zg,phiwriteSTD(:,:,iff-6),iff)
2026 1443 CALL histwrite_phy(o_hus,qwriteSTD(:,:,iff-6),iff)
2027 1443 CALL histwrite_phy(o_hur,rhwriteSTD(:,:,iff-6),iff)
2028 1443 CALL histwrite_phy(o_ua,uwriteSTD(:,:,iff-6),iff)
2029 1443 CALL histwrite_phy(o_va,vwriteSTD(:,:,iff-6),iff)
2030 1443 CALL histwrite_phy(o_wap,wwriteSTD(:,:,iff-6),iff)
2031
2/2
✓ Branch 0 taken 1440 times.
✓ Branch 1 taken 3 times.
1443 IF (vars_defined) THEN
2032
2/2
✓ Branch 0 taken 24480 times.
✓ Branch 1 taken 1440 times.
25920 DO k=1, nlevSTD
2033
2/2
✓ Branch 0 taken 24333120 times.
✓ Branch 1 taken 24480 times.
24359040 DO i=1, klon
2034
1/2
✓ Branch 0 taken 24333120 times.
✗ Branch 1 not taken.
24357600 IF (tnondef(i,k,iff-6).NE.missing_val) THEN
2035
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 24333120 times.
24333120 IF (freq_outNMC(iff-6).LT.0) THEN
2036 freq_moyNMC(iff-6)=(mth_len*un_jour)/freq_calNMC(iff-6)
2037 ELSE
2038 24333120 freq_moyNMC(iff-6)=freq_outNMC(iff-6)/freq_calNMC(iff-6)
2039 ENDIF
2040 24333120 zx_tmp_fi3d_STD(i,k) = (100.*tnondef(i,k,iff-6))/freq_moyNMC(iff-6)
2041 ELSE
2042 zx_tmp_fi3d_STD(i,k) = missing_val
2043 ENDIF
2044 ENDDO
2045 ENDDO
2046 ENDIF
2047 1443 CALL histwrite_phy(o_psbg,zx_tmp_fi3d_STD,iff)
2048
2/2
✓ Branch 0 taken 1440 times.
✓ Branch 1 taken 3 times.
1443 IF (vars_defined) THEN
2049
2/2
✓ Branch 0 taken 24480 times.
✓ Branch 1 taken 1440 times.
25920 DO k=1, nlevSTD
2050
2/2
✓ Branch 0 taken 24333120 times.
✓ Branch 1 taken 24480 times.
24359040 DO i=1, klon
2051
2/2
✓ Branch 0 taken 24298898 times.
✓ Branch 1 taken 34222 times.
24357600 IF (O3sumSTD(i,k,iff-6).NE.missing_val) THEN
2052 24298898 zx_tmp_fi3d_STD(i,k) = O3sumSTD(i,k,iff-6) * 1.e+9
2053 ELSE
2054 34222 zx_tmp_fi3d_STD(i,k) = missing_val
2055 ENDIF
2056 ENDDO
2057 ENDDO !k=1, nlevSTD
2058 ENDIF
2059 1443 CALL histwrite_phy(o_tro3,zx_tmp_fi3d_STD,iff)
2060
1/2
✗ Branch 0 not taken.
✓ Branch 1 taken 1443 times.
1443 IF (read_climoz == 2) THEN
2061 IF (vars_defined) THEN
2062 DO k=1, nlevSTD
2063 DO i=1, klon
2064 IF (O3daysumSTD(i,k,iff-6).NE.missing_val) THEN
2065 zx_tmp_fi3d_STD(i,k) = O3daysumSTD(i,k,iff-6) * 1.e+9
2066 ELSE
2067 zx_tmp_fi3d_STD(i,k) = missing_val
2068 ENDIF
2069 ENDDO
2070 ENDDO !k=1, nlevSTD
2071 ENDIF
2072