Changeset 2194

Timestamp:

Feb 5, 2015, 12:29:37 PM (9 years ago)

Author:

fhourdin

Message:

Use of add_phys_tend to add the tendencies from radiation.
In order to control temperature range after radiation.

Utilisation de add_phys_tend pour ajouter les tendances
du rayonnement, de façon à contrôler les températures après
cet ajout.

Location:

LMDZ5/trunk/libf/phylmd

Files:

• phys_local_var_mod.F90 (modified) (3 diffs)
• phys_output_write_mod.F90 (modified) (4 diffs)
• physiq.F90 (modified) (1 diff)

LMDZ5/trunk/libf/phylmd/phys_local_var_mod.F90

@@ -35 +35 @@
       REAL, SAVE, ALLOCATABLE :: d_t_lsc(:,:),d_q_lsc(:,:),d_ql_lsc(:,:),d_qi_lsc(:,:)
       !$OMP THREADPRIVATE(d_t_lsc,d_q_lsc,d_ql_lsc,d_qi_lsc)
+      REAL, SAVE, ALLOCATABLE :: d_t_lwr(:,:),d_t_lw0(:,:),d_t_swr(:,:),d_t_sw0(:,:)
+      !$OMP THREADPRIVATE(d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0)
       REAL, SAVE, ALLOCATABLE :: d_t_ajsb(:,:), d_q_ajsb(:,:)
       !$OMP THREADPRIVATE(d_t_ajsb, d_q_ajsb)
@@ -393 +395 @@
       allocate(d_t_wake(klon,klev),d_q_wake(klon,klev))
       allocate(d_t_lsc(klon,klev),d_q_lsc(klon,klev))
+      allocate(d_t_lwr(klon,klev),d_t_lw0(klon,klev))
+      allocate(d_t_swr(klon,klev),d_t_sw0(klon,klev))
       allocate(d_ql_lsc(klon,klev),d_qi_lsc(klon,klev))
       allocate(d_t_ajsb(klon,klev),d_q_ajsb(klon,klev))
@@ -599 +603 @@
       deallocate(d_t_wake,d_q_wake)
       deallocate(d_t_lsc,d_q_lsc)
+      deallocate(d_t_lwr,d_t_lw0)
+      deallocate(d_t_swr,d_t_sw0)
       deallocate(d_ql_lsc,d_qi_lsc)
       deallocate(d_t_ajsb,d_q_ajsb)

LMDZ5/trunk/libf/phylmd/phys_output_write_mod.F90

@@ -166 +166 @@
          wake_deltaq, ftd, fqd, ale_bl_trig, albsol1, &
          rnebcon, wo, falb1, albsol2, coefh, clwcon0, &
-         ratqs, entr_therm, zqasc, detr_therm, f0, heat, &
-         heat0, cool, cool0, lwup, lwdn, lwup0, coefm, &
+         ratqs, entr_therm, zqasc, detr_therm, f0, &
+         lwup, lwdn, lwup0, coefm, &
          swupp, lwupp, swup0p, lwup0p, swdnp, lwdnp, &
          swdn0p, lwdn0p, tnondef, O3sumSTD, uvsumSTD, &
@@ -215 +215 @@
          d_u_ajs, d_v_ajs, &
          d_u_con, d_v_con, d_q_con, d_q_ajs, d_t_lsc, &
+         d_t_lwr,d_t_lw0,d_t_swr,d_t_sw0, &
          d_t_eva, d_q_lsc, beta_prec, d_t_lscth, &
          d_t_lscst, d_q_lscth, d_q_lscst, plul_th, &
@@ -1127 +1128 @@
        IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_q_ajsb(1:klon,1:klev)/pdtphys
        CALL histwrite_phy(o_dqajs, zx_tmp_fi3d)
-       IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY
+       IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys
        CALL histwrite_phy(o_dtswr, zx_tmp_fi3d)
-       IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=heat0(1:klon,1:klev)/RDAY
+       IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_sw0(1:klon,1:klev)/pdtphys
        CALL histwrite_phy(o_dtsw0, zx_tmp_fi3d)
-       IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=-1.*cool(1:klon,1:klev)/RDAY
+       IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lwr(1:klon,1:klev)/pdtphys
        CALL histwrite_phy(o_dtlwr, zx_tmp_fi3d)
-       IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=-1.*cool0(1:klon,1:klev)/RDAY
+       IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_lw0(1:klon,1:klev)/pdtphys
        CALL histwrite_phy(o_dtlw0, zx_tmp_fi3d)
        IF(vars_defined) zx_tmp_fi3d(1:klon,1:klev)=d_t_ec(1:klon,1:klev)/pdtphys
@@ -1193 +1194 @@
        CALL histwrite_phy(o_tnt, zx_tmp_fi3d)
        IF(vars_defined) THEN
-          zx_tmp_fi3d(1:klon,1:klev)=heat(1:klon,1:klev)/RDAY - &
-               cool(1:klon,1:klev)/RDAY
+          zx_tmp_fi3d(1:klon,1:klev)=d_t_swr(1:klon,1:klev)/pdtphys + &
+               d_t_lwr(1:klon,1:klev)/pdtphys
        ENDIF
        CALL histwrite_phy(o_tntr, zx_tmp_fi3d)

LMDZ5/trunk/libf/phylmd/physiq.F90

@@ -3455 +3455 @@
   ! Ajouter la tendance des rayonnements (tous les pas)
   !
-  DO k = 1, klev
-     DO i = 1, klon
-        t_seri(i,k) = t_seri(i,k) &
-             + (heat(i,k)-cool(i,k)) * dtime/RDAY
-     ENDDO
-  ENDDO
+  d_t_swr(:,:)=heat(:,:)*dtime/RDAY
+  d_t_lwr(:,:)=-cool(:,:)*dtime/RDAY
+  d_t_sw0(:,:)=heat0(:,:)*dtime/RDAY
+  d_t_lw0(:,:)=-cool0(:,:)*dtime/RDAY
+  CALL add_phys_tend(du0,dv0,d_t_swr,dq0,dql0,dqi0,paprs,'SW')
+  CALL add_phys_tend(du0,dv0,d_t_lwr,dq0,dql0,dqi0,paprs,'LW')
+
   !
   if (mydebug) then