- Timestamp:
- Jul 23, 2024, 5:57:06 PM (2 months ago)
- File:
-
- 1 moved
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/phylmd/Dust/minmaxsource.f90
r5103 r5104 1 SUBROUTINE minmaxsource(zq,qmin,qmax,comment)1 SUBROUTINE minmaxsource(zq, qmin, qmax, comment) 2 2 3 4 3 USE dimphy 4 USE infotrac 5 5 6 6 INCLUDE "dimensions.h" 7 7 8 !character*20 comment9 character*(*)comment10 real qmin,qmax11 real zq(klon,nbtr)8 ! character*20 comment 9 character(len = *) :: comment 10 real :: qmin, qmax 11 real :: zq(klon, nbtr) 12 12 13 integerijmin, lmin, ijlmin14 integerijmax, lmax, ijlmax13 integer :: ijmin, lmin, ijlmin 14 integer :: ijmax, lmax, ijlmax 15 15 16 integer ismin,ismax16 integer :: ismin, ismax 17 17 18 ijlmin=ismin(klon*nbtr,zq,1)19 lmin=(ijlmin-1)/klon+120 ijmin=ijlmin-(lmin-1)*klon21 zqmin=zq(ijmin,lmin)18 ijlmin = ismin(klon * nbtr, zq, 1) 19 lmin = (ijlmin - 1) / klon + 1 20 ijmin = ijlmin - (lmin - 1) * klon 21 zqmin = zq(ijmin, lmin) 22 22 23 ijlmax=ismax(klon*nbtr,zq,1)24 lmax=(ijlmax-1)/klon+125 ijmax=ijlmax-(lmax-1)*klon26 zqmax=zq(ijmax,lmax)23 ijlmax = ismax(klon * nbtr, zq, 1) 24 lmax = (ijlmax - 1) / klon + 1 25 ijmax = ijlmax - (lmax - 1) * klon 26 zqmax = zq(ijmax, lmax) 27 27 28 if(zqmin<qmin.or.zqmax>qmax)29 s write(*,9999) comment,30 s ijmin,lmin,zqmin,ijmax,lmax,zqmax28 if(zqmin<qmin.or.zqmax>qmax) & 29 write(*, 9999) comment, & 30 ijmin, lmin, zqmin, ijmax, lmax, zqmax 31 31 32 33 9999 format(a20,2(' q(',i4,',',i2,')=',e12.5))34 end 32 return 33 9999 format(a20, 2(' q(', i4, ',', i2, ')=', e12.5)) 34 end subroutine minmaxsource
Note: See TracChangeset
for help on using the changeset viewer.