/[lmdze]/trunk/libf/dyn3d/prather.f
ViewVC logotype

Diff of /trunk/libf/dyn3d/prather.f

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 3 by guez, Wed Feb 27 13:16:39 2008 UTC revision 39 by guez, Tue Jan 25 15:11:05 2011 UTC
# Line 7  Line 7 
7        use comconst        use comconst
8        use comvert        use comvert
9        use comgeom        use comgeom
10          USE nr_util, ONLY : pi
11        IMPLICIT NONE        IMPLICIT NONE
12    
13  c=======================================================================  c=======================================================================
# Line 27  c======================================= Line 28  c=======================================
28  c   Arguments:  c   Arguments:
29  c   ----------  c   ----------
30        INTEGER iq,nt        INTEGER iq,nt
31        REAL pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm )        REAL, intent(in):: pbaru( ip1jmp1,llm ),pbarv( ip1jm,llm )
32        REAL masse(iip1,jjp1,llm)        REAL masse(iip1,jjp1,llm)
33        REAL q( iip1,jjp1,llm,0:9)        REAL q( iip1,jjp1,llm,0:9)
34        REAL w( ip1jmp1,llm )        REAL w( ip1jmp1,llm )
# Line 248  c   filtrages aux poles Line 249  c   filtrages aux poles
249            q( i,1,llm+1-l,0)=qpn            q( i,1,llm+1-l,0)=qpn
250            q( i,jjp1,llm+1-l,0)=qps            q( i,jjp1,llm+1-l,0)=qps
251           enddo           enddo
 c       enddo  
 c         print*,'qpn',qpn,'qps',qps  
 c          print*,'dqzpn',dqzpn,'dqzps',dqzps  
 c       enddo  
252             dyn1=0.             dyn1=0.
253             dys1=0.             dys1=0.
254             dyn2=0.             dyn2=0.
# Line 336  c   bouclage en longitude Line 333  c   bouclage en longitude
333           PRINT*,'SY(',i,j,l,')=',q(i,j,l,2),           PRINT*,'SY(',i,j,l,')=',q(i,j,l,2),
334       $   q(i,j-1,l,2)         $   q(i,j-1,l,2)  
335           PRINT*,'SZ(',i,j,l,')=',q(i,j,l,3)           PRINT*,'SZ(',i,j,l,')=',q(i,j,l,3)
 c                    PRINT*,' PBL EN SORTIE D'' ADVZP'  
336                       q(i,j,l,0)=0.                       q(i,j,l,0)=0.
 c                  STOP  
337                 ENDIF                 ENDIF
338             ENDDO             ENDDO
339           ENDDO           ENDDO
Legend:
Removed from v.3  
changed lines
  Added in v.39
  ViewVC Help
Powered by ViewVC 1.1.21