[524] | 1 | ! |
---|
| 2 | ! $Header$ |
---|
| 3 | ! |
---|
[5246] | 4 | SUBROUTINE gr_dyn_fi(nfield,im,jm,ngrid,pdyn,pfi) |
---|
| 5 | IMPLICIT NONE |
---|
| 6 | !======================================================================= |
---|
| 7 | ! passage d'un champ de la grille scalaire a la grille physique |
---|
| 8 | !======================================================================= |
---|
[524] | 9 | |
---|
[5246] | 10 | !----------------------------------------------------------------------- |
---|
| 11 | ! declarations: |
---|
| 12 | ! ------------- |
---|
[524] | 13 | |
---|
[5246] | 14 | INTEGER :: im,jm,ngrid,nfield |
---|
| 15 | REAL :: pdyn(im,jm,nfield) |
---|
| 16 | REAL :: pfi(ngrid,nfield) |
---|
[524] | 17 | |
---|
[5246] | 18 | INTEGER :: j,ifield,ig |
---|
[524] | 19 | |
---|
[5246] | 20 | !----------------------------------------------------------------------- |
---|
| 21 | ! calcul: |
---|
| 22 | ! ------- |
---|
[524] | 23 | |
---|
[5246] | 24 | IF (ngrid.NE.2+(jm-2)*(im-1)) then |
---|
| 25 | call abort_gcm("gr_dyn_fi", 'probleme de dim', 1) |
---|
| 26 | end if |
---|
| 27 | ! traitement des poles |
---|
| 28 | CALL SCOPY(nfield,pdyn,im*jm,pfi,ngrid) |
---|
| 29 | CALL SCOPY(nfield,pdyn(1,jm,1),im*jm,pfi(ngrid,1),ngrid) |
---|
[524] | 30 | |
---|
[5246] | 31 | ! traitement des point normaux |
---|
| 32 | DO ifield=1,nfield |
---|
| 33 | DO j=2,jm-1 |
---|
| 34 | ig=2+(j-2)*(im-1) |
---|
| 35 | CALL SCOPY(im-1,pdyn(1,j,ifield),1,pfi(ig,ifield),1) |
---|
| 36 | ENDDO |
---|
| 37 | ENDDO |
---|
[524] | 38 | |
---|
[5246] | 39 | RETURN |
---|
| 40 | END SUBROUTINE gr_dyn_fi |
---|