9 #include "dimensions.h"
21 c re-inversion des niveaux
22 c le programme
ppm3d travaille avec une 3è e coordonnée inversée par rapport
c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
c On passe donc des niveaux de Lin à ceux du LMDZ
do l=1,llm
do j=1,jjp1
do i=1,iim
q(i,j,l)=qppm(i,j,llm-l+1)
enddo
enddo
enddo
c BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D
do l=1,llm
do j=1,jjp1
q(iip1,j,l)=q(1,j,l)
enddo
enddo
return
end
me coordonné inversée par rapport
c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
c On passe donc des niveaux de Lin à ceux du LMDZ
do l=1,llm
do j=1,jjp1
do i=1,iim
q(i,j,l)=qppm(i,j,llm-l+1)
enddo
enddo
enddo
c BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D
do l=1,llm
do j=1,jjp1
q(iip1,j,l)=q(1,j,l)
enddo
enddo
return
end
e inversé par rapport
c de celle du LMDZ: z=1<=>niveau max, z=llm+1<=>surface
c On passe donc des niveaux de Lin à ceux du LMDZ
do l=1,llm
do j=1,jjp1
do i=1,iim
q(i,j,l)=qppm(i,j,llm-l+1)
enddo
enddo
enddo
c BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D
do l=1,llm
do j=1,jjp1
q(iip1,j,l)=q(1,j,l)
enddo
enddo
return
end
e par rapport
23 c
de celle du lmdz: z=1<=>niveau max, z=
llm+1<=>surface
24 c on passe donc des niveaux
de lin à ceux du LMDZ
do l=1,llm
do j=1,jjp1
do i=1,iim
q(i,j,l)=qppm(i,j,llm-l+1)
enddo
enddo
enddo
c BOUCLAGE EN LONGITUDE PAS EFFECTUE DANS PPM3D
do l=1,llm
do j=1,jjp1
q(iip1,j,l)=q(1,j,l)
enddo
enddo
return
end
ceux du lmdz
29 q(i,j,l)=qppm(i,j,
llm-l+1)
34 c bouclage en longitude pas effectue dans
ppm3d
subroutine interpost(q, qppm)
!$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
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)
c c zjulian c cym CALL iim cym klev iim