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

Last change on this file since 5159 was 5159, checked in by abarral, 7 weeks ago

Put dimensions.h and paramet.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.1 KB
Line 
1! $Header$
2
3SUBROUTINE covnat(klevel, ucov, vcov, unat, vnat)
4  USE lmdz_comgeom
5
6  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
7  USE lmdz_paramet
8  IMPLICIT NONE
9
10  !=======================================================================
11
12  !   Auteur:  F Hourdin Phu LeVan
13  !   -------
14
15  !   Objet:
16  !   ------
17
18  !  *********************************************************************
19  !    calcul des compos. naturelles a partir des comp.covariantes
20  !  ********************************************************************
21
22  !=======================================================================
23
24
25
26
27  INTEGER :: klevel
28  REAL :: ucov(ip1jmp1, klevel), vcov(ip1jm, klevel)
29  REAL :: unat(ip1jmp1, klevel), vnat(ip1jm, klevel)
30  INTEGER :: l, ij
31
32  DO l = 1, klevel
33    DO ij = 1, iip1
34      unat (ij, l) = 0.
35    END DO
36
37    DO ij = iip2, ip1jm
38      unat(ij, l) = ucov(ij, l) / cu(ij)
39    ENDDO
40    DO ij = ip1jm + 1, ip1jmp1
41      unat (ij, l) = 0.
42    END DO
43
44    DO ij = 1, ip1jm
45      vnat(ij, l) = vcov(ij, l) / cv(ij)
46    ENDDO
47
48  ENDDO
49
50END SUBROUTINE covnat
Note: See TracBrowser for help on using the repository browser.