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

Photochemistry architecture for Venus mesoscale and LES model

File:
1 edited

Legend:

Unmodified
Added
Removed
  • 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      , &
Note: See TracChangeset for help on using the changeset viewer.