Ignore:
Timestamp:
Dec 16, 2019, 8:55:01 AM (5 years ago)
Author:
flefevre
Message:

correction de bugs dans le supercycling de la chimie

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.VENUS/libf/phyvenus/phytrac_chimie.F

    r2197 r2200  
    270270
    271271      do iq = 1,nqmax - nmicro
    272          ztrac(:,:,iq) = ztrac(:,:,iq)*m_tr(iq)/mmean(:,:)
     272         ztrac(:,:,iq) = max(ztrac(:,:,iq)*m_tr(iq)/mmean(:,:),
     273     $                       1.e-30)
    273274         d_tr_chem(:,:,iq) = (ztrac(:,:,iq) - trac(:,:,iq))/pdtphys
    274275      end do
    … …  
    277278
    278279      if (ok_cloud .and. cl_scheme == 1) then
    279          ztrac(:,:,i_h2so4liq) = ztrac(:,:,i_h2so4liq)*m_tr(i_h2so4liq)
    280      $                           /mmean(:,:)
    281          ztrac(:,:,i_h2oliq) = ztrac(:,:,i_h2oliq)*m_tr(i_h2oliq)
    282      $                         /mmean(:,:)
     280         ztrac(:,:,i_h2so4liq) = max(ztrac(:,:,i_h2so4liq)
     281     $                               *m_tr(i_h2so4liq)/mmean(:,:),
     282     $                               1.e-30)
     283         ztrac(:,:,i_h2oliq)   = max(ztrac(:,:,i_h2oliq)
     284     $                               *m_tr(i_h2oliq)/mmean(:,:),
     285     $                               1.e-30)
    283286         d_tr_chem(:,:,i_h2so4liq) = (ztrac(:,:,i_h2so4liq)
    284287     $                              - trac(:,:,i_h2so4liq))/pdtphys
Note: See TracChangeset for help on using the changeset viewer.