source: LMDZ6/branches/Amaury_dev/libf/dyn3d/covnat.F90 @ 5136

Last change on this file since 5136 was 5136, checked in by abarral, 8 weeks ago

Put comgeom.h, comgeom2.h into modules

  • 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.0 KB
Line 
1! $Header$
2
3SUBROUTINE covnat(klevel, ucov, vcov, unat, vnat)
4  USE lmdz_comgeom
5
6  IMPLICIT NONE
7
8  !=======================================================================
9  !
10  !   Auteur:  F Hourdin Phu LeVan
11  !   -------
12  !
13  !   Objet:
14  !   ------
15  !
16  !  *********************************************************************
17  !    calcul des compos. naturelles a partir des comp.covariantes
18  !  ********************************************************************
19  !
20  !=======================================================================
21
22 INCLUDE "dimensions.h"
23 INCLUDE "paramet.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  DO l = 1, klevel
31    DO ij = 1, iip1
32      unat (ij, l) = 0.
33    END DO
34
35    DO ij = iip2, ip1jm
36      unat(ij, l) = ucov(ij, l) / cu(ij)
37    ENDDO
38    DO ij = ip1jm + 1, ip1jmp1
39      unat (ij, l) = 0.
40    END DO
41
42    DO ij = 1, ip1jm
43      vnat(ij, l) = vcov(ij, l) / cv(ij)
44    ENDDO
45
46  ENDDO
47
48END SUBROUTINE covnat
Note: See TracBrowser for help on using the repository browser.