Ignore:
Timestamp:
Jun 6, 2016, 4:04:57 PM (8 years ago)
Author:
Laurent Fairhead
Message:

Merged trunk changes r2487:2541 into testing branch

Location:
LMDZ5/branches/testing
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/branches/testing

  • LMDZ5/branches/testing/libf/phylmd/phys_output_write_mod.F90

    r2488 r2542  
    1717  SUBROUTINE phys_output_write(itap, pdtphys, paprs, pphis, &
    1818       pplay, lmax_th, aerosol_couple,         &
    19        ok_ade, ok_aie, ivap, new_aod, ok_sync, &
     19       ok_ade, ok_aie, ivap, iliq, isol, new_aod, ok_sync, &
    2020       ptconv, read_climoz, clevSTD, ptconvth, &
    2121       d_t, qx, d_qx, zmasse, flag_aerosol, flag_aerosol_strat, ok_cdnc)
    … …  
    3636         o_wind10m, o_wind10max, o_gusts, o_sicf, &
    3737         o_psol, o_mass, o_qsurf, o_qsol, &
    38          o_precip, o_ndayrain, o_plul, o_pluc, &
     38         o_precip, o_ndayrain, o_plul, o_pluc, o_plun, &
    3939         o_snow, o_msnow, o_fsnow, o_evap, o_ep,o_epmax_diag, & ! epmax_cape
    4040         o_tops, o_tops0, o_topl, o_topl0, &
    … …  
    5050         o_bils_ec,o_bils_ech, o_bils_tke, o_bils_kinetic, &
    5151         o_bils_latent, o_bils_enthalp, o_sens, &
    52          o_fder, o_ffonte, o_fqcalving, o_fqfonte, &
     52         o_fder, o_ffonte, o_fqcalving, o_fqfonte, o_runofflic, &
    5353         o_taux, o_tauy, o_snowsrf, o_qsnow, &
    5454         o_snowhgt, o_toice, o_sissnow, o_runoff, &
    … …  
    6767         o_wbeff, o_cape_max, o_upwd, o_ep,o_epmax_diag, o_Ma, &
    6868         o_dnwd, o_dnwd0, o_ftime_con, o_mc, &
    69          o_prw, o_s_pblh, o_s_pblt, o_s_lcl, &
     69         o_prw, o_prlw, o_prsw, o_s_pblh, o_s_pblt, o_s_lcl, &
    7070         o_s_therm, o_uSTDlevs, o_vSTDlevs, &
    7171         o_wSTDlevs, o_zSTDlevs, o_qSTDlevs, &
    … …  
    112112         o_swtoacf_zero, o_swsrfcf_zero, &
    113113         o_topswai, o_solswai, o_scdnc, &
    114          o_cldncl, o_reffclws, o_reffclwc, &
     114         o_cldncl, o_reffclws, o_reffclwc, o_stratomask,&
    115115         o_cldnvi, o_lcc, o_lcc3d, o_lcc3dcon, &
    116116         o_lcc3dstra, o_reffclwtop, o_ec550aer, &
    … …  
    120120         o_zfull, o_zhalf, o_rneb, o_rnebjn, o_rnebcon, &
    121121         o_rnebls, o_rhum, o_ozone, o_ozone_light, &
    122          o_dtphy, o_dqphy, o_dqphy2d, o_albe_srf, o_z0m_srf, o_z0h_srf, &
     122         o_dtphy, o_dqphy, o_dqphy2d, o_dqlphy, o_dqlphy2d, &
     123         o_dqsphy, o_dqsphy2d, o_albe_srf, o_z0m_srf, o_z0h_srf, &
    123124         o_ages_srf, o_snow_srf, o_alb1, o_alb2, o_tke, &
    124125         o_tke_max, o_kz, o_kz_max, o_clwcon, &
    125          o_dtdyn, o_dqdyn, o_dqdyn2d, o_dudyn, o_dvdyn, &
     126         o_dtdyn, o_dqdyn, o_dqdyn2d, o_dqldyn, o_dqldyn2d, &
     127         o_dqsdyn, o_dqsdyn2d, o_dudyn, o_dvdyn, &
    126128         o_dtcon, o_tntc, o_ducon, o_dvcon, &
    127129         o_dqcon, o_dqcon2d, o_tnhusc, o_tnhusc, o_dtlsc, &
    … …  
    161163         o_dtr_sat, o_dtr_uscav, o_trac_cum, o_du_gwd_rando, o_dv_gwd_rando, &
    162164         o_ustr_gwd_hines,o_vstr_gwd_hines,o_ustr_gwd_rando,o_vstr_gwd_rando, &
    163          o_ustr_gwd_front,o_vstr_gwd_front
     165         o_ustr_gwd_front,o_vstr_gwd_front, &
     166         o_sens_prec_liq_oce, o_sens_prec_liq_sic, &
     167         o_sens_prec_sol_oce, o_sens_prec_sol_sic, &
     168         o_lat_prec_liq_oce, o_lat_prec_liq_sic, &
     169         o_lat_prec_sol_oce, o_lat_prec_sol_sic
    164170
    165171    USE phys_state_var_mod, only: pctsrf, paire_ter, rain_fall, snow_fall, &
    … …  
    193199         t2m_min_mon, t2m_max_mon, evap, &
    194200         zu10m, zv10m, zq2m, zustar, zxqsurf, &
    195          rain_lsc, snow_lsc, bils, sens, fder, &
    196          zxffonte, zxfqcalving, zxfqfonte, fluxu, &
     201         rain_lsc, rain_num, snow_lsc, bils, sens, fder, &
     202         zxffonte, zxfqcalving, zxfqfonte, zxrunofflic, fluxu, &
    197203         fluxv, zxsnow, qsnow, snowhgt, to_ice, &
    198204         sissnow, runoff, albsol3_lic, evap_pot, &
    … …  
    201207         cldh, cldt, JrNt, cldljn, cldmjn, cldhjn, &
    202208         cldtjn, cldq, flwp, fiwp, ue, ve, uq, vq, &
    203          plcl, plfc, wbeff, upwd, dnwd, dnwd0, prw, &
     209         plcl, plfc, wbeff, upwd, dnwd, dnwd0, prw, prlw, prsw, &
    204210         s_pblh, s_pblt, s_lcl, s_therm, uwriteSTD, &
    205211         vwriteSTD, wwriteSTD, phiwriteSTD, qwriteSTD, &
    … …  
    228234         toplwad_aero, toplwad0_aero, sollwad_aero, &
    229235         sollwad0_aero, toplwai_aero, sollwai_aero, &
    230          scdnc, cldncl, reffclws, reffclwc, cldnvi, &
     236         scdnc, cldncl, reffclws, reffclwc, cldnvi, stratomask,&
    231237         lcc, lcc3d, lcc3dcon, lcc3dstra, reffclwtop, &
    232238         ec550aer, flwc, fiwc, t_seri, theta, q_seri, &
    233 !jyg<
    234 !!         ql_seri, zphi, u_seri, v_seri, omega, cldfra, &
    235239         ql_seri, tr_seri, &
    236240         zphi, u_seri, v_seri, omega, cldfra, &
    237 !>jyg
    238          rneb, rnebjn, zx_rh, d_t_dyn, d_q_dyn, &
     241         rneb, rnebjn, zx_rh, d_t_dyn,  &
     242         d_q_dyn,  d_ql_dyn, d_qs_dyn, &
     243         d_q_dyn2d,  d_ql_dyn2d, d_qs_dyn2d, &
    239244         d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &
    240245         d_u_ajs, d_v_ajs, &
    … …  
    256261         itau_con, nfiles, clef_files, nid_files, &
    257262         zustr_gwd_hines, zvstr_gwd_hines,zustr_gwd_rando, zvstr_gwd_rando, &
    258          zustr_gwd_front, zvstr_gwd_front                                   
     263         zustr_gwd_front, zvstr_gwd_front,     &
     264         sens_prec_liq_o, sens_prec_sol_o, lat_prec_liq_o, lat_prec_sol_o
     265 
     266
    259267    USE ocean_slab_mod, only: tslab, slab_bils, slab_bilg, tice, seaice
    260268    USE pbl_surface_mod, only: snow
    … …  
    288296
    289297    ! Input
    290     INTEGER :: itap, ivap, read_climoz
     298    INTEGER :: itap, ivap, iliq, isol, read_climoz
    291299    INTEGER, DIMENSION(klon) :: lmax_th
    292300    LOGICAL :: aerosol_couple, ok_sync
    … …  
    301309    REAL, DIMENSION(klon,klev,nqtot) :: qx, d_qx
    302310    REAL, DIMENSION(klon, klev) :: zmasse
    303     LOGICAL :: flag_aerosol_strat
     311    INTEGER :: flag_aerosol_strat
    304312    INTEGER :: flag_aerosol
    305313    LOGICAL :: ok_cdnc
    … …  
    440448       ENDIF
    441449       CALL histwrite_phy(o_plul, zx_tmp_fi2d)
     450       CALL histwrite_phy(o_plun, rain_num)
    442451
    443452       IF (vars_defined) THEN
    … …  
    545554       CALL histwrite_phy(o_fqcalving, zxfqcalving)
    546555       CALL histwrite_phy(o_fqfonte, zxfqfonte)
     556       CALL histwrite_phy(o_runofflic, zxrunofflic)
    547557       IF (vars_defined) THEN
    548558          zx_tmp_fi2d=0.
    … …  
    618628
    619629       ENDDO
     630
     631       IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_liq_o(1 : klon, 1)
     632       CALL histwrite_phy(o_sens_prec_liq_oce, zx_tmp_fi2d)       
     633       IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_liq_o(1 : klon, 2)
     634       CALL histwrite_phy(o_sens_prec_liq_sic, zx_tmp_fi2d)       
     635       IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_sol_o(1 : klon, 1)
     636       CALL histwrite_phy(o_sens_prec_sol_oce, zx_tmp_fi2d)       
     637       IF (vars_defined) zx_tmp_fi2d(1 : klon) = sens_prec_sol_o(1 : klon, 2)
     638       CALL histwrite_phy(o_sens_prec_sol_sic, zx_tmp_fi2d)       
     639
     640       IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_liq_o(1 : klon, 1)
     641       CALL histwrite_phy(o_lat_prec_liq_oce, zx_tmp_fi2d)       
     642       IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_liq_o(1 : klon, 2)
     643       CALL histwrite_phy(o_lat_prec_liq_sic, zx_tmp_fi2d)       
     644       IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_sol_o(1 : klon, 1)
     645       CALL histwrite_phy(o_lat_prec_sol_oce, zx_tmp_fi2d)       
     646       IF (vars_defined) zx_tmp_fi2d(1 : klon) = lat_prec_sol_o(1 : klon, 2)
     647       CALL histwrite_phy(o_lat_prec_sol_sic, zx_tmp_fi2d)       
     648
    620649       DO nsrf=1,nbsrf+1
    621650          CALL histwrite_phy(o_wstar(nsrf), wstar(1 : klon, nsrf))
    … …  
    671700       ENDIF !iflag_con .GE. 3
    672701       CALL histwrite_phy(o_prw, prw)
     702       CALL histwrite_phy(o_prlw, prlw)
     703       CALL histwrite_phy(o_prsw, prsw)
    673704       CALL histwrite_phy(o_s_pblh, s_pblh)
    674705       CALL histwrite_phy(o_s_pblt, s_pblt)
    … …  
    938969             !--STRAT AER
    939970          ENDIF
    940           IF (flag_aerosol.GT.0.OR.flag_aerosol_strat) THEN
     971          IF (flag_aerosol.GT.0.OR.flag_aerosol_strat.GT.0) THEN
    941972!             DO naero = 1, naero_spc
    942973!--correction mini bug OB
    … …  
    946977             END DO
    947978          ENDIF
    948           IF (flag_aerosol_strat) THEN
     979          IF (flag_aerosol_strat.GT.0) THEN
    949980             CALL histwrite_phy(o_tausumaero_lw, &
    950981                  tausum_aero(:,6,id_STRAT_phy) )
    … …  
    10321063       CALL histwrite_phy(o_zfull,zphi/RG)
    10331064
     1065       IF (flag_aerosol_strat.EQ.2) THEN
     1066         CALL histwrite_phy(o_stratomask, stratomask)
     1067       ENDIF
     1068     
    10341069       IF (vars_defined)  THEN
    10351070        zx_tmp_fi3d(:,1)= pphis(:)/RG
    … …  
    10641099
    10651100       CALL histwrite_phy(o_dtphy, d_t)
     1101
    10661102       CALL histwrite_phy(o_dqphy,  d_qx(:,:,ivap))
    10671103       CALL water_int(klon,klev,d_qx(:,:,ivap),zmasse,zx_tmp_fi2d)
    10681104       CALL histwrite_phy(o_dqphy2d,  zx_tmp_fi2d)
     1105
     1106       CALL histwrite_phy(o_dqlphy,  d_qx(:,:,iliq))
     1107       CALL water_int(klon,klev,d_qx(:,:,iliq),zmasse,zx_tmp_fi2d)
     1108       CALL histwrite_phy(o_dqlphy2d,  zx_tmp_fi2d)
     1109
     1110       IF (nqo.EQ.3) THEN
     1111       CALL histwrite_phy(o_dqsphy,  d_qx(:,:,isol))
     1112       CALL water_int(klon,klev,d_qx(:,:,isol),zmasse,zx_tmp_fi2d)
     1113       CALL histwrite_phy(o_dqsphy2d,  zx_tmp_fi2d)
     1114       ELSE
     1115       zx_tmp_fi3d=0.0
     1116       CALL histwrite_phy(o_dqsphy,  zx_tmp_fi3d)
     1117       zx_tmp_fi2d=0.0
     1118       CALL histwrite_phy(o_dqsphy2d,  zx_tmp_fi2d)
     1119       ENDIF
     1120
    10691121       DO nsrf=1, nbsrf
    10701122          IF (vars_defined) zx_tmp_fi2d(1 : klon) = falb1( 1 : klon, nsrf)
    … …  
    11031155       CALL histwrite_phy(o_clwcon, clwcon0)
    11041156       CALL histwrite_phy(o_dtdyn, d_t_dyn)
     1157
    11051158       CALL histwrite_phy(o_dqdyn, d_q_dyn)
    1106        CALL water_int(klon,klev,d_q_dyn,zmasse,zx_tmp_fi2d)
    1107        CALL histwrite_phy(o_dqdyn2d,zx_tmp_fi2d)
     1159
     1160       CALL histwrite_phy(o_dqdyn2d,d_q_dyn2d)
     1161
     1162       CALL histwrite_phy(o_dqldyn, d_ql_dyn)
     1163
     1164       CALL histwrite_phy(o_dqldyn2d, d_ql_dyn2d)
     1165
     1166       CALL histwrite_phy(o_dqsdyn, d_qs_dyn)
     1167
     1168       CALL histwrite_phy(o_dqsdyn2d, d_qs_dyn2d)
     1169
    11081170       CALL histwrite_phy(o_dudyn, d_u_dyn)
    11091171       CALL histwrite_phy(o_dvdyn, d_v_dyn)
Note: See TracChangeset for help on using the changeset viewer.