source: LMDZ5/trunk/libf/dyn3dmem/gr_fi_dyn_p.F @ 1985

Last change on this file since 1985 was 1907, checked in by lguez, 11 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
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

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.

  • 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
File size: 1.3 KB
RevLine 
[1632]1!
[1673]2! $Id$
[1632]3!
4      SUBROUTINE gr_fi_dyn_p(nfield,ngrid,im,jm,pfi,pdyn)
[1673]5#ifdef CPP_PHYS
[1632]6! Interface with parallel physics,
7      USE mod_interface_dyn_phys
8      USE dimphy
[1823]9      USE parallel_lmdz
[1632]10      IMPLICIT NONE
11c=======================================================================
12c   passage d'un champ de la grille scalaire a la grille physique
13c=======================================================================
14
15c-----------------------------------------------------------------------
16c   declarations:
17c   -------------
18
19      INTEGER im,jm,ngrid,nfield
20      REAL pdyn(im,jm,nfield)
21      REAL pfi(ngrid,nfield)
22
23      INTEGER i,j,ifield,ig
24
25c-----------------------------------------------------------------------
26c   calcul:
27c   -------
28c$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
29      DO ifield=1,nfield
30
31        do ig=1,klon
32          i=index_i(ig)
33          j=index_j(ig)
34          pdyn(i,j,ifield)=pfi(ig,ifield)
35          if (i==1) pdyn(im,j,ifield)=pdyn(i,j,ifield)
36        enddo
37
38c   traitement des poles
39      if (pole_nord) then
40        do i=1,im
41          pdyn(i,1,ifield)=pdyn(1,1,ifield)
42        enddo
43      endif
44       
45      if (pole_sud) then
46        do i=1,im
47          pdyn(i,jm,ifield)=pdyn(1,jm,ifield)
48        enddo
49      endif
50     
51      ENDDO
52c$OMP END DO NOWAIT
53#endif
[1673]54! of #ifdef CPP_PHYS
[1632]55      RETURN
56      END
Note: See TracBrowser for help on using the repository browser.