Ignore:
Timestamp:
Apr 16, 2014, 7:16:58 AM (10 years ago)
Author:
fhourdin
Message:

Passage au format libre pour inclure les ficheirs du 1D dans lmdz1d.F90

File:
1 edited

Legend:

Unmodified
Added
Removed

  • LMDZ5/trunk/libf/phylmd/1D_nudge_sandu_astex.h

    r2017 r2019  
    2020!     print *,'l dt  relax dtadv',l,dt_phys(l),relax_thl(l),d_th_adv(l)
    2121      enddo
    22         u(1:mxcalc)=u(1:mxcalc) + timestep*(
    23      .              du_phys(1:mxcalc) - relax_u(1:mxcalc))
    24         v(1:mxcalc)=v(1:mxcalc) + timestep*(
    25      .              dv_phys(1:mxcalc) - relax_v(1:mxcalc))
    26 c       q(1:mxcalc,:)=q(1:mxcalc,:)+timestep*(
    27 c    .              dq(1:mxcalc,:) - relax_q(1:mxcalc,:)+
    28 c    .              d_q_adv(1:mxcalc,:))
    29         q(1:mxcalc,1)=q(1:mxcalc,1)+timestep*(
    30      .      dq(1:mxcalc,1) - relax_q(1:mxcalc,1)+d_q_adv(1:mxcalc,1))
    31         q(1:mxcalc,2)=q(1:mxcalc,2)+timestep*(
    32      .      dq(1:mxcalc,2) - relax_q(1:mxcalc,2)+d_q_adv(1:mxcalc,2))
    33         temp(1:mxcalc)=temp(1:mxcalc)+timestep*(
    34      .      dt_phys(1:mxcalc)-relax_thl(1:mxcalc)+d_th_adv(1:mxcalc))
     22        u(1:mxcalc)=u(1:mxcalc) + timestep*(                                &         
     23     &              du_phys(1:mxcalc) - relax_u(1:mxcalc))
     24        v(1:mxcalc)=v(1:mxcalc) + timestep*(                                &
     25     &               dv_phys(1:mxcalc) - relax_v(1:mxcalc))
     26!       q(1:mxcalc,:)=q(1:mxcalc,:)+timestep*(
     27!    .              dq(1:mxcalc,:) - relax_q(1:mxcalc,:)+
     28!    .              d_q_adv(1:mxcalc,:))
     29        q(1:mxcalc,1)=q(1:mxcalc,1)+timestep*(                              &
     30     &      dq(1:mxcalc,1) - relax_q(1:mxcalc,1)+d_q_adv(1:mxcalc,1))
     31        q(1:mxcalc,2)=q(1:mxcalc,2)+timestep*(                              &
     32     &      dq(1:mxcalc,2) - relax_q(1:mxcalc,2)+d_q_adv(1:mxcalc,2))
     33        temp(1:mxcalc)=temp(1:mxcalc)+timestep*(                            &
     34     &      dt_phys(1:mxcalc)-relax_thl(1:mxcalc)+d_th_adv(1:mxcalc))
Note: See TracChangeset for help on using the changeset viewer.