Changeset 1943 for trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90
- Timestamp:
- Jun 4, 2018, 11:59:42 AM (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TabularUnified trunk/LMDZ.TITAN/libf/phytitan/physiq_mod.F90 ¶
r1927 r1943 17 17 use radcommon_h, only: sigma, glat, grav, BWNV 18 18 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 20 20 use comdiurn_h, only: coslat, sinlat, coslon, sinlon 21 21 use comsaison_h, only: mu0, fract, dist_star, declin, right_ascen … … 1593 1593 ! Chemical tracers 1594 1594 if (callchim) then 1595 1596 ! Advected fields 1595 1597 do iq=1,nkim 1596 CALL send_xios_field( cnames(iq),zq(:,:,iq+nmicro)/rat_mmol(iq+nmicro)) !mol/mol1598 CALL send_xios_field(trim(cnames(iq)),zq(:,:,iq+nmicro)/rat_mmol(iq+nmicro)) ! kg/kg -> mol/mol 1597 1599 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' 1599 1623 1600 1624 ! Microphysical tracers
Note: See TracChangeset
for help on using the changeset viewer.