cvltr_spl Subroutine

subroutine cvltr_spl(pdtime, da, phi, phi2, d1a, dam, mpIN, epIN, sigd, sij, wght_cvfd, clw, elij, epmlmMm, eplaMm, pmflxrIN, pmflxsIN, ev, te, wdtrainA, wdtrainM, paprs, it, tr, upd, dnd, inb, icb, kk, henry, zrho, ccntrAA_spla, ccntrENV_spla, coefcoli_spla, id_prec, id_fine, id_coss, id_codu, id_scdu, dtrcv, trsptd, dtrSscav, dtrsat, dtrUscav, qDi, qPr, qPa, qMel, qTrdi, dtrcvMA, Mint, zmfd1a, zmfphi2, zmfdam)

impaction sur la surface de la colonne de la descente insaturee On prend la moyenne des precip entre le niveau i+1 et i I=3/4 (P(1+1)+P(i))/2 / (sigdr*rho_l) 1000kg/m3= densité de l'eau 0.75e-3 = 3/4 /1000 Par la suite, I est tout le temps multiplié par sig_d pour avoir l'impaction sur la surface de la maille on le néglige ici pour simplifier le code

rho(i,j)=pplay(i,j)/(rd*te(i,j))

!!! quantite de traceur dans la descente d'air insaturee : 4 juin 2012 !!!

    zmfa(i,j,it)=da(i,j)*(tr(i,1,it)-tr(i,j,it))                     ! da

jyg (20130119) correction pour le sommet du nuage if(j.ge.inb(i)) then !au-dessus du nuage, sommet inclu print test descente insaturee

     'zmfdam',zmfdam(i,j,it),'zmfpsi',zmfphi2(i,j,it),&


      'zmfa',zmfa(i,j,it),'zmfp',zmfp(i,j,it),&

      'Pa',Pa(i,j),'eplaMm',eplaMm(i,j)
     'zmfd1a=da1*qa',zmfd1a(i,j,it),'Pa*qPa',wdtrainA(i,j)*qPa(i,j,it),'da1',d1a(i,j)

JE adds

Je end

!!!!

Arguments