Ignore:
Timestamp:
Nov 19, 2013, 1:14:17 PM (11 years ago)
Author:
musat
Message:

Correction sorties dvlif, dtlif seulement si ok_orolf (comme pour dulif).
Ajout sorties tendances duthe, dvthe.
MPL/JY/FC
IM

Location:
LMDZ5/trunk/libf/phylmd
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

    r1872 r1903  
    4141      REAL, SAVE, ALLOCATABLE :: d_u_ajs(:,:), d_v_ajs(:,:)
    4242      !$OMP THREADPRIVATE(d_u_ajs, d_v_ajs)
     43      REAL, SAVE, ALLOCATABLE :: d_u_ajsb(:,:), d_v_ajsb(:,:)
     44      !$OMP THREADPRIVATE(d_u_ajsb, d_v_ajsb)
    4345      REAL, SAVE, ALLOCATABLE :: d_t_eva(:,:),d_q_eva(:,:)
    4446      !$OMP THREADPRIVATE(d_t_eva,d_q_eva)
     
    297299      allocate(d_t_ajs(klon,klev),d_q_ajs(klon,klev))
    298300      allocate(d_u_ajs(klon,klev),d_v_ajs(klon,klev))
     301      allocate(d_u_ajsb(klon,klev),d_v_ajsb(klon,klev))
    299302      allocate(d_t_eva(klon,klev),d_q_eva(klon,klev))
    300303      allocate(d_t_lscst(klon,klev),d_q_lscst(klon,klev))
     
    447450      deallocate(d_t_ajs,d_q_ajs)
    448451      deallocate(d_u_ajs,d_v_ajs)
     452      deallocate(d_u_ajsb,d_v_ajsb)
    449453      deallocate(d_t_eva,d_q_eva)
    450454      deallocate(d_t_lscst,d_q_lscst)
  • LMDZ5/trunk/libf/phylmd/phys_output_ctrlout_mod.F90

    r1895 r1903  
    10121012  TYPE(ctrl_out), SAVE :: o_dtthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
    10131013    'dtthe', 'Thermal dT', 'K/s', (/ ('', i=1, 9) /))
     1014  TYPE(ctrl_out), SAVE :: o_duthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
     1015    'duthe', 'Thermal du', 'm/s2', (/ ('', i=1, 9) /))
     1016  TYPE(ctrl_out), SAVE :: o_dvthe = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
     1017    'dvthe', 'Thermal dv', 'm/s2', (/ ('', i=1, 9) /))
    10141018  TYPE(ctrl_out), SAVE :: o_f_th = ctrl_out((/ 4, 10, 10, 10, 10, 10, 11, 11, 11 /), &
    10151019    'f_th', 'Thermal plume mass flux', 'kg/(m2*s)', (/ ('', i=1, 9) /))
  • LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

    r1898 r1903  
    119119         o_ptconvth, o_lmaxth, o_dtvdf, &
    120120         o_dtdis, o_dqvdf, o_dteva, o_dqeva, &
    121          o_ptconv, o_ratqs, o_dtthe, o_ftime_th, &
     121         o_ptconv, o_ratqs, o_dtthe, &
     122         o_duthe, o_dvthe, o_ftime_th, &
    122123         o_f_th, o_e_th, o_w_th, o_q_th, &
    123124         o_a_th, o_d_th, o_f0_th, o_zmax_th, &
     
    202203         rneb, zx_rh, frugs, agesno, d_t_dyn, d_q_dyn, &
    203204         d_u_dyn, d_v_dyn, d_t_con, d_t_ajsb, d_t_ajs, &
     205         d_u_ajsb, d_u_ajs, d_v_ajsb, d_v_ajs, &
    204206         d_u_con, d_v_con, d_q_con, d_q_ajs, d_t_lsc, &
    205207         d_t_eva, d_q_lsc, beta_prec, d_t_lscth, &
     
    10261028       ENDIF
    10271029       CALL histwrite_phy(o_dtthe, zx_tmp_fi3d)
     1030       IF (vars_defined) THEN
     1031          zx_tmp_fi3d(1:klon,1:klev)=d_u_ajs(1:klon,1:klev)/pdtphys - &
     1032               d_u_ajsb(1:klon,1:klev)/pdtphys
     1033       ENDIF
     1034       CALL histwrite_phy(o_duthe, zx_tmp_fi3d)
     1035       IF (vars_defined) THEN
     1036          zx_tmp_fi3d(1:klon,1:klev)=d_v_ajs(1:klon,1:klev)/pdtphys - &
     1037               d_v_ajsb(1:klon,1:klev)/pdtphys
     1038       ENDIF
     1039       CALL histwrite_phy(o_dvthe, zx_tmp_fi3d)
     1040
    10281041       IF (iflag_thermals>=1) THEN
    10291042          ! Pour l instant 0 a y reflichir pour les thermiques
     
    10731086          IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_u_lif(1:klon,1:klev)/pdtphys
    10741087          CALL histwrite_phy(o_dulif, zx_tmp_fi3d)
    1075        ENDIF
    1076        IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
    1077        CALL histwrite_phy(o_dvlif, zx_tmp_fi3d)
    1078 
    1079        IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys
    1080        CALL histwrite_phy(o_dtlif, zx_tmp_fi3d)
     1088
     1089          IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_v_lif(1:klon,1:klev)/pdtphys
     1090          CALL histwrite_phy(o_dvlif, zx_tmp_fi3d)
     1091
     1092          IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lif(1:klon,1:klev)/pdtphys
     1093          CALL histwrite_phy(o_dtlif, zx_tmp_fi3d)
     1094       ENDIF
    10811095
    10821096       IF (ok_hines) THEN
Note: See TracChangeset for help on using the changeset viewer.