source: trunk/LMDZ.COMMON/libf/dyn3dpar/covcont_p.F @ 3553

Last change on this file since 3553 was 1441, checked in by emillour, 10 years ago

Updates in common dynamics (seq and ) to keep up with updates
in LMDZ5 (up to LMDZ5 trunk, rev 2250):

  • compilation:
  • added test in grid/dimension/makdim to check that # of longitudes is a multiple of 8
  • dyn3d_common:

Bug correction concerning zoom (cf LMDZ5 rev 2218)

  • coefpoly.F becomes coefpoly_m.F90 (in misc)
  • fxhyp.F => fxhyp_m.F90 , fyhyp.F => fyhyp_m.F90
  • new routines for zoom: invert_zoom_x_m.F90 and principal_cshift_m.F90
  • inigeom.F adapted to new zoom definition routines
  • fluxstokenc.F : got rid of calls to initial0()
  • dyn3d:
  • advtrac.F90 : got rid of calls to initial0()
  • conf_gcm.F90 : cosmetic changes and change in default dzoomx,dzoomy values
  • guide_mod.F90 : followed updates from Earth Model
  • gcm.F is now gcm.F90
  • dyn3dpar:
  • advtrac_p.F90, covcont_p.F90, mod_hallo.F90 : cosmetic changes
  • conf_gcm.F90 : cosmetic and changed in default dzoomx,dzoomy values
  • parallel_lmdz.F90 : updates to keep up with Earth model
  • misc:
  • arth.F90 becomes arth_m.F90
  • wxios.F90 updated wrt Earth model changes
  • nrtype.F90 and coefpoly_m.F90 added
  • ran1.F, sort.F, minmax.F, minmax2.F, juldate.F moved over from dyn3d_common

EM

File size: 1.5 KB
Line 
1      SUBROUTINE covcont_p (klevel,ucov, vcov, ucont, vcont )
2      USE parallel_lmdz, ONLY: ij_begin,ij_end,OMP_CHUNK,
3     &                         pole_nord, pole_sud
4      IMPLICIT NONE
5
6c=======================================================================
7c
8c   Auteur:  P. Le Van
9c   -------
10c
11c   Objet:
12c   ------
13c
14c  *********************************************************************
15c    calcul des compos. contravariantes a partir des comp.covariantes
16c  ********************************************************************
17c
18c=======================================================================
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 ucont( ip1jmp1,klevel ), vcont( ip1jm,klevel )
27      INTEGER   l,ij
28      INTEGER ijb_u,ijb_v,ije_u,ije_v
29
30     
31      ijb_u=ij_begin-iip1
32      ijb_v=ij_begin-iip1
33      ije_u=ij_end+iip1
34      ije_v=ij_end+iip1
35     
36      if (pole_nord) then
37        ijb_u=ij_begin+iip1
38        ijb_v=ij_begin
39      endif
40     
41      if (pole_sud) then
42        ije_u=ij_end-iip1
43        ije_v=ij_end-iip1
44      endif
45
46c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)   
47      DO 10 l = 1,klevel
48
49      DO 2  ij = ijb_u,ije_u
50      ucont( ij,l ) = ucov( ij,l ) * unscu2( ij )
51   2  CONTINUE
52
53      DO 4 ij = ijb_v,ije_v
54      vcont( ij,l ) = vcov( ij,l ) * unscv2( ij )
55   4  CONTINUE
56
57  10  CONTINUE
58c$OMP END DO NOWAIT
59      RETURN
60      END
Note: See TracBrowser for help on using the repository browser.