source: LMDZ6/trunk/libf/dyn3dpar/covnat.F @ 3770

Last change on this file since 3770 was 1907, checked in by lguez, 11 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.

  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
File size: 1.1 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat )
5      IMPLICIT NONE
6
7c=======================================================================
8c
9c   Auteur:  F Hourdin Phu LeVan
10c   -------
11c
12c   Objet:
13c   ------
14c
15c  *********************************************************************
16c    calcul des compos. naturelles a partir des comp.covariantes
17c  ********************************************************************
18c
19c=======================================================================
20
21#include "dimensions.h"
22#include "paramet.h"
23#include "comgeom.h"
24
25      INTEGER klevel
26      REAL ucov( ip1jmp1,klevel ),  vcov( ip1jm,klevel )
27      REAL unat( ip1jmp1,klevel ), vnat( ip1jm,klevel )
28      INTEGER   l,ij
29
30
31      DO l = 1,klevel
32         DO ij = 1, iip1
33            unat (ij,l) =0.
34         END DO
35
36         DO ij = iip2, ip1jm
37            unat( ij,l ) = ucov( ij,l ) / cu(ij)
38         ENDDO
39         DO ij = ip1jm+1, ip1jmp1 
40            unat (ij,l) =0.
41         END DO
42
43         DO ij = 1,ip1jm
44            vnat( ij,l ) = vcov( ij,l ) / cv(ij)
45         ENDDO
46
47      ENDDO
48      RETURN
49      END
Note: See TracBrowser for help on using the repository browser.