source: codes/icosagcm/branches/SATURN_DYNAMICO/LMDZ.COMMON/libf/dyn3dpar/calfis_icosa.f90 @ 222

Last change on this file since 222 was 222, checked in by ymipsl, 10 years ago

Creating temporary dynamico/lmdz/saturn branche

YM
File size: 2.9 KB
Line 
1SUBROUTINE calfis_icosa(dtphy, lafin, day, time, presnivs,      &
2                        plev, play, phi, phis, ulon, ulat, T, q, wflux, &
3                        du, dv, dT, dq, dps) 
4  USE infotrac
5  USE dimphy
6  USE control_mod
7  USE mod_phys_lmdz_para
8  IMPLICIT NONE
9  INCLUDE "comconst.h"
10 
11  REAL, INTENT(IN)    :: dtphy
12  LOGICAL, INTENT(IN) :: lafin
13  REAL, INTENT(IN)    :: day
14  REAL, INTENT(IN)    :: time
15  REAL,INTENT(IN) :: presnivs(klev)
16  REAL,INTENT(IN) :: plev(klon_mpi,klev+1)
17  REAL,INTENT(IN) :: play(klon_mpi,klev)
18  REAL,INTENT(IN) :: phi(klon_mpi,klev)
19  REAL,INTENT(IN) :: phis(klon_mpi)
20  REAL,INTENT(IN) :: ulon(klon_mpi,klev) 
21  REAL,INTENT(IN) :: ulat(klon_mpi,klev)
22  REAL,INTENT(IN) :: T(klon_mpi,klev)
23  REAL,INTENT(IN) :: q(klon_mpi,klev,nqtot)
24  REAL,INTENT(IN) :: wflux(klon_mpi,klev)
25  REAL,INTENT(OUT) :: du(klon_mpi,klev)
26  REAL,INTENT(OUT) :: dv(klon_mpi,klev)
27  REAL,INTENT(OUT) :: dT(klon_mpi,klev)
28  REAL,INTENT(OUT) :: dq(klon_mpi,klev,nqtot)
29  REAL,INTENT(OUT) :: dps(klon_mpi)
30
31  REAL :: ducov(klon_mpi,klev)
32  REAL :: plev_omp(klon_omp,klev+1)
33  REAL :: play_omp(klon_omp,klev)
34  REAL :: phi_omp(klon_omp,klev)
35  REAL :: ulon_omp(klon_omp,klev)
36  REAL :: ulat_omp(klon_omp,klev)
37  REAL :: T_omp(klon_omp,klev)
38  REAL :: q_omp(klon_omp,klev,nqtot)
39  CHARACTER(len=255) :: t_name(nqtot)
40  REAL :: wflux_omp(klon_omp,klev)
41  REAL :: du_omp(klon_omp,klev)
42  REAL :: dv_omp(klon_omp,klev)
43  REAL :: dt_omp(klon_omp,klev)
44  REAL :: dq_omp(klon_omp,klev,nqtot)
45  REAL :: dps_omp(klon_omp)
46  LOGICAL :: tracerdyn
47   
48  LOGICAL,SAVE :: debut=.TRUE.
49  INTEGER :: k,offset
50 
51   
52    offset=klon_omp_begin-1
53    klon=klon_omp
54    DO k=1,klon_omp
55      plev_omp(k,:) = plev(k+offset,:)
56      play_omp(k,:) = play(k+offset,:)
57      phi_omp(k,:)  = phi(k+offset,:)
58      ulon_omp(k,:) = ulon(k+offset,:)
59      ulat_omp(k,:) = ulat(k+offset,:)
60      T_omp(k,:)    = T(k+offset,:)
61      q_omp(k,:,:)  = q(k+offset,:,:)
62      wflux_omp(k,:)= wflux(k+offset,:)
63    ENDDO   
64 
65    CALL physiq  (klon,      &
66                  klev,      &
67                  nqtot,     &
68                  t_name,    &
69                  debut,     &
70                  lafin,     &
71                  day,       &
72                  time,      &
73                  dtphy,     &
74                  plev_omp,  &
75                  play_omp,  &
76                  phi_omp,   &
77                  ulon_omp,  &
78                  ulat_omp,  &
79                  T_omp,     &
80                  q_omp,     &
81                  wflux_omp, &
82                  du_omp,    &
83                  dv_omp,    &
84                  dt_omp,    &
85                  dq_omp,    &
86                  dps_omp,   &
87                  tracerdyn)
88 
89    DO k=1,klon_omp
90      du(k+offset,:) = du_omp(k,:)
91      dv(k+offset,:) = dv_omp(k,:)
92      dq(k+offset,:,:) = dq_omp(k,:,:)
93      dT(k+offset,:) = dT_omp(k,:)
94      dps(k+offset) = dps_omp(k)
95   ENDDO
96  debut=.FALSE.
97 
98END SUBROUTINE calfis_icosa
Note: See TracBrowser for help on using the repository browser.