source: LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/gr_dyn_fi.f90 @ 5409

Last change on this file since 5409 was 5128, checked in by abarral, 5 months ago

Correct bug in vlspltqs_loc.f90 from r2270 where we call SSUM with incorrect arguments.
Merge the three different versions of abort_gcm into one
Fix seq, para 3D compilation broken from r5107 onwards
(lint) usual + Remove uneeded fixed-form continuations

  • 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.1 KB
Line 
1! $Header$
2
3SUBROUTINE gr_dyn_fi(nfield, im, jm, ngrid, pdyn, pfi)
4  USE lmdz_ssum_scopy, ONLY: scopy
5
6
7  IMPLICIT NONE
8  !=======================================================================
9  !   passage d'un champ de la grille scalaire a la grille physique
10  !=======================================================================
11
12  !-----------------------------------------------------------------------
13  !   declarations:
14  !   -------------
15
16  INTEGER :: im, jm, ngrid, nfield
17  REAL :: pdyn(im, jm, nfield)
18  REAL :: pfi(ngrid, nfield)
19
20  INTEGER :: j, ifield, ig
21
22  !-----------------------------------------------------------------------
23  !   calcul:
24  !   -------
25
26  IF (ngrid/=2 + (jm - 2) * (im - 1)) THEN
27    CALL abort_gcm("gr_dyn_fi", 'probleme de dim', 1)
28  end if
29  !   traitement des poles
30  CALL SCOPY(nfield, pdyn, im * jm, pfi, ngrid)
31  CALL SCOPY(nfield, pdyn(1, jm, 1), im * jm, pfi(ngrid, 1), ngrid)
32
33  !   traitement des point normaux
34  DO ifield = 1, nfield
35    DO j = 2, jm - 1
36      ig = 2 + (j - 2) * (im - 1)
37      CALL SCOPY(im - 1, pdyn(1, j, ifield), 1, pfi(ig, ifield), 1)
38    ENDDO
39  ENDDO
40
41  RETURN
42END SUBROUTINE gr_dyn_fi
Note: See TracBrowser for help on using the repository browser.