Changeset 2021 for trunk


Ignore:
Timestamp:
Oct 19, 2018, 6:16:41 PM (6 years ago)
Author:
mlefevre
Message:

Photochemistry architecture for Venus mesoscale and LES model

Location:
trunk/MESOSCALE/LMD_MM_MARS/SRC
Files:
10 edited

Legend:

Unmodified
Added
Removed

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM

    r2018 r2021  
    147147####MARS MARS MARS MARS MARS
    148148####MARS MARS MARS MARS MARS
    149 
     149###VENUS TRACER
     150state  real  qco2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO2"       "CO2 number mixing ratio"   "kg kg-1"
     151state  real  qco       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO"        "CO number mixing ratio"   "kg kg-1"
     152state  real  qh2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2"        "H2 number mixing ratio"   "kg kg-1"
     153state  real  qo1d      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO1D"        "O1d number mixing ratio"   "kg kg-1"
     154state  real  qo        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO"        "O number mixing ratio"   "kg kg-1"
     155state  real  qo2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2"        "O2 number mixing ratio"   "kg kg-1"
     156state  real  qo2dg     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2dg"        "O2dg number mixing ratio"   "kg kg-1"
     157state  real  qo3       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO3"        "O3 number mixing ratio"   "kg kg-1"
     158state  real  qh        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH"        "H number mixing ratio"   "kg kg-1"
     159state  real  qoh       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOH"        "OH number mixing ratio"   "kg kg-1"
     160state  real  qho2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHO2"        "HO2 number mixing ratio"   "kg kg-1"
     161state  real  qh2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2O2"        "H2O2 number mixing ratio"   "kg kg-1"
     162state  real  qcl       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL"        "Cl number mixing ratio"   "kg kg-1"
     163state  real  qclo      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLO"        "ClO number mixing ratio"   "kg kg-1"
     164state  real  qcl2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL2"        "Cl2 number mixing ratio"   "kg kg-1"
     165state  real  qhcl      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHCL"        "HCl number mixing ratio"   "kg kg-1"
     166state  real  qhocl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHOCL"        "HOCl number mixing ratio"   "kg kg-1"
     167state  real  qclco     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO"        "ClCO number mixing ratio"   "kg kg-1"
     168state  real  qclco3    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO3"        "ClCO3 number mixing ratio"   "kg kg-1"
     169state  real  qcocl2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCOCL2"        "COCl2 number mixing ratio"   "kg kg-1"
     170state  real  qss       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS"        "S number mixing ratio"   "kg kg-1"
     171state  real  qso       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO"        "SO number mixing ratio"   "kg kg-1"
     172state  real  qso2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO2"        "SO2 number mixing ratio"   "kg kg-1"
     173state  real  qso3      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO3"        "SO3 number mixing ratio"   "kg kg-1"
     174state  real  qs2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2O2"        "S2O2 number mixing ratio"   "kg kg-1"
     175state  real  qocs      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOCS"        "OCS number mixing ratio"   "kg kg-1"
     176state  real  qhso3     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHSO3"        "HSO3 number mixing ratio"   "kg kg-1"
     177state  real  qh2so4    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4"        "H2SO4 number mixing ratio"   "kg kg-1"
     178state  real  qs2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2"        "S2 number mixing ratio"   "kg kg-1"
     179state  real  qclso2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLSO2"        "ClSO2 number mixing ratio"   "kg kg-1"
     180state  real  qoscl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOSCL"        "OSCl number mixing ratio"   "kg kg-1"
     181state  real  qh2oliq   ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2OLIQ"        "H2O liq number mixing ratio"   "kg kg-1"
     182state  real  qh2so4liq ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4LIQ"        "H2SO4 liq number mixing ratio"   "kg kg-1"
    150183
    151184
    … …  
    15171550package   radioac      mars==20                     -              scalar:qtrac1
    15181551package   radioac2     mars==21                     -              scalar:upward,downward
     1552package   VenusPC      mars==34                     -              scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq
    15191553##### MARS OPTIONS
    15201554##### MARS OPTIONS

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/Registry.EM.newphys

    r2018 r2021  
    144144####MARS MARS MARS MARS MARS
    145145####MARS MARS MARS MARS MARS
     146###VENUS TRACER
     147state  real  qco2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO2"       "CO2 number mixing ratio"   "kg kg-1"
     148state  real  qco       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO"        "CO number mixing ratio"   "kg kg-1"
     149state  real  qh2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2"        "H2 number mixing ratio"   "kg kg-1"
     150state  real  qo1d      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO1D"        "O1d number mixing ratio"   "kg kg-1"
     151state  real  qo        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO"        "O number mixing ratio"   "kg kg-1"
     152state  real  qo2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2"        "O2 number mixing ratio"   "kg kg-1"
     153state  real  qo2dg     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2dg"        "O2dg number mixing ratio"   "kg kg-1"
     154state  real  qo3       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO3"        "O3 number mixing ratio"   "kg kg-1"
     155state  real  qh        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH"        "H number mixing ratio"   "kg kg-1"
     156state  real  qoh       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOH"        "OH number mixing ratio"   "kg kg-1"
     157state  real  qho2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHO2"        "HO2 number mixing ratio"   "kg kg-1"
     158state  real  qh2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2O2"        "H2O2 number mixing ratio"   "kg kg-1"
     159state  real  qcl       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL"        "Cl number mixing ratio"   "kg kg-1"
     160state  real  qclo      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLO"        "ClO number mixing ratio"   "kg kg-1"
     161state  real  qcl2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL2"        "Cl2 number mixing ratio"   "kg kg-1"
     162state  real  qhcl      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHCL"        "HCl number mixing ratio"   "kg kg-1"
     163state  real  qhocl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHOCL"        "HOCl number mixing ratio"   "kg kg-1"
     164state  real  qclco     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO"        "ClCO number mixing ratio"   "kg kg-1"
     165state  real  qclco3    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO3"        "ClCO3 number mixing ratio"   "kg kg-1"
     166state  real  qcocl2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCOCL2"        "COCl2 number mixing ratio"   "kg kg-1"
     167state  real  qss       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS"        "S number mixing ratio"   "kg kg-1"
     168state  real  qso       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO"        "SO number mixing ratio"   "kg kg-1"
     169state  real  qso2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO2"        "SO2 number mixing ratio"   "kg kg-1"
     170state  real  qso3      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO3"        "SO3 number mixing ratio"   "kg kg-1"
     171state  real  qs2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2O2"        "S2O2 number mixing ratio"   "kg kg-1"
     172state  real  qocs      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOCS"        "OCS number mixing ratio"   "kg kg-1"
     173state  real  qhso3     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHSO3"        "HSO3 number mixing ratio"   "kg kg-1"
     174state  real  qh2so4    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4"        "H2SO4 number mixing ratio"   "kg kg-1"
     175state  real  qs2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2"        "S2 number mixing ratio"   "kg kg-1"
     176state  real  qclso2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLSO2"        "ClSO2 number mixing ratio"   "kg kg-1"
     177state  real  qoscl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOSCL"        "OSCl number mixing ratio"   "kg kg-1"
     178state  real  qh2oliq   ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2OLIQ"        "H2O liq number mixing ratio"   "kg kg-1"
     179state  real  qh2so4liq ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4LIQ"        "H2SO4 liq number mixing ratio"   "kg kg-1"
     180
    146181
    147182
    … …  
    15161551package   newwatersca  mars==12                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn
    15171552package   radioac      mars==20                     -              scalar:qtrac1
     1553package   VenusPC      mars==34                     -              scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq
    15181554##### MARS OPTIONS
    15191555##### MARS OPTIONS

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/LES/modif_mars/module_initialize_les.F

    r2020 r2021  
    8787   INTEGER, PARAMETER :: nl_max = 1000
    8888   REAL, DIMENSION(nl_max) :: zk, p_in, theta, rho, u, v, qv, pd_in
     89   REAL*8, DIMENSION(nl_max) :: pd_in8
    8990   INTEGER :: nl_in
    90 
    9191
    9292   INTEGER :: icm,jcm, ii, im1, jj, jm1, loop, error, fid, nxc, nyc
    9393   REAL    :: u_mean,v_mean, f0, p_surf, p_level, qvf, z_at_v, z_at_u
     94   REAL*8 :: p_level8
    9495   REAL    :: xrad, yrad, zrad, rad, delt, cof1, cof2
    9596!   REAL, EXTERNAL :: interp_0
    … …  
    121122 REAL, DIMENSION(nl_max) :: venus_hrdyn
    122123 REAL, DIMENSION(nl_max) :: altitude
     124 REAL*8, DIMENSION(nl_max) :: trac
    123125!!MARS
    124126
    … …  
    969971    ENDIF
    970972
     973    IF ((planet.eq."venus") .AND. ( config_flags%mars == 34 )) Then
     974       pd_in8(:)=pd_in(:)
     975       do i = 1,34
     976         call read_tracer(trac,num_scalar,i,nl_in)
     977         DO k=1,kte-1
     978           p_level = grid%znu(k)*(pd_surf - grid%p_top) + grid%p_top
     979           p_level8=p_level
     980           scalar(its:ite,k,jts:jte,i+1) =  interp_0_log2( trac, pd_in8, p_level8, nl_in )
     981         ENDDO
     982       ENDDO
     983       !close(unit=22)
     984    ENDIF
    971985
    972986    open(unit=21,file="altitude.txt",action="write")
    … …  
    13581372      end subroutine read_lsf
    13591373
     1374      subroutine read_tracer(trace,nq,qn,n)
     1375      implicit none
     1376      integer n,qn,nq ! qn : number of the tracer
     1377      real*8 tra(nq-1,n+1)
     1378      real*8 trace(n+1) !output
     1379      logical end_of_file
     1380
     1381      integer k,j
     1382
     1383! first element is the surface
     1384      open(unit=14,file='input_tracer',form='formatted',status='old')
     1385      rewind(14)
     1386      end_of_file = .false.
     1387        DO k=1,n
     1388          read(14,*) tra(:,k)
     1389          write(*,*) k,tra(qn,k)
     1390      ENDDO
     1391
     1392      close(14)
     1393      trace(:)=tra(qn,:)
     1394      end subroutine read_tracer
     1395
     1396
    13601397END MODULE module_initialize_ideal

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/PREP_MARS/Venus/readmeteo.F90

    r1756 r2021  
    252252allocate(aps(altlen),bps(altlen),levels(altlen))
    253253#ifdef PHOTOCHEM
    254 nchemtrac = 14
     254nchemtrac = 34
    255255allocate(wtnom(nchemtrac))
    256 wtnom(1)  = "c_co2"
    257 wtnom(2)  = "c_co"
    258 wtnom(3)  = "c_o"
    259 wtnom(4)  = "c_o1d"
    260 wtnom(5)  = "c_o2"
    261 wtnom(6)  = "c_o3"
    262 wtnom(7)  = "c_h"
    263 wtnom(8)  = "c_h2"
    264 wtnom(9)  = "c_oh"
    265 wtnom(10) = "c_ho2"
    266 wtnom(11) = "c_h2o2"
    267 wtnom(12) = "c_ch4"
    268 wtnom(13) = "c_n2"
    269 wtnom(14) = "c_ar"
     256print*,'PHOTOCHEM2.1'
     257wtnom(1)  = "co2"
     258wtnom(2)  = "co"
     259wtnom(3)  = "h2"
     260wtnom(4) = "h2o"
     261wtnom(5)  = "o1d"
     262wtnom(6)  = "o"
     263wtnom(7)  = "o2"
     264wtnom(8)  = "o2dg"
     265wtnom(9)  = "o3"
     266wtnom(10)  = "h"
     267wtnom(11)  = "oh"
     268wtnom(12)  = "ho2"
     269wtnom(13) = "h2o2"
     270wtnom(14)  = "cl"
     271wtnom(15)  = "clo"
     272wtnom(16)  = "cl2"
     273wtnom(17)  = "hcl"
     274wtnom(18)  = "hocl"
     275wtnom(19)  = "clco"
     276wtnom(20)  = "clco3"
     277wtnom(21)  = "cocl2"
     278wtnom(22)  = "s"
     279wtnom(23)  = "so"
     280wtnom(24)  = "so2"
     281wtnom(25)  = "so3"
     282wtnom(26)  = "s2o2"
     283wtnom(27)  = "ocs"
     284wtnom(28)  = "hso3"
     285wtnom(29)  = "h2so4"
     286wtnom(30)  = "s2"
     287wtnom(31)  = "clso2"
     288wtnom(32)  = "oscl"
     289wtnom(33)  = "h2oliq"
     290wtnom(34)  = "h2so4liq"
    270291allocate(chemtrac(lonlen,latlen,altlen,timelen,nchemtrac))
    271292chemtrac(:,:,:,:,:)=0
    … …  
    772793  vfile(j,i,:,:)=vfile(lonlen+1-j,latlen+1-i,:,:)
    773794  tsoilfile(j,i,:,:)=tsoilfile(lonlen+1-j,latlen+1-i,:,:)
     795#ifdef PHOTOCHEM
     796 chemtrac(j,i,:,:,:)=chemtrac(lonlen+1-j,latlen+1-i,:,:,:)
     797#endif
    774798ENDDO
    775799ENDDO 
    … …  
    13431367!------------------------!
    13441368#ifdef PHOTOCHEM
    1345     DO i=1,nchemtrac
    1346        FIELD=wtnom(i)
    1347        UNITS='units'
    1348        DESC='desc'
    1349        XLVL=200100.
    1350        SLAB=chemtrac(:,:,1,time_out(l),i)
    1351        ! And now put everything in the destination file
    1352        ! ... Header
    1353        write(1) IFV
    1354        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
    1355        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
    1356        ! ... Data
    1357        write(1) SLAB
    1358     ENDDO
     1369        FIELD='CO2'
     1370        UNITS='kg/kg'
     1371        DESC='CO2 mixing ratio'
     1372        XLVL=levels(k)
     1373        SLAB=chemtrac(:,:,1,time_out(l),1)
     1374                ! And now put everything in the destination file
     1375                ! ... Header
     1376        write(1) IFV
     1377        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1378        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1379                ! ... Data
     1380        write(1) SLAB
     1381
     1382
     1383        FIELD='CO'
     1384        UNITS='kg/kg'
     1385        DESC='CO mixing ratio'
     1386        XLVL=levels(k)
     1387        SLAB=chemtrac(:,:,1,time_out(l),2)
     1388                ! And now put everything in the destination file
     1389                ! ... Header
     1390        write(1) IFV
     1391        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1392        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1393                ! ... Data
     1394        write(1) SLAB
     1395
     1396
     1397        FIELD='H2'
     1398        UNITS='kg/kg'
     1399        DESC='H2 mixing ratio'
     1400        XLVL=levels(k)
     1401        SLAB=chemtrac(:,:,1,time_out(l),3)
     1402                ! And now put everything in the destination file
     1403                ! ... Header
     1404        write(1) IFV
     1405        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1406        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1407                ! ... Data
     1408        write(1) SLAB
     1409
     1410
     1411        FIELD='H2O'
     1412        UNITS='kg/kg'
     1413        DESC='H2O mixing ratio'
     1414        XLVL=levels(k)
     1415        SLAB=chemtrac(:,:,1,time_out(l),4)
     1416                ! And now put everything in the destination file
     1417                ! ... Header
     1418        write(1) IFV
     1419        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1420        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1421                ! ... Data
     1422        write(1) SLAB
     1423
     1424
     1425        FIELD='O1D'
     1426        UNITS='kg/kg'
     1427        DESC='O1D mixing ratio'
     1428        XLVL=levels(k)
     1429        SLAB=chemtrac(:,:,1,time_out(l),5)
     1430                ! And now put everything in the destination file
     1431                ! ... Header
     1432        write(1) IFV
     1433        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1434        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1435                ! ... Data
     1436        write(1) SLAB
     1437
     1438
     1439        FIELD='O'
     1440        UNITS='kg/kg'
     1441        DESC='O mixing ratio'
     1442        XLVL=levels(k)
     1443        SLAB=chemtrac(:,:,1,time_out(l),6)
     1444                ! And now put everything in the destination file
     1445                ! ... Header
     1446        write(1) IFV
     1447        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1448        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1449                ! ... Data
     1450        write(1) SLAB
     1451
     1452
     1453        FIELD='O2'
     1454        UNITS='kg/kg'
     1455        DESC='O2 mixing ratio'
     1456        XLVL=levels(k)
     1457        SLAB=chemtrac(:,:,1,time_out(l),7)
     1458                ! And now put everything in the destination file
     1459                ! ... Header
     1460        write(1) IFV
     1461        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1462        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1463                ! ... Data
     1464        write(1) SLAB
     1465
     1466
     1467        FIELD='O2DG'
     1468        UNITS='kg/kg'
     1469        DESC='O2DG mixing ratio'
     1470        XLVL=levels(k)
     1471        SLAB=chemtrac(:,:,1,time_out(l),8)
     1472                ! And now put everything in the destination file
     1473                ! ... Header
     1474        write(1) IFV
     1475        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1476        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1477                ! ... Data
     1478        write(1) SLAB
     1479
     1480
     1481        FIELD='O3'
     1482        UNITS='kg/kg'
     1483        DESC='O3 mixing ratio'
     1484        XLVL=levels(k)
     1485        SLAB=chemtrac(:,:,1,time_out(l),9)
     1486                ! And now put everything in the destination file
     1487                ! ... Header
     1488        write(1) IFV
     1489        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1490        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1491                ! ... Data
     1492        write(1) SLAB
     1493
     1494
     1495        FIELD='H'
     1496        UNITS='kg/kg'
     1497        DESC='H mixing ratio'
     1498        XLVL=levels(k)
     1499        SLAB=chemtrac(:,:,1,time_out(l),10)
     1500                ! And now put everything in the destination file
     1501                ! ... Header
     1502        write(1) IFV
     1503        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1504        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1505                ! ... Data
     1506        write(1) SLAB
     1507
     1508
     1509        FIELD='OH'
     1510        UNITS='kg/kg'
     1511        DESC='OH mixing ratio'
     1512        XLVL=levels(k)
     1513        SLAB=chemtrac(:,:,1,time_out(l),11)
     1514                ! And now put everything in the destination file
     1515                ! ... Header
     1516        write(1) IFV
     1517        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1518        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1519                ! ... Data
     1520        write(1) SLAB
     1521
     1522
     1523        FIELD='HO2'
     1524        UNITS='kg/kg'
     1525        DESC='hO2 mixing ratio'
     1526        XLVL=levels(k)
     1527        SLAB=chemtrac(:,:,1,time_out(l),12)
     1528                ! And now put everything in the destination file
     1529                ! ... Header
     1530        write(1) IFV
     1531        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1532        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1533                ! ... Data
     1534        write(1) SLAB
     1535
     1536
     1537        FIELD='H2O2'
     1538        UNITS='kg/kg'
     1539        DESC='H2O2 mixing ratio'
     1540        XLVL=levels(k)
     1541        SLAB=chemtrac(:,:,1,time_out(l),13)
     1542                ! And now put everything in the destination file
     1543                ! ... Header
     1544        write(1) IFV
     1545        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1546        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1547                ! ... Data
     1548        write(1) SLAB
     1549
     1550
     1551        FIELD='Cl'
     1552        UNITS='kg/kg'
     1553        DESC='Cl mixing ratio'
     1554        XLVL=levels(k)
     1555        SLAB=chemtrac(:,:,1,time_out(l),14)
     1556                ! And now put everything in the destination file
     1557                ! ... Header
     1558        write(1) IFV
     1559        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1560        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1561                ! ... Data
     1562        write(1) SLAB
     1563
     1564
     1565        FIELD='ClO'
     1566        UNITS='kg/kg'
     1567        DESC='ClO mixing ratio'
     1568        XLVL=levels(k)
     1569        SLAB=chemtrac(:,:,1,time_out(l),15)
     1570                ! And now put everything in the destination file
     1571                ! ... Header
     1572        write(1) IFV
     1573        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1574        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1575                ! ... Data
     1576        write(1) SLAB
     1577
     1578
     1579        FIELD='Cl2'
     1580        UNITS='kg/kg'
     1581        DESC='Cl2 mixing ratio'
     1582        XLVL=levels(k)
     1583        SLAB=chemtrac(:,:,1,time_out(l),16)
     1584                ! And now put everything in the destination file
     1585                ! ... Header
     1586        write(1) IFV
     1587        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1588        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1589                ! ... Data
     1590        write(1) SLAB
     1591
     1592
     1593        FIELD='HCl'
     1594        UNITS='kg/kg'
     1595        DESC='HCl mixing ratio'
     1596        XLVL=levels(k)
     1597        SLAB=chemtrac(:,:,1,time_out(l),17)
     1598                ! And now put everything in the destination file
     1599                ! ... Header
     1600        write(1) IFV
     1601        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1602        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1603                ! ... Data
     1604        write(1) SLAB
     1605
     1606
     1607        FIELD='HOCl'
     1608        UNITS='kg/kg'
     1609        DESC='HOCl mixing ratio'
     1610        XLVL=levels(k)
     1611        SLAB=chemtrac(:,:,1,time_out(l),18)
     1612                ! And now put everything in the destination file
     1613                ! ... Header
     1614        write(1) IFV
     1615        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1616        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1617                ! ... Data
     1618        write(1) SLAB
     1619
     1620
     1621        FIELD='ClCO'
     1622        UNITS='kg/kg'
     1623        DESC='ClCO mixing ratio'
     1624        XLVL=levels(k)
     1625        SLAB=chemtrac(:,:,1,time_out(l),19)
     1626                ! And now put everything in the destination file
     1627                ! ... Header
     1628        write(1) IFV
     1629        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1630        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1631                ! ... Data
     1632        write(1) SLAB
     1633
     1634
     1635        FIELD='ClCO3'
     1636        UNITS='kg/kg'
     1637        DESC='ClCO3 mixing ratio'
     1638        XLVL=levels(k)
     1639        SLAB=chemtrac(:,:,1,time_out(l),20)
     1640                ! And now put everything in the destination file
     1641                ! ... Header
     1642        write(1) IFV
     1643        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1644        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1645                ! ... Data
     1646        write(1) SLAB
     1647
     1648
     1649        FIELD='COCl2'
     1650        UNITS='kg/kg'
     1651        DESC='COClC2 mixing ratio'
     1652        XLVL=levels(k)
     1653        SLAB=chemtrac(:,:,1,time_out(l),21)
     1654                ! And now put everything in the destination file
     1655                ! ... Header
     1656        write(1) IFV
     1657        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1658        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1659                ! ... Data
     1660        write(1) SLAB
     1661
     1662
     1663        FIELD='S'
     1664        UNITS='kg/kg'
     1665        DESC='S mixing ratio'
     1666        XLVL=levels(k)
     1667        SLAB=chemtrac(:,:,1,time_out(l),22)
     1668                ! And now put everything in the destination file
     1669                ! ... Header
     1670        write(1) IFV
     1671        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1672        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1673                ! ... Data
     1674        write(1) SLAB
     1675
     1676
     1677        FIELD='SO'
     1678        UNITS='kg/kg'
     1679        DESC='SO mixing ratio'
     1680        XLVL=levels(k)
     1681        SLAB=chemtrac(:,:,1,time_out(l),23)
     1682                ! And now put everything in the destination file
     1683                ! ... Header
     1684        write(1) IFV
     1685        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1686        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1687                ! ... Data
     1688        write(1) SLAB
     1689
     1690
     1691        FIELD='SO2'
     1692        UNITS='kg/kg'
     1693        DESC='SO2 mixing ratio'
     1694        XLVL=levels(k)
     1695        SLAB=chemtrac(:,:,1,time_out(l),24)
     1696                ! And now put everything in the destination file
     1697                ! ... Header
     1698        write(1) IFV
     1699        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1700        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1701                ! ... Data
     1702        write(1) SLAB
     1703
     1704
     1705        FIELD='SO3'
     1706        UNITS='kg/kg'
     1707        DESC='SO3 mixing ratio'
     1708        XLVL=levels(k)
     1709        SLAB=chemtrac(:,:,1,time_out(l),25)
     1710                ! And now put everything in the destination file
     1711                ! ... Header
     1712        write(1) IFV
     1713        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1714        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1715                ! ... Data
     1716        write(1) SLAB
     1717
     1718
     1719        FIELD='S2O2'
     1720        UNITS='kg/kg'
     1721        DESC='S2O2 mixing ratio'
     1722        XLVL=levels(k)
     1723        SLAB=chemtrac(:,:,1,time_out(l),26)
     1724                ! And now put everything in the destination file
     1725                ! ... Header
     1726        write(1) IFV
     1727        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1728        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1729                ! ... Data
     1730        write(1) SLAB
     1731
     1732
     1733        FIELD='OCS'
     1734        UNITS='kg/kg'
     1735        DESC='OCS mixing ratio'
     1736        XLVL=levels(k)
     1737        SLAB=chemtrac(:,:,1,time_out(l),27)
     1738                ! And now put everything in the destination file
     1739                ! ... Header
     1740        write(1) IFV
     1741        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1742        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1743                ! ... Data
     1744        write(1) SLAB
     1745
     1746
     1747        FIELD='HSO3'
     1748        UNITS='kg/kg'
     1749        DESC='HSO3 mixing ratio'
     1750        XLVL=levels(k)
     1751        SLAB=chemtrac(:,:,1,time_out(l),28)
     1752                ! And now put everything in the destination file
     1753                ! ... Header
     1754        write(1) IFV
     1755        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1756        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1757                ! ... Data
     1758        write(1) SLAB
     1759
     1760
     1761        FIELD='H2SO4'
     1762        UNITS='kg/kg'
     1763        DESC='H2SO4 mixing ratio'
     1764        XLVL=levels(k)
     1765        SLAB=chemtrac(:,:,1,time_out(l),29)
     1766                ! And now put everything in the destination file
     1767                ! ... Header
     1768        write(1) IFV
     1769        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1770        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1771                ! ... Data
     1772        write(1) SLAB
     1773
     1774
     1775        FIELD='S2'
     1776        UNITS='kg/kg'
     1777        DESC='S2 mixing ratio'
     1778        XLVL=levels(k)
     1779        SLAB=chemtrac(:,:,1,time_out(l),30)
     1780                ! And now put everything in the destination file
     1781                ! ... Header
     1782        write(1) IFV
     1783        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1784        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1785                ! ... Data
     1786        write(1) SLAB
     1787
     1788
     1789        FIELD='ClSO2'
     1790        UNITS='kg/kg'
     1791        DESC='ClSO2 mixing ratio'
     1792        XLVL=levels(k)
     1793        SLAB=chemtrac(:,:,1,time_out(l),31)
     1794                ! And now put everything in the destination file
     1795                ! ... Header
     1796        write(1) IFV
     1797        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1798        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1799                ! ... Data
     1800        write(1) SLAB
     1801
     1802
     1803        FIELD='OSCl'
     1804        UNITS='kg/kg'
     1805        DESC='OSCl mixing ratio'
     1806        XLVL=levels(k)
     1807        SLAB=chemtrac(:,:,1,time_out(l),32)
     1808                ! And now put everything in the destination file
     1809                ! ... Header
     1810        write(1) IFV
     1811        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1812        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1813                ! ... Data
     1814        write(1) SLAB
     1815
     1816
     1817        FIELD='H2Oliq'
     1818        UNITS='kg/kg'
     1819        DESC='H2Oliq mixing ratio'
     1820        XLVL=levels(k)
     1821        SLAB=chemtrac(:,:,1,time_out(l),33)
     1822                ! And now put everything in the destination file
     1823                ! ... Header
     1824        write(1) IFV
     1825        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1826        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1827                ! ... Data
     1828        write(1) SLAB
     1829
     1830
     1831        FIELD='H2SO4liq'
     1832        UNITS='kg/kg'
     1833        DESC='H2SO4liq mixing ratio'
     1834        XLVL=levels(k)
     1835        SLAB=chemtrac(:,:,1,time_out(l),34)
     1836                ! And now put everything in the destination file
     1837                ! ... Header
     1838        write(1) IFV
     1839        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     1840        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     1841                ! ... Data
     1842        write(1) SLAB
    13591843#endif
    13601844
    … …  
    16862170!------------------------!
    16872171#ifdef PHOTOCHEM
    1688     DO i=1,nchemtrac
    1689        FIELD=wtnom(i)
    1690        UNITS='units'
    1691        DESC='desc'
    1692        DO k = 1,altlen
    1693          XLVL=levels(k)
    1694          SLAB=chemtrac(:,:,k,time_out(l),i)
    1695          ! And now put everything in the destination file
    1696          ! ... Header
    1697          write(1) IFV
    1698          write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
    1699          write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
    1700          ! ... Data
    1701          write(1) SLAB
    1702        END DO
    1703     ENDDO
     2172FIELD='CO2'
     2173UNITS='kg/kg'
     2174DESC='CO2 mixing ratio'
     2175DO k = 1,altlen
     2176        XLVL=levels(k)
     2177        SLAB=chemtrac(:,:,k,time_out(l),1)
     2178                ! And now put everything in the destination file
     2179                ! ... Header
     2180        write(1) IFV
     2181        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2182        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2183                ! ... Data
     2184                write(1) SLAB
     2185END DO
     2186
     2187FIELD='CO'
     2188UNITS='kg/kg'
     2189DESC='CO mixing ratio'
     2190DO k = 1,altlen
     2191        XLVL=levels(k)
     2192        SLAB=chemtrac(:,:,k,time_out(l),2)
     2193                ! And now put everything in the destination file
     2194                ! ... Header
     2195        write(1) IFV
     2196        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2197        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2198                ! ... Data
     2199                write(1) SLAB
     2200END DO
     2201
     2202FIELD='H2'
     2203UNITS='kg/kg'
     2204DESC='H2 mixing ratio'
     2205DO k = 1,altlen
     2206        XLVL=levels(k)
     2207        SLAB=chemtrac(:,:,k,time_out(l),3)
     2208                ! And now put everything in the destination file
     2209                ! ... Header
     2210        write(1) IFV
     2211        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2212        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2213                ! ... Data
     2214                write(1) SLAB
     2215END DO
     2216
     2217FIELD='H2O'
     2218UNITS='kg/kg'
     2219DESC='H2O mixing ratio'
     2220DO k = 1,altlen
     2221        XLVL=levels(k)
     2222        SLAB=chemtrac(:,:,k,time_out(l),4)
     2223                ! And now put everything in the destination file
     2224                ! ... Header
     2225        write(1) IFV
     2226        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2227        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2228                ! ... Data
     2229                write(1) SLAB
     2230END DO
     2231
     2232FIELD='O1D'
     2233UNITS='kg/kg'
     2234DESC='O1D mixing ratio'
     2235DO k = 1,altlen
     2236        XLVL=levels(k)
     2237        SLAB=chemtrac(:,:,k,time_out(l),5)
     2238                ! And now put everything in the destination file
     2239                ! ... Header
     2240        write(1) IFV
     2241        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2242        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2243                ! ... Data
     2244                write(1) SLAB
     2245END DO
     2246
     2247FIELD='O'
     2248UNITS='kg/kg'
     2249DESC='O mixing ratio'
     2250DO k = 1,altlen
     2251        XLVL=levels(k)
     2252        SLAB=chemtrac(:,:,k,time_out(l),6)
     2253                ! And now put everything in the destination file
     2254                ! ... Header
     2255        write(1) IFV
     2256        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2257        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2258                ! ... Data
     2259                write(1) SLAB
     2260END DO
     2261
     2262FIELD='O2'
     2263UNITS='kg/kg'
     2264DESC='O2 mixing ratio'
     2265DO k = 1,altlen
     2266        XLVL=levels(k)
     2267        SLAB=chemtrac(:,:,k,time_out(l),7)
     2268                ! And now put everything in the destination file
     2269                ! ... Header
     2270        write(1) IFV
     2271        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2272        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2273                ! ... Data
     2274                write(1) SLAB
     2275END DO
     2276
     2277FIELD='O2DG'
     2278UNITS='kg/kg'
     2279DESC='O2DG mixing ratio'
     2280DO k = 1,altlen
     2281        XLVL=levels(k)
     2282        SLAB=chemtrac(:,:,k,time_out(l),8)
     2283                ! And now put everything in the destination file
     2284                ! ... Header
     2285        write(1) IFV
     2286        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2287        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2288                ! ... Data
     2289                write(1) SLAB
     2290END DO
     2291
     2292FIELD='O3'
     2293UNITS='kg/kg'
     2294DESC='O3 mixing ratio'
     2295DO k = 1,altlen
     2296        XLVL=levels(k)
     2297        SLAB=chemtrac(:,:,k,time_out(l),9)
     2298                ! And now put everything in the destination file
     2299                ! ... Header
     2300        write(1) IFV
     2301        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2302        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2303                ! ... Data
     2304                write(1) SLAB
     2305END DO
     2306
     2307FIELD='H'
     2308UNITS='kg/kg'
     2309DESC='H mixing ratio'
     2310DO k = 1,altlen
     2311        XLVL=levels(k)
     2312        SLAB=chemtrac(:,:,k,time_out(l),10)
     2313                ! And now put everything in the destination file
     2314                ! ... Header
     2315        write(1) IFV
     2316        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2317        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2318                ! ... Data
     2319                write(1) SLAB
     2320END DO
     2321
     2322FIELD='OH'
     2323UNITS='kg/kg'
     2324DESC='OH mixing ratio'
     2325DO k = 1,altlen
     2326        XLVL=levels(k)
     2327        SLAB=chemtrac(:,:,k,time_out(l),11)
     2328                ! And now put everything in the destination file
     2329                ! ... Header
     2330        write(1) IFV
     2331        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2332        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2333                ! ... Data
     2334                write(1) SLAB
     2335END DO
     2336
     2337FIELD='HO2'
     2338UNITS='kg/kg'
     2339DESC='hO2 mixing ratio'
     2340DO k = 1,altlen
     2341        XLVL=levels(k)
     2342        SLAB=chemtrac(:,:,k,time_out(l),12)
     2343                ! And now put everything in the destination file
     2344                ! ... Header
     2345        write(1) IFV
     2346        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2347        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2348                ! ... Data
     2349                write(1) SLAB
     2350END DO
     2351
     2352FIELD='H2O2'
     2353UNITS='kg/kg'
     2354DESC='H2O2 mixing ratio'
     2355DO k = 1,altlen
     2356        XLVL=levels(k)
     2357        SLAB=chemtrac(:,:,k,time_out(l),13)
     2358                ! And now put everything in the destination file
     2359                ! ... Header
     2360        write(1) IFV
     2361        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2362        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2363                ! ... Data
     2364                write(1) SLAB
     2365END DO
     2366
     2367FIELD='Cl'
     2368UNITS='kg/kg'
     2369DESC='Cl mixing ratio'
     2370DO k = 1,altlen
     2371        XLVL=levels(k)
     2372        SLAB=chemtrac(:,:,k,time_out(l),14)
     2373                ! And now put everything in the destination file
     2374                ! ... Header
     2375        write(1) IFV
     2376        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2377        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2378                ! ... Data
     2379                write(1) SLAB
     2380END DO
     2381
     2382FIELD='ClO'
     2383UNITS='kg/kg'
     2384DESC='ClO mixing ratio'
     2385DO k = 1,altlen
     2386        XLVL=levels(k)
     2387        SLAB=chemtrac(:,:,k,time_out(l),15)
     2388                ! And now put everything in the destination file
     2389                ! ... Header
     2390        write(1) IFV
     2391        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2392        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2393                ! ... Data
     2394                write(1) SLAB
     2395END DO
     2396
     2397FIELD='Cl2'
     2398UNITS='kg/kg'
     2399DESC='Cl2 mixing ratio'
     2400DO k = 1,altlen
     2401        XLVL=levels(k)
     2402        SLAB=chemtrac(:,:,k,time_out(l),16)
     2403                ! And now put everything in the destination file
     2404                ! ... Header
     2405        write(1) IFV
     2406        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2407        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2408                ! ... Data
     2409                write(1) SLAB
     2410END DO
     2411
     2412FIELD='HCl'
     2413UNITS='kg/kg'
     2414DESC='HCl mixing ratio'
     2415DO k = 1,altlen
     2416        XLVL=levels(k)
     2417        SLAB=chemtrac(:,:,k,time_out(l),17)
     2418                ! And now put everything in the destination file
     2419                ! ... Header
     2420        write(1) IFV
     2421        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2422        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2423                ! ... Data
     2424                write(1) SLAB
     2425END DO
     2426
     2427FIELD='HOCl'
     2428UNITS='kg/kg'
     2429DESC='HOCl mixing ratio'
     2430DO k = 1,altlen
     2431        XLVL=levels(k)
     2432        SLAB=chemtrac(:,:,k,time_out(l),18)
     2433                ! And now put everything in the destination file
     2434                ! ... Header
     2435        write(1) IFV
     2436        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2437        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2438                ! ... Data
     2439                write(1) SLAB
     2440END DO
     2441
     2442FIELD='ClCO'
     2443UNITS='kg/kg'
     2444DESC='ClCO mixing ratio'
     2445DO k = 1,altlen
     2446        XLVL=levels(k)
     2447        SLAB=chemtrac(:,:,k,time_out(l),19)
     2448                ! And now put everything in the destination file
     2449                ! ... Header
     2450        write(1) IFV
     2451        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2452        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2453                ! ... Data
     2454                write(1) SLAB
     2455END DO
     2456
     2457FIELD='ClCO3'
     2458UNITS='kg/kg'
     2459DESC='ClCO3 mixing ratio'
     2460DO k = 1,altlen
     2461        XLVL=levels(k)
     2462        SLAB=chemtrac(:,:,k,time_out(l),20)
     2463                ! And now put everything in the destination file
     2464                ! ... Header
     2465        write(1) IFV
     2466        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2467        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2468                ! ... Data
     2469                write(1) SLAB
     2470END DO
     2471
     2472FIELD='COCl2'
     2473UNITS='kg/kg'
     2474DESC='COClC2 mixing ratio'
     2475DO k = 1,altlen
     2476        XLVL=levels(k)
     2477        SLAB=chemtrac(:,:,k,time_out(l),21)
     2478                ! And now put everything in the destination file
     2479                ! ... Header
     2480        write(1) IFV
     2481        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2482        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2483                ! ... Data
     2484                write(1) SLAB
     2485END DO
     2486
     2487FIELD='S'
     2488UNITS='kg/kg'
     2489DESC='S mixing ratio'
     2490DO k = 1,altlen
     2491        XLVL=levels(k)
     2492        SLAB=chemtrac(:,:,k,time_out(l),22)
     2493                ! And now put everything in the destination file
     2494                ! ... Header
     2495        write(1) IFV
     2496        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2497        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2498                ! ... Data
     2499                write(1) SLAB
     2500END DO
     2501
     2502FIELD='SO'
     2503UNITS='kg/kg'
     2504DESC='SO mixing ratio'
     2505DO k = 1,altlen
     2506        XLVL=levels(k)
     2507        SLAB=chemtrac(:,:,k,time_out(l),23)
     2508                ! And now put everything in the destination file
     2509                ! ... Header
     2510        write(1) IFV
     2511        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2512        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2513                ! ... Data
     2514                write(1) SLAB
     2515END DO
     2516
     2517FIELD='SO2'
     2518UNITS='kg/kg'
     2519DESC='SO2 mixing ratio'
     2520DO k = 1,altlen
     2521        XLVL=levels(k)
     2522        SLAB=chemtrac(:,:,k,time_out(l),24)
     2523                ! And now put everything in the destination file
     2524                ! ... Header
     2525        write(1) IFV
     2526        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2527        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2528                ! ... Data
     2529                write(1) SLAB
     2530END DO
     2531
     2532FIELD='SO3'
     2533UNITS='kg/kg'
     2534DESC='SO3 mixing ratio'
     2535DO k = 1,altlen
     2536        XLVL=levels(k)
     2537        SLAB=chemtrac(:,:,k,time_out(l),25)
     2538                ! And now put everything in the destination file
     2539                ! ... Header
     2540        write(1) IFV
     2541        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2542        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2543                ! ... Data
     2544                write(1) SLAB
     2545END DO
     2546
     2547FIELD='S2O2'
     2548UNITS='kg/kg'
     2549DESC='S2O2 mixing ratio'
     2550DO k = 1,altlen
     2551        XLVL=levels(k)
     2552        SLAB=chemtrac(:,:,k,time_out(l),26)
     2553                ! And now put everything in the destination file
     2554                ! ... Header
     2555        write(1) IFV
     2556        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2557        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2558                ! ... Data
     2559                write(1) SLAB
     2560END DO
     2561
     2562FIELD='OCS'
     2563UNITS='kg/kg'
     2564DESC='OCS mixing ratio'
     2565DO k = 1,altlen
     2566        XLVL=levels(k)
     2567        SLAB=chemtrac(:,:,k,time_out(l),27)
     2568                ! And now put everything in the destination file
     2569                ! ... Header
     2570        write(1) IFV
     2571        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2572        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2573                ! ... Data
     2574                write(1) SLAB
     2575END DO
     2576
     2577FIELD='HSO3'
     2578UNITS='kg/kg'
     2579DESC='HSO3 mixing ratio'
     2580DO k = 1,altlen
     2581        XLVL=levels(k)
     2582        SLAB=chemtrac(:,:,k,time_out(l),28)
     2583                ! And now put everything in the destination file
     2584                ! ... Header
     2585        write(1) IFV
     2586        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2587        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2588                ! ... Data
     2589                write(1) SLAB
     2590END DO
     2591
     2592FIELD='H2SO4'
     2593UNITS='kg/kg'
     2594DESC='H2SO4 mixing ratio'
     2595DO k = 1,altlen
     2596        XLVL=levels(k)
     2597        SLAB=chemtrac(:,:,k,time_out(l),29)
     2598                ! And now put everything in the destination file
     2599                ! ... Header
     2600        write(1) IFV
     2601        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2602        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2603                ! ... Data
     2604                write(1) SLAB
     2605END DO
     2606
     2607FIELD='S2'
     2608UNITS='kg/kg'
     2609DESC='S2 mixing ratio'
     2610DO k = 1,altlen
     2611        XLVL=levels(k)
     2612        SLAB=chemtrac(:,:,k,time_out(l),30)
     2613                ! And now put everything in the destination file
     2614                ! ... Header
     2615        write(1) IFV
     2616        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2617        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2618                ! ... Data
     2619                write(1) SLAB
     2620END DO
     2621
     2622FIELD='ClSO2'
     2623UNITS='kg/kg'
     2624DESC='ClSO2 mixing ratio'
     2625DO k = 1,altlen
     2626        XLVL=levels(k)
     2627        SLAB=chemtrac(:,:,k,time_out(l),31)
     2628                ! And now put everything in the destination file
     2629                ! ... Header
     2630        write(1) IFV
     2631        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2632        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2633                ! ... Data
     2634                write(1) SLAB
     2635END DO
     2636
     2637FIELD='OSCl'
     2638UNITS='kg/kg'
     2639DESC='OSCl mixing ratio'
     2640DO k = 1,altlen
     2641        XLVL=levels(k)
     2642        SLAB=chemtrac(:,:,k,time_out(l),32)
     2643                ! And now put everything in the destination file
     2644                ! ... Header
     2645        write(1) IFV
     2646        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2647        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2648                ! ... Data
     2649                write(1) SLAB
     2650END DO
     2651
     2652FIELD='H2Oliq'
     2653UNITS='kg/kg'
     2654DESC='H2Oliq mixing ratio'
     2655DO k = 1,altlen
     2656        XLVL=levels(k)
     2657        SLAB=chemtrac(:,:,k,time_out(l),33)
     2658                ! And now put everything in the destination file
     2659                ! ... Header
     2660        write(1) IFV
     2661        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2662        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2663                ! ... Data
     2664                write(1) SLAB
     2665END DO
     2666
     2667FIELD='H2SO4liq'
     2668UNITS='kg/kg'
     2669DESC='H2SO4liq mixing ratio'
     2670DO k = 1,altlen
     2671        XLVL=levels(k)
     2672        SLAB=chemtrac(:,:,k,time_out(l),34)
     2673                ! And now put everything in the destination file
     2674                ! ... Header
     2675        write(1) IFV
     2676        write(1) HDATE,XFCST,SOURCE,FIELD,UNITS,DESC,XLVL,NX,NY,IPROJ
     2677        write(1) STARTLOC,STARTLAT,STARTLON,DELTALAT,DELTALON
     2678                ! ... Data
     2679                write(1) SLAB
     2680END DO
    17042681#endif
    17052682

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM

    r2018 r2021  
    171171####MARS MARS MARS MARS MARS
    172172####MARS MARS MARS MARS MARS
     173###VENUS TRACER
     174state  real  qco       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO"        "CO number mixing ratio"   "kg kg-1"
     175state  real  qh2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2"        "H2 number mixing ratio"   "kg kg-1"
     176state  real  qo1d      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO1D"        "O1d number mixing ratio"   "kg kg-1"
     177state  real  qo        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO"        "O number mixing ratio"   "kg kg-1"
     178state  real  qo2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2"        "O2 number mixing ratio"   "kg kg-1"
     179state  real  qo2dg     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2dg"        "O2dg number mixing ratio"   "kg kg-1"
     180state  real  qo3       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO3"        "O3 number mixing ratio"   "kg kg-1"
     181state  real  qh        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH"        "H number mixing ratio"   "kg kg-1"
     182state  real  qoh       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOH"        "OH number mixing ratio"   "kg kg-1"
     183state  real  qho2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHO2"        "HO2 number mixing ratio"   "kg kg-1"
     184state  real  qh2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2O2"        "H2O2 number mixing ratio"   "kg kg-1"
     185state  real  qcl       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL"        "Cl number mixing ratio"   "kg kg-1"
     186state  real  qclo      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLO"        "ClO number mixing ratio"   "kg kg-1"
     187state  real  qcl2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL2"        "Cl2 number mixing ratio"   "kg kg-1"
     188state  real  qhcl      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHCL"        "HCl number mixing ratio"   "kg kg-1"
     189state  real  qhocl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHOCL"        "HOCl number mixing ratio"   "kg kg-1"
     190state  real  qclco     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO"        "ClCO number mixing ratio"   "kg kg-1"
     191state  real  qclco3    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO3"        "ClCO3 number mixing ratio"   "kg kg-1"
     192state  real  qcocl2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCOCL2"        "COCl2 number mixing ratio"   "kg kg-1"
     193state  real  qss       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS"        "S number mixing ratio"   "kg kg-1"
     194state  real  qso       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO"        "SO number mixing ratio"   "kg kg-1"
     195state  real  qso2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO2"        "SO2 number mixing ratio"   "kg kg-1"
     196state  real  qso3      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO3"        "SO3 number mixing ratio"   "kg kg-1"
     197state  real  qs2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2O2"        "S2O2 number mixing ratio"   "kg kg-1"
     198state  real  qocs      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOCS"        "OCS number mixing ratio"   "kg kg-1"
     199state  real  qhso3     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHSO3"        "HSO3 number mixing ratio"   "kg kg-1"
     200state  real  qh2so4    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4"        "H2SO4 number mixing ratio"   "kg kg-1"
     201state  real  qs2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2"        "S2 number mixing ratio"   "kg kg-1"
     202state  real  qclso2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLSO2"        "ClSO2 number mixing ratio"   "kg kg-1"
     203state  real  qoscl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOSCL"        "OSCl number mixing ratio"   "kg kg-1"
     204state  real  qh2oliq   ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2OLIQ"        "H2O liq number mixing ratio"   "kg kg-1"
     205state  real  qh2so4liq ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4LIQ"        "H2SO4 liq number mixing ratio"   "kg kg-1"
    173206
    174207
    … …  
    13921425package   newwater     mars==11                     -              scalar:qh2o,qh2o_ice,qdust,qdustn
    13931426package   newwatersca  mars==12                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn
     1427package   VenusPC      mars==34                     -              scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq
    13941428package   photochem    mars==42                     -              scalar:qco2,chem_co,chem_o,chem_o1d,chem_o2,chem_o3,chem_h,chem_h2,chem_oh,chem_ho2,chem_h2o2,chem_ch4,chem_n2,chem_ar,qh2o_ice,qh2o,qdust,qdustn
    13951429##### MARS OPTIONS

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/Registry/Registry.EM.newphys

    r2018 r2021  
    156156####MARS MARS MARS MARS MARS
    157157####MARS MARS MARS MARS MARS
    158 
    159 
     158###VENUS TRACER
     159state  real  qco       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCO"        "CO number mixing ratio"   "kg kg-1"
     160state  real  qh2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2"        "H2 number mixing ratio"   "kg kg-1"
     161state  real  qo1d      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO1D"        "O1d number mixing ratio"   "kg kg-1"
     162state  real  qo        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO"        "O number mixing ratio"   "kg kg-1"
     163state  real  qo2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2"        "O2 number mixing ratio"   "kg kg-1"
     164state  real  qo2dg     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO2dg"        "O2dg number mixing ratio"   "kg kg-1"
     165state  real  qo3       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QO3"        "O3 number mixing ratio"   "kg kg-1"
     166state  real  qh        ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH"        "H number mixing ratio"   "kg kg-1"
     167state  real  qoh       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOH"        "OH number mixing ratio"   "kg kg-1"
     168state  real  qho2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHO2"        "HO2 number mixing ratio"   "kg kg-1"
     169state  real  qh2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2O2"        "H2O2 number mixing ratio"   "kg kg-1"
     170state  real  qcl       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL"        "Cl number mixing ratio"   "kg kg-1"
     171state  real  qclo      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLO"        "ClO number mixing ratio"   "kg kg-1"
     172state  real  qcl2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCL2"        "Cl2 number mixing ratio"   "kg kg-1"
     173state  real  qhcl      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHCL"        "HCl number mixing ratio"   "kg kg-1"
     174state  real  qhocl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHOCL"        "HOCl number mixing ratio"   "kg kg-1"
     175state  real  qclco     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO"        "ClCO number mixing ratio"   "kg kg-1"
     176state  real  qclco3    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLCO3"        "ClCO3 number mixing ratio"   "kg kg-1"
     177state  real  qcocl2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCOCL2"        "COCl2 number mixing ratio"   "kg kg-1"
     178state  real  qss       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS"        "S number mixing ratio"   "kg kg-1"
     179state  real  qso       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO"        "SO number mixing ratio"   "kg kg-1"
     180state  real  qso2      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO2"        "SO2 number mixing ratio"   "kg kg-1"
     181state  real  qso3      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QSO3"        "SO3 number mixing ratio"   "kg kg-1"
     182state  real  qs2o2     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2O2"        "S2O2 number mixing ratio"   "kg kg-1"
     183state  real  qocs      ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOCS"        "OCS number mixing ratio"   "kg kg-1"
     184state  real  qhso3     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QHSO3"        "HSO3 number mixing ratio"   "kg kg-1"
     185state  real  qh2so4    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4"        "H2SO4 number mixing ratio"   "kg kg-1"
     186state  real  qs2       ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QS2"        "S2 number mixing ratio"   "kg kg-1"
     187state  real  qclso2    ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QCLSO2"        "ClSO2 number mixing ratio"   "kg kg-1"
     188state  real  qoscl     ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QOSCL"        "OSCl number mixing ratio"   "kg kg-1"
     189state  real  qh2oliq   ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2OLIQ"        "H2O liq number mixing ratio"   "kg kg-1"
     190state  real  qh2so4liq ikjftb  scalar  1  -  i01rhusdf=(bdy_interp:dt)  "QH2SO4LIQ"        "H2SO4 liq number mixing ratio"   "kg kg-1"
    160191
    161192################################################################################
    … …  
    13781409package   newwatersca  mars==12                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn
    13791410package   cloudco2     mars==32                     -              scalar:qh2o,qh2o_ice,qdust,qdustn,qccn,qccnn,qco2,qco2_ice,qccn_co2,qccnn_co2
     1411package   VenusPC      mars==34                     -              scalar:qco2,qco,qh2,qh2o,qo1d,qo,qo2,qo2dg,qo3,qh,qoh,qho2,qh2o2,qcl,qclo,qcl2,qhcl,qhocl,qclco,qclco3,qcocl2,qss,qso,qso2,qso3,qs2o2,qocs,qhso3,qh2so4,qs2,qclso2,qoscl,qh2oliq,qh2so4liq
    13801412package   photochem    mars==42                     -              scalar:qco2,chem_co,chem_o,chem_o1d,chem_o2,chem_o3,chem_h,chem_h2,chem_oh,chem_ho2,chem_h2o2,chem_ch4,chem_n2,chem_ar,qh2o_ice,qh2o,qdust,qdustn
    13811413##### MARS OPTIONS

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/module_initialize_real.F

    r1750 r2021  
    14181418#endif
    14191419
     1420!VENUS photochemistry
     1421if ( config_flags%mars == 34 ) then
     1422       print*,'grid%em_qco2_gc',grid%em_qco2_gc(0,:,0)
     1423       CALL vert_interp_old ( grid%em_qco2_gc , grid%em_pd_gc , scalar(:,:,:,2) , grid%em_pb , &
     1424                                    num_metgrid_levels , 'Q' , &
     1425                                    interp_type , lagrange_order , lowest_lev_from_sfc , &
     1426                                    zap_close_levels , force_sfc_in_vinterp , &
     1427                                    ids , ide , jds , jde , kds , kde , &
     1428                                    ims , ime , jms , jme , kms , kme , &
     1429                                    its , ite , jts , jte , kts , kte )
     1430
     1431       CALL vert_interp_old ( grid%em_qco_gc , grid%em_pd_gc , scalar(:,:,:,3), grid%em_pb , &
     1432                                    num_metgrid_levels , 'Q' , &
     1433                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
     1434                                    zap_close_levels , force_sfc_in_vinterp , &
     1435                                    ids , ide , jds , jde , kds , kde , &
     1436                                    ims , ime , jms , jme , kms , kme , &
     1437                                    its , ite , jts , jte , kts , kte )
     1438
     1439       CALL vert_interp_old ( grid%em_qh2_gc , grid%em_pd_gc , scalar(:,:,:,4),grid%em_pb , &
     1440                                    num_metgrid_levels , 'Q' , &
     1441                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1442                                    zap_close_levels , force_sfc_in_vinterp , &
     1443                                    ids , ide , jds , jde , kds , kde , &
     1444                                    ims , ime , jms , jme , kms , kme , &
     1445                                    its , ite , jts , jte , kts , kte )
     1446
     1447       CALL vert_interp_old ( grid%em_qh2o_gc , grid%em_pd_gc , scalar(:,:,:,5),grid%em_pb , &
     1448                                    num_metgrid_levels , 'Q' , &
     1449                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1450                                    zap_close_levels , force_sfc_in_vinterp , &
     1451                                    ids , ide , jds , jde , kds , kde , &
     1452                                    ims , ime , jms , jme , kms , kme , &
     1453                                    its , ite , jts , jte , kts , kte )
     1454
     1455       CALL vert_interp_old ( grid%em_qo1d_gc , grid%em_pd_gc , scalar(:,:,:,6),grid%em_pb , &
     1456                                    num_metgrid_levels , 'Q' , &
     1457                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1458                                    zap_close_levels , force_sfc_in_vinterp , &
     1459                                    ids , ide , jds , jde , kds , kde , &
     1460                                    ims , ime , jms , jme , kms , kme , &
     1461                                    its , ite , jts , jte , kts , kte )
     1462
     1463       CALL vert_interp_old ( grid%em_qo_gc , grid%em_pd_gc , scalar(:,:,:,7),grid%em_pb , &
     1464                                    num_metgrid_levels , 'Q' , &
     1465                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1466                                    zap_close_levels , force_sfc_in_vinterp , &
     1467                                    ids , ide , jds , jde , kds , kde , &
     1468                                    ims , ime , jms , jme , kms , kme , &
     1469                                    its , ite , jts , jte , kts , kte )
     1470
     1471       CALL vert_interp_old ( grid%em_qo2_gc , grid%em_pd_gc , scalar(:,:,:,8),grid%em_pb , &
     1472                                    num_metgrid_levels , 'Q' , &
     1473                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1474                                    zap_close_levels , force_sfc_in_vinterp , &
     1475                                    ids , ide , jds , jde , kds , kde , &
     1476                                    ims , ime , jms , jme , kms , kme , &
     1477                                    its , ite , jts , jte , kts , kte )
     1478
     1479       CALL vert_interp_old ( grid%em_qo2dg_gc , grid%em_pd_gc , scalar(:,:,:,9),grid%em_pb , &
     1480                                    num_metgrid_levels , 'Q' , &
     1481                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1482                                    zap_close_levels , force_sfc_in_vinterp , &
     1483                                    ids , ide , jds , jde , kds , kde , &
     1484                                    ims , ime , jms , jme , kms , kme , &
     1485                                    its , ite , jts , jte , kts , kte )
     1486
     1487       CALL vert_interp_old ( grid%em_qo3_gc , grid%em_pd_gc , scalar(:,:,:,10),grid%em_pb , &
     1488                                    num_metgrid_levels , 'Q' , &
     1489                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1490                                    zap_close_levels , force_sfc_in_vinterp , &
     1491                                    ids , ide , jds , jde , kds , kde , &
     1492                                    ims , ime , jms , jme , kms , kme , &
     1493                                    its , ite , jts , jte , kts , kte )
     1494
     1495       CALL vert_interp_old ( grid%em_qh_gc , grid%em_pd_gc , scalar(:,:,:,11),grid%em_pb , &
     1496                                    num_metgrid_levels , 'Q' , &
     1497                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1498                                    zap_close_levels , force_sfc_in_vinterp , &
     1499                                    ids , ide , jds , jde , kds , kde , &
     1500                                    ims , ime , jms , jme , kms , kme , &
     1501                                    its , ite , jts , jte , kts , kte )
     1502
     1503       CALL vert_interp_old ( grid%em_qoh_gc , grid%em_pd_gc , scalar(:,:,:,12), grid%em_pb , &
     1504                                    num_metgrid_levels , 'Q' , &
     1505                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1506                                    zap_close_levels , force_sfc_in_vinterp , &
     1507                                    ids , ide , jds , jde , kds , kde , &
     1508                                    ims , ime , jms , jme , kms , kme , &
     1509                                    its , ite , jts , jte , kts , kte )
     1510
     1511       CALL vert_interp_old ( grid%em_qho2_gc , grid%em_pd_gc , scalar(:,:,:,13),grid%em_pb , &
     1512                                    num_metgrid_levels , 'Q' , &
     1513                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1514                                    zap_close_levels , force_sfc_in_vinterp , &
     1515                                    ids , ide , jds , jde , kds , kde , &
     1516                                    ims , ime , jms , jme , kms , kme , &
     1517                                    its , ite , jts , jte , kts , kte )
     1518
     1519       CALL vert_interp_old ( grid%em_qh2o2_gc , grid%em_pd_gc , scalar(:,:,:,14),grid%em_pb , &
     1520                                    num_metgrid_levels , 'Q' , &
     1521                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1522                                    zap_close_levels , force_sfc_in_vinterp , &
     1523                                    ids , ide , jds , jde , kds , kde , &
     1524                                    ims , ime , jms , jme , kms , kme , &
     1525                                    its , ite , jts , jte , kts , kte )
     1526
     1527       CALL vert_interp_old ( grid%em_qcl_gc , grid%em_pd_gc , scalar(:,:,:,15),grid%em_pb , &
     1528                                    num_metgrid_levels , 'Q' , &
     1529                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
     1530                                    zap_close_levels , force_sfc_in_vinterp , &
     1531                                    ids , ide , jds , jde , kds , kde , &
     1532                                    ims , ime , jms , jme , kms , kme , &
     1533                                    its , ite , jts , jte , kts , kte )
     1534
     1535       CALL vert_interp_old ( grid%em_qclo_gc , grid%em_pd_gc , scalar(:,:,:,16),grid%em_pb , &
     1536                                    num_metgrid_levels , 'Q' , &
     1537                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1538                                    zap_close_levels , force_sfc_in_vinterp , &
     1539                                    ids , ide , jds , jde , kds , kde , &
     1540                                    ims , ime , jms , jme , kms , kme , &
     1541                                    its , ite , jts , jte , kts , kte )
     1542
     1543       CALL vert_interp_old ( grid%em_qcl2_gc , grid%em_pd_gc , scalar(:,:,:,17),grid%em_pb , &
     1544                                    num_metgrid_levels , 'Q' , &
     1545                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1546                                    zap_close_levels , force_sfc_in_vinterp , &
     1547                                    ids , ide , jds , jde , kds , kde , &
     1548                                    ims , ime , jms , jme , kms , kme , &
     1549                                    its , ite , jts , jte , kts , kte )
     1550
     1551       CALL vert_interp_old ( grid%em_qhcl_gc , grid%em_pd_gc , scalar(:,:,:,18),grid%em_pb , &
     1552                                    num_metgrid_levels , 'Q' , &
     1553                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1554                                    zap_close_levels , force_sfc_in_vinterp , &
     1555                                    ids , ide , jds , jde , kds , kde , &
     1556                                    ims , ime , jms , jme , kms , kme , &
     1557                                    its , ite , jts , jte , kts , kte )
     1558
     1559       CALL vert_interp_old ( grid%em_qhocl_gc , grid%em_pd_gc , scalar(:,:,:,19),grid%em_pb , &
     1560                                    num_metgrid_levels , 'Q' , &
     1561                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1562                                    zap_close_levels , force_sfc_in_vinterp , &
     1563                                    ids , ide , jds , jde , kds , kde , &
     1564                                    ims , ime , jms , jme , kms , kme , &
     1565                                    its , ite , jts , jte , kts , kte )
     1566
     1567       CALL vert_interp_old ( grid%em_qclco_gc , grid%em_pd_gc , scalar(:,:,:,20),grid%em_pb , &
     1568                                    num_metgrid_levels , 'Q' , &
     1569                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1570                                    zap_close_levels , force_sfc_in_vinterp , &
     1571                                    ids , ide , jds , jde , kds , kde , &
     1572                                    ims , ime , jms , jme , kms , kme , &
     1573                                    its , ite , jts , jte , kts , kte )
     1574
     1575       CALL vert_interp_old ( grid%em_qclco3_gc , grid%em_pd_gc , scalar(:,:,:,21),grid%em_pb , &
     1576                                    num_metgrid_levels , 'Q' , &
     1577                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1578                                    zap_close_levels , force_sfc_in_vinterp , &
     1579                                    ids , ide , jds , jde , kds , kde , &
     1580                                    ims , ime , jms , jme , kms , kme , &
     1581                                    its , ite , jts , jte , kts , kte )
     1582
     1583       CALL vert_interp_old ( grid%em_qcocl2_gc , grid%em_pd_gc , scalar(:,:,:,22),grid%em_pb , &
     1584                                    num_metgrid_levels , 'Q' , &
     1585                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1586                                    zap_close_levels , force_sfc_in_vinterp , &
     1587                                    ids , ide , jds , jde , kds , kde , &
     1588                                    ims , ime , jms , jme , kms , kme , &
     1589                                    its , ite , jts , jte , kts , kte )
     1590
     1591       CALL vert_interp_old ( grid%em_qs_gc , grid%em_pd_gc , scalar(:,:,:,23), grid%em_pb , &
     1592                                    num_metgrid_levels , 'Q' , &
     1593                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
     1594                                    zap_close_levels , force_sfc_in_vinterp , &
     1595                                    ids , ide , jds , jde , kds , kde , &
     1596                                    ims , ime , jms , jme , kms , kme , &
     1597                                    its , ite , jts , jte , kts , kte )
     1598
     1599       CALL vert_interp_old ( grid%em_qso_gc , grid%em_pd_gc , scalar(:,:,:,24),grid%em_pb , &
     1600                                    num_metgrid_levels , 'Q' , &
     1601                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1602                                    zap_close_levels , force_sfc_in_vinterp , &
     1603                                    ids , ide , jds , jde , kds , kde , &
     1604                                    ims , ime , jms , jme , kms , kme , &
     1605                                    its , ite , jts , jte , kts , kte )
     1606
     1607       CALL vert_interp_old ( grid%em_qso2_gc , grid%em_pd_gc , scalar(:,:,:,25), grid%em_pb , &
     1608                                    num_metgrid_levels , 'Q' , &
     1609                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
     1610                                    zap_close_levels , force_sfc_in_vinterp , &
     1611                                    ids , ide , jds , jde , kds , kde , &
     1612                                    ims , ime , jms , jme , kms , kme , &
     1613                                    its , ite , jts , jte , kts , kte )
     1614
     1615       CALL vert_interp_old ( grid%em_qo3_gc , grid%em_pd_gc , scalar(:,:,:,26),grid%em_pb , &
     1616                                    num_metgrid_levels , 'Q' , &
     1617                                    interp_type , lagrange_order ,lowest_lev_from_sfc , &
     1618                                    zap_close_levels , force_sfc_in_vinterp , &
     1619                                    ids , ide , jds , jde , kds , kde , &
     1620                                    ims , ime , jms , jme , kms , kme , &
     1621                                    its , ite , jts , jte , kts , kte )
     1622
     1623       CALL vert_interp_old ( grid%em_qs2o2_gc , grid%em_pd_gc , scalar(:,:,:,27),grid%em_pb , &
     1624                                    num_metgrid_levels , 'Q' , &
     1625                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1626                                    zap_close_levels , force_sfc_in_vinterp , &
     1627                                    ids , ide , jds , jde , kds , kde , &
     1628                                    ims , ime , jms , jme , kms , kme , &
     1629                                    its , ite , jts , jte , kts , kte )
     1630
     1631       CALL vert_interp_old ( grid%em_qocs_gc , grid%em_pd_gc , scalar(:,:,:,28),grid%em_pb , &
     1632                                    num_metgrid_levels , 'Q' , &
     1633                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1634                                    zap_close_levels , force_sfc_in_vinterp , &
     1635                                    ids , ide , jds , jde , kds , kde , &
     1636                                    ims , ime , jms , jme , kms , kme , &
     1637                                    its , ite , jts , jte , kts , kte )
     1638
     1639       CALL vert_interp_old ( grid%em_qhso3_gc , grid%em_pd_gc , scalar(:,:,:,29),grid%em_pb , &
     1640                                    num_metgrid_levels , 'Q' , &
     1641                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1642                                    zap_close_levels , force_sfc_in_vinterp , &
     1643                                    ids , ide , jds , jde , kds , kde , &
     1644                                    ims , ime , jms , jme , kms , kme , &
     1645                                    its , ite , jts , jte , kts , kte )
     1646
     1647       CALL vert_interp_old ( grid%em_qh2so4_gc , grid%em_pd_gc , scalar(:,:,:,30),grid%em_pb , &
     1648                                    num_metgrid_levels , 'Q' , &
     1649                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1650                                    zap_close_levels , force_sfc_in_vinterp , &
     1651                                    ids , ide , jds , jde , kds , kde , &
     1652                                    ims , ime , jms , jme , kms , kme , &
     1653                                    its , ite , jts , jte , kts , kte )
     1654
     1655       CALL vert_interp_old ( grid%em_qs2_gc , grid%em_pd_gc , scalar(:,:,:,31),grid%em_pb , &
     1656                                    num_metgrid_levels , 'Q' , &
     1657                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1658                                    zap_close_levels , force_sfc_in_vinterp , &
     1659                                    ids , ide , jds , jde , kds , kde , &
     1660                                    ims , ime , jms , jme , kms , kme , &
     1661                                    its , ite , jts , jte , kts , kte )
     1662
     1663       CALL vert_interp_old ( grid%em_qclso2_gc , grid%em_pd_gc , scalar(:,:,:,32),grid%em_pb , &
     1664                                    num_metgrid_levels , 'Q' , &
     1665                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1666                                    zap_close_levels , force_sfc_in_vinterp , &
     1667                                    ids , ide , jds , jde , kds , kde , &
     1668                                    ims , ime , jms , jme , kms , kme , &
     1669                                    its , ite , jts , jte , kts , kte )
     1670
     1671       CALL vert_interp_old ( grid%em_qoscl_gc , grid%em_pd_gc , scalar(:,:,:,33),grid%em_pb , &
     1672                                    num_metgrid_levels , 'Q' , &
     1673                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1674                                    zap_close_levels , force_sfc_in_vinterp , &
     1675                                    ids , ide , jds , jde , kds , kde , &
     1676                                    ims , ime , jms , jme , kms , kme , &
     1677                                    its , ite , jts , jte , kts , kte )
     1678
     1679       CALL vert_interp_old ( grid%em_qh2oliq_gc , grid%em_pd_gc , scalar(:,:,:,34),grid%em_pb , &
     1680                                    num_metgrid_levels , 'Q' , &
     1681                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1682                                    zap_close_levels , force_sfc_in_vinterp , &
     1683                                    ids , ide , jds , jde , kds , kde , &
     1684                                    ims , ime , jms , jme , kms , kme , &
     1685                                    its , ite , jts , jte , kts , kte )
     1686
     1687       CALL vert_interp_old ( grid%em_qh2so4liq_gc , grid%em_pd_gc , scalar(:,:,:,35),grid%em_pb , &
     1688                                    num_metgrid_levels , 'Q' , &
     1689                                    interp_type , lagrange_order,lowest_lev_from_sfc , &
     1690                                    zap_close_levels , force_sfc_in_vinterp , &
     1691                                    ids , ide , jds , jde , kds , kde , &
     1692                                    ims , ime , jms , jme , kms , kme , &
     1693                                    its , ite , jts , jte , kts , kte )
     1694
     1695endif
     1696
    14201697!!! we want any scalar (i.e. tracer) to be positive
    14211698!!! and because of interpolation it is possible that negative values occur...

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/dyn_em/solve_em.F

    r2018 r2021  
    27512751 .OR. (  config_flags%mars .eq. 11                    ) &    !! 11  pass EVERYTHING
    27522752 .OR. (  config_flags%mars .eq. 12                    ) &    !! 12  pass EVERYTHING
     2753 .OR. (  config_flags%mars .eq. 34                    ) &
    27532754 .OR. config_flags%nested ) THEN                             !! *   pass all tracers if nested
    27542755         CALL relax_bdy_scalar ( scalar_tend(ims,kms,jms,is),            &

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/main/real_em.F

    r802 r2021  
    568568   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q5bdy3dtemp1
    569569   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q6bdy3dtemp1
     570   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q7bdy3dtemp1
     571   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q8bdy3dtemp1
     572   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q9bdy3dtemp1
     573   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q10bdy3dtemp1
     574   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q11bdy3dtemp1
     575   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q12bdy3dtemp1
     576   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q13bdy3dtemp1
     577   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q14bdy3dtemp1
     578   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q15bdy3dtemp1
     579   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q16bdy3dtemp1
     580   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q17bdy3dtemp1
     581   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q18bdy3dtemp1
     582   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q19bdy3dtemp1
     583   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q20bdy3dtemp1
     584   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q21bdy3dtemp1
     585   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q22bdy3dtemp1
     586   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q23bdy3dtemp1
     587   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q24bdy3dtemp1
     588   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q25bdy3dtemp1
     589   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q26bdy3dtemp1
     590   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q27bdy3dtemp1
     591   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q28bdy3dtemp1
     592   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q29bdy3dtemp1
     593   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q30bdy3dtemp1
     594   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q31bdy3dtemp1
     595   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q32bdy3dtemp1
     596   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q33bdy3dtemp1
     597   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q34bdy3dtemp1
     598   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q35bdy3dtemp1
    570599!!!!***MARS <<
    571600   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: mbdy2dtemp1
    … …  
    577606   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q5bdy3dtemp2
    578607   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q6bdy3dtemp2
     608   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q7bdy3dtemp2
     609   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q8bdy3dtemp2
     610   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q9bdy3dtemp2
     611   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q10bdy3dtemp2
     612   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q11bdy3dtemp2
     613   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q12bdy3dtemp2
     614   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q13bdy3dtemp2
     615   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q14bdy3dtemp2
     616   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q15bdy3dtemp2
     617   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q16bdy3dtemp2
     618   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q17bdy3dtemp2
     619   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q18bdy3dtemp2
     620   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q19bdy3dtemp2
     621   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q20bdy3dtemp2
     622   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q21bdy3dtemp2
     623   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q22bdy3dtemp2
     624   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q23bdy3dtemp2
     625   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q24bdy3dtemp2
     626   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q25bdy3dtemp2
     627   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q26bdy3dtemp2
     628   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q27bdy3dtemp2
     629   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q28bdy3dtemp2
     630   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q29bdy3dtemp2
     631   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q30bdy3dtemp2
     632   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q31bdy3dtemp2
     633   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q32bdy3dtemp2
     634   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q33bdy3dtemp2
     635   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q34bdy3dtemp2
     636   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: q35bdy3dtemp2
    579637!!!!***MARS <<
    580638   REAL , DIMENSION(:,:,:) , ALLOCATABLE , SAVE :: mbdy2dtemp2
    … …  
    631689      IF ( ALLOCATED ( q5bdy3dtemp1 ) ) DEALLOCATE ( q5bdy3dtemp1 )
    632690      IF ( ALLOCATED ( q6bdy3dtemp1 ) ) DEALLOCATE ( q6bdy3dtemp1 )
     691      IF ( ALLOCATED ( q7bdy3dtemp1 ) ) DEALLOCATE ( q7bdy3dtemp1 )
     692      IF ( ALLOCATED ( q8bdy3dtemp1 ) ) DEALLOCATE ( q8bdy3dtemp1 )
     693      IF ( ALLOCATED ( q9bdy3dtemp1 ) ) DEALLOCATE ( q9bdy3dtemp1 )
     694      IF ( ALLOCATED ( q10bdy3dtemp1 ) ) DEALLOCATE ( q10bdy3dtemp1 )
     695      IF ( ALLOCATED ( q11bdy3dtemp1 ) ) DEALLOCATE ( q11bdy3dtemp1 )
     696      IF ( ALLOCATED ( q12bdy3dtemp1 ) ) DEALLOCATE ( q12bdy3dtemp1 )
     697      IF ( ALLOCATED ( q13bdy3dtemp1 ) ) DEALLOCATE ( q13bdy3dtemp1 )
     698      IF ( ALLOCATED ( q14bdy3dtemp1 ) ) DEALLOCATE ( q14bdy3dtemp1 )
     699      IF ( ALLOCATED ( q15bdy3dtemp1 ) ) DEALLOCATE ( q15bdy3dtemp1 )
     700      IF ( ALLOCATED ( q16bdy3dtemp1 ) ) DEALLOCATE ( q16bdy3dtemp1 )
     701      IF ( ALLOCATED ( q17bdy3dtemp1 ) ) DEALLOCATE ( q17bdy3dtemp1 )
     702      IF ( ALLOCATED ( q18bdy3dtemp1 ) ) DEALLOCATE ( q18bdy3dtemp1 )
     703      IF ( ALLOCATED ( q19bdy3dtemp1 ) ) DEALLOCATE ( q19bdy3dtemp1 )
     704      IF ( ALLOCATED ( q20bdy3dtemp1 ) ) DEALLOCATE ( q20bdy3dtemp1 )
     705      IF ( ALLOCATED ( q21bdy3dtemp1 ) ) DEALLOCATE ( q21bdy3dtemp1 )
     706      IF ( ALLOCATED ( q22bdy3dtemp1 ) ) DEALLOCATE ( q22bdy3dtemp1 )
     707      IF ( ALLOCATED ( q23bdy3dtemp1 ) ) DEALLOCATE ( q23bdy3dtemp1 )
     708      IF ( ALLOCATED ( q24bdy3dtemp1 ) ) DEALLOCATE ( q24bdy3dtemp1 )
     709      IF ( ALLOCATED ( q25bdy3dtemp1 ) ) DEALLOCATE ( q25bdy3dtemp1 )
     710      IF ( ALLOCATED ( q26bdy3dtemp1 ) ) DEALLOCATE ( q26bdy3dtemp1 )
     711      IF ( ALLOCATED ( q27bdy3dtemp1 ) ) DEALLOCATE ( q27bdy3dtemp1 )
     712      IF ( ALLOCATED ( q28bdy3dtemp1 ) ) DEALLOCATE ( q28bdy3dtemp1 )
     713      IF ( ALLOCATED ( q29bdy3dtemp1 ) ) DEALLOCATE ( q29bdy3dtemp1 )
     714      IF ( ALLOCATED ( q30bdy3dtemp1 ) ) DEALLOCATE ( q30bdy3dtemp1 )
     715      IF ( ALLOCATED ( q31bdy3dtemp1 ) ) DEALLOCATE ( q31bdy3dtemp1 )
     716      IF ( ALLOCATED ( q32bdy3dtemp1 ) ) DEALLOCATE ( q32bdy3dtemp1 )
     717      IF ( ALLOCATED ( q33bdy3dtemp1 ) ) DEALLOCATE ( q33bdy3dtemp1 )
     718      IF ( ALLOCATED ( q34bdy3dtemp1 ) ) DEALLOCATE ( q34bdy3dtemp1 )
     719      IF ( ALLOCATED ( q35bdy3dtemp1 ) ) DEALLOCATE ( q35bdy3dtemp1 )
    633720!!!!***MARS <<
    634721      IF ( ALLOCATED ( mbdy2dtemp1 ) ) DEALLOCATE ( mbdy2dtemp1 )
    … …  
    644731      IF ( ALLOCATED ( q5bdy3dtemp2 ) ) DEALLOCATE ( q5bdy3dtemp2 )
    645732      IF ( ALLOCATED ( q6bdy3dtemp2 ) ) DEALLOCATE ( q6bdy3dtemp2 )
     733      IF ( ALLOCATED ( q7bdy3dtemp2 ) ) DEALLOCATE ( q7bdy3dtemp2 )
     734      IF ( ALLOCATED ( q8bdy3dtemp2 ) ) DEALLOCATE ( q8bdy3dtemp2 )
     735      IF ( ALLOCATED ( q9bdy3dtemp2 ) ) DEALLOCATE ( q9bdy3dtemp2 )
     736      IF ( ALLOCATED ( q10bdy3dtemp2 ) ) DEALLOCATE ( q10bdy3dtemp2 )
     737      IF ( ALLOCATED ( q11bdy3dtemp2 ) ) DEALLOCATE ( q11bdy3dtemp2 )
     738      IF ( ALLOCATED ( q12bdy3dtemp2 ) ) DEALLOCATE ( q12bdy3dtemp2 )
     739      IF ( ALLOCATED ( q13bdy3dtemp2 ) ) DEALLOCATE ( q13bdy3dtemp2 )
     740      IF ( ALLOCATED ( q14bdy3dtemp2 ) ) DEALLOCATE ( q14bdy3dtemp2 )
     741      IF ( ALLOCATED ( q15bdy3dtemp2 ) ) DEALLOCATE ( q15bdy3dtemp2 )
     742      IF ( ALLOCATED ( q16bdy3dtemp2 ) ) DEALLOCATE ( q16bdy3dtemp2 )
     743      IF ( ALLOCATED ( q17bdy3dtemp2 ) ) DEALLOCATE ( q17bdy3dtemp2 )
     744      IF ( ALLOCATED ( q18bdy3dtemp2 ) ) DEALLOCATE ( q18bdy3dtemp2 )
     745      IF ( ALLOCATED ( q19bdy3dtemp2 ) ) DEALLOCATE ( q19bdy3dtemp2 )
     746      IF ( ALLOCATED ( q20bdy3dtemp2 ) ) DEALLOCATE ( q20bdy3dtemp2 )
     747      IF ( ALLOCATED ( q21bdy3dtemp2 ) ) DEALLOCATE ( q21bdy3dtemp2 )
     748      IF ( ALLOCATED ( q22bdy3dtemp2 ) ) DEALLOCATE ( q22bdy3dtemp2 )
     749      IF ( ALLOCATED ( q23bdy3dtemp2 ) ) DEALLOCATE ( q23bdy3dtemp2 )
     750      IF ( ALLOCATED ( q24bdy3dtemp2 ) ) DEALLOCATE ( q24bdy3dtemp2 )
     751      IF ( ALLOCATED ( q25bdy3dtemp2 ) ) DEALLOCATE ( q25bdy3dtemp2 )
     752      IF ( ALLOCATED ( q26bdy3dtemp2 ) ) DEALLOCATE ( q26bdy3dtemp2 )
     753      IF ( ALLOCATED ( q27bdy3dtemp2 ) ) DEALLOCATE ( q27bdy3dtemp2 )
     754      IF ( ALLOCATED ( q28bdy3dtemp2 ) ) DEALLOCATE ( q28bdy3dtemp2 )
     755      IF ( ALLOCATED ( q29bdy3dtemp2 ) ) DEALLOCATE ( q29bdy3dtemp2 )
     756      IF ( ALLOCATED ( q30bdy3dtemp2 ) ) DEALLOCATE ( q30bdy3dtemp2 )
     757      IF ( ALLOCATED ( q31bdy3dtemp2 ) ) DEALLOCATE ( q31bdy3dtemp2 )
     758      IF ( ALLOCATED ( q32bdy3dtemp2 ) ) DEALLOCATE ( q32bdy3dtemp2 )
     759      IF ( ALLOCATED ( q33bdy3dtemp2 ) ) DEALLOCATE ( q33bdy3dtemp2 )
     760      IF ( ALLOCATED ( q34bdy3dtemp2 ) ) DEALLOCATE ( q34bdy3dtemp2 )
     761      IF ( ALLOCATED ( q35bdy3dtemp2 ) ) DEALLOCATE ( q35bdy3dtemp2 )
    646762!!!!***MARS <<
    647763      IF ( ALLOCATED ( mbdy2dtemp2 ) ) DEALLOCATE ( mbdy2dtemp2 )
    … …  
    658774      ALLOCATE ( q5bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
    659775      ALLOCATE ( q6bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     776      ALLOCATE ( q7bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     777      ALLOCATE ( q8bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     778      ALLOCATE ( q9bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     779      ALLOCATE ( q10bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     780      ALLOCATE ( q11bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     781      ALLOCATE ( q12bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     782      ALLOCATE ( q13bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     783      ALLOCATE ( q14bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     784      ALLOCATE ( q15bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     785      ALLOCATE ( q16bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     786      ALLOCATE ( q17bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     787      ALLOCATE ( q18bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     788      ALLOCATE ( q19bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     789      ALLOCATE ( q20bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     790      ALLOCATE ( q21bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     791      ALLOCATE ( q22bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     792      ALLOCATE ( q23bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     793      ALLOCATE ( q24bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     794      ALLOCATE ( q25bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     795      ALLOCATE ( q26bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     796      ALLOCATE ( q27bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     797      ALLOCATE ( q28bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     798      ALLOCATE ( q29bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     799      ALLOCATE ( q30bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     800      ALLOCATE ( q31bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     801      ALLOCATE ( q32bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     802      ALLOCATE ( q33bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     803      ALLOCATE ( q34bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
     804      ALLOCATE ( q35bdy3dtemp1(ims:ime,kms:kme,jms:jme) )
    660805!!!!***MARS <<
    661806      ALLOCATE ( mbdy2dtemp1(ims:ime,1:1,    jms:jme) )
    … …  
    671816      ALLOCATE ( q5bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
    672817      ALLOCATE ( q6bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     818      ALLOCATE ( q7bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     819      ALLOCATE ( q8bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     820      ALLOCATE ( q9bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     821      ALLOCATE ( q10bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     822      ALLOCATE ( q11bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     823      ALLOCATE ( q12bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     824      ALLOCATE ( q13bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     825      ALLOCATE ( q14bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     826      ALLOCATE ( q15bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     827      ALLOCATE ( q16bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     828      ALLOCATE ( q17bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     829      ALLOCATE ( q18bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     830      ALLOCATE ( q19bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     831      ALLOCATE ( q20bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     832      ALLOCATE ( q21bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     833      ALLOCATE ( q22bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     834      ALLOCATE ( q23bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     835      ALLOCATE ( q24bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     836      ALLOCATE ( q25bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     837      ALLOCATE ( q26bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     838      ALLOCATE ( q27bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     839      ALLOCATE ( q28bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     840      ALLOCATE ( q29bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     841      ALLOCATE ( q30bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     842      ALLOCATE ( q31bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     843      ALLOCATE ( q32bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     844      ALLOCATE ( q33bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     845      ALLOCATE ( q34bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
     846      ALLOCATE ( q35bdy3dtemp2(ims:ime,kms:kme,jms:jme) )
    673847!!!!***MARS <<
    674848      ALLOCATE ( mbdy2dtemp2(ims:ime,1:1,    jms:jme) )
    … …  
    741915      CALL couple ( grid%em_mu_2 , grid%em_mub , q6bdy3dtemp1 , grid%scalar(:,:,:,7) , 't' , grid%msft , &
    742916                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )
     917ENDIF
     918
     919!VENUS
     920IF (config_flags%mars .eq. 34) THEN
     921      CALL couple ( grid%em_mu_2 , grid%em_mub , qbdy3dtemp1 ,grid%scalar(:,:,:,2) , 't' , grid%msft , &
     922                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme,ips, ipe, jps, jpe, kps, kpe )
     923      CALL couple ( grid%em_mu_2 , grid%em_mub , q2bdy3dtemp1,grid%scalar(:,:,:,3) , 't' , grid%msft , &
     924                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     925      CALL couple ( grid%em_mu_2 , grid%em_mub , q3bdy3dtemp1,grid%scalar(:,:,:,4) , 't' , grid%msft , &
     926                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     927      CALL couple ( grid%em_mu_2 , grid%em_mub , q4bdy3dtemp1,grid%scalar(:,:,:,5) , 't' , grid%msft , &
     928                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     929      CALL couple ( grid%em_mu_2 , grid%em_mub , q5bdy3dtemp1,grid%scalar(:,:,:,6) , 't' , grid%msft , &
     930                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     931      CALL couple ( grid%em_mu_2 , grid%em_mub , q6bdy3dtemp1,grid%scalar(:,:,:,7) , 't' , grid%msft , &
     932                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     933      CALL couple ( grid%em_mu_2 , grid%em_mub , q7bdy3dtemp1,grid%scalar(:,:,:,8) , 't' , grid%msft , &
     934                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     935      CALL couple ( grid%em_mu_2 , grid%em_mub , q8bdy3dtemp1,grid%scalar(:,:,:,9) , 't' , grid%msft , &
     936                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     937      CALL couple ( grid%em_mu_2 , grid%em_mub , q9bdy3dtemp1,grid%scalar(:,:,:,10) , 't' , grid%msft , &
     938                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     939      CALL couple ( grid%em_mu_2 , grid%em_mub , q10bdy3dtemp1,grid%scalar(:,:,:,11) , 't' , grid%msft , &
     940                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     941      CALL couple ( grid%em_mu_2 , grid%em_mub , q11bdy3dtemp1,grid%scalar(:,:,:,12) , 't' , grid%msft , &
     942                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     943      CALL couple ( grid%em_mu_2 , grid%em_mub , q12bdy3dtemp1,grid%scalar(:,:,:,13) , 't' , grid%msft , &
     944                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     945      CALL couple ( grid%em_mu_2 , grid%em_mub , q13bdy3dtemp1,grid%scalar(:,:,:,14) , 't' , grid%msft , &
     946                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     947      CALL couple ( grid%em_mu_2 , grid%em_mub , q14bdy3dtemp1,grid%scalar(:,:,:,15) , 't' , grid%msft , &
     948                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     949      CALL couple ( grid%em_mu_2 , grid%em_mub , q15bdy3dtemp1,grid%scalar(:,:,:,16) , 't' , grid%msft , &
     950                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     951      CALL couple ( grid%em_mu_2 , grid%em_mub , q16bdy3dtemp1,grid%scalar(:,:,:,17) , 't' , grid%msft , &
     952                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     953      CALL couple ( grid%em_mu_2 , grid%em_mub , q17bdy3dtemp1,grid%scalar(:,:,:,18) , 't' , grid%msft , &
     954                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     955      CALL couple ( grid%em_mu_2 , grid%em_mub , q18bdy3dtemp1,grid%scalar(:,:,:,19) , 't' , grid%msft , &
     956                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     957      CALL couple ( grid%em_mu_2 , grid%em_mub , q19bdy3dtemp1,grid%scalar(:,:,:,20) , 't' , grid%msft , &
     958                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     959      CALL couple ( grid%em_mu_2 , grid%em_mub , q20bdy3dtemp1,grid%scalar(:,:,:,21) , 't' , grid%msft , &
     960                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     961      CALL couple ( grid%em_mu_2 , grid%em_mub , q21bdy3dtemp1,grid%scalar(:,:,:,22) , 't' , grid%msft , &
     962                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     963      CALL couple ( grid%em_mu_2 , grid%em_mub , q22bdy3dtemp1,grid%scalar(:,:,:,23) , 't' , grid%msft , &
     964                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     965      CALL couple ( grid%em_mu_2 , grid%em_mub , q23bdy3dtemp1,grid%scalar(:,:,:,24) , 't' , grid%msft , &
     966                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     967      CALL couple ( grid%em_mu_2 , grid%em_mub , q24bdy3dtemp1,grid%scalar(:,:,:,25) , 't' , grid%msft , &
     968                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     969      CALL couple ( grid%em_mu_2 , grid%em_mub , q25bdy3dtemp1,grid%scalar(:,:,:,26) , 't' , grid%msft , &
     970                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     971      CALL couple ( grid%em_mu_2 , grid%em_mub , q26bdy3dtemp1,grid%scalar(:,:,:,27) , 't' , grid%msft , &
     972                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     973      CALL couple ( grid%em_mu_2 , grid%em_mub , q27bdy3dtemp1,grid%scalar(:,:,:,28) , 't' , grid%msft , &
     974                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     975      CALL couple ( grid%em_mu_2 , grid%em_mub , q28bdy3dtemp1,grid%scalar(:,:,:,29) , 't' , grid%msft , &
     976                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     977      CALL couple ( grid%em_mu_2 , grid%em_mub , q29bdy3dtemp1,grid%scalar(:,:,:,30) , 't' , grid%msft , &
     978                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     979      CALL couple ( grid%em_mu_2 , grid%em_mub , q30bdy3dtemp1,grid%scalar(:,:,:,31) , 't' , grid%msft , &
     980                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     981      CALL couple ( grid%em_mu_2 , grid%em_mub , q31bdy3dtemp1,grid%scalar(:,:,:,32) , 't' , grid%msft , &
     982                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     983      CALL couple ( grid%em_mu_2 , grid%em_mub , q32bdy3dtemp1,grid%scalar(:,:,:,33) , 't' , grid%msft , &
     984                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     985      CALL couple ( grid%em_mu_2 , grid%em_mub , q33bdy3dtemp1,grid%scalar(:,:,:,34) , 't' , grid%msft , &
     986                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     987      CALL couple ( grid%em_mu_2 , grid%em_mub , q34bdy3dtemp1,grid%scalar(:,:,:,35) , 't' , grid%msft , &
     988                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
    743989ENDIF
    744990!!!!!MARS!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
    … …  
    8331079                                                                 ips , ipe , jps , jpe , kps , kpe )
    8341080ENDIF
     1081
     1082!VENUS
     1083IF (config_flags%mars .eq. 34) THEN
     1084      CALL stuff_bdy     ( qbdy3dtemp1 , grid%scalar_b(:,:,:,:,2)   , 'T' , ijds, ijde , spec_bdy_width      , &
     1085                                                                 ids , ide , jds, jde , kds , kde , &
     1086                                                                 ims , ime , jms, jme , kms , kme , &
     1087                                                                 ips , ipe , jps, jpe , kps , kpe )
     1088      CALL stuff_bdy     ( q2bdy3dtemp1 , grid%scalar_b(:,:,:,:,3)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1089                                                                 ids , ide , jds, jde , kds , kde , &
     1090                                                                 ims , ime , jms, jme , kms , kme , &
     1091                                                                 ips , ipe , jps, jpe , kps , kpe )
     1092      CALL stuff_bdy     ( q3bdy3dtemp1 , grid%scalar_b(:,:,:,:,4)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1093                                                                 ids , ide , jds, jde , kds , kde , &
     1094                                                                 ims , ime , jms, jme , kms , kme , &
     1095                                                                 ips , ipe , jps, jpe , kps , kpe )
     1096      CALL stuff_bdy     ( q4bdy3dtemp1 , grid%scalar_b(:,:,:,:,5)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1097                                                                 ids , ide , jds, jde , kds , kde , &
     1098                                                                 ims , ime , jms, jme , kms , kme , &
     1099                                                                 ips , ipe , jps, jpe , kps , kpe )
     1100      CALL stuff_bdy     ( q5bdy3dtemp1 , grid%scalar_b(:,:,:,:,6)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1101                                                                 ids , ide , jds, jde , kds , kde , &
     1102                                                                 ims , ime , jms, jme , kms , kme , &
     1103                                                                 ips , ipe , jps, jpe , kps , kpe )
     1104      CALL stuff_bdy     ( q6bdy3dtemp1 , grid%scalar_b(:,:,:,:,7)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1105                                                                 ids , ide , jds, jde , kds , kde , &
     1106                                                                 ims , ime , jms, jme , kms , kme , &
     1107                                                                 ips , ipe , jps, jpe , kps , kpe )
     1108      CALL stuff_bdy     ( q7bdy3dtemp1 , grid%scalar_b(:,:,:,:,8)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1109                                                                 ids , ide , jds, jde , kds , kde , &
     1110                                                                 ims , ime , jms, jme , kms , kme , &
     1111                                                                 ips , ipe , jps, jpe , kps , kpe )
     1112      CALL stuff_bdy     ( q8bdy3dtemp1 , grid%scalar_b(:,:,:,:,9)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1113                                                                 ids , ide , jds, jde , kds , kde , &
     1114                                                                 ims , ime , jms, jme , kms , kme , &
     1115                                                                 ips , ipe , jps, jpe , kps , kpe )
     1116      CALL stuff_bdy     ( q9bdy3dtemp1 , grid%scalar_b(:,:,:,:,10)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1117                                                                 ids , ide , jds, jde , kds , kde , &
     1118                                                                 ims , ime , jms, jme , kms , kme , &
     1119                                                                 ips , ipe , jps, jpe , kps , kpe )
     1120      CALL stuff_bdy     ( q10bdy3dtemp1 , grid%scalar_b(:,:,:,:,11)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1121                                                                 ids , ide , jds, jde , kds , kde , &
     1122                                                                 ims , ime , jms, jme , kms , kme , &
     1123                                                                 ips , ipe , jps, jpe , kps , kpe )
     1124      CALL stuff_bdy     ( q11bdy3dtemp1 , grid%scalar_b(:,:,:,:,12)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1125                                                                 ids , ide , jds, jde , kds , kde , &
     1126                                                                 ims , ime , jms, jme , kms , kme , &
     1127                                                                 ips , ipe , jps, jpe , kps , kpe )
     1128      CALL stuff_bdy     ( q12bdy3dtemp1 , grid%scalar_b(:,:,:,:,13)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1129                                                                 ids , ide , jds, jde , kds , kde , &
     1130                                                                 ims , ime , jms, jme , kms , kme , &
     1131                                                                 ips , ipe , jps, jpe , kps , kpe )
     1132      CALL stuff_bdy     ( q13bdy3dtemp1 , grid%scalar_b(:,:,:,:,14)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1133                                                                 ids , ide , jds, jde , kds , kde , &
     1134                                                                 ims , ime , jms, jme , kms , kme , &
     1135                                                                 ips , ipe , jps, jpe , kps , kpe )
     1136      CALL stuff_bdy     ( q14bdy3dtemp1 , grid%scalar_b(:,:,:,:,15)   , 'T' ,ijds , ijde , spec_bdy_width      , &                                                                 
     1137                                                                 ids , ide , jds, jde , kds , kde , &
     1138                                                                 ims , ime , jms, jme , kms , kme , &
     1139                                                                 ips , ipe , jps, jpe , kps , kpe )
     1140      CALL stuff_bdy     ( q15bdy3dtemp1 , grid%scalar_b(:,:,:,:,16)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1141                                                                 ids , ide , jds, jde , kds , kde , &
     1142                                                                 ims , ime , jms, jme , kms , kme , &
     1143                                                                 ips , ipe , jps, jpe , kps , kpe )
     1144      CALL stuff_bdy     ( q16bdy3dtemp1 , grid%scalar_b(:,:,:,:,17)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1145                                                                 ids , ide , jds, jde , kds , kde , &
     1146                                                                 ims , ime , jms, jme , kms , kme , &
     1147                                                                 ips , ipe , jps, jpe , kps , kpe )
     1148      CALL stuff_bdy     ( q17bdy3dtemp1 , grid%scalar_b(:,:,:,:,18)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1149                                                                 ids , ide , jds, jde , kds , kde , &
     1150                                                                 ims , ime , jms, jme , kms , kme , &
     1151                                                                 ips , ipe , jps, jpe , kps , kpe )
     1152      CALL stuff_bdy     ( q18bdy3dtemp1 , grid%scalar_b(:,:,:,:,19)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1153                                                                 ids , ide , jds, jde , kds , kde , &
     1154                                                                 ims , ime , jms, jme , kms , kme , &
     1155                                                                 ips , ipe , jps, jpe , kps , kpe )
     1156      CALL stuff_bdy     ( q19bdy3dtemp1 , grid%scalar_b(:,:,:,:,20)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1157                                                                 ids , ide , jds, jde , kds , kde , &
     1158                                                                 ims , ime , jms, jme , kms , kme , &
     1159                                                                 ips , ipe , jps, jpe , kps , kpe )
     1160      CALL stuff_bdy     ( q20bdy3dtemp1 , grid%scalar_b(:,:,:,:,21)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1161                                                                 ids , ide , jds, jde , kds , kde , &
     1162                                                                 ims , ime , jms, jme , kms , kme , &
     1163                                                                 ips , ipe , jps, jpe , kps , kpe )
     1164      CALL stuff_bdy     ( q21bdy3dtemp1 , grid%scalar_b(:,:,:,:,22)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1165                                                                 ids , ide , jds, jde , kds , kde , &
     1166                                                                 ims , ime , jms, jme , kms , kme , &
     1167                                                                 ips , ipe , jps, jpe , kps , kpe )
     1168      CALL stuff_bdy     ( q22bdy3dtemp1 , grid%scalar_b(:,:,:,:,23)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1169                                                                 ids , ide , jds, jde , kds , kde , &
     1170                                                                 ims , ime , jms, jme , kms , kme , &
     1171                                                                 ips , ipe , jps, jpe , kps , kpe )
     1172      CALL stuff_bdy     ( q23bdy3dtemp1 , grid%scalar_b(:,:,:,:,24)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1173                                                                 ids , ide , jds, jde , kds , kde , &
     1174                                                                 ims , ime , jms, jme , kms , kme , &
     1175                                                                 ips , ipe , jps, jpe , kps , kpe )
     1176      CALL stuff_bdy     ( q24bdy3dtemp1 , grid%scalar_b(:,:,:,:,25)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1177                                                                 ids , ide , jds, jde , kds , kde , &
     1178                                                                 ims , ime , jms, jme , kms , kme , &
     1179                                                                 ips , ipe , jps, jpe , kps , kpe )
     1180      CALL stuff_bdy     ( q25bdy3dtemp1 , grid%scalar_b(:,:,:,:,26)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1181                                                                 ids , ide , jds, jde , kds , kde , &
     1182                                                                 ims , ime , jms, jme , kms , kme , &
     1183                                                                 ips , ipe , jps, jpe , kps , kpe )
     1184      CALL stuff_bdy     ( q26bdy3dtemp1 , grid%scalar_b(:,:,:,:,27)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1185                                                                 ids , ide , jds, jde , kds , kde , &
     1186                                                                 ims , ime , jms, jme , kms , kme , &
     1187                                                                 ips , ipe , jps, jpe , kps , kpe )
     1188      CALL stuff_bdy     ( q27bdy3dtemp1 , grid%scalar_b(:,:,:,:,28)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1189                                                                 ids , ide , jds, jde , kds , kde , &
     1190                                                                 ims , ime , jms, jme , kms , kme , &
     1191                                                                 ips , ipe , jps, jpe , kps , kpe )
     1192      CALL stuff_bdy     ( q28bdy3dtemp1 , grid%scalar_b(:,:,:,:,29)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1193                                                                 ids , ide , jds, jde , kds , kde , &
     1194                                                                 ims , ime , jms, jme , kms , kme , &
     1195                                                                 ips , ipe , jps, jpe , kps , kpe )
     1196      CALL stuff_bdy     ( q29bdy3dtemp1 , grid%scalar_b(:,:,:,:,30)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1197                                                                 ids , ide , jds, jde , kds , kde , &
     1198                                                                 ims , ime , jms, jme , kms , kme , &
     1199                                                                 ips , ipe , jps, jpe , kps , kpe )
     1200      CALL stuff_bdy     ( q30bdy3dtemp1 , grid%scalar_b(:,:,:,:,31)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1201                                                                 ids , ide , jds, jde , kds , kde , &
     1202                                                                 ims , ime , jms, jme , kms , kme , &
     1203                                                                 ips , ipe , jps, jpe , kps , kpe )
     1204      CALL stuff_bdy     ( q31bdy3dtemp1 , grid%scalar_b(:,:,:,:,33)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1205                                                                 ids , ide , jds, jde , kds , kde , &
     1206                                                                 ims , ime , jms, jme , kms , kme , &
     1207                                                                 ips , ipe , jps, jpe , kps , kpe )
     1208      CALL stuff_bdy     ( q32bdy3dtemp1 , grid%scalar_b(:,:,:,:,33)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1209                                                                 ids , ide , jds, jde , kds , kde , &
     1210                                                                 ims , ime , jms, jme , kms , kme , &
     1211                                                                 ips , ipe , jps, jpe , kps , kpe )
     1212      CALL stuff_bdy     ( q33bdy3dtemp1 , grid%scalar_b(:,:,:,:,34)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1213                                                                 ids , ide , jds, jde , kds , kde , &
     1214                                                                 ims , ime , jms, jme , kms , kme , &
     1215                                                                 ips , ipe , jps, jpe , kps , kpe )
     1216      CALL stuff_bdy     ( q34bdy3dtemp1 , grid%scalar_b(:,:,:,:,35)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1217                                                                 ids , ide , jds, jde , kds , kde , &
     1218                                                                 ims , ime , jms, jme , kms , kme , &
     1219                                                                 ips , ipe , jps, jpe , kps , kpe )
     1220ENDIF
    8351221!!!!MARS <<
    8361222
    … …  
    9111297                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme, ips, ipe, jps, jpe, kps, kpe )
    9121298ENDIF
     1299
     1300!VENUS
     1301IF (config_flags%mars .eq. 34) THEN
     1302      CALL couple ( grid%em_mu_2 , grid%em_mub , qbdy3dtemp2 ,grid%scalar(:,:,:,2) , 't' , grid%msft , &
     1303                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms, kme,ips, ipe, jps, jpe, kps, kpe )
     1304      CALL couple ( grid%em_mu_2 , grid%em_mub , q2bdy3dtemp2,grid%scalar(:,:,:,3) , 't' , grid%msft , &
     1305                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1306      CALL couple ( grid%em_mu_2 , grid%em_mub , q3bdy3dtemp2,grid%scalar(:,:,:,4) , 't' , grid%msft , &
     1307                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1308      CALL couple ( grid%em_mu_2 , grid%em_mub , q4bdy3dtemp2,grid%scalar(:,:,:,5) , 't' , grid%msft , &
     1309                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1310      CALL couple ( grid%em_mu_2 , grid%em_mub , q5bdy3dtemp2,grid%scalar(:,:,:,6) , 't' , grid%msft , &
     1311                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1312      CALL couple ( grid%em_mu_2 , grid%em_mub , q6bdy3dtemp2,grid%scalar(:,:,:,7) , 't' , grid%msft , &
     1313                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1314      CALL couple ( grid%em_mu_2 , grid%em_mub , q7bdy3dtemp2,grid%scalar(:,:,:,8) , 't' , grid%msft , &
     1315                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1316      CALL couple ( grid%em_mu_2 , grid%em_mub , q8bdy3dtemp2,grid%scalar(:,:,:,9) , 't' , grid%msft , &
     1317                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1318      CALL couple ( grid%em_mu_2 , grid%em_mub , q9bdy3dtemp2,grid%scalar(:,:,:,10) , 't' , grid%msft , &
     1319                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1320      CALL couple ( grid%em_mu_2 , grid%em_mub , q10bdy3dtemp2,grid%scalar(:,:,:,11) , 't' , grid%msft , &
     1321                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1322      CALL couple ( grid%em_mu_2 , grid%em_mub , q11bdy3dtemp2,grid%scalar(:,:,:,12) , 't' , grid%msft , &
     1323                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1324      CALL couple ( grid%em_mu_2 , grid%em_mub , q12bdy3dtemp2,grid%scalar(:,:,:,13) , 't' , grid%msft , &
     1325                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1326      CALL couple ( grid%em_mu_2 , grid%em_mub , q13bdy3dtemp2,grid%scalar(:,:,:,14) , 't' , grid%msft , &
     1327                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1328      CALL couple ( grid%em_mu_2 , grid%em_mub , q14bdy3dtemp2,grid%scalar(:,:,:,15) , 't' , grid%msft , &
     1329                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1330      CALL couple ( grid%em_mu_2 , grid%em_mub , q15bdy3dtemp2,grid%scalar(:,:,:,16) , 't' , grid%msft , &
     1331                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1332      CALL couple ( grid%em_mu_2 , grid%em_mub , q16bdy3dtemp2,grid%scalar(:,:,:,17) , 't' , grid%msft , &
     1333                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1334      CALL couple ( grid%em_mu_2 , grid%em_mub , q17bdy3dtemp2,grid%scalar(:,:,:,18) , 't' , grid%msft , &
     1335                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1336      CALL couple ( grid%em_mu_2 , grid%em_mub , q18bdy3dtemp2,grid%scalar(:,:,:,19) , 't' , grid%msft , &
     1337                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1338      CALL couple ( grid%em_mu_2 , grid%em_mub , q19bdy3dtemp2,grid%scalar(:,:,:,20) , 't' , grid%msft , &
     1339                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1340      CALL couple ( grid%em_mu_2 , grid%em_mub , q20bdy3dtemp2,grid%scalar(:,:,:,21) , 't' , grid%msft , &
     1341                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1342      CALL couple ( grid%em_mu_2 , grid%em_mub , q21bdy3dtemp2,grid%scalar(:,:,:,22) , 't' , grid%msft , &
     1343                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1344      CALL couple ( grid%em_mu_2 , grid%em_mub , q22bdy3dtemp2,grid%scalar(:,:,:,23) , 't' , grid%msft , &
     1345                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1346      CALL couple ( grid%em_mu_2 , grid%em_mub , q23bdy3dtemp2,grid%scalar(:,:,:,24) , 't' , grid%msft , &
     1347                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1348      CALL couple ( grid%em_mu_2 , grid%em_mub , q24bdy3dtemp2,grid%scalar(:,:,:,25) , 't' , grid%msft , &
     1349                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1350      CALL couple ( grid%em_mu_2 , grid%em_mub , q25bdy3dtemp2,grid%scalar(:,:,:,26) , 't' , grid%msft , &
     1351                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1352      CALL couple ( grid%em_mu_2 , grid%em_mub , q26bdy3dtemp2,grid%scalar(:,:,:,27) , 't' , grid%msft , &
     1353                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1354      CALL couple ( grid%em_mu_2 , grid%em_mub , q27bdy3dtemp2,grid%scalar(:,:,:,28) , 't' , grid%msft , &
     1355                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1356      CALL couple ( grid%em_mu_2 , grid%em_mub , q28bdy3dtemp2,grid%scalar(:,:,:,29) , 't' , grid%msft , &
     1357                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1358      CALL couple ( grid%em_mu_2 , grid%em_mub , q29bdy3dtemp2,grid%scalar(:,:,:,30) , 't' , grid%msft , &
     1359                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1360      CALL couple ( grid%em_mu_2 , grid%em_mub , q30bdy3dtemp2,grid%scalar(:,:,:,31) , 't' , grid%msft , &
     1361                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1362      CALL couple ( grid%em_mu_2 , grid%em_mub , q31bdy3dtemp2,grid%scalar(:,:,:,32) , 't' , grid%msft , &
     1363                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1364      CALL couple ( grid%em_mu_2 , grid%em_mub , q32bdy3dtemp2,grid%scalar(:,:,:,33) , 't' , grid%msft , &
     1365                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1366      CALL couple ( grid%em_mu_2 , grid%em_mub , q33bdy3dtemp2,grid%scalar(:,:,:,34) , 't' , grid%msft , &
     1367                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1368      CALL couple ( grid%em_mu_2 , grid%em_mub , q34bdy3dtemp2,grid%scalar(:,:,:,35) , 't' , grid%msft , &
     1369                    ids, ide, jds, jde, kds, kde, ims, ime, jms, jme, kms,kme,ips, ipe, jps, jpe, kps, kpe )
     1370ENDIF
    9131371!!!!!MARS <<
    9141372
    … …  
    10101468                                                            ims , ime , jms , jme , kms , kme , &
    10111469                                                            ips , ipe , jps , jpe , kps , kpe )
     1470ENDIF
     1471
     1472!VENUS
     1473IF (config_flags%mars .eq. 34) THEN
     1474      CALL stuff_bdytend ( qbdy3dtemp2 , qbdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,2) , 'T' , &
     1475                                                            ijds , ijde , spec_bdy_width      , &
     1476                                                            ids , ide , jds , jde , kds , kde , &
     1477                                                            ims , ime , jms , jme , kms , kme , &
     1478                                                            ips , ipe , jps , jpe , kps , kpe )
     1479      CALL stuff_bdytend ( q2bdy3dtemp2 , q2bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,3) , 'T' , &
     1480                                                            ijds , ijde , spec_bdy_width      , &
     1481                                                            ids , ide , jds , jde , kds , kde , &
     1482                                                            ims , ime , jms , jme , kms , kme , &
     1483                                                            ips , ipe , jps , jpe , kps , kpe )
     1484      CALL stuff_bdytend ( q3bdy3dtemp2 , q3bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,4) , 'T' , &
     1485                                                            ijds , ijde , spec_bdy_width      , &
     1486                                                            ids , ide , jds , jde , kds , kde , &
     1487                                                            ims , ime , jms , jme , kms , kme , &
     1488                                                            ips , ipe , jps , jpe , kps , kpe )
     1489      CALL stuff_bdytend ( q4bdy3dtemp2 , q4bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,5) , 'T' , &
     1490                                                            ijds , ijde , spec_bdy_width      , &
     1491                                                            ids , ide , jds , jde , kds , kde , &
     1492                                                            ims , ime , jms , jme , kms , kme , &
     1493                                                            ips , ipe , jps , jpe , kps , kpe )
     1494      CALL stuff_bdytend ( q5bdy3dtemp2 , q5bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,6) , 'T' , &
     1495                                                            ijds , ijde , spec_bdy_width      , &
     1496                                                            ids , ide , jds , jde , kds , kde , &
     1497                                                            ims , ime , jms , jme , kms , kme , &
     1498                                                            ips , ipe , jps , jpe , kps , kpe )
     1499      CALL stuff_bdytend ( q6bdy3dtemp2 , q6bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,7) , 'T' , &
     1500                                                            ijds , ijde , spec_bdy_width      , &
     1501                                                            ids , ide , jds , jde , kds , kde , &
     1502                                                            ims , ime , jms , jme , kms , kme , &
     1503                                                            ips , ipe , jps , jpe , kps , kpe )
     1504      CALL stuff_bdytend ( q7bdy3dtemp2 , q7bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,8) , 'T' , &
     1505                                                            ijds , ijde , spec_bdy_width      , &
     1506                                                            ids , ide , jds , jde , kds , kde , &
     1507                                                            ims , ime , jms , jme , kms , kme , &
     1508                                                            ips , ipe , jps , jpe , kps , kpe )
     1509      CALL stuff_bdytend ( q8bdy3dtemp2 , q8bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,9) , 'T' , &
     1510                                                            ijds , ijde , spec_bdy_width      , &
     1511                                                            ids , ide , jds , jde , kds , kde , &
     1512                                                            ims , ime , jms , jme , kms , kme , &
     1513                                                            ips , ipe , jps , jpe , kps , kpe )
     1514      CALL stuff_bdytend ( q9bdy3dtemp2 , q9bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,10) , 'T' , &
     1515                                                            ijds , ijde , spec_bdy_width      , &
     1516                                                            ids , ide , jds , jde , kds , kde , &
     1517                                                            ims , ime , jms , jme , kms , kme , &
     1518                                                            ips , ipe , jps , jpe , kps , kpe )
     1519      CALL stuff_bdytend ( q10bdy3dtemp2 , q10bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,11) , 'T' , &
     1520                                                            ijds , ijde , spec_bdy_width      , &
     1521                                                            ids , ide , jds , jde , kds , kde , &
     1522                                                            ims , ime , jms , jme , kms , kme , &
     1523                                                            ips , ipe , jps , jpe , kps , kpe )
     1524      CALL stuff_bdytend ( q11bdy3dtemp2 , q11bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,12) , 'T' , &
     1525                                                            ijds , ijde , spec_bdy_width      , &
     1526                                                            ids , ide , jds , jde , kds , kde , &
     1527                                                            ims , ime , jms , jme , kms , kme , &
     1528                                                            ips , ipe , jps , jpe , kps , kpe )
     1529      CALL stuff_bdytend ( q12bdy3dtemp2 , q12bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,13) , 'T' , &
     1530                                                            ijds , ijde , spec_bdy_width      , &
     1531                                                            ids , ide , jds , jde , kds , kde , &
     1532                                                            ims , ime , jms , jme , kms , kme , &
     1533                                                            ips , ipe , jps , jpe , kps , kpe )
     1534      CALL stuff_bdytend ( q13bdy3dtemp2 , q13bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,14) , 'T' , &
     1535                                                            ijds , ijde , spec_bdy_width      , &
     1536                                                            ids , ide , jds , jde , kds , kde , &
     1537                                                            ims , ime , jms , jme , kms , kme , &
     1538                                                            ips , ipe , jps , jpe , kps , kpe )
     1539      CALL stuff_bdytend ( q14bdy3dtemp2 , q14bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,15) , 'T' , &
     1540                                                            ijds , ijde , spec_bdy_width      , &
     1541                                                            ids , ide , jds , jde , kds , kde , &
     1542                                                            ims , ime , jms , jme , kms , kme , &
     1543                                                            ips , ipe , jps , jpe , kps , kpe )
     1544      CALL stuff_bdytend ( q15bdy3dtemp2 , q15bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,16) , 'T' , &
     1545                                                            ijds , ijde , spec_bdy_width      , &
     1546                                                            ids , ide , jds , jde , kds , kde , &
     1547                                                            ims , ime , jms , jme , kms , kme , &
     1548                                                            ips , ipe , jps , jpe , kps , kpe )
     1549      CALL stuff_bdytend ( q16bdy3dtemp2 , q16bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,17) , 'T' , &
     1550                                                            ijds , ijde , spec_bdy_width      , &
     1551                                                            ids , ide , jds , jde , kds , kde , &
     1552                                                            ims , ime , jms , jme , kms , kme , &
     1553                                                            ips , ipe , jps , jpe , kps , kpe )
     1554      CALL stuff_bdytend ( q17bdy3dtemp2 , q17bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,18) , 'T' , &
     1555                                                            ijds , ijde , spec_bdy_width      , &
     1556                                                            ids , ide , jds , jde , kds , kde , &
     1557                                                            ims , ime , jms , jme , kms , kme , &
     1558                                                            ips , ipe , jps , jpe , kps , kpe )
     1559      CALL stuff_bdytend ( q18bdy3dtemp2 , q18bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,19) , 'T' , &
     1560                                                            ijds , ijde , spec_bdy_width      , &
     1561                                                            ids , ide , jds , jde , kds , kde , &
     1562                                                            ims , ime , jms , jme , kms , kme , &
     1563                                                            ips , ipe , jps , jpe , kps , kpe )
     1564      CALL stuff_bdytend ( q19bdy3dtemp2 , q19bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,20) , 'T' , &
     1565                                                            ijds , ijde , spec_bdy_width      , &
     1566                                                            ids , ide , jds , jde , kds , kde , &
     1567                                                            ims , ime , jms , jme , kms , kme , &
     1568                                                            ips , ipe , jps , jpe , kps , kpe )
     1569      CALL stuff_bdytend ( q20bdy3dtemp2 , q20bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,21) , 'T' , &
     1570                                                            ijds , ijde , spec_bdy_width      , &
     1571                                                            ids , ide , jds , jde , kds , kde , &
     1572                                                            ims , ime , jms , jme , kms , kme , &
     1573                                                            ips , ipe , jps , jpe , kps , kpe )
     1574      CALL stuff_bdytend ( q21bdy3dtemp2 , q21bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,22) , 'T' , &
     1575                                                            ijds , ijde , spec_bdy_width      , &
     1576                                                            ids , ide , jds , jde , kds , kde , &
     1577                                                            ims , ime , jms , jme , kms , kme , &
     1578                                                            ips , ipe , jps , jpe , kps , kpe )
     1579      CALL stuff_bdytend ( q22bdy3dtemp2 , q22bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,23) , 'T' , &
     1580                                                            ijds , ijde , spec_bdy_width      , &
     1581                                                            ids , ide , jds , jde , kds , kde , &
     1582                                                            ims , ime , jms , jme , kms , kme , &
     1583                                                            ips , ipe , jps , jpe , kps , kpe )
     1584      CALL stuff_bdytend ( q23bdy3dtemp2 , q23bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,24) , 'T' , &
     1585                                                            ijds , ijde , spec_bdy_width      , &
     1586                                                            ids , ide , jds , jde , kds , kde , &
     1587                                                            ims , ime , jms , jme , kms , kme , &
     1588                                                            ips , ipe , jps , jpe , kps , kpe )
     1589      CALL stuff_bdytend ( q24bdy3dtemp2 , q24bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,25) , 'T' , &
     1590                                                            ijds , ijde , spec_bdy_width      , &
     1591                                                            ids , ide , jds , jde , kds , kde , &
     1592                                                            ims , ime , jms , jme , kms , kme , &
     1593                                                            ips , ipe , jps , jpe , kps , kpe )
     1594      CALL stuff_bdytend ( q25bdy3dtemp2 , q25bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,26) , 'T' , &
     1595                                                            ijds , ijde , spec_bdy_width      , &
     1596                                                            ids , ide , jds , jde , kds , kde , &
     1597                                                            ims , ime , jms , jme , kms , kme , &
     1598                                                            ips , ipe , jps , jpe , kps , kpe )
     1599      CALL stuff_bdytend ( q26bdy3dtemp2 , q26bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,27) , 'T' , &
     1600                                                            ijds , ijde , spec_bdy_width      , &
     1601                                                            ids , ide , jds , jde , kds , kde , &
     1602                                                            ims , ime , jms , jme , kms , kme , &
     1603                                                            ips , ipe , jps , jpe , kps , kpe )
     1604      CALL stuff_bdytend ( q27bdy3dtemp2 , q27bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,28) , 'T' , &
     1605                                                            ijds , ijde , spec_bdy_width      , &
     1606                                                            ids , ide , jds , jde , kds , kde , &
     1607                                                            ims , ime , jms , jme , kms , kme , &
     1608                                                            ips , ipe , jps , jpe , kps , kpe )
     1609      CALL stuff_bdytend ( q28bdy3dtemp2 , q28bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,29) , 'T' , &
     1610                                                            ijds , ijde , spec_bdy_width      , &
     1611                                                            ids , ide , jds , jde , kds , kde , &
     1612                                                            ims , ime , jms , jme , kms , kme , &
     1613                                                            ips , ipe , jps , jpe , kps , kpe )
     1614      CALL stuff_bdytend ( q29bdy3dtemp2 , q29bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,30) , 'T' , &
     1615                                                            ijds , ijde , spec_bdy_width      , &
     1616                                                            ids , ide , jds , jde , kds , kde , &
     1617                                                            ims , ime , jms , jme , kms , kme , &
     1618                                                            ips , ipe , jps , jpe , kps , kpe )
     1619      CALL stuff_bdytend ( q30bdy3dtemp2 , q30bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,31) , 'T' , &
     1620                                                            ijds , ijde , spec_bdy_width      , &
     1621                                                            ids , ide , jds , jde , kds , kde , &
     1622                                                            ims , ime , jms , jme , kms , kme , &
     1623                                                            ips , ipe , jps , jpe , kps , kpe )
     1624      CALL stuff_bdytend ( q31bdy3dtemp2 , q31bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,32) , 'T' , &
     1625                                                            ijds , ijde , spec_bdy_width      , &
     1626                                                            ids , ide , jds , jde , kds , kde , &
     1627                                                            ims , ime , jms , jme , kms , kme , &
     1628                                                            ips , ipe , jps , jpe , kps , kpe )
     1629      CALL stuff_bdytend ( q32bdy3dtemp2 , q32bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,33) , 'T' , &
     1630                                                            ijds , ijde , spec_bdy_width      , &
     1631                                                            ids , ide , jds , jde , kds , kde , &
     1632                                                            ims , ime , jms , jme , kms , kme , &
     1633                                                            ips , ipe , jps , jpe , kps , kpe )
     1634      CALL stuff_bdytend ( q33bdy3dtemp2 , q33bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,34) , 'T' , &
     1635                                                            ijds , ijde , spec_bdy_width      , &
     1636                                                            ids , ide , jds , jde , kds , kde , &
     1637                                                            ims , ime , jms , jme , kms , kme , &
     1638                                                            ips , ipe , jps , jpe , kps , kpe )
     1639      CALL stuff_bdytend ( q34bdy3dtemp2 , q34bdy3dtemp1 , REAL(interval_seconds) , grid%scalar_bt(:,:,:,:,35) , 'T' , &
     1640                                                            ijds , ijde , spec_bdy_width      , &
     1641                                                            ids , ide , jds , jde , kds , kde , &
     1642                                                            ims , ime , jms , jme , kms , kme , &
     1643                                                            ips , ipe , jps , jpe , kps , kpe )
     1644
    10121645ENDIF
    10131646!!!!!MARS <<
    … …  
    11011734                  q5bdy3dtemp1(i,k,j) = q5bdy3dtemp2(i,k,j)
    11021735                  q6bdy3dtemp1(i,k,j) = q6bdy3dtemp2(i,k,j)
     1736                  q7bdy3dtemp1(i,k,j) = q7bdy3dtemp2(i,k,j)
     1737                  q8bdy3dtemp1(i,k,j) = q8bdy3dtemp2(i,k,j)
     1738                  q9bdy3dtemp1(i,k,j) = q9bdy3dtemp2(i,k,j)
     1739                  q10bdy3dtemp1(i,k,j) = q10bdy3dtemp2(i,k,j)
     1740                  q11bdy3dtemp1(i,k,j) = q11bdy3dtemp2(i,k,j)
     1741                  q12bdy3dtemp1(i,k,j) = q12bdy3dtemp2(i,k,j)
     1742                  q13bdy3dtemp1(i,k,j) = q13bdy3dtemp2(i,k,j)
     1743                  q14bdy3dtemp1(i,k,j) = q14bdy3dtemp2(i,k,j)
     1744                  q15bdy3dtemp1(i,k,j) = q15bdy3dtemp2(i,k,j)
     1745                  q16bdy3dtemp1(i,k,j) = q16bdy3dtemp2(i,k,j)
     1746                  q17bdy3dtemp1(i,k,j) = q17bdy3dtemp2(i,k,j)
     1747                  q18bdy3dtemp1(i,k,j) = q18bdy3dtemp2(i,k,j)
     1748                  q19bdy3dtemp1(i,k,j) = q19bdy3dtemp2(i,k,j)
     1749                  q20bdy3dtemp1(i,k,j) = q20bdy3dtemp2(i,k,j)
     1750                  q21bdy3dtemp1(i,k,j) = q21bdy3dtemp2(i,k,j)
     1751                  q22bdy3dtemp1(i,k,j) = q22bdy3dtemp2(i,k,j)
     1752                  q23bdy3dtemp1(i,k,j) = q23bdy3dtemp2(i,k,j)
     1753                  q24bdy3dtemp1(i,k,j) = q24bdy3dtemp2(i,k,j)
     1754                  q25bdy3dtemp1(i,k,j) = q25bdy3dtemp2(i,k,j)
     1755                  q26bdy3dtemp1(i,k,j) = q26bdy3dtemp2(i,k,j)
     1756                  q27bdy3dtemp1(i,k,j) = q27bdy3dtemp2(i,k,j)
     1757                  q28bdy3dtemp1(i,k,j) = q28bdy3dtemp2(i,k,j)
     1758                  q29bdy3dtemp1(i,k,j) = q29bdy3dtemp2(i,k,j)
     1759                  q30bdy3dtemp1(i,k,j) = q30bdy3dtemp2(i,k,j)
     1760                  q31bdy3dtemp1(i,k,j) = q31bdy3dtemp2(i,k,j)
     1761                  q32bdy3dtemp1(i,k,j) = q32bdy3dtemp2(i,k,j)
     1762                  q33bdy3dtemp1(i,k,j) = q33bdy3dtemp2(i,k,j)
     1763                  q34bdy3dtemp1(i,k,j) = q34bdy3dtemp2(i,k,j)                 
     1764                  q35bdy3dtemp1(i,k,j) = q35bdy3dtemp2(i,k,j)
    11031765!!!!MARS <<
    11041766               END DO
    … …  
    11921854                                                                    ips , ipe , jps , jpe , kps , kpe )
    11931855ENDIF
     1856
     1857!VENUS
     1858IF (config_flags%mars .eq. 34) THEN
     1859      CALL stuff_bdy     ( qbdy3dtemp1 , grid%scalar_b(:,:,:,:,2)   , 'T' , ijds, ijde , spec_bdy_width      , &
     1860                                                                 ids , ide , jds, jde , kds , kde , &
     1861                                                                 ims , ime , jms, jme , kms , kme , &
     1862                                                                 ips , ipe , jps, jpe , kps , kpe )
     1863      CALL stuff_bdy     ( q2bdy3dtemp1 , grid%scalar_b(:,:,:,:,3)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1864                                                                 ids , ide , jds, jde , kds , kde , &
     1865                                                                 ims , ime , jms, jme , kms , kme , &
     1866                                                                 ips , ipe , jps, jpe , kps , kpe )
     1867      CALL stuff_bdy     ( q3bdy3dtemp1 , grid%scalar_b(:,:,:,:,4)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1868                                                                 ids , ide , jds, jde , kds , kde , &
     1869                                                                 ims , ime , jms, jme , kms , kme , &
     1870                                                                 ips , ipe , jps, jpe , kps , kpe )
     1871      CALL stuff_bdy     ( q4bdy3dtemp1 , grid%scalar_b(:,:,:,:,5)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1872                                                                 ids , ide , jds, jde , kds , kde , &
     1873                                                                 ims , ime , jms, jme , kms , kme , &
     1874                                                                 ips , ipe , jps, jpe , kps , kpe )
     1875      CALL stuff_bdy     ( q5bdy3dtemp1 , grid%scalar_b(:,:,:,:,6)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1876                                                                 ids , ide , jds, jde , kds , kde , &
     1877                                                                 ims , ime , jms, jme , kms , kme , &
     1878                                                                 ips , ipe , jps, jpe , kps , kpe )
     1879      CALL stuff_bdy     ( q6bdy3dtemp1 , grid%scalar_b(:,:,:,:,7)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1880                                                                 ids , ide , jds, jde , kds , kde , &
     1881                                                                 ims , ime , jms, jme , kms , kme , &
     1882                                                                 ips , ipe , jps, jpe , kps , kpe )
     1883      CALL stuff_bdy     ( q7bdy3dtemp1 , grid%scalar_b(:,:,:,:,8)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1884                                                                 ids , ide , jds, jde , kds , kde , &
     1885                                                                 ims , ime , jms, jme , kms , kme , &
     1886                                                                 ips , ipe , jps, jpe , kps , kpe )
     1887      CALL stuff_bdy     ( q8bdy3dtemp1 , grid%scalar_b(:,:,:,:,9)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1888                                                                 ids , ide , jds, jde , kds , kde , &
     1889                                                                 ims , ime , jms, jme , kms , kme , &
     1890                                                                 ips , ipe , jps, jpe , kps , kpe )
     1891      CALL stuff_bdy     ( q9bdy3dtemp1 , grid%scalar_b(:,:,:,:,10)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1892                                                                 ids , ide , jds, jde , kds , kde , &
     1893                                                                 ims , ime , jms, jme , kms , kme , &
     1894                                                                 ips , ipe , jps, jpe , kps , kpe )
     1895      CALL stuff_bdy     ( q10bdy3dtemp1 , grid%scalar_b(:,:,:,:,11)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1896                                                                 ids , ide , jds, jde , kds , kde , &
     1897                                                                 ims , ime , jms, jme , kms , kme , &
     1898                                                                 ips , ipe , jps, jpe , kps , kpe )
     1899      CALL stuff_bdy     ( q11bdy3dtemp1 , grid%scalar_b(:,:,:,:,12)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1900                                                                 ids , ide , jds, jde , kds , kde , &
     1901                                                                 ims , ime , jms, jme , kms , kme , &
     1902                                                                 ips , ipe , jps, jpe , kps , kpe )
     1903      CALL stuff_bdy     ( q12bdy3dtemp1 , grid%scalar_b(:,:,:,:,13)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1904                                                                 ids , ide , jds, jde , kds , kde , &
     1905                                                                 ims , ime , jms, jme , kms , kme , &
     1906                                                                 ips , ipe , jps, jpe , kps , kpe )
     1907      CALL stuff_bdy     ( q13bdy3dtemp1 , grid%scalar_b(:,:,:,:,14)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1908                                                                 ids , ide , jds, jde , kds , kde , &
     1909                                                                 ims , ime , jms, jme , kms , kme , &
     1910                                                                 ips , ipe , jps, jpe , kps , kpe )
     1911      CALL stuff_bdy     ( q14bdy3dtemp1 , grid%scalar_b(:,:,:,:,15)   , 'T' ,ijds , ijde , spec_bdy_width      , &                                                                 
     1912                                                                 ids , ide , jds, jde , kds , kde , &
     1913                                                                 ims , ime , jms, jme , kms , kme , &
     1914                                                                 ips , ipe , jps, jpe , kps , kpe )
     1915      CALL stuff_bdy     ( q15bdy3dtemp1 , grid%scalar_b(:,:,:,:,16)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1916                                                                 ids , ide , jds, jde , kds , kde , &
     1917                                                                 ims , ime , jms, jme , kms , kme , &
     1918                                                                 ips , ipe , jps, jpe , kps , kpe )
     1919      CALL stuff_bdy     ( q16bdy3dtemp1 , grid%scalar_b(:,:,:,:,17)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1920                                                                 ids , ide , jds, jde , kds , kde , &
     1921                                                                 ims , ime , jms, jme , kms , kme , &
     1922                                                                 ips , ipe , jps, jpe , kps , kpe )
     1923      CALL stuff_bdy     ( q17bdy3dtemp1 , grid%scalar_b(:,:,:,:,18)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1924                                                                 ids , ide , jds, jde , kds , kde , &
     1925                                                                 ims , ime , jms, jme , kms , kme , &
     1926                                                                 ips , ipe , jps, jpe , kps , kpe )
     1927      CALL stuff_bdy     ( q18bdy3dtemp1 , grid%scalar_b(:,:,:,:,19)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1928                                                                 ids , ide , jds, jde , kds , kde , &
     1929                                                                 ims , ime , jms, jme , kms , kme , &
     1930                                                                 ips , ipe , jps, jpe , kps , kpe )
     1931      CALL stuff_bdy     ( q19bdy3dtemp1 , grid%scalar_b(:,:,:,:,20)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1932                                                                 ids , ide , jds, jde , kds , kde , &
     1933                                                                 ims , ime , jms, jme , kms , kme , &
     1934                                                                 ips , ipe , jps, jpe , kps , kpe )
     1935      CALL stuff_bdy     ( q20bdy3dtemp1 , grid%scalar_b(:,:,:,:,21)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1936                                                                 ids , ide , jds, jde , kds , kde , &
     1937                                                                 ims , ime , jms, jme , kms , kme , &
     1938                                                                 ips , ipe , jps, jpe , kps , kpe )
     1939      CALL stuff_bdy     ( q21bdy3dtemp1 , grid%scalar_b(:,:,:,:,22)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1940                                                                 ids , ide , jds, jde , kds , kde , &
     1941                                                                 ims , ime , jms, jme , kms , kme , &
     1942                                                                 ips , ipe , jps, jpe , kps , kpe )
     1943      CALL stuff_bdy     ( q22bdy3dtemp1 , grid%scalar_b(:,:,:,:,23)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1944                                                                 ids , ide , jds, jde , kds , kde , &
     1945                                                                 ims , ime , jms, jme , kms , kme , &
     1946                                                                 ips , ipe , jps, jpe , kps , kpe )
     1947      CALL stuff_bdy     ( q23bdy3dtemp1 , grid%scalar_b(:,:,:,:,24)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1948                                                                 ids , ide , jds, jde , kds , kde , &
     1949                                                                 ims , ime , jms, jme , kms , kme , &
     1950                                                                 ips , ipe , jps, jpe , kps , kpe )
     1951      CALL stuff_bdy     ( q24bdy3dtemp1 , grid%scalar_b(:,:,:,:,25)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1952                                                                 ids , ide , jds, jde , kds , kde , &
     1953                                                                 ims , ime , jms, jme , kms , kme , &
     1954                                                                 ips , ipe , jps, jpe , kps , kpe )
     1955      CALL stuff_bdy     ( q25bdy3dtemp1 , grid%scalar_b(:,:,:,:,26)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1956                                                                 ids , ide , jds, jde , kds , kde , &
     1957                                                                 ims , ime , jms, jme , kms , kme , &
     1958                                                                 ips , ipe , jps, jpe , kps , kpe )
     1959      CALL stuff_bdy     ( q26bdy3dtemp1 , grid%scalar_b(:,:,:,:,27)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1960                                                                 ids , ide , jds, jde , kds , kde , &
     1961                                                                 ims , ime , jms, jme , kms , kme , &
     1962                                                                 ips , ipe , jps, jpe , kps , kpe )
     1963      CALL stuff_bdy     ( q27bdy3dtemp1 , grid%scalar_b(:,:,:,:,28)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1964                                                                 ids , ide , jds, jde , kds , kde , &
     1965                                                                 ims , ime , jms, jme , kms , kme , &
     1966                                                                 ips , ipe , jps, jpe , kps , kpe )
     1967      CALL stuff_bdy     ( q28bdy3dtemp1 , grid%scalar_b(:,:,:,:,29)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1968                                                                 ids , ide , jds, jde , kds , kde , &
     1969                                                                 ims , ime , jms, jme , kms , kme , &
     1970                                                                 ips , ipe , jps, jpe , kps , kpe )
     1971      CALL stuff_bdy     ( q29bdy3dtemp1 , grid%scalar_b(:,:,:,:,30)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1972                                                                 ids , ide , jds, jde , kds , kde , &
     1973                                                                 ims , ime , jms, jme , kms , kme , &
     1974                                                                 ips , ipe , jps, jpe , kps , kpe )
     1975      CALL stuff_bdy     ( q30bdy3dtemp1 , grid%scalar_b(:,:,:,:,31)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1976                                                                 ids , ide , jds, jde , kds , kde , &
     1977                                                                 ims , ime , jms, jme , kms , kme , &
     1978                                                                 ips , ipe , jps, jpe , kps , kpe )
     1979      CALL stuff_bdy     ( q31bdy3dtemp1 , grid%scalar_b(:,:,:,:,33)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1980                                                                 ids , ide , jds, jde , kds , kde , &
     1981                                                                 ims , ime , jms, jme , kms , kme , &
     1982                                                                 ips , ipe , jps, jpe , kps , kpe )
     1983      CALL stuff_bdy     ( q32bdy3dtemp1 , grid%scalar_b(:,:,:,:,33)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1984                                                                 ids , ide , jds, jde , kds , kde , &
     1985                                                                 ims , ime , jms, jme , kms , kme , &
     1986                                                                 ips , ipe , jps, jpe , kps , kpe )
     1987      CALL stuff_bdy     ( q33bdy3dtemp1 , grid%scalar_b(:,:,:,:,34)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1988                                                                 ids , ide , jds, jde , kds , kde , &
     1989                                                                 ims , ime , jms, jme , kms , kme , &
     1990                                                                 ips , ipe , jps, jpe , kps , kpe )
     1991      CALL stuff_bdy     ( q34bdy3dtemp1 , grid%scalar_b(:,:,:,:,35)   , 'T' ,ijds , ijde , spec_bdy_width      , &
     1992                                                                 ids , ide , jds, jde , kds , kde , &
     1993                                                                 ims , ime , jms, jme , kms , kme , &
     1994                                                                 ips , ipe , jps, jpe , kps , kpe )
     1995ENDIF
    11941996!!!!MARS <<
    11951997         CALL stuff_bdy     ( mbdy2dtemp1 , grid%em_mu_b    , 'M' , ijds , ijde , spec_bdy_width      , &

  • trunk/MESOSCALE/LMD_MM_MARS/SRC/WRFV2/phys/dynphy_wrf_venus_lmd_new/update_inputs_physiq_mod.F

    r1755 r2021  
    7070SUBROUTINE update_inputs_physiq_tracers(nq,MARS_MODE)
    7171
    72   !use infotrac
    73  
     72  use infotrac, only: nqtot,tname
     73  use logic_mod, only: iflag_trac
    7474  implicit none
    7575  INTEGER, INTENT(IN) :: nq,MARS_MODE
    7676
    77   ALLOCATE(traceurs(nq))
    78   traceurs(:)='zolbxs'
    79   !call infotrac_init ! Need traceur.def
    80 
     77  IF (.not.ALLOCATED(tname)) ALLOCATE(tname(nq))
     78  nqtot=nq
     79
     80  IF (MARS_MODE .EQ. 34) THEN
     81    tname(1)="co2"
     82    tname(2)="co"
     83    tname(3)="h2"
     84    tname(4)="h2o"
     85    tname(5)="o1d"
     86    tname(6)="o"
     87    tname(7)="o2"
     88    tname(8)="o2dg"
     89    tname(9)="o3"
     90    tname(10)="h"
     91    tname(11)="oh"
     92    tname(12)="ho2"
     93    tname(13)="h2o2"
     94    tname(14)="cl"
     95    tname(15)="clo"
     96    tname(16)="cl2"
     97    tname(17)="hcl"
     98    tname(18)="hocl"
     99    tname(19)="clco"
     100    tname(20)="clco3"
     101    tname(21)="cocl2"
     102    tname(22)="s"
     103    tname(23)="so"
     104    tname(24)="so2"
     105    tname(25)="so3"
     106    tname(26)="s2o2"
     107    tname(27)="ocs"
     108    tname(28)="hso3"
     109    tname(29)="h2so4"
     110    tname(30)="s2"
     111    tname(31)="clso2"
     112    tname(32)="oscl"
     113    tname(33)="h2oliq"
     114    tname(34)="h2so4liq"
     115    iflag_trac=1
     116  ELSE
     117    tname(:)='zolbxs'
     118  ENDIF
    81119
    82120END SUBROUTINE update_inputs_physiq_tracers
Note: See TracChangeset for help on using the changeset viewer.