*=*=*=*= SUAER.html =*=*=*=*
SUBROUTINE SUAER

SUBROUTINE SUAER


      SUBROUTINE SUAER
      implicit none
C
C**** *SUAER*   - INITIALIZE COMMON YOMAER
C
C     PURPOSE.
C     --------
C           INITIALIZE YOMAER, THE COMMON THAT CONTAINS THE
C           RADIATIVE CHARACTERISTICS OF THE AEROSOLS
c
c
c
C
C**   INTERFACE.
C     ----------
C        *CALL* *SUAER*
C
C        EXPLICIT ARGUMENTS :
C        --------------------
C        NONE
C
C        IMPLICIT ARGUMENTS :
C        --------------------
C        COMMON YOMAER
C
C     METHOD.
C     -------
C        SEE DOCUMENTATION
C
C     EXTERNALS.
C     ----------
C
C     REFERENCE.
C     ----------
C        ECMWF RESEARCH DEPARTMENT DOCUMENTATION OF THE
C     "IN CORE MODEL"
C
C     AUTHOR.
C     -------
C        JEAN-JACQUES MORCRETTE *ECMWF*
C
C     MODIFICATIONS.
C     --------------
C        ORIGINAL : 88-02-15
c        Frederic Hourdin
c        Richard Fournier (1996) Francois Forget (1996)
C     ------------------------------------------------------------------
C
C-----------------------------------------------------------------------
C
#include "dimensions.h"
#include "dimphys.h"
#include "dimradmars.h"
#include "yomaer.h"
C
C-----------------------------------------------------------------------
c
      INTEGER isun,iir,itra
      REAL tsun
      REAL longref,epref
      REAL longsun(nsun+1)
      REAL tsol
      REAL longir(nir+1)
      REAL epir(nir),omegir(nir),gir(nir)
      REAL t9s9

c     Aerosol Spectral properties :
#include "aerdata.h"

c   quelques initialisations a 0

      do itra=1,ntra
         do isun=1,nsun
            cga(isun,itra)=0.
            piza(isun,itra)=0.
            taua(isun,itra)=0.
         enddo
         do iir=1,nir
            caer(iir,itra)=0.
         enddo
      enddo
c
C-----------------------------------------------------------------------
C
C      ----------------------------------------------------------------
C
C*       1.    SHORTWAVE COEFFICIENTS
C              ----------------------
C
c
c on utilise le premier traceur uniquement
c
c._. a peine comprehensible mais :
c._.
c    1) on commence par calculer les proprietes optiques
c       moyennes sur les deux bandes solaires.
c       les epaisseurs optiques moyennes sont "en fait" des
c       $\overbar{\tau} \over \tau_{0.67micron}$
c._.
c
      tsun=6000.E+0
      longsun(1)=long1vis
      longsun(2)=long2vis
      longsun(3)=long3vis
      longref=longrefvis
      epref=1.E+0
      CALL aerave ( ndatavis,
     & longdatavis,epdatavis,omegdatavis,gdatavis,
     & longref,epref,tsun
     &             ,nsun,longsun,taua,piza,cga              )
c
C      ----------------------------------------------------------------
C
C*       2.    LONGWAVE COEFFICIENTS
C              ---------------------
C
c
c._.
c    2) puis on calcule les proprietes optiques moyennes
c       sur la bande IRTM a 9 micron.
c       l' epaisseur optique moyenne est "en fait"
c       $\overbar{\tau} \over \tau_{9micron}$
c._.
c
      tsol=215.D+0
      longir(1)=8.3E-6
      longir(2)=9.7E-6
      longref=longrefir
      epref=1.E+0
      CALL aerave ( ndatair,
     & longdatair,epdatair,omegdatair,gdatair,
     & longref,epref,tsol
     &             ,1,longir,epir,omegir,gir )
      t9s9=epir(1)
c
c._.
c    3) puis on calcule les proprietes optiques moyennes
c       sur les trois bandes IR.
c       les epaisseurs optiques moyennes sont "en fait"
c       et "dans un premier temps" des
c       $\overbar{\tau} \over \tau_{9micron}$
c._.
c
      tsol=215.D+0
      longir(1)=long1ir
      longir(2)=long1co2
      longir(3)=long2co2
      longir(4)=long2ir
      longref=longrefir
      epref=1.E+0
      CALL aerave ( ndatair,
     & longdatair,epdatair,omegdatair,gdatair,
     & longref,epref,tsol
     &             ,nir-1,longir,epir,omegir,gir )
c
c._.
c    4) enfin il reste a ecrire les epaisseurs optiques moyennes
c       en
c       $\overbar{\tau} \over \tau_{0.67micron}$
c._.
c
      DO iir=1,nir-1
        epir(iir)=epir(iir)/t9s9/solsir
      ENDDO
c
c
c
c
c co2
c
      qextaer(1)=epir(2)
      qextaer(2)=epir(2)
      omegaer(1)=omegir(2)
      omegaer(2)=omegir(2)
      gaer(1)=gir(2)
      gaer(2)=gir(2)
c
c bande 9 micron
c
      qextaer(3)=epir(1)
      omegaer(3)=omegir(1)
      gaer(3)=gir(1)
c
c grandes longueurs d'ondes
c
      DO iir=4,nir
        qextaer(iir)=epir(iir-1)
        omegaer(iir)=omegir(iir-1)
        gaer(iir)=gir(iir-1)
      ENDDO
c
      PRINT*,'Rapport Solaire/IR :',solsir
      PRINT*

      PRINT *,'Les donnees spectrales :'
      PRINT *,'Solaire ---->'
      PRINT *,'Qext / Qext0.67 ; omega ; g'
      DO isun=1,nsun
c._.
c            taua(isun,1)=1.E+0
c            piza(isun,1)=0.86E+0
c            cga(isun,1)=0.79E+0
c._.
        PRINT *,taua(isun,1),piza(isun,1),cga(isun,1)
      ENDDO
      PRINT *,'IR ---->'
      PRINT *,'Qext / Qext0.67 ; omega ; g'
      DO iir=1,nir
        PRINT *,qextaer(iir),omegaer(iir),gaer(iir)
      ENDDO
c



          DO iir=1,5
            caer(iir,1)=qextaer(1)*(1.E+0 - omegaer(1))
          ENDDO
          print *,'dans le co2 on prend Qabs : caer=',caer(1,1)



      RETURN
      END