source: trunk/LMDZ.COMMON/libf/dynphy_lonlat/gr_dyn_fi.F @ 3365

Last change on this file since 3365 was 1403, checked in by emillour, 10 years ago

All models: Reorganizing the physics/dynamics interface.

  • makelmdz and makelmdz_fcm scripts adapted to handle the new directory settings
  • misc: (replaces what was the "bibio" directory)
  • Should only contain extremely generic (and non physics or dynamics-specific) routines
  • Therefore moved initdynav.F90, initfluxsto.F, inithist.F, writedynav.F90, write_field.F90, writehist.F to "dyn3d_common"
  • dynlonlat_phylonlat: (new interface directory)
  • This directory contains routines relevent to physics/dynamics grid interactions, e.g. routines gr_dyn_fi or gr_fi_dyn and calfis
  • Moreover the dynlonlat_phylonlat contains directories "phy*" corresponding to each physics package "phy*" to be used. These subdirectories should only contain specific interfaces (e.g. iniphysiq) or main programs (e.g. newstart)
  • phy*/dyn1d: this subdirectory contains the 1D model using physics from phy*

EM

File size: 1.0 KB
Line 
1!
2! $Header$
3!
4      SUBROUTINE gr_dyn_fi(nfield,im,jm,ngrid,pdyn,pfi)
5      IMPLICIT NONE
6c=======================================================================
7c   passage d'un champ de la grille scalaire a la grille physique
8c=======================================================================
9
10c-----------------------------------------------------------------------
11c   declarations:
12c   -------------
13
14      INTEGER im,jm,ngrid,nfield
15      REAL pdyn(im,jm,nfield)
16      REAL pfi(ngrid,nfield)
17
18      INTEGER j,ifield,ig
19
20c-----------------------------------------------------------------------
21c   calcul:
22c   -------
23
24      IF(ngrid.NE.2+(jm-2)*(im-1)) STOP 'probleme de dim'
25c   traitement des poles
26      CALL SCOPY(nfield,pdyn,im*jm,pfi,ngrid)
27      CALL SCOPY(nfield,pdyn(1,jm,1),im*jm,pfi(ngrid,1),ngrid)
28
29c   traitement des point normaux
30      DO ifield=1,nfield
31         DO j=2,jm-1
32            ig=2+(j-2)*(im-1)
33            CALL SCOPY(im-1,pdyn(1,j,ifield),1,pfi(ig,ifield),1)
34         ENDDO
35      ENDDO
36
37      RETURN
38      END
Note: See TracBrowser for help on using the repository browser.