subroutine maxmin(u,char,ktau,fact,il,kl) ! has more general format & scaling factor ! has entry maxmin0 for single level fields !include 'newmpar.h' integer il,jl,kl,ifull character*2 char dimension u(il,6*il,kl) & ,umin(kl),umax(kl),iumax(kl),jumax(kl),iumin(kl),jumin(kl) jl=6*il ifull=il*jl kup=kl 2 do k=1,kup umax(k)=u(1,1,k) umin(k)=u(1,1,k) jumax(k)=1 iumax(k)=1 jumin(k)=1 iumin(k)=1 do j=1,jl do i=1,il if(umax(k).lt.u(i,j,k))then umax(k)=u(i,j,k) jumax(k)=j iumax(k)=i endif if(umin(k).gt.u(i,j,k))then c print *,'k,i,j,iu,ju,umin,u: ', c . k,i,j,iumin(k),jumin(k),umin(k),u(i,j,k) umin(k)=u(i,j,k) jumin(k)=j iumin(k)=i endif enddo ! i loop enddo ! j loop umax(k)=fact*umax(k) umin(k)=fact*umin(k) enddo ! k loop if(umax(kup).gt.30.)then ! format for T, & usually u,v print 971,ktau,char,(umax(k),k=1,kup) 971 format(i7,1x,a2,'max ',10f7.2/(14x,10f7.2)/(14x,10f7.2)) print 977,ktau,(iumax(k),jumax(k),k=1,kup) print 972,ktau,char,(umin(k),k=1,kup) 972 format(i7,1x,a2,'min ',10f7.2/(14x,10f7.2)/(14x,10f7.2)) print 977,ktau,(iumin(k),jumin(k),k=1,kup) 977 format(i7,' posij',10(i3,i4)/(14x,10(i3,i4))/(14x,10(i3,i4))) else ! for qg & sd print 981,ktau,char,(umax(k),k=1,kup) 981 format(i7,1x,a2,'max ',10f7.3/(14x,10f7.3)/(14x,10f7.3)) print 977,ktau,(iumax(k),jumax(k),k=1,kup) print 982,ktau,char,(umin(k),k=1,kup) 982 format(i7,1x,a2,'min ',10f7.3/(14x,10f7.3)/(14x,10f7.3)) print 977,ktau,(iumin(k),jumin(k),k=1,kup) endif return entry maxmin0(u,char,ktau,fact,il,kl) kup=1 jl=6*il ifull=il*jl go to 2 end