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

!
! $Header$
!
SUBROUTINE covnat (klevel,ucov, vcov, unat, vnat )
  USE comgeom_mod_h
  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
  USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
  IMPLICIT NONE

  !=======================================================================
  !
  !   Auteur:  F Hourdin Phu LeVan
  !   -------
  !
  !   Objet:
  !   ------
  !
  !  *********************************************************************
  !    calcul des compos. naturelles a partir des comp.covariantes
  !  ********************************************************************
  !
  !=======================================================================

  INTEGER :: klevel
  REAL :: ucov( ip1jmp1,klevel ),  vcov( ip1jm,klevel )
  REAL :: unat( ip1jmp1,klevel ), vnat( ip1jm,klevel )
  INTEGER :: l,ij


  DO l = 1,klevel
     DO ij = 1, iip1
        unat (ij,l) =0.
     END DO

     DO ij = iip2, ip1jm
        unat( ij,l ) = ucov( ij,l ) / cu(ij)
     ENDDO
     DO ij = ip1jm+1, ip1jmp1
        unat (ij,l) =0.
     END DO

     DO ij = 1,ip1jm
        vnat( ij,l ) = vcov( ij,l ) / cv(ij)
     ENDDO

  ENDDO
  RETURN
END SUBROUTINE covnat