23 #include<file_version.h> 51 subroutine mpp_array_global_min_max(in_array, tmask,isd,jsd,isc,iec,jsc,jec,nk, g_min, g_max, &
52 geo_x,geo_y,geo_z, xgmin, ygmin, zgmin, xgmax, ygmax, zgmax)
56 real,
dimension(isd:,jsd:,:),
intent(in) :: in_array
57 real,
dimension(isd:,jsd:,:),
intent(in) :: tmask
58 integer,
intent(in) :: isd,jsd,isc,iec,jsc,jec,nk
59 real,
intent(out):: g_min, g_max
60 real,
dimension(isd:,jsd:),
intent(in) :: geo_x,geo_y
61 real,
dimension(:),
intent(in) :: geo_z
62 real,
intent(out):: xgmin, ygmin, zgmin, xgmax, ygmax, zgmax
66 real :: tmax, tmin, tmax0, tmin0
67 integer :: itmax, jtmax, ktmax, itmin, jtmin, ktmin
68 integer :: igmax, jgmax, kgmax, igmin, jgmin, kgmin
72 integer :: iminarr(3),imaxarr(3)
74 g_min=-88888888888.0 ; g_max=-999999999.0
76 tmax=-1.e10;tmin=1.e10
77 itmax=0;jtmax=0;ktmax=0
78 itmin=0;jtmin=0;ktmin=0
80 if(any(tmask(isc:iec,jsc:jec,:) > 0.))
then 81 iminarr=minloc(in_array(isc:iec,jsc:jec,:),tmask(isc:iec,jsc:jec,:) > 0.)
82 imaxarr=maxloc(in_array(isc:iec,jsc:jec,:),tmask(isc:iec,jsc:jec,:) > 0.)
83 itmin=iminarr(1)+isc-1
84 jtmin=iminarr(2)+jsc-1
86 itmax=imaxarr(1)+isc-1
87 jtmax=imaxarr(2)+jsc-1
89 tmin=in_array(itmin,jtmin,ktmin)
90 tmax=in_array(itmax,jtmax,ktmax)
94 fudge = 1.0 + 1.e-12*
real(mpp_pe() )
98 tmax = tmax + 1.e-12*
real(mpp_pe() )
101 tmin = tmin + 1.e-12*
real(mpp_pe() )
105 tmax0=tmax;tmin0=tmin
122 xgmax=-999.; ygmax=-999.; zgmax=-999.
123 xgmin=-999.; ygmin=-999.; zgmin=-999.
127 if (tmax0 == tmax)
then 128 xgmax=geo_x(itmax,jtmax)
129 ygmax=geo_y(itmax,jtmax)
137 if (tmin0 == tmin)
then 138 xgmin=geo_x(itmin,jtmin)
139 ygmin=geo_y(itmin,jtmin)
subroutine, public mpp_array_global_min_max(in_array, tmask, isd, jsd, isc, iec, jsc, jec, nk, g_min, g_max, geo_x, geo_y, geo_z, xgmin, ygmin, zgmin, xgmax, ygmax, zgmax)