*=*=*=*= inirappel.html =*=*=*=*
SUBROUTINE inirappel

SUBROUTINE inirappel


      SUBROUTINE inirappel(teq)
c%W%    %G%
c=======================================================================
c
c   Author:    Frederic Hourdin      original: 15/01/93
c   -------
c
c   Subject:
c   ------
c
c   Method:
c   --------
c
c   Interface:
c   ----------
c
c      Input:
c      ------
c
c      Output:
c      -------
c
c=======================================================================
      IMPLICIT NONE
c-----------------------------------------------------------------------
c   Declararations:
c   ---------------

#include "dimensions.h"
#include "paramet.h"
#include "comvert.h"
#include "comconst.h"
#include "comgeom.h"

c   Arguments:
c   ----------

      REAL teq(iip1,jjp1,llm)

c   Local:
c   ------

      REAL dsin,teqj,teql
      INTEGER i,j,l

c-----------------------------------------------------------------------

      DO l=1,llm
         teql=1./8.*(-log(s(l))/kappa-.5)
         DO j=1,jjp1
            dsin=sin(rlatu(j))-sin(pi/20.)
            teqj=cpp/1000.**kappa*
     s		300.*(1+1./18.*(1.-3.*dsin*dsin)+teql)
      		PRINT*,j,l,teqj
            DO i=1,iip1
               teq(i,j,l)=teqj
            ENDDO
         ENDDO
      ENDDO
      RETURN
      END