Changeset 308 for trunk/LMDZ.MARS/libf


Ignore:
Timestamp:
Sep 26, 2011, 7:53:28 PM (13 years ago)
Author:
aslmd
Message:

MESOSCALE: follow-up to previous commit.

Location:
trunk/LMDZ.MARS/libf/phymars
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • trunk/LMDZ.MARS/libf/phymars/meso_inc/meso_inc_var.F

    r285 r308  
    1818      INTEGER tracerset    !!! this corresponds to config%mars
    1919      CHARACTER (len=20) :: wtnom(nqmx) ! tracer name
     20
     21      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! JF
     22      REAL qsurfice_dust(ngridmx)     ! useful for dust diagnostics
     23      REAL TAU_lay(ngridmx)           ! true opacity (it's not e reference life tauref)
     24      REAL dsodust(ngridmx,nlayermx)
     25      REAL zdqsdif_diag(ngridmx)      ! useful for lifting diagnostics                               
     26      REAL zdqssed_diag(ngridmx)      ! useful for sedimentation diagnostics                         
     27      REAL pdq_diag(ngridmx)          ! useful for dust perturbation diagnostics   
     28      REAL dustot(ngridmx)            ! Total mass of dust integrated along vertical axe (kg/m2)
     29      REAL zdqnorm(ngridmx,nlayermx,2)
     30      !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! JF

  • trunk/LMDZ.MARS/libf/phymars/physiq.F

    r299 r308  
    14261426      wco2ice(1:ngrid) = co2ice(1:ngrid)  !! co2 ice
    14271427      mtot(1:ngrid) = mtot(1:ngrid) * 1.e6 / rho_ice
     1428      !! JF
     1429      TAU_lay(:)=tau(:,1)!!true opacity (not a reference like tauref)
     1430      IF (igcm_dust_mass .ne. 0) THEN
     1431        qsurfice_dust(1:ngrid) = qsurf(1:ngrid,igcm_dust_mass)
     1432      ENDIF
    14281433      IF (igcm_h2o_ice .ne. 0) THEN     
    14291434        qsurfice(1:ngrid) = qsurf(1:ngrid,igcm_h2o_ice)
    … …  
    14311436     .           * mugaz / mmol(igcm_h2o_ice)
    14321437      ENDIF
     1438      !! Dust quantity integration along the vertical axe
     1439      dustot(:)=0
     1440      do ig=1,ngrid
     1441       do l=1,nlayermx
     1442        dustot(ig) = dustot(ig) +
     1443     &               zq(ig,l,igcm_dust_mass)
     1444     &               *  (pplev(ig,l) - pplev(ig,l+1)) / g
     1445       enddo
     1446      enddo
    14331447c AUTOMATICALLY GENERATED FROM REGISTRY
    14341448#include "fill_save.inc"
Note: See TracChangeset for help on using the changeset viewer.