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

Last change on this file since 5118 was 5106, checked in by abarral, 4 months ago

Turn coefils.h into lmdz_coefils.f90
Put filtreg.F90 inside lmdz_filtreg.F90
Turn mod_filtreg_p.F90 into lmdz_filtreg_p.F90
Delete obsolete parafilt.h*
(lint) remove spaces between routine name and args

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