source: LMDZ.3.3/trunk/libf/dyn3d/inirappelF @ 61

Last change on this file since 61 was 2, checked in by lmdz, 25 years ago

Initial revision

  • Property svn:eol-style set to native
  • Property svn:executable set to *
  • Property svn:keywords set to Author Date Id Revision
File size: 1.2 KB
Line 
1      SUBROUTINE inirappel(teq)
2c%W%    %G%
3c=======================================================================
4c
5c   Author:    Frederic Hourdin      original: 15/01/93
6c   -------
7c
8c   Subject:
9c   ------
10c
11c   Method:
12c   --------
13c
14c   Interface:
15c   ----------
16c
17c      Input:
18c      ------
19c
20c      Output:
21c      -------
22c
23c=======================================================================
24      IMPLICIT NONE
25c-----------------------------------------------------------------------
26c   Declararations:
27c   ---------------
28
29#include "dimensions.h"
30#include "paramet.h"
31#include "comvert.h"
32#include "comconst.h"
33#include "comgeom.h"
34
35c   Arguments:
36c   ----------
37
38      REAL teq(iip1,jjp1,llm)
39
40c   Local:
41c   ------
42
43      REAL dsin,teqj,teql
44      INTEGER i,j,l
45
46c-----------------------------------------------------------------------
47
48      DO l=1,llm
49         teql=1./8.*(-log(s(l))/kappa-.5)
50         DO j=1,jjp1
51            dsin=sin(rlatu(j))-sin(pi/20.)
52            teqj=cpp/1000.**kappa*
53     s          300.*(1+1./18.*(1.-3.*dsin*dsin)+teql)
54                        PRINT*,j,l,teqj
55            DO i=1,iip1
56               teq(i,j,l)=teqj
57            ENDDO
58         ENDDO
59      ENDDO
60      RETURN
61      END
Note: See TracBrowser for help on using the repository browser.