- 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/minmaxqfi2.f90
r5103 r5104 1 SUBROUTINE minmaxqfi2(zq,qmin,qmax,comment)2 c 3 4 5 1 SUBROUTINE minmaxqfi2(zq, qmin, qmax, comment) 2 ! 3 USE dimphy 4 USE infotrac 5 INCLUDE "dimensions.h" 6 6 7 !character*20 comment8 character*(*)comment9 real qmin,qmax10 real zq(klon,klev)7 ! character*20 comment 8 character(len = *) :: comment 9 real :: qmin, qmax 10 real :: zq(klon, klev) 11 11 12 integerijmin, lmin, ijlmin13 integerijmax, lmax, ijlmax12 integer :: ijmin, lmin, ijlmin 13 integer :: ijmax, lmax, ijlmax 14 14 15 integer ismin,ismax15 integer :: ismin, ismax 16 16 17 ijlmin=ismin(klon*klev,zq,1)18 lmin=(ijlmin-1)/klon+119 ijmin=ijlmin-(lmin-1)*klon20 zqmin=zq(ijmin,lmin)17 ijlmin = ismin(klon * klev, zq, 1) 18 lmin = (ijlmin - 1) / klon + 1 19 ijmin = ijlmin - (lmin - 1) * klon 20 zqmin = zq(ijmin, lmin) 21 21 22 ijlmax=ismax(klon*klev,zq,1) 23 lmax=(ijlmax-1)/klon+1 24 ijmax=ijlmax-(lmax-1)*klon 25 zqmax=zq(ijmax,lmax) 26 27 if(zqmin<qmin.or.zqmax>qmax) 28 s write(*,9999) comment, 29 s ijmin,lmin,zqmin,ijmax,lmax,zqmax 22 ijlmax = ismax(klon * klev, zq, 1) 23 lmax = (ijlmax - 1) / klon + 1 24 ijmax = ijlmax - (lmax - 1) * klon 25 zqmax = zq(ijmax, lmax) 30 26 31 return 32 9999 format(a20,2(' q(',i4,',',i2,')=',e12.5)) 33 end 27 if(zqmin<qmin.or.zqmax>qmax) & 28 write(*, 9999) comment, & 29 ijmin, lmin, zqmin, ijmax, lmax, zqmax 30 31 return 32 9999 format(a20, 2(' q(', i4, ',', i2, ')=', e12.5)) 33 end subroutine minmaxqfi2
Note: See TracChangeset
for help on using the changeset viewer.