source: LMDZ6/trunk/libf/dyn3d/covnat.F90 @ 5273

Last change on this file since 5273 was 5272, checked in by abarral, 2 days ago

Turn paramet.h into a module
  • 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.2 KB
Line 
1!
2! $Header$
3!
4SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat )
5  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
6  USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
7          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
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  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
49END SUBROUTINE covnat
Note: See TracBrowser for help on using the repository browser.