Ignore:
Timestamp:
Jun 29, 2018, 12:31:11 PM (7 years ago)
Author:
Laurent Fairhead
Message:

First attempt at merging with trunk

Location:
LMDZ6/branches/DYNAMICO-conv
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ6/branches/DYNAMICO-conv

  • LMDZ6/branches/DYNAMICO-conv/libf/phylmd/phys_output_write_mod.F90

    r3355 r3356  
    1414CONTAINS
    1515
    16   ! ug Routine pour définir (los du premier passageà) ET sortir les variables
     16  ! ug Routine pour définir (lors du premier passage) ET sortir les variables
    1717  SUBROUTINE phys_output_write(itap, pdtphys, paprs, pphis, &
    1818       pplay, lmax_th, aerosol_couple,         &
    … …  
    3131    USE vertical_layers_mod, ONLY : ap, bp, aps, bps
    3232    USE phys_output_ctrlout_mod, ONLY: o_phis, o_aire, is_ter, is_lic, is_oce, &
     33         o_longitude, o_latitude, &
    3334         o_Ahyb, o_Bhyb,o_Ahyb_inter, o_Bhyb_inter, &
    3435         is_ave, is_sic, o_contfracATM, o_contfracOR, &
    … …  
    3940         o_wind10m, o_wind10max, o_gusts, o_sicf, &
    4041         o_psol, o_mass, o_qsurf, o_qsol, &
    41          o_precip, o_ndayrain, o_plul, o_pluc, o_plun, &
     42         o_precip, o_rain_fall, o_rain_con, o_ndayrain, o_plul, o_pluc, o_plun, &
    4243         o_snow, o_msnow, o_fsnow, o_evap, o_ep,o_epmax_diag, & ! epmax_cape
    4344         o_tops, o_tops0, o_topl, o_topl0, &
    44          o_SWupTOA, o_SWupTOAclr, o_SWdnTOA, &
     45         o_SWupTOA, o_SWupTOAclr, o_SWupTOAcleanclr, o_SWdnTOA, &
    4546         o_SWdnTOAclr, o_nettop, o_SWup200, &
    4647         o_SWup200clr, o_SWdn200, o_SWdn200clr, &
    … …  
    4849         o_LWdn200clr, o_sols, o_sols0, &
    4950         o_soll, o_radsol, o_soll0, o_SWupSFC, &
    50          o_SWupSFCclr, o_SWdnSFC, o_SWdnSFCclr, &
     51         o_SWupSFCclr, o_SWupSFCcleanclr, o_SWdnSFC, o_SWdnSFCclr, o_SWdnSFCcleanclr, &
    5152         o_LWupSFC, o_LWdnSFC, o_LWupSFCclr, &
    52          o_LWdnSFCclr, o_bils, o_bils_diss, &
     53         o_LWdnSFCclr, o_LWupTOAcleanclr, o_LWdnSFCcleanclr, o_bils, o_bils_diss, &
    5354         o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &
    5455         o_bils_latent, o_bils_enthalp, o_sens, &
    55          o_fder, o_ffonte, o_fqcalving, o_fqfonte, o_runofflic, &
     56         o_fder, o_ffonte, o_fqcalving, o_fqfonte, o_mrroli, o_runofflic, &
    5657         o_taux, o_tauy, o_snowsrf, o_qsnow, &
    5758         o_snowhgt, o_toice, o_sissnow, o_runoff, &
    … …  
    6970         o_cldhjn, o_cldtjn, o_cldq, o_lwp, o_iwp, &
    7071         o_ue, o_ve, o_uq, o_vq, o_cape, o_pbase, &
     72         o_uwat, o_vwat, &
    7173         o_ptop, o_fbase, o_plcl, o_plfc, &
    7274         o_wbeff, o_convoccur, o_cape_max, o_upwd, o_ep,o_epmax_diag, o_Ma, &
    73          o_dnwd, o_dnwd0, o_ftime_con, o_mc, &
     75         o_dnwd, o_dnwd0, o_ftime_deepcv, o_ftime_con, o_mc, &
    7476         o_prw, o_prlw, o_prsw, o_s_pblh, o_s_pblt, o_s_lcl, &
    7577         o_s_therm, o_uSTDlevs, o_vSTDlevs, &
    … …  
    8284         o_cdragh_x   , o_cdragh_w   , o_cdragm_x   , o_cdragm_w   , &
    8385         o_kh         , o_kh_x       , o_kh_w       , &
    84          o_ale, o_alp, o_cin, o_WAPE, o_wake_h, &
     86         o_ale, o_alp, o_cin, o_WAPE, o_wake_h, o_cv_gen, o_wake_dens, &
    8587         o_wake_s, o_wake_deltat, o_wake_deltaq, &
    8688         o_wake_omg, o_dtwak, o_dqwak, o_dqwak2d, o_Vprecip, &
    … …  
    9799         o_cldemi, o_pr_con_l, o_pr_con_i, &
    98100         o_pr_lsc_l, o_pr_lsc_i, o_re, o_fl, &
    99          o_rh2m, o_rh2m_min, o_rh2m_max, &
     101         o_rh2m, &
    100102         o_qsat2m, o_tpot, o_tpote, o_SWnetOR, &
    101          o_SWdownOR, o_LWdownOR, o_snowl, &
     103         o_LWdownOR, o_snowl, &
    102104         o_solldown, o_dtsvdfo, o_dtsvdft, &
    103105         o_dtsvdfg, o_dtsvdfi, o_z0m, o_z0h,  o_od443aer, o_od550aer, &
    … …  
    119121         o_swtoacf_zero, o_swsrfcf_zero, &
    120122         o_topswai, o_solswai, o_toplwai, o_sollwai, o_scdnc, &
    121          o_cldncl, o_reffclws, o_reffclwc, o_stratomask,&
     123         o_cldncl, o_reffclws, o_reffclwc, o_solbnd, o_stratomask,&
    122124         o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
    123          o_lcc3dstra, o_reffclwtop, o_ec550aer, &
     125         o_lcc3dstra, o_icc3dcon, o_icc3dstra, &
     126         o_cldicemxrat, o_cldwatmxrat, o_reffclwtop, o_ec550aer, &
    124127         o_lwcon, o_iwcon, o_temp, o_theta, &
    125128         o_ovapinit, o_ovap, o_oliq, o_ocond, o_geop, &
    … …  
    153156         o_du_gwd_front, o_dv_gwd_front, &
    154157         o_east_gwstress, o_west_gwstress, &
    155          o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, &
     158         o_rsd, o_rlu, o_rld, o_rsucs, o_rsdcs, o_rsucsaf, o_rsdcsaf, &
    156159         o_rlucs, o_rldcs, o_tnt, o_tntr, &
    157160         o_tntscpbl, o_tnhus, o_tnhusscpbl, &
    … …  
    166169         o_uxv, o_vxq, o_vxT, o_wxq, o_vxphi, &
    167170         o_wxT, o_uxu, o_vxv, o_TxT, o_trac, &
     171#ifdef REPROBUS
     172         o_nas, &
     173#endif
    168174         o_dtr_vdf, o_dtr_the, o_dtr_con, &
    169175         o_dtr_lessi_impa, o_dtr_lessi_nucl, &
    … …  
    177183         o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
    178184         o_lat_prec_sol_oce, o_lat_prec_sol_sic, &
     185         o_sza, &
    179186! Marine
    180187         o_map_prop_hc, o_map_prop_hist, o_map_emis_hc, o_map_iwp_hc, &
    … …  
    197204         o_budg_sed_part, o_R2SO4, o_OCS_lifetime, o_SO2_lifetime, &
    198205         o_budg_3D_backgr_ocs, o_budg_3D_backgr_so2, &
     206         o_budg_dep_dry_ocs, o_budg_dep_wet_ocs, &
     207         o_budg_dep_dry_so2, o_budg_dep_wet_so2, &
    199208         o_budg_dep_dry_h2so4, o_budg_dep_wet_h2so4, &
    200209         o_budg_dep_dry_part, o_budg_dep_wet_part, &
    … …  
    209218         nday_rain, rain_con, snow_con, &
    210219         topsw, toplw, toplw0, swup, swdn, &
    211          topsw0, swup0, swdn0, SWup200, SWup200clr, &
     220         topsw0, swupc0, swdnc0, swup0, swdn0, SWup200, SWup200clr, &
    212221         SWdn200, SWdn200clr, LWup200, LWup200clr, &
    213222         LWdn200, LWdn200clr, solsw, solsw0, sollw, &
    214223         radsol, swradcorr, sollw0, sollwdown, sollw, gustiness, &
    215          sollwdownclr, lwdn0, ftsol, ustar, u10m, &
     224         sollwdownclr, lwdnc0, lwdn0, ftsol, ustar, u10m, &
    216225         v10m, pbl_tke, wake_delta_pbl_TKE, &
    217226         delta_tsurf, &
    218227         wstar, cape, ema_pcb, ema_pct, &
    219228         ema_cbmf, Ma, fm_therm, ale_bl, alp_bl, ale, &
    220          alp, cin, wake_pe, wake_s, wake_deltat, &
     229         alp, cin, wake_pe, wake_dens, wake_s, wake_deltat, &
    221230         wake_deltaq, ftd, fqd, ale_bl_trig, albsol1, &
     231         ale_wake, ale_bl_stat, &
    222232         rnebcon, wo, falb1, albsol2, coefh, clwcon0, &
    223233         ratqs, entr_therm, zqasc, detr_therm, f0, &
    224          lwup, lwdn, lwup0, coefm, &
    225          swupp, lwupp, swup0p, lwup0p, swdnp, lwdnp, &
    226          swdn0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
     234         lwup, lwdn, lwupc0, lwup0, coefm, &
     235         swupp, lwupp, swupc0p, swup0p, lwupc0p, lwup0p, swdnp, lwdnp, &
     236         swdnc0p, swdn0p, lwdnc0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
    227237         vqsumSTD, vTsumSTD, O3daysumSTD, wqsumSTD, &
    228238         vphisumSTD, wTsumSTD, u2sumSTD, v2sumSTD, &
    … …  
    245255         cldh, cldt, JrNt, cldljn, cldmjn, cldhjn, &
    246256         cldtjn, cldq, flwp, fiwp, ue, ve, uq, vq, &
     257         uwat, vwat, &
    247258         plcl, plfc, wbeff, convoccur, upwd, dnwd, dnwd0, prw, prlw, prsw, &
    248259         s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
    249260         vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &
    250          twriteSTD, ale_wake, alp_wake, &
    251          dtvdf_x    ,dtvdf_w    ,dqvdf_x    ,dqvdf_w    , &
     261         twriteSTD, alp_wake, &
     262!!         dtvdf_x    ,dtvdf_w    ,dqvdf_x    ,dqvdf_w    , &
     263         d_t_vdf_x    ,d_t_vdf_w    ,d_q_vdf_x    ,d_q_vdf_w    , &
    252264         sens_x     ,sens_w     ,zxfluxlat_x,zxfluxlat_w, &
    253265         cdragh_x   ,cdragh_w   ,cdragm_x   ,cdragm_w   , &
    254266         kh         ,kh_x       ,kh_w       , &
    255          wake_h, &
     267         cv_gen, wake_h, &
    256268         wake_omg, d_t_wake, d_q_wake, Vprecip, &
    257269         wdtrainA, wdtrainM, n2, s2, proba_notrig, &
    258          random_notrig, ale_bl_stat, &
     270         random_notrig, &
    259271         alp_bl_det, alp_bl_fluct_m, alp_bl_conv, &
    260272         alp_bl_stat, alp_bl_fluct_tke, slab_wfbils, &
    … …  
    273285         sollwad0_aero, toplwai_aero, sollwai_aero, &
    274286         scdnc, cldncl, reffclws, reffclwc, cldnvi, stratomask,&
    275          lcc, lcc3d, lcc3dcon, lcc3dstra, reffclwtop, &
     287         lcc, lcc3d, lcc3dcon, lcc3dstra, &
     288         icc3dcon, icc3dstra, zfice, reffclwtop, &
    276289         ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
    277290         ql_seri, qs_seri, tr_seri, &
    … …  
    301314         budg_sed_part, R2SO4, OCS_lifetime, SO2_lifetime, &
    302315         budg_3D_backgr_ocs, budg_3D_backgr_so2, &
     316         budg_dep_dry_ocs, budg_dep_wet_ocs, &
     317         budg_dep_dry_so2, budg_dep_wet_so2, &
    303318         budg_dep_dry_h2so4, budg_dep_wet_h2so4, &
    304319         budg_dep_dry_part, budg_dep_wet_part, &
    … …  
    309324#endif
    310325
     326#ifdef REPROBUS
     327    USE CHEM_REP, ONLY : nas, nbnas, tnamenas, ttextnas
     328#endif
     329
    311330    USE phys_output_var_mod, ONLY: vars_defined, snow_o, zfra_o, bils_diss, &
    312331         bils_ec,bils_ech, bils_tke, bils_kinetic, bils_latent, bils_enthalp, &
    313332         itau_con, nfiles, clef_files, nid_files, dryaod_diag, &
    314333         zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
    315          zustr_gwd_front, zvstr_gwd_front,     &
     334         zustr_gwd_front, zvstr_gwd_front, sza_o,    &
    316335         sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o, &
    317336         cloudth_sth,cloudth_senv,cloudth_sigmath,cloudth_sigmaenv, &
    … …  
    336355    USE indice_sol_mod, ONLY: nbsrf
    337356    USE infotrac_phy, ONLY: nqtot, nqo, type_trac, tname, niadv
    338     USE geometry_mod, ONLY: cell_area
     357    USE geometry_mod, ONLY: cell_area, latitude_deg, longitude_deg
    339358    USE surface_data, ONLY: type_ocean, version_ocean, ok_veget, ok_snow
    340359    USE aero_mod, ONLY: naero_tot, id_STRAT_phy
    … …  
    386405    INTEGER :: itau_w
    387406    INTEGER :: i, iinit, iinitend=1, iff, iq, iiq, nsrf, k, ll, naero
    388     REAL, DIMENSION (klon) :: zx_tmp_fi2d
     407    REAL, DIMENSION (klon) :: zx_tmp_fi2d, zpt_conv2d
    389408    REAL, DIMENSION (klon,klev) :: zx_tmp_fi3d, zpt_conv
    390409    REAL, DIMENSION (klon,klev+1) :: zx_tmp_fi3d1
     410    REAL, DIMENSION (klon,NSW) :: zx_tmp_fi3dsp
    391411    CHARACTER (LEN=4)              :: bb2
    392412    INTEGER, DIMENSION(nbp_lon*nbp_lat)  :: ndex2d
    … …  
    550570       CALL histwrite_phy(o_Ahyb_inter, Ahyb_bounds)
    551571       CALL histwrite_phy(o_Bhyb_inter, Bhyb_bounds)
     572       CALL histwrite_phy(o_longitude, longitude_deg)
     573       CALL histwrite_phy(o_latitude, latitude_deg)
    552574!
    553575#ifdef CPP_RRTM
    … …  
    609631     ENDIF
    610632
     633       CALL histwrite_phy(o_sza, sza_o)
    611634       CALL histwrite_phy(o_flat, zxfluxlat)
    612635       CALL histwrite_phy(o_ptstar, ptstar)
    … …  
    667690
    668691       CALL histwrite_phy(o_precip, zx_tmp_fi2d)
     692       CALL histwrite_phy(o_rain_fall, rain_fall)
    669693       CALL histwrite_phy(o_ndayrain, nday_rain)
    670694
    … …  
    687711       ENDIF
    688712       CALL histwrite_phy(o_pluc, zx_tmp_fi2d)
     713       CALL histwrite_phy(o_rain_con, rain_con)
    689714       CALL histwrite_phy(o_snow, snow_fall)
    690715       CALL histwrite_phy(o_msnow, zxsnow)
    … …  
    714739       ENDIF
    715740       CALL histwrite_phy(o_SWupTOAclr, zx_tmp_fi2d)
     741
     742       IF (vars_defined) THEN
     743          zx_tmp_fi2d(:) = swupc0(:,klevp1)*swradcorr(:)
     744       ENDIF
     745       CALL histwrite_phy(o_SWupTOAcleanclr, zx_tmp_fi2d)
    716746
    717747       IF (vars_defined) THEN
    … …  
    781811
    782812       IF (vars_defined) THEN
     813          zx_tmp_fi2d(:) = swupc0(:,1)*swradcorr(:)
     814       ENDIF
     815       CALL histwrite_phy(o_SWupSFCcleanclr, zx_tmp_fi2d)
     816
     817       IF (vars_defined) THEN
    783818          zx_tmp_fi2d(:) = swdn(:,1)*swradcorr(:)
    784819       ENDIF
    … …  
    789824       ENDIF
    790825       CALL histwrite_phy(o_SWdnSFCclr, zx_tmp_fi2d)
     826
     827       IF (vars_defined) THEN
     828          zx_tmp_fi2d(:) = swdnc0(:,1)*swradcorr(:)
     829       ENDIF
     830       CALL histwrite_phy(o_SWdnSFCcleanclr, zx_tmp_fi2d)
    791831
    792832       IF (vars_defined) THEN
    … …  
    802842       CALL histwrite_phy(o_LWupSFCclr, zx_tmp_fi2d)
    803843       CALL histwrite_phy(o_LWdnSFCclr, sollwdownclr)
     844
     845       IF (vars_defined) THEN
     846          zx_tmp_fi2d(:) = lwupc0(:,klevp1)
     847       ENDIF
     848       CALL histwrite_phy(o_LWupTOAcleanclr, zx_tmp_fi2d)
     849       IF (vars_defined) THEN
     850          zx_tmp_fi2d(:) = -1.*lwdnc0(:,1)
     851       ENDIF
     852       CALL histwrite_phy(o_LWdnSFCcleanclr, zx_tmp_fi2d)
     853
    804854       CALL histwrite_phy(o_bils, bils)
    805855       CALL histwrite_phy(o_bils_diss, bils_diss)
    … …  
    819869       CALL histwrite_phy(o_fqcalving, zxfqcalving)
    820870       CALL histwrite_phy(o_fqfonte, zxfqfonte)
     871       IF (vars_defined) THEN
     872          zx_tmp_fi2d(1:klon)=(zxfqfonte(1:klon)+rain_fall(1:klon))*pctsrf(1:klon,is_lic)
     873       ENDIF
     874       CALL histwrite_phy(o_mrroli, zx_tmp_fi2d)
    821875       CALL histwrite_phy(o_runofflic, zxrunofflic)
    822876       IF (vars_defined) THEN
    … …  
    9581012       CALL histwrite_phy(o_uq, uq)
    9591013       CALL histwrite_phy(o_vq, vq)
     1014       CALL histwrite_phy(o_uwat, uwat)
     1015       CALL histwrite_phy(o_vwat, vwat)
    9601016       IF (iflag_con.GE.3) THEN ! sb
    9611017          CALL histwrite_phy(o_cape, cape)
    … …  
    9761032          CALL histwrite_phy(o_dnwd, dnwd)
    9771033          CALL histwrite_phy(o_dnwd0, dnwd0)
    978           IF (vars_defined)         zx_tmp_fi2d=float(itau_con)/float(itap)
    979           CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
     1034          !! The part relative to the frequency of occurence of convection
     1035          !! is now grouped with the part relative to thermals and shallow
     1036          !! convection (output of the 3 fields: ftime_deepcv, ftime_th and
     1037          !!  ftime_con).
    9801038          IF (vars_defined) THEN
    9811039             IF (iflag_thermals>=1)THEN
    982                 zx_tmp_fi3d=dnwd+dnwd0+upwd+fm_therm(:,1:klev)
     1040                zx_tmp_fi3d=-dnwd+dnwd0+upwd+fm_therm(:,1:klev)
    9831041             ELSE
    984                 zx_tmp_fi3d=dnwd+dnwd0+upwd
     1042                zx_tmp_fi3d=-dnwd+dnwd0+upwd
    9851043             ENDIF
    9861044          ENDIF
    … …  
    10991157             CALL histwrite_phy(o_alp_wk, alp_wake)
    11001158             IF (iflag_pbl_split>=1) THEN
    1101                !ym warning : dtvdf_x, dtvdf_w, dqvdf_x, dqvdf_w are not defined by model
    1102                !ym => init to 0 in phys_local_var_mod
    1103                IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_x(1:klon,1:klev)/pdtphys
    1104                CALL histwrite_phy(o_dtvdf_x    ,zx_tmp_fi3d)
    1105                IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dtvdf_w(1:klon,1:klev)/pdtphys
    1106                CALL histwrite_phy(o_dtvdf_w    ,zx_tmp_fi3d)
    1107                IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_x(1:klon,1:klev)/pdtphys
    1108                CALL histwrite_phy(o_dqvdf_x    ,zx_tmp_fi3d)
    1109                IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=dqvdf_w(1:klon,1:klev)/pdtphys
     1159               IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf_x(1:klon,1:klev)/pdtphys
     1160               CALL histwrite_phy(o_dtvdf_x    ,zx_tmp_fi3d)
     1161               IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf_w(1:klon,1:klev)/pdtphys
     1162               CALL histwrite_phy(o_dtvdf_w    ,zx_tmp_fi3d)
     1163               IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf_x(1:klon,1:klev)/pdtphys
     1164               CALL histwrite_phy(o_dqvdf_x    ,zx_tmp_fi3d)
     1165               IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_vdf_w(1:klon,1:klev)/pdtphys
     1166
     1167!
    11101168               CALL histwrite_phy(o_dqvdf_w    ,zx_tmp_fi3d)
    11111169               CALL histwrite_phy(o_sens_x     ,sens_x     )
    … …  
    11261184             CALL histwrite_phy(o_cin, cin)
    11271185             CALL histwrite_phy(o_WAPE, wake_pe)
     1186             CALL histwrite_phy(o_cv_gen, cv_gen)
    11281187             CALL histwrite_phy(o_wake_h, wake_h)
     1188             CALL histwrite_phy(o_wake_dens, wake_dens)
    11291189             CALL histwrite_phy(o_wake_s, wake_s)
    11301190             CALL histwrite_phy(o_wake_deltat, wake_deltat)
    … …  
    12211281       CALL histwrite_phy(o_rh2m, zx_tmp_fi2d)
    12221282
    1223        IF (vars_defined) THEN
    1224           DO i=1, klon
    1225              zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
    1226           ENDDO
    1227        ENDIF
    1228        CALL histwrite_phy(o_rh2m_min, zx_tmp_fi2d)
    1229 
    1230        IF (vars_defined) THEN
    1231           DO i=1, klon
    1232              zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
    1233           ENDDO
    1234        ENDIF
    1235        CALL histwrite_phy(o_rh2m_max, zx_tmp_fi2d)
     1283!       IF (vars_defined) THEN
     1284!          DO i=1, klon
     1285!             zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
     1286!          ENDDO
     1287!       ENDIF
     1288!       CALL histwrite_phy(o_rh2m_min, zx_tmp_fi2d)
     1289
     1290!       IF (vars_defined) THEN
     1291!          DO i=1, klon
     1292!             zx_tmp_fi2d(i)=MIN(100.,rh2m(i)*100.)
     1293!          ENDDO
     1294!       ENDIF
     1295!       CALL histwrite_phy(o_rh2m_max, zx_tmp_fi2d)
    12361296
    12371297       CALL histwrite_phy(o_qsat2m, qsat2m)
    … …  
    12401300       IF (vars_defined) zx_tmp_fi2d(1 : klon) = fsolsw( 1 : klon, is_ter)
    12411301       CALL histwrite_phy(o_SWnetOR,  zx_tmp_fi2d)
    1242        IF (vars_defined) zx_tmp_fi2d(1:klon) = solsw(1:klon)/(1.-albsol1(1:klon))
    1243        CALL histwrite_phy(o_SWdownOR,  zx_tmp_fi2d)
    12441302       CALL histwrite_phy(o_LWdownOR, sollwdown)
    12451303       CALL histwrite_phy(o_snowl, snow_lsc)
    … …  
    13301388          CALL histwrite_phy(o_budg_3D_nucl,         budg_3D_nucl)
    13311389          CALL histwrite_phy(o_budg_3D_cond_evap,    budg_3D_cond_evap)
     1390          CALL histwrite_phy(o_budg_dep_dry_ocs,     budg_dep_dry_ocs)
     1391          CALL histwrite_phy(o_budg_dep_wet_ocs,     budg_dep_wet_ocs)
     1392          CALL histwrite_phy(o_budg_dep_dry_so2,     budg_dep_dry_so2)
     1393          CALL histwrite_phy(o_budg_dep_wet_so2,     budg_dep_wet_so2)
    13321394          CALL histwrite_phy(o_budg_dep_dry_h2so4,   budg_dep_dry_h2so4)
    13331395          CALL histwrite_phy(o_budg_dep_wet_h2so4,   budg_dep_wet_h2so4)
    … …  
    14351497          CALL histwrite_phy(o_lcc3dcon, lcc3dcon)
    14361498          CALL histwrite_phy(o_lcc3dstra, lcc3dstra)
     1499          CALL histwrite_phy(o_icc3dcon, icc3dcon)
     1500          CALL histwrite_phy(o_icc3dstra, icc3dstra)
     1501          CALL histwrite_phy(o_cldicemxrat, zfice)
     1502          zx_tmp_fi3d(:,:)=1-zfice(:,:)
     1503          CALL histwrite_phy(o_cldwatmxrat, zx_tmp_fi3d)
    14371504          CALL histwrite_phy(o_reffclwtop, reffclwtop)
    14381505       ENDIF
    … …  
    14621529       CALL histwrite_phy(o_zfull,zx_tmp_fi3d)
    14631530
     1531#ifdef CPP_XIOS
     1532!solbnd begin
     1533#ifdef CPP_RRTM
     1534      IF (iflag_rrtm.EQ.1) THEN
     1535       IF (vars_defined) THEN
     1536        DO ISW=1, NSW
     1537          zx_tmp_fi3dsp(:,ISW) = swdn(:,klevp1)*swradcorr(:)*RSUN(ISW)
     1538        ENDDO
     1539        CALL histwrite_phy(o_solbnd, zx_tmp_fi3dsp)
     1540       ENDIF
     1541      ENDIF
     1542#endif
     1543!solbnd end
     1544#endif
     1545
    14641546       IF (flag_aerosol_strat.EQ.2) THEN
    14651547         CALL histwrite_phy(o_stratomask, stratomask)
    … …  
    14701552        DO k = 2, klev
    14711553         DO i = 1, klon
    1472             zx_tmp_fi3d(i,k) = zphi(i,k)/RG + &
     1554            zx_tmp_fi3d(i,k) = zphi(i,k-1)/RG + &
    14731555                          (zphi(i,k)-zphi(i,k-1))/RG * &
    1474                           (paprs(i,k)-pplay(i,k))/(pplay(i,k)-pplay(i,k-1))
     1556                          (paprs(i,k)-pplay(i,k-1))/(pplay(i,k)-pplay(i,k-1))
    14751557         ENDDO
    14761558        ENDDO
    … …  
    16421724          CALL histwrite_phy(o_plulst, plul_st)
    16431725          IF (vars_defined) THEN
     1726             do i=1,klon
     1727                zx_tmp_fi2d(1:klon)=lmax_th(:)
     1728             enddo
     1729          ENDIF
     1730          CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d)
     1731          IF (vars_defined) THEN
    16441732             DO k=1,klev
    16451733                DO i=1,klon
    … …  
    16531741          ENDIF
    16541742          CALL histwrite_phy(o_ptconvth, zx_tmp_fi3d)
    1655           IF (vars_defined) THEN
    1656              do i=1,klon
    1657                 zx_tmp_fi2d(1:klon)=lmax_th(:)
    1658              enddo
    1659           ENDIF
    1660           CALL histwrite_phy(o_lmaxth, zx_tmp_fi2d)
    16611743       ENDIF ! iflag_thermals>=1
     1744!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
     1745       zpt_conv = 0.
     1746       WHERE (ptconv) zpt_conv = 1.
     1747       CALL histwrite_phy(o_ptconv, zpt_conv)
     1748!!       IF (vars_defined)         zx_tmp_fi2d=float(itau_con)/float(itap)
     1749!!       CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
     1750       IF (vars_defined) THEN
     1751          zpt_conv2d(:) = 0.
     1752          DO k=1,klev
     1753            WHERE (ptconv(:,k)) zpt_conv2d(:) = 1.
     1754          ENDDO
     1755       ENDIF
     1756       CALL histwrite_phy(o_ftime_deepcv, zpt_conv2d)
     1757       IF (vars_defined) THEN
     1758          zx_tmp_fi2d(:) = 0.
     1759          DO k=1,klev
     1760            WHERE (ptconvth(:,k)) zx_tmp_fi2d(:) = 1.
     1761          ENDDO
     1762       ENDIF
     1763       CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
     1764       IF (vars_defined) THEN
     1765           zx_tmp_fi2d(:) = max(zx_tmp_fi2d(:),zpt_conv2d(:))
     1766       ENDIF
     1767       CALL histwrite_phy(o_ftime_con, zx_tmp_fi2d)
    16621768!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    16631769       IF (vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_vdf(1:klon,1:klev)/pdtphys
    … …  
    16751781       IF (vars_defined) CALL water_int(klon,klev,zx_tmp_fi3d,zmasse,zx_tmp_fi2d)
    16761782       CALL histwrite_phy(o_dqeva2d, zx_tmp_fi2d)
    1677        zpt_conv = 0.
    1678        WHERE (ptconv) zpt_conv = 1.
    1679        CALL histwrite_phy(o_ptconv, zpt_conv)
    16801783       CALL histwrite_phy(o_ratqs, ratqs)
    16811784       IF (vars_defined) THEN
    … …  
    16951798       IF (iflag_thermals>=1) THEN
    16961799          ! Pour l instant 0 a y reflichir pour les thermiques
    1697           zx_tmp_fi2d=0.
    1698           CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
     1800          ! regroupe avec ftime_deepcv et ftime_con
     1801          !!zx_tmp_fi2d=0.
     1802          !!CALL histwrite_phy(o_ftime_th, zx_tmp_fi2d)
    16991803          CALL histwrite_phy(o_f_th, fm_therm)
    17001804          CALL histwrite_phy(o_e_th, entr_therm)
    … …  
    18281932       ENDIF
    18291933       CALL histwrite_phy(o_rsdcs, zx_tmp_fi3d1)
     1934       DO k=1, klevp1
     1935         zx_tmp_fi3d1(:,k)=swdnc0(:,k)*swradcorr(:)
     1936       ENDDO
     1937       CALL histwrite_phy(o_rsdcsaf, zx_tmp_fi3d1)
    18301938
    18311939       CALL histwrite_phy(o_rlu, lwup)
    … …  
    20622170            IF (type_trac == 'lmdz' .OR. type_trac == 'repr' .OR. type_trac == 'coag') THEN
    20632171             !--3D fields
    2064 !             CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
    2065 !             CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
    2066 !             CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
    2067 !             CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))
    2068 !             CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
    2069 !             CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
    2070 !             CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))
    2071 !             CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))
    2072 !             CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))
    2073 !             CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
    2074 !             CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
    2075 !             CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
    2076 !             CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
    2077 !             CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))
     2172             CALL histwrite_phy(o_trac(iq-nqo), tr_seri(:,:,iq-nqo))
     2173             CALL histwrite_phy(o_dtr_vdf(iq-nqo),d_tr_cl(:,:,iq-nqo))
     2174             CALL histwrite_phy(o_dtr_the(iq-nqo),d_tr_th(:,:,iq-nqo))
     2175             CALL histwrite_phy(o_dtr_con(iq-nqo),d_tr_cv(:,:,iq-nqo))
     2176             CALL histwrite_phy(o_dtr_lessi_impa(iq-nqo),d_tr_lessi_impa(:,:,iq-nqo))
     2177             CALL histwrite_phy(o_dtr_lessi_nucl(iq-nqo),d_tr_lessi_nucl(:,:,iq-nqo))
     2178             CALL histwrite_phy(o_dtr_insc(iq-nqo),d_tr_insc(:,:,iq-nqo))
     2179             CALL histwrite_phy(o_dtr_bcscav(iq-nqo),d_tr_bcscav(:,:,iq-nqo))
     2180             CALL histwrite_phy(o_dtr_evapls(iq-nqo),d_tr_evapls(:,:,iq-nqo))
     2181             CALL histwrite_phy(o_dtr_ls(iq-nqo),d_tr_ls(:,:,iq-nqo))
     2182             CALL histwrite_phy(o_dtr_trsp(iq-nqo),d_tr_trsp(:,:,iq-nqo))
     2183             CALL histwrite_phy(o_dtr_sscav(iq-nqo),d_tr_sscav(:,:,iq-nqo))
     2184             CALL histwrite_phy(o_dtr_sat(iq-nqo),d_tr_sat(:,:,iq-nqo))
     2185             CALL histwrite_phy(o_dtr_uscav(iq-nqo),d_tr_uscav(:,:,iq-nqo))
    20782186             !--2D fields
    2079 !             CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo))
     2187             CALL histwrite_phy(o_dtr_dry(iq-nqo), flux_tr_dry(:,iq-nqo))
    20802188             zx_tmp_fi2d=0.
    20812189             IF (vars_defined) THEN
    … …  
    20842192                ENDDO
    20852193             ENDIF
    2086 !             CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
     2194#ifndef REPROBUS
     2195             CALL histwrite_phy(o_trac_cum(iq-nqo), zx_tmp_fi2d)
     2196#endif
    20872197            ENDIF
    20882198          ENDDO
    20892199       ENDIF
     2200
     2201       IF (type_trac == 'repr') THEN
     2202#ifdef REPROBUS
     2203           DO iq=1,nbnas
     2204             CALL histwrite_phy(o_nas(iq), nas(:,:,iq))
     2205           ENDDO
     2206#endif
     2207       ENDIF
     2208
    20902209       ENDIF   !(iflag_phytrac==1)
    20912210
Note: See TracChangeset for help on using the changeset viewer.