source: LMDZ6/trunk/libf/dyn3d_common/covcont.f90 @ 5272

Last change on this file since 5272 was 5272, checked in by abarral, 23 hours ago

Turn paramet.h into a module

File size: 1.4 KB
Line 
1SUBROUTINE covcont (klevel,ucov, vcov, ucont, vcont )
2!
3!-------------------------------------------------------------------------------
4! Author: P. Le Van
5!-------------------------------------------------------------------------------
6! Purpose: Compute contravariant components from covariant components.
7!-------------------------------------------------------------------------------
8  USE dimensions_mod, ONLY: iim, jjm, llm, ndm
9USE paramet_mod_h, ONLY: iip1, iip2, iip3, jjp1, llmp1, llmp2, llmm1, kftd, ip1jm, ip1jmp1, &
10          ip1jmi1, ijp1llm, ijmllm, mvar, jcfil, jcfllm
11IMPLICIT NONE
12
13
14  include "comgeom.h"
15!===============================================================================
16! Arguments:
17  INTEGER, INTENT(IN)  :: klevel                    !--- VERTICAL LEVELS NUMBER
18  REAL,    INTENT(IN)  :: ucov ( ip1jmp1,klevel )   !--- U COVARIANT WIND
19  REAL,    INTENT(IN)  :: vcov ( ip1jm  ,klevel )   !--- V COVARIANT WIND
20  REAL,    INTENT(OUT) :: ucont( ip1jmp1,klevel )   !--- U CONTRAVAR WIND
21  REAL,    INTENT(OUT) :: vcont( ip1jm  ,klevel )   !--- V CONTRAVAR WIND
22!===============================================================================
23!   Local variables:
24  INTEGER :: l
25!===============================================================================
26  DO l=1,klevel
27    ucont(iip2:ip1jm,l)=ucov(iip2:ip1jm,l) * unscu2(iip2:ip1jm)
28    vcont(   1:ip1jm,l)=vcov(   1:ip1jm,l) * unscv2(   1:ip1jm)
29  END DO
30
31END SUBROUTINE covcont
32
Note: See TracBrowser for help on using the repository browser.