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
RevLine 
[524]1! $Header$
[5099]2
[5106]3SUBROUTINE covnat(klevel, ucov, vcov, unat, vnat)
[5136]4  USE lmdz_comgeom
5
[5103]6  IMPLICIT NONE
[524]7
[5103]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  !=======================================================================
[524]21
[5134]22 INCLUDE "dimensions.h"
23 INCLUDE "paramet.h"
[524]24
[5103]25  INTEGER :: klevel
26  REAL :: ucov(ip1jmp1, klevel), vcov(ip1jm, klevel)
27  REAL :: unat(ip1jmp1, klevel), vnat(ip1jm, klevel)
28  INTEGER :: l, ij
[524]29
[5103]30  DO l = 1, klevel
31    DO ij = 1, iip1
32      unat (ij, l) = 0.
33    END DO
[524]34
[5103]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
[524]41
[5103]42    DO ij = 1, ip1jm
43      vnat(ij, l) = vcov(ij, l) / cv(ij)
44    ENDDO
[524]45
[5103]46  ENDDO
[5105]47
[5103]48END SUBROUTINE covnat
Note: See TracBrowser for help on using the repository browser.