source: codes/icosagcm/trunk/src/pression.f90 @ 156

Last change on this file since 156 was 151, checked in by ymipsl, 11 years ago

Implementation of mixte parallelism MPI/OpenMP into src directory

YM
File size: 967 bytes
Line 
1MODULE pression_mod
2
3CONTAINS
4
5  SUBROUTINE pression(f_ps,f_p)
6  USE icosa
7  IMPLICIT NONE
8    TYPE(t_field), POINTER :: f_ps(:)
9    TYPE(t_field), POINTER :: f_p(:)
10 
11    REAL(rstd), POINTER :: ps(:)
12    REAL(rstd), POINTER :: p(:,:)
13    INTEGER :: ind
14
15    DO ind=1,ndomain
16      CALL swap_dimensions(ind)
17      CALL swap_geometry(ind)
18      ps=f_ps(ind)
19      p=f_p(ind)
20      CALL compute_pression(ps, p,0)
21    ENDDO
22 
23  END SUBROUTINE pression
24 
25  SUBROUTINE compute_pression(ps,p,offset)
26  USE icosa
27  USE disvert_mod
28  IMPLICIT NONE
29    REAL(rstd),INTENT(IN) :: ps(iim*jjm)
30    REAL(rstd),INTENT(OUT) :: p(iim*jjm,llm+1)
31    INTEGER,INTENT(IN) :: offset
32    INTEGER :: i,j,ij,l
33
34    DO    l    = 1, llm+1
35      DO j=jj_begin-offset,jj_end+offset
36        DO i=ii_begin-offset,ii_end+offset
37          ij=(j-1)*iim+i
38          p(ij,l) = ap(l) + bp(l) * ps(ij)
39        ENDDO
40      ENDDO
41    ENDDO
42   
43  END SUBROUTINE compute_pression
44
45END MODULE pression_mod
Note: See TracBrowser for help on using the repository browser.