source: LMDZ5/trunk/libf/dyn3dmem/gr_u_scal_loc.F @ 1632

Last change on this file since 1632 was 1632, checked in by Laurent Fairhead, 12 years ago

Import initial du répertoire dyn3dmem

Attention! ceci n'est qu'une version préliminaire du code "basse mémoire":
le code contenu dans ce répertoire est basé sur la r1320 et a donc besoin
d'être mis à jour par rapport à la dynamique parallèle d'aujourd'hui.
Ce code est toutefois mis à disposition pour circonvenir à des problèmes
de mémoire que certaines configurations du modèle pourraient rencontrer.
Dans l'état, il compile et tourne sur vargas et au CCRT


Initial import of dyn3dmem

Warning! this is just a preliminary version of the memory light code:
it is based on r1320 of the code and thus needs to be updated before
it can replace the present dyn3dpar code. It is nevertheless put at your
disposal to circumvent some memory problems some LMDZ configurations may
encounter. In its present state, it will compile and run on vargas and CCRT

File size: 1.4 KB
RevLine 
[1632]1!
2! $Header$
3!
4      SUBROUTINE gr_u_scal_loc(nx,x_u,x_scal)
5c%W%    %G%
6c=======================================================================
7c
8c   Author:    Frederic Hourdin      original: 11/11/92
9c   -------
10c
11c   Subject:
12c   ------
13c
14c   Method:
15c   --------
16c
17c   Interface:
18c   ----------
19c
20c      Input:
21c      ------
22c
23c      Output:
24c      -------
25c
26c=======================================================================
27      USE parallel
28      IMPLICIT NONE
29c-----------------------------------------------------------------------
30c   Declararations:
31c   ---------------
32
33#include "dimensions.h"
34#include "paramet.h"
35#include "comgeom.h"
36
37c   Arguments:
38c   ----------
39
40      INTEGER nx
41      REAL x_u(ijb_u:ije_u,nx),x_scal(ijb_u:ije_u,nx)
42
43c   Local:
44c   ------
45
46      INTEGER l,ij
47      INTEGER :: ijb,ije
48
49c-----------------------------------------------------------------------
50      ijb=ij_begin
51      ije=ij_end
52
53!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)       
54      DO l=1,nx
55         DO ij=ijb+1,ije
56            x_scal(ij,l)=
57     s      (aireu(ij)*x_u(ij,l)+aireu(ij-1)*x_u(ij-1,l))
58     s      /(aireu(ij)+aireu(ij-1))
59         ENDDO
60      ENDDO
61!$OMP ENDDO NOWAIT
62
63      ijb=ij_begin
64      ije=ij_end
65
66!$OMP DO SCHEDULE(STATIC,OMP_CHUNK)       
67      DO l=1,nx
68         DO ij=ijb,ije-iip1+1,iip1
69           x_scal(ij,l)=x_scal(ij+iip1-1,l)
70         ENDDO
71      ENDDO
72!$OMP ENDDO NOWAIT
73      RETURN
74     
75      END
Note: See TracBrowser for help on using the repository browser.