source: LMDZ6/trunk/libf/dynphy_lonlat/gr_fi_dyn_p.F @ 3981

Last change on this file since 3981 was 2239, checked in by Ehouarn Millour, 9 years ago

Reorganizing physics/dynamics interface:

  • what is related to dynamics-physics interface is now in a seperate directory: dynlmdz_phy* for physics in phy*
  • 1d model and related dependencies (including a couple from "dynamics", set up as symbolic links) is now in subdirectory "dyn1d" of phy*.
  • "bibio" directory is now "misc" and should only contain autonomous utilities.
  • "cosp" is now a subdirectory of phylmd.

EM
  • 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
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 1.3 KB
Line 
1!
2! $Id: gr_fi_dyn_p.F 2239 2015-03-23 07:27:30Z fairhead $
3!
4      SUBROUTINE gr_fi_dyn_p(nfield,ngrid,im,jm,pfi,pdyn)
5#ifdef CPP_PARA
6! Interface with parallel physics,
7      USE mod_interface_dyn_phys
8      USE dimphy
9      USE parallel_lmdz
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
54! of #ifdef CPP_PARA
55      RETURN
56      END
Note: See TracBrowser for help on using the repository browser.