source: trunk/LMDZ.MARS/libf/dynphy_lonlat/gr_fi_dyn.F @ 2322

Last change on this file since 2322 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      SUBROUTINE gr_fi_dyn(nfield,ngrid,im,jm,pfi,pdyn)
2      IMPLICIT NONE
3c=======================================================================
4c   passage d'un champ de la grille scalaire a la grille physique
5c=======================================================================
6
7c-----------------------------------------------------------------------
8c   declarations:
9c   -------------
10
11      INTEGER im,jm,ngrid,nfield
12      REAL pdyn(im,jm,nfield)
13      REAL pfi(ngrid,nfield)
14
15      INTEGER i,j,ifield,ig
16      EXTERNAL SCOPY
17
18c-----------------------------------------------------------------------
19c   calcul:
20c   -------
21
22      DO ifield=1,nfield
23c   traitement des poles
24         DO i=1,im
25            pdyn(i,1,ifield)=pfi(1,ifield)
26            pdyn(i,jm,ifield)=pfi(ngrid,ifield)
27         ENDDO
28
29c   traitement des point normaux
30         DO j=2,jm-1
31            ig=2+(j-2)*(im-1)
32            CALL SCOPY(im-1,pfi(ig,ifield),1,pdyn(1,j,ifield),1)
33            pdyn(im,j,ifield)=pdyn(1,j,ifield)
34         ENDDO
35      ENDDO
36
37      RETURN
38      END
Note: See TracBrowser for help on using the repository browser.