2 SUBROUTINE ana_diag (ng, tile, model)
23 integer,
intent(in) :: ng, tile, model
27 character (len=*),
parameter :: MyFile = &
32 CALL ana_diag_tile (ng, tile, model, &
33 & lbi, ubi, lbj, ubj, &
34 & imins, imaxs, jmins, jmaxs, &
56 SUBROUTINE ana_diag_tile (ng, tile, model, &
57 & LBi, UBi, LBj, UBj, &
58 & IminS, ImaxS, JminS, JmaxS, &
74 integer,
intent(in) :: ng, tile, model
75 integer,
intent(in) :: LBi, UBi, LBj, UBj
76 integer,
intent(in) :: IminS, ImaxS, JminS, JmaxS
80 real(r8),
intent(in) :: u(LBi:,LBj:,:,:)
81 real(r8),
intent(in) :: v(LBi:,LBj:,:,:)
83 real(r8),
intent(in) :: ubar(LBi:,LBj:,:)
84 real(r8),
intent(in) :: vbar(LBi:,LBj:,:)
87 real(r8),
intent(in) :: u(LBi:UBi,LBj:UBj,N(ng),2)
88 real(r8),
intent(in) :: v(LBi:UBi,LBj:UBj,N(ng),2)
90 real(r8),
intent(in) :: ubar(LBi:UBi,LBj:UBj,:)
91 real(r8),
intent(in) :: vbar(LBi:UBi,LBj:UBj,:)
96 integer :: i, io_error, j, k
98 real(r8) :: umax, ubarmax, vmax, vbarmax
100 character (len=256) :: io_errmsg
102#include "set_bounds.h"
113 OPEN (
usrout,file=
usrname,form=
'formatted',status=
'unknown', &
114 & iostat=io_err, iomsg=io_errmsg)
115 IF (io_err.ne.0)
THEN
119 10
FORMAT (
' ANA_DIAG - unable to open output file: ',a, &
133 umax=max(umax,u(i,j,k,
nnew(ng)))
138 vmax=max(vmax,v(i,j,k,
nnew(ng)))
144 ubarmax=max(ubarmax,ubar(i,j,
knew(ng)))
149 vbarmax=max(vbarmax,vbar(i,j,
knew(ng)))
155 WRITE (
usrout,20)
tdays(ng), ubarmax, vbarmax, umax, vmax
156 20
FORMAT (2x,f13.6,2x,1pe13.6,2x,1pe13.6,2x,1pe13.6,2x,1pe13.6)
subroutine ana_diag_tile(ng, tile, model, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs, u, v, ubar, vbar)