Ignore:
Timestamp:
Jun 21, 2022, 11:05:31 AM (3 years ago)
Author:
aslmd
Message:

fixed update of tendencies for generic tracer

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/libf/phystd/physiq_mod.F90

    r2690 r2692  
    15461546  !   VI.3. Aerosol particles
    15471547  ! -------------------------
    1548          ! Generic condensation and evaporation
    1549          ! i need the tracers and to check wether it ends with _vap
     1548         ! Generic tracers condensation and evaporation
    15501549         do iq=1,nq
    1551             ! print*,"starting the generic looping"
    1552             ! write(*,*)"for iq = ",iq, "is_generic = ",is_generic(iq)
    1553             ! write(*,*)"noms(iq) = ",noms(iq)
    1554             ! write(*,*),"noms(iq)(len(noms(iq)-2:)=",noms(iq)(len(noms(iq))-2:)
    1555             !loop on tracers and check if is_generic and _vap in noms
     1550            !loop on tracers and check if is_generic and  if "vap" in noms
    15561551            if((is_generic(iq)==1) .and. (index(noms(iq),"vap") .ne. 0)) then
    15571552               ! call to largescale_generic
     
    15601555                                       pt,pq,pdt,pdq,noms(iq),dtlscale_generic, &
    15611556                                       dqvaplscale_generic,dqcldlscale_generic)
    1562                !question sur la tendance de temperature qui sort. On l'ajoute dans une variable temp
    1563                ! et apres la fin de la boucle sur les tracers on l'ajouteà la tendance du pas de temps physique ?
    1564                !temp_dt = temp_dt+dtlscale_generic ?
     1557               pdt(1:ngrid,1:nlayer) = pdt(1:ngrid,1:nlayer)+dtlscale_generic(1:ngrid,1:nlayer)
    15651558               pdq(1:ngrid,1:nlayer,ig) = pdq(1:ngrid,1:nlayer,ig)+dtlscale_generic(1:ngrid,1:nlayer)
    15661559               pdq(1:ngrid,1:nlayer,ig+1) = pdq(1:ngrid,1:nlayer,ig+1)+dqcldlscale_generic(1:ngrid,1:nlayer)
     
    15681561            endif ! (is_generic(iq)==1 .and. (index(noms(iq),"vap")==1))
    15691562         enddo ! iq=1,nq
    1570          !should I add the pdt = pdt + temp_dt ?
    15711563
    15721564         ! Sedimentation.
Note: See TracChangeset for help on using the changeset viewer.