*=*=*=*= test_period.html =*=*=*=*
SUBROUTINE test_period ( ucov, vcov, h, q, pext, phis ) c c Auteur : P. Le Van c --------- c .... Cette routine teste la periodicite en longitude des champs ucov, c h, q , pext et phis .......... c c IMPLICIT NONE c #include "dimensions.h" #include "paramet.h" c c ...... Arguments ...... c REAL ucov(ip1jmp1,llm), vcov(ip1jm,llm), h(ip1jmp1,llm) , , q(ip1jmp1,llm,nqmx), pext(ip1jmp1), phis(ip1jmp1) c c ..... Variables locales ..... c INTEGER ij,l,nq c DO l = 1, llm DO ij = 1, ip1jmp1, iip1 IF( ucov(ij,l).NE.ucov(ij+iim,l) ) THEN PRINT *,'STOP dans test_period car --- UCOV --- n est pas', , ' periodique en longitude ! ' PRINT *,' l, ij = ', l, ij, ij+iim STOP ENDIF IF( h(ij,l).NE.h(ij+iim,l) ) THEN PRINT *,'STOP dans test_period car --- H --- n est pas', , ' periodique en longitude ! ' PRINT *,' l, ij = ', l, ij, ij+iim STOP ENDIF ENDDO ENDDO c DO l = 1, llm DO ij = 1, ip1jm, iip1 IF( vcov(ij,l).NE.vcov(ij+iim,l) ) THEN PRINT *,'STOP dans test_period car --- VCOV --- n est pas', , ' periodique en longitude !' PRINT *,' l, ij = ', l, ij, ij+iim STOP ENDIF ENDDO ENDDO c DO nq =1, nqmx DO l =1, llm DO ij = 1, ip1jmp1, iip1 IF( q(ij,l,nq).NE.q(ij+iim,l,nq) ) THEN PRINT *,'STOP dans test_period car --- Q --- n est pas ', , 'periodique en longitude !' PRINT *,' nq , l, ij = ', nq, l, ij, ij+iim STOP ENDIF ENDDO ENDDO ENDDO c DO ij = 1, ip1jmp1, iip1 IF( pext(ij).NE.pext(ij+iim) ) THEN PRINT *,'STOP dans test_period car --- PEXT --- n est pas', , ' periodique en longitude !' PRINT *,' ij = ', ij, ij+iim STOP ENDIF IF( phis(ij).NE.phis(ij+iim) ) THEN PRINT *,'STOP dans test_period car --- PHIS --- n est pas', , ' periodique en longitude ! l, IJ = ', l, ij,ij+iim PRINT *,' ij = ', ij, ij+iim STOP ENDIF ENDDO c c RETURN END