Changeset 5116 for LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigrads.f90
- Timestamp:
- Jul 24, 2024, 2:54:37 PM (2 months ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
LMDZ6/branches/Amaury_dev/libf/dyn3d_common/inigrads.f90
r5113 r5116 7 7 IMPLICIT NONE 8 8 9 integer:: if, im, jm, lm, i, j, l10 real:: x(im), y(jm), z(lm), fx, fy, fz, dt11 real:: xmin, xmax, ymin, ymax9 INTEGER :: if, im, jm, lm, i, j, l 10 REAL :: x(im), y(jm), z(lm), fx, fy, fz, dt 11 REAL :: xmin, xmax, ymin, ymax 12 12 13 character(len= *), intent(in) :: file14 character(len= *), intent(in) :: titlel13 CHARACTER(LEN = *), intent(in) :: file 14 CHARACTER(LEN = *), intent(in) :: titlel 15 15 16 16 INCLUDE "gradsdef.h" 17 17 18 18 ! data unit/66,32,34,36,38,40,42,44,46,48/ 19 integer:: nf19 INTEGER :: nf 20 20 save nf 21 21 data nf/0/ … … 49 49 do i = 1, im 50 50 xd(i, if) = x(i) * fx 51 if(xd(i, if)<xmin) iid(if) = i + 152 if(xd(i, if)<=xmax) ifd(if) = i51 IF(xd(i, if)<xmin) iid(if) = i + 1 52 IF(xd(i, if)<=xmax) ifd(if) = i 53 53 enddo 54 54 PRINT*, 'On stoke du point ', iid(if), ' a ', ifd(if), ' en x' … … 59 59 do j = 1, jm 60 60 yd(j, if) = y(j) * fy 61 if(yd(j, if)>ymax) jid(if) = j + 162 if(yd(j, if)>=ymin) jfd(if) = j61 IF(yd(j, if)>ymax) jid(if) = j + 1 62 IF(yd(j, if)>=ymin) jfd(if) = j 63 63 enddo 64 64 PRINT*, 'On stoke du point ', jid(if), ' a ', jfd(if), ' en y' … … 89 89 90 90 91 end subroutineinigrads91 END SUBROUTINE inigrads
Note: See TracChangeset
for help on using the changeset viewer.