LMDZ
advtrac_p.F90
Go to the documentation of this file.
1 ! $Id: advtrac_p.F90 2239 2015-03-23 07:27:30Z emillour $
2 
3 SUBROUTINE advtrac_p(pbaru,pbarv , p, masse,q,iapptrac,teta, flxw, pk)
4 
5  ! Auteur : F. Hourdin
6  !
7  ! Modif. P. Le Van (20/12/97)
8  ! F. Codron (10/99)
9  ! D. Le Croller (07/2001)
10  ! M.A Filiberti (04/2002)
11  !
12  USE parallel_lmdz, ONLY: ij_begin,ij_end,omp_chunk,pole_nord,pole_sud,&
13  setdistrib
14  USE write_field_p, ONLY: writefield_p
15  USE bands, ONLY: jj_nb_caldyn,jj_nb_vanleer
16  USE mod_hallo
17  USE vampir
18  USE times
19  USE infotrac, ONLY: nqtot, iadv
20  USE control_mod, ONLY: iapp_tracvl, day_step, planet_type
21  IMPLICIT NONE
22  !
23  include "dimensions.h"
24  include "paramet.h"
25  include "comconst.h"
26  include "comvert.h"
27  include "comdissip.h"
28  include "comgeom2.h"
29  include "logic.h"
30  include "temps.h"
31  include "ener.h"
32  include "description.h"
33 
34  !-------------------------------------------------------------------
35  ! Arguments
36  !-------------------------------------------------------------------
37  INTEGER,INTENT(OUT) :: iapptrac
38  REAL,INTENT(IN) :: pbaru(ip1jmp1,llm)
39  REAL,INTENT(IN) :: pbarv(ip1jm,llm)
40  REAL,INTENT(INOUT) :: q(ip1jmp1,llm,nqtot)
41  REAL,INTENT(IN) :: masse(ip1jmp1,llm)
42  REAL,INTENT(IN) :: p( ip1jmp1,llmp1 )
43  REAL,INTENT(IN) :: teta(ip1jmp1,llm)
44  REAL,INTENT(IN) :: pk(ip1jmp1,llm)
45  REAL,INTENT(OUT) :: flxw(ip1jmp1,llm)
46  !-------------------------------------------------------------------
47  ! Ajout PPM
48  !--------------------------------------------------------
49  REAL massebx(ip1jmp1,llm),masseby(ip1jm,llm)
50  !-------------------------------------------------------------
51  ! Variables locales
52  !-------------------------------------------------------------
53 
54  REAL pbaruc(ip1jmp1,llm),pbarvc(ip1jm,llm)
55  REAL massem(ip1jmp1,llm),zdp(ip1jmp1)
56  REAL,SAVE::pbarug(ip1jmp1,llm),pbarvg(ip1jm,llm),wg(ip1jmp1,llm)
57  REAL (kind=kind(1.d0)) :: t_initial, t_final, tps_cpu
58  INTEGER iadvtr
59  INTEGER ij,l,iq,iiq
60  REAL zdpmin, zdpmax
61  SAVE iadvtr, massem, pbaruc, pbarvc
62  DATA iadvtr/0/
63  !$OMP THREADPRIVATE(iadvtr)
64  !----------------------------------------------------------
65  ! Rajouts pour PPM
66  !----------------------------------------------------------
67  INTEGER indice,n
68  REAL dtbon ! Pas de temps adaptatif pour que CFL<1
69  REAL CFLmaxz,aaa,bbb ! CFL maximum
70  REAL psppm(iim,jjp1) ! pression au sol
71  REAL unatppm(iim,jjp1,llm),vnatppm(iim,jjp1,llm)
72  REAL qppm(iim*jjp1,llm,nqtot)
73  REAL fluxwppm(iim,jjp1,llm)
74  REAL apppm(llmp1), bpppm(llmp1)
75  LOGICAL dum,fill
76  DATA fill/.true./
77  DATA dum/.true./
78  REAL,SAVE :: finmasse(ip1jmp1,llm)
79  integer ijb,ije,ijb_u,ijb_v,ije_u,ije_v,j
80  type(request) :: Request_vanleer
81  REAL,SAVE :: p_tmp( ip1jmp1,llmp1 )
82  REAL,SAVE :: teta_tmp(ip1jmp1,llm)
83  REAL,SAVE :: pk_tmp(ip1jmp1,llm)
84 
85  ijb_u=ij_begin
86  ije_u=ij_end
87 
88  ijb_v=ij_begin-iip1
89  ije_v=ij_end
90  if (pole_nord) ijb_v=ij_begin
91  if (pole_sud) ije_v=ij_end-iip1
92 
93  IF(iadvtr.EQ.0) THEN
94  ! CALL initial0(ijp1llm,pbaruc)
95  ! CALL initial0(ijmllm,pbarvc)
96  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
97  DO l=1,llm
98  pbaruc(ijb_u:ije_u,l)=0.
99  pbarvc(ijb_v:ije_v,l)=0.
100  ENDDO
101  !$OMP END DO NOWAIT
102  ENDIF
103 
104  ! accumulation des flux de masse horizontaux
105  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
106  DO l=1,llm
107  DO ij = ijb_u,ije_u
108  pbaruc(ij,l) = pbaruc(ij,l) + pbaru(ij,l)
109  ENDDO
110  DO ij = ijb_v,ije_v
111  pbarvc(ij,l) = pbarvc(ij,l) + pbarv(ij,l)
112  ENDDO
113  ENDDO
114  !$OMP END DO NOWAIT
115 
116  ! selection de la masse instantannee des mailles avant le transport.
117  IF(iadvtr.EQ.0) THEN
118 
119  ! CALL SCOPY(ip1jmp1*llm,masse,1,massem,1)
120  ijb=ij_begin
121  ije=ij_end
122 
123  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
124  DO l=1,llm
125  massem(ijb:ije,l)=masse(ijb:ije,l)
126  ENDDO
127  !$OMP END DO NOWAIT
128 
129  !cc CALL filtreg ( massem ,jjp1, llm,-2, 2, .TRUE., 1 )
130  !
131  ENDIF ! of IF(iadvtr.EQ.0)
132 
133  iadvtr = iadvtr+1
134 
135  !$OMP MASTER
136  iapptrac = iadvtr
137  !$OMP END MASTER
138 
139  ! Test pour savoir si on advecte a ce pas de temps
140 
141  IF ( iadvtr.EQ.iapp_tracvl ) THEN
142  !$OMP MASTER
143  call suspend_timer(timer_caldyn)
144  !$OMP END MASTER
145 
146  ijb=ij_begin
147  ije=ij_end
148 
149 
150  !c .. Modif P.Le Van ( 20/12/97 ) ....
151  !c
152 
153  ! traitement des flux de masse avant advection.
154  ! 1. calcul de w
155  ! 2. groupement des mailles pres du pole.
156 
157  CALL groupe_p( massem, pbaruc,pbarvc, pbarug,pbarvg,wg )
158 
159  !$OMP BARRIER
160 
161  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
162  DO l=1,llmp1
163  p_tmp(ijb:ije,l)=p(ijb:ije,l)
164  ENDDO
165  !$OMP END DO NOWAIT
166 
167  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
168  DO l=1,llm
169  pk_tmp(ijb:ije,l)=pk(ijb:ije,l)
170  teta_tmp(ijb:ije,l)=teta(ijb:ije,l)
171  ENDDO
172  !$OMP END DO NOWAIT
173 
174  !$OMP MASTER
175  call vtb(vthallo)
176  !$OMP END MASTER
177 
178  call register_swapfieldhallo(pbarug,pbarug,ip1jmp1,llm, &
179  jj_nb_vanleer,0,0,request_vanleer)
180  call register_swapfieldhallo(pbarvg,pbarvg,ip1jm,llm, &
181  jj_nb_vanleer,1,0,request_vanleer)
182  call register_swapfieldhallo(massem,massem,ip1jmp1,llm, &
183  jj_nb_vanleer,0,0,request_vanleer)
184  call register_swapfieldhallo(wg,wg,ip1jmp1,llm, &
185  jj_nb_vanleer,0,0,request_vanleer)
186  call register_swapfieldhallo(teta_tmp,teta_tmp,ip1jmp1,llm, &
187  jj_nb_vanleer,1,1,request_vanleer)
188  call register_swapfieldhallo(p_tmp,p_tmp,ip1jmp1,llmp1, &
189  jj_nb_vanleer,1,1,request_vanleer)
190  call register_swapfieldhallo(pk_tmp,pk_tmp,ip1jmp1,llm, &
191  jj_nb_vanleer,1,1,request_vanleer)
192  do j=1,nqtot
193  call register_swapfieldhallo(q(1,1,j),q(1,1,j),ip1jmp1,llm, &
194  jj_nb_vanleer,0,0,request_vanleer)
195  enddo
196 
197  call sendrequest(request_vanleer)
198  !$OMP BARRIER
199  call waitrequest(request_vanleer)
200 
201 
202  !$OMP BARRIER
203  !$OMP MASTER
204  call setdistrib(jj_nb_vanleer)
205  call vte(vthallo)
206  call vtb(vtadvection)
207  call start_timer(timer_vanleer)
208  !$OMP END MASTER
209  !$OMP BARRIER
210 
211  ! ... Flux de masse diaganostiques traceurs
212  ijb=ij_begin
213  ije=ij_end
214  flxw(ijb:ije,1:llm)=wg(ijb:ije,1:llm)/REAL(iapp_tracvl)
215 
216  ! test sur l'eventuelle creation de valeurs negatives de la masse
217  ijb=ij_begin
218  ije=ij_end
219  if (pole_nord) ijb=ij_begin+iip1
220  if (pole_sud) ije=ij_end-iip1
221 
222  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
223  DO l=1,llm-1
224  DO ij = ijb+1,ije
225  zdp(ij) = pbarug(ij-1,l) - pbarug(ij,l) &
226  - pbarvg(ij-iip1,l) + pbarvg(ij,l) &
227  + wg(ij,l+1) - wg(ij,l)
228  ENDDO
229 
230  ! CALL SCOPY( jjm -1 ,zdp(iip1+iip1),iip1,zdp(iip2),iip1 )
231  ! ym ---> pourquoi jjm-1 et non jjm ? a cause du pole ?
232 
233  do ij=ijb,ije-iip1+1,iip1
234  zdp(ij)=zdp(ij+iip1-1)
235  enddo
236 
237  DO ij = ijb,ije
238  zdp(ij)= zdp(ij)*dtvr/ massem(ij,l)
239  ENDDO
240 
241 
242  ! CALL minmax ( ip1jm-iip1, zdp(iip2), zdpmin,zdpmax )
243  ! ym ---> eventuellement a revoir
244  CALL minmax ( ije-ijb+1, zdp(ijb), zdpmin,zdpmax )
245 
246  IF(max(abs(zdpmin),abs(zdpmax)).GT.0.5) THEN
247  print*,'WARNING DP/P l=',l,' MIN:',zdpmin, &
248  ' MAX:', zdpmax
249  ENDIF
250 
251  ENDDO
252  !$OMP END DO NOWAIT
253 
254  !-------------------------------------------------------------------
255  ! Advection proprement dite (Modification Le Croller (07/2001)
256  !-------------------------------------------------------------------
257 
258  !----------------------------------------------------
259  ! Calcul des moyennes basées sur la masse
260  !----------------------------------------------------
261 
262  !ym ----> Normalement, inutile pour les schémas classiques
263  !ym ----> Revérifier lors de la parallélisation des autres schemas
264 
265  !ym call massbar_p(massem,massebx,masseby)
266 
267  call vlspltgen_p( q,iadv, 2., massem, wg , &
268  pbarug,pbarvg,dtvr,p_tmp,pk_tmp,teta_tmp )
269 
270 
271  GOTO 1234
272  !-----------------------------------------------------------
273  ! Appel des sous programmes d'advection
274  !-----------------------------------------------------------
275  do iq=1,nqtot
276  ! call clock(t_initial)
277  if(iadv(iq) == 0) cycle
278  ! ----------------------------------------------------------------
279  ! Schema de Van Leer I MUSCL
280  ! ----------------------------------------------------------------
281  if(iadv(iq).eq.10) THEN
282 
283  call vlsplt_p(q(1,1,iq),2.,massem,wg,pbarug,pbarvg,dtvr)
284 
285  ! ----------------------------------------------------------------
286  ! Schema "pseudo amont" + test sur humidite specifique
287  ! pour la vapeur d'eau. F. Codron
288  ! ----------------------------------------------------------------
289  else if(iadv(iq).eq.14) then
290  !
291  !ym stop 'advtrac : appel à vlspltqs :schema non parallelise'
292  CALL vlspltqs_p( q(1,1,1), 2., massem, wg , &
293  pbarug,pbarvg,dtvr,p_tmp,pk_tmp,teta_tmp )
294  ! ----------------------------------------------------------------
295  ! Schema de Frederic Hourdin
296  ! ----------------------------------------------------------------
297  else if(iadv(iq).eq.12) then
298  stop 'advtrac : schema non parallelise'
299  ! Pas de temps adaptatif
300  call adaptdt(iadv(iq),dtbon,n,pbarug,massem)
301  if (n.GT.1) then
302  write(*,*) 'WARNING horizontal dt=',dtbon,'dtvr=', &
303  dtvr,'n=',n
304  endif
305  do indice=1,n
306  call advn(q(1,1,iq),massem,wg,pbarug,pbarvg,dtbon,1)
307  end do
308  else if(iadv(iq).eq.13) then
309  stop 'advtrac : schema non parallelise'
310  ! Pas de temps adaptatif
311  call adaptdt(iadv(iq),dtbon,n,pbarug,massem)
312  if (n.GT.1) then
313  write(*,*) 'WARNING horizontal dt=',dtbon,'dtvr=', &
314  dtvr,'n=',n
315  endif
316  do indice=1,n
317  call advn(q(1,1,iq),massem,wg,pbarug,pbarvg,dtbon,2)
318  end do
319  ! ----------------------------------------------------------------
320  ! Schema de pente SLOPES
321  ! ----------------------------------------------------------------
322  else if (iadv(iq).eq.20) then
323  stop 'advtrac : schema non parallelise'
324 
325  call pentes_ini (q(1,1,iq),wg,massem,pbarug,pbarvg,0)
326 
327  ! ----------------------------------------------------------------
328  ! Schema de Prather
329  ! ----------------------------------------------------------------
330  else if (iadv(iq).eq.30) then
331  stop 'advtrac : schema non parallelise'
332  ! Pas de temps adaptatif
333  call adaptdt(iadv(iq),dtbon,n,pbarug,massem)
334  if (n.GT.1) then
335  write(*,*) 'WARNING horizontal dt=',dtbon,'dtvr=', &
336  dtvr,'n=',n
337  endif
338  call prather(q(1,1,iq),wg,massem,pbarug,pbarvg, &
339  n,dtbon)
340  ! ----------------------------------------------------------------
341  ! Schemas PPM Lin et Rood
342  ! ----------------------------------------------------------------
343  else if (iadv(iq).eq.11.OR.(iadv(iq).GE.16.AND. &
344  iadv(iq).LE.18)) then
345 
346  stop 'advtrac : schema non parallelise'
347 
348  ! Test sur le flux horizontal
349  ! Pas de temps adaptatif
350  call adaptdt(iadv(iq),dtbon,n,pbarug,massem)
351  if (n.GT.1) then
352  write(*,*) 'WARNING horizontal dt=',dtbon,'dtvr=', &
353  dtvr,'n=',n
354  endif
355  ! Test sur le flux vertical
356  cflmaxz=0.
357  do l=2,llm
358  do ij=iip2,ip1jm
359  aaa=wg(ij,l)*dtvr/massem(ij,l)
360  cflmaxz=max(cflmaxz,aaa)
361  bbb=-wg(ij,l)*dtvr/massem(ij,l-1)
362  cflmaxz=max(cflmaxz,bbb)
363  enddo
364  enddo
365  if (cflmaxz.GE.1) then
366  write(*,*) 'WARNING vertical','CFLmaxz=', cflmaxz
367  endif
368 
369  !-----------------------------------------------------------
370  ! Ss-prg interface LMDZ.4->PPM3d
371  !-----------------------------------------------------------
372 
373  call interpre(q(1,1,iq),qppm(1,1,iq),wg,fluxwppm,massem, &
374  apppm,bpppm,massebx,masseby,pbarug,pbarvg, &
375  unatppm,vnatppm,psppm)
376 
377  do indice=1,n
378  !----------------------------------------------------------------
379  ! VL (version PPM) horiz. et PPM vert.
380  !----------------------------------------------------------------
381  if (iadv(iq).eq.11) then
382  ! Ss-prg PPM3d de Lin
383  call ppm3d(1,qppm(1,1,iq), &
384  psppm,psppm, &
385  unatppm,vnatppm,fluxwppm,dtbon,2,2,2,1, &
386  iim,jjp1,2,llm,apppm,bpppm,0.01,6400000, &
387  fill,dum,220.)
388 
389  !-------------------------------------------------------------
390  ! Monotonic PPM
391  !-------------------------------------------------------------
392  else if (iadv(iq).eq.16) then
393  ! Ss-prg PPM3d de Lin
394  call ppm3d(1,qppm(1,1,iq), &
395  psppm,psppm, &
396  unatppm,vnatppm,fluxwppm,dtbon,3,3,3,1, &
397  iim,jjp1,2,llm,apppm,bpppm,0.01,6400000, &
398  fill,dum,220.)
399  !-------------------------------------------------------------
400 
401  !-------------------------------------------------------------
402  ! Semi Monotonic PPM
403  !-------------------------------------------------------------
404  else if (iadv(iq).eq.17) then
405  ! Ss-prg PPM3d de Lin
406  call ppm3d(1,qppm(1,1,iq), &
407  psppm,psppm, &
408  unatppm,vnatppm,fluxwppm,dtbon,4,4,4,1, &
409  iim,jjp1,2,llm,apppm,bpppm,0.01,6400000, &
410  fill,dum,220.)
411  !-------------------------------------------------------------
412 
413  !-------------------------------------------------------------
414  ! Positive Definite PPM
415  !-------------------------------------------------------------
416  else if (iadv(iq).eq.18) then
417  ! Ss-prg PPM3d de Lin
418  call ppm3d(1,qppm(1,1,iq), &
419  psppm,psppm, &
420  unatppm,vnatppm,fluxwppm,dtbon,5,5,5,1, &
421  iim,jjp1,2,llm,apppm,bpppm,0.01,6400000, &
422  fill,dum,220.)
423  !-------------------------------------------------------------
424  endif
425  enddo
426  !-----------------------------------------------------------------
427  ! Ss-prg interface PPM3d-LMDZ.4
428  !-----------------------------------------------------------------
429  call interpost(q(1,1,iq),qppm(1,1,iq))
430  endif
431  !----------------------------------------------------------------------
432 
433  !-----------------------------------------------------------------
434  ! On impose une seule valeur du traceur au pôle Sud j=jjm+1=jjp1
435  ! et Nord j=1
436  !-----------------------------------------------------------------
437 
438  ! call traceurpole(q(1,1,iq),massem)
439 
440  ! calcul du temps cpu pour un schema donne
441 
442  ! call clock(t_final)
443  !ym tps_cpu=t_final-t_initial
444  !ym cpuadv(iq)=cpuadv(iq)+tps_cpu
445 
446  end DO
447 
448 1234 CONTINUE
449  !$OMP BARRIER
450 
451  if (planet_type=="earth") then
452 
453  ijb=ij_begin
454  ije=ij_end
455 
456  !$OMP DO SCHEDULE(STATIC,OMP_CHUNK)
457  DO l = 1, llm
458  DO ij = ijb, ije
459  finmasse(ij,l) = p(ij,l) - p(ij,l+1)
460  ENDDO
461  ENDDO
462  !$OMP END DO
463 
464  CALL qminimum_p( q, 2, finmasse )
465  endif ! of if (planet_type=="earth")
466 
467  !------------------------------------------------------------------
468  ! on reinitialise a zero les flux de masse cumules
469  !---------------------------------------------------
470  ! iadvtr=0
471 
472  !$OMP MASTER
473  call vte(vtadvection)
474  call stop_timer(timer_vanleer)
475  call vtb(vthallo)
476  !$OMP END MASTER
477 
478  do j=1,nqtot
479  call register_swapfieldhallo(q(1,1,j),q(1,1,j),ip1jmp1,llm, &
480  jj_nb_caldyn,0,0,request_vanleer)
481  enddo
482 
483  call register_swapfieldhallo(flxw,flxw,ip1jmp1,llm, &
484  jj_nb_caldyn,0,0,request_vanleer)
485 
486  call sendrequest(request_vanleer)
487  !$OMP BARRIER
488  call waitrequest(request_vanleer)
489 
490  !$OMP BARRIER
491  !$OMP MASTER
492  call setdistrib(jj_nb_caldyn)
493  call vte(vthallo)
494  call resume_timer(timer_caldyn)
495  !$OMP END MASTER
496  !$OMP BARRIER
497  iadvtr=0
498  ENDIF ! if iadvtr.EQ.iapp_tracvl
499 
500 END SUBROUTINE advtrac_p
bands
Definition: bands.F90:4
pentes_ini
subroutine pentes_ini(q, w, masse, pbaru, pbarv, mode)
Definition: pentes_ini.F:5
control_mod::iapp_tracvl
integer, save iapp_tracvl
Definition: control_mod.F90:17
iip2
!$Header iip2
Definition: paramet.h:14
ip1jmp1
!$Header llmm1 INTEGER ip1jmp1
Definition: paramet.h:14
interpost
subroutine interpost(q, qppm)
Definition: interpost.F:5
groupe_p
subroutine groupe_p(pext, pbaru, pbarv, pbarum, pbarvm, wm)
Definition: groupe_p.F:2
bands::jj_nb_caldyn
integer, dimension(:), allocatable jj_nb_caldyn
Definition: bands.F90:10
qminimum_p
subroutine qminimum_p(q, nq, deltap)
Definition: qminimum_p.F:2
times::stop_timer
subroutine stop_timer(no_timer)
Definition: times.F90:103
adaptdt
subroutine adaptdt(nadv, dtbon, n, pbaru, masse)
Definition: adaptdt.F:6
llmp1
!$Header llmp1
Definition: paramet.h:14
vampir
Definition: vampir.F90:1
advtrac_p
subroutine advtrac_p(pbaru, pbarv, p, masse, q, iapptrac, teta, flxw, pk)
Definition: advtrac_p.F90:4
parallel_lmdz::ij_end
integer, save ij_end
Definition: parallel_lmdz.F90:23
parallel_lmdz::pole_sud
logical, save pole_sud
Definition: parallel_lmdz.F90:25
vlspltqs_p
subroutine vlspltqs_p(q, pente_max, masse, w, pbaru, pbarv, pdt, p, pk, teta)
Definition: vlspltqs_p.F:6
vampir::vtb
subroutine vtb(number)
Definition: vampir.F90:52
control_mod::planet_type
character(len=10), save planet_type
Definition: control_mod.F90:32
vlsplt_p
subroutine vlsplt_p(q, pente_max, masse, w, pbaru, pbarv, pdt)
Definition: vlsplt_p.F:6
llm
!$Id Turb_fcg_gcssold get_uvd hqturb_gcssold endif!large scale llm day day1 day day1 *dt_toga endif!time annee_ref dt_toga u_toga vq_toga w_prof vq_prof llm day day1 day day1 *dt_dice endif!time annee_ref dt_dice swup_dice vg_dice omega_dice tg_prof vg_profd w_profd omega_profd!do llm!print llm l llm
Definition: 1D_interp_cases.h:103
vlspltgen_p
subroutine vlspltgen_p(q, iadv, pente_max, masse, w, pbaru, pbarv, pdt,
Definition: vlspltgen_p.F:5
advn
subroutine advn(q, masse, w, pbaru, pbarv, pdt, mode)
Definition: advn.F:5
control_mod::day_step
integer, save day_step
Definition: control_mod.F90:15
times::timer_caldyn
integer, parameter timer_caldyn
Definition: times.F90:7
times::resume_timer
subroutine resume_timer(no_timer)
Definition: times.F90:87
infotrac::nqtot
integer, save nqtot
Definition: infotrac.F90:6
ip1jm
!$Header llmm1 INTEGER ip1jm
Definition: paramet.h:14
write_field_p::writefield_p
Definition: write_field_p.F90:4
write_field_p
Definition: write_field_p.F90:1
mod_hallo::register_swapfieldhallo
subroutine register_swapfieldhallo(FieldS, FieldR, ij, ll, jj_Nb_New, Up, Down, a_request)
Definition: mod_hallo.F90:302
parallel_lmdz::pole_nord
logical, save pole_nord
Definition: parallel_lmdz.F90:24
parallel_lmdz
Definition: parallel_lmdz.F90:4
jjp1
!$Header jjp1
Definition: paramet.h:14
minmax
subroutine minmax(imax, xi, zmin, zmax)
Definition: minmax.F:5
vampir::vthallo
integer, parameter vthallo
Definition: vampir.F90:7
times
Definition: times.F90:1
ppm3d
subroutine ppm3d(IGD, Q, PS1, PS2, U, V, W, NDT, IORD, JORD, KORD, NC, IMR, JNP, j1, NLAY, AP, BP, PT, AE, fill, dum, Umax)
Definition: ppm3d.F:67
prather
subroutine prather(q, w, masse, pbaru, pbarv, nt, dt)
Definition: prather.F:5
mod_hallo::sendrequest
subroutine sendrequest(a_Request)
Definition: mod_hallo.F90:1072
true
!$Id itapm1 ENDIF!IM on interpole les champs sur les niveaux STD de pression!IM a chaque pas de temps de la physique c!positionnement de l argument logique a false c!pour ne pas recalculer deux fois la meme chose!c!a cet effet un appel a plevel_new a ete deplace c!a la fin de la serie d appels c!la boucle DO nlevSTD a ete internalisee c!dans d ou la creation de cette routine c c!CALL ulevSTD CALL &zphi philevSTD CALL &zx_rh rhlevSTD!DO klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon klev DO klon du jour ou toutes les read_climoz CALL true
Definition: calcul_STDlev.h:139
parallel_lmdz::ij_begin
integer, save ij_begin
Definition: parallel_lmdz.F90:22
parallel_lmdz::setdistrib
subroutine setdistrib(jj_Nb_New)
Definition: parallel_lmdz.F90:206
dtvr
!$Id mode_top_bound COMMON comconstr dtvr
Definition: comconst.h:7
parallel_lmdz::omp_chunk
integer, save omp_chunk
Definition: parallel_lmdz.F90:46
vampir::vte
subroutine vte(number)
Definition: vampir.F90:69
times::suspend_timer
subroutine suspend_timer(no_timer)
Definition: times.F90:70
vampir::vtadvection
integer, parameter vtadvection
Definition: vampir.F90:5
iim
c c zjulian c cym CALL iim cym klev iim
Definition: ini_bilKP_ave.h:24
times::start_timer
subroutine start_timer(no_timer)
Definition: times.F90:51
mod_hallo::request
Definition: mod_hallo.F90:38
times::timer_vanleer
integer, parameter timer_vanleer
Definition: times.F90:8
bands::jj_nb_vanleer
integer, dimension(:), allocatable jj_nb_vanleer
Definition: bands.F90:11
mod_hallo
Definition: mod_hallo.F90:1
mod_hallo::waitrequest
subroutine waitrequest(a_Request)
Definition: mod_hallo.F90:1196
control_mod
Definition: control_mod.F90:5
infotrac::iadv
integer, dimension(:), allocatable, save iadv
Definition: infotrac.F90:22