Ignore:
Timestamp:
Jun 4, 2018, 11:59:42 AM (7 years ago)
Author:
jvatant
Message:

Add XIOS outputs for chemistry including "*_tot" fields (concat GCM+upper atm : 0->1300km )
Only in pseudo-pressure axis for now, pseudo-altitude TBD.
--JVO

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TabularUnified trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90

    r1927 r1943  
    1717      use radcommon_h, only: sigma, glat, grav, BWNV
    1818      use surfdat_h, only: phisfi, zmea, zstd, zsig, zgam, zthe
    19       use comchem_h, only: nkim, cnames
     19      use comchem_h, only: nkim, cnames, nlaykim_up, ykim_up, ykim_tot
    2020      use comdiurn_h, only: coslat, sinlat, coslon, sinlon
    2121      use comsaison_h, only: mu0, fract, dist_star, declin, right_ascen
     
    15931593      ! Chemical tracers
    15941594      if (callchim) then
     1595     
     1596        ! Advected fields
    15951597        do iq=1,nkim
    1596           CALL send_xios_field(cnames(iq),zq(:,:,iq+nmicro)/rat_mmol(iq+nmicro)) ! mol/mol
     1598          CALL send_xios_field(trim(cnames(iq)),zq(:,:,iq+nmicro)/rat_mmol(iq+nmicro)) ! kg/kg -> mol/mol
    15971599        enddo
    1598       endif
     1600       
     1601        ! Upper chemistry fields
     1602        do iq=1,nkim
     1603          CALL send_xios_field(trim(cnames(iq))//"_up",ykim_up(iq,:,:)) ! mol/mol
     1604        enddo
     1605       
     1606        ! Append fields in ykim_tot for output on the total vertical grid (0->1300km)
     1607        do iq=1,nkim
     1608         
     1609          ! GCM levels
     1610          do l=1,nlayer
     1611            ykim_tot(iq,:,l) = zq(:,l,iq+nmicro)/rat_mmol(iq+nmicro)
     1612          enddo
     1613          ! Upper levels
     1614          do l=1,nlaykim_up
     1615            ykim_tot(iq,:,nlayer+l) = ykim_up(iq,:,l)
     1616          enddo
     1617         
     1618          CALL send_xios_field(trim(cnames(iq))//"_tot",ykim_tot(iq,:,:)) ! mol/mol
     1619         
     1620        enddo
     1621       
     1622      endif ! of 'if callchim'
    15991623
    16001624      ! Microphysical tracers
Note: See TracChangeset for help on using the changeset viewer.