Changeset 3223 for LMDZ6/trunk/libf/phylmd/dyn1d/1D_nudge_sandu_astex.h

Timestamp:

Feb 17, 2018, 6:02:23 PM (6 years ago)

Author:

fhourdin

Message:

Corrections forcages cas 1D sandu

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/1D_nudge_sandu_astex.h (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/1D_nudge_sandu_astex.h

@@ -20 +20 @@
 !     print *,'l dt  relax dtadv',l,dt_phys(l),relax_thl(l),d_t_adv(l)
       enddo
-        u(1:mxcalc)=u(1:mxcalc) + timestep*(                                &          
+
+        u(1:mxcalc)=u(1:mxcalc) + timestep*( du_age(1:mxcalc)+             &          
      &              du_phys(1:mxcalc) - relax_u(1:mxcalc))
-        v(1:mxcalc)=v(1:mxcalc) + timestep*(                                &
+        v(1:mxcalc)=v(1:mxcalc) + timestep*( dv_age(1:mxcalc)+              &
      &               dv_phys(1:mxcalc) - relax_v(1:mxcalc))
-!       q(1:mxcalc,:)=q(1:mxcalc,:)+timestep*(
-!    .              dq(1:mxcalc,:) - relax_q(1:mxcalc,:)+
-!    .              d_q_adv(1:mxcalc,:))
         q(1:mxcalc,1)=q(1:mxcalc,1)+timestep*(                              &
      &      dq(1:mxcalc,1) - relax_q(1:mxcalc,1)+d_q_adv(1:mxcalc,1))
@@ -34 +32 @@
      &      dt_phys(1:mxcalc)-relax_thl(1:mxcalc)+d_t_adv(1:mxcalc))
 
+