Changeset 10087 for vendors/AGRIF/CMEMS_2020/AGRIF_FILES/modupdatebasic.F90
- Timestamp:
- 2018-09-05T15:33:44+02:00 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
vendors/AGRIF/CMEMS_2020/AGRIF_FILES/modupdatebasic.F90
r5656 r10087 49 49 integer, intent(in) :: np !< Length of parent array 50 50 integer, intent(in) :: nc !< Length of child array 51 real ,intent(in) :: s_parent !< Parent grid position (s_root = 0)52 real ,intent(in) :: s_child !< Child grid position (s_root = 0)53 real ,intent(in) :: ds_parent !< Parent grid dx (ds_root = 1)54 real ,intent(in) :: ds_child !< Child grid dx (ds_root = 1)51 real(kind=8), intent(in) :: s_parent !< Parent grid position (s_root = 0) 52 real(kind=8), intent(in) :: s_child !< Child grid position (s_root = 0) 53 real(kind=8), intent(in) :: ds_parent !< Parent grid dx (ds_root = 1) 54 real(kind=8), intent(in) :: ds_child !< Child grid dx (ds_root = 1) 55 55 !--------------------------------------------------------------------------------------------------- 56 56 integer :: i, locind_child_left, coeffraf … … 84 84 integer, intent(in) :: np !< Length of parent array 85 85 integer, intent(in) :: nc !< Length of child array 86 real ,intent(in) :: s_parent !< Parent grid position (s_root = 0)87 real ,intent(in) :: s_child !< Child grid position (s_root = 0)88 real ,intent(in) :: ds_parent !< Parent grid dx (ds_root = 1)89 real ,intent(in) :: ds_child !< Child grid dx (ds_root = 1)86 real(kind=8), intent(in) :: s_parent !< Parent grid position (s_root = 0) 87 real(kind=8), intent(in) :: s_child !< Child grid position (s_root = 0) 88 real(kind=8), intent(in) :: ds_parent !< Parent grid dx (ds_root = 1) 89 real(kind=8), intent(in) :: ds_child !< Child grid dx (ds_root = 1) 90 90 integer, intent(in) :: dir !< Direction 91 91 !--------------------------------------------------------------------------------------------------- … … 157 157 REAL, DIMENSION(nc), intent(in) :: y 158 158 INTEGER, intent(in) :: np,nc 159 REAL ,intent(in) :: s_parent, s_child160 REAL ,intent(in) :: ds_parent, ds_child159 REAL(kind=8), intent(in) :: s_parent, s_child 160 REAL(kind=8), intent(in) :: ds_parent, ds_child 161 161 ! 162 162 INTEGER :: i, ii, locind_child_left, coeffraf 163 REAL :: xpos, invcoeffraf 163 REAL(kind=8) :: xpos 164 REAL :: invcoeffraf 164 165 INTEGER :: nbnonnuls 165 166 INTEGER :: diffmod … … 229 230 !--------------------------------------------------------------------------------------------------- 230 231 INTEGER, intent(in) :: nc2, np, nc 231 REAL , intent(in) :: s_parent, s_child232 REAL , intent(in) :: ds_parent, ds_child232 REAL(kind=8), intent(in) :: s_parent, s_child 233 REAL(kind=8), intent(in) :: ds_parent, ds_child 233 234 INTEGER, intent(in) :: dir 234 235 ! 235 236 INTEGER, DIMENSION(:,:), ALLOCATABLE :: indchildaverage_tmp 236 237 INTEGER :: i, locind_child_left, coeffraf 237 REAL :: xpos238 REAL(kind=8) :: xpos 238 239 INTEGER :: diffmod 239 240 ! … … 281 282 REAL, DIMENSION(nc), intent(in) :: y 282 283 INTEGER, intent(in) :: np, nc 283 REAL , intent(in) :: s_parent, s_child284 REAL , intent(in) :: ds_parent, ds_child284 REAL(kind=8), intent(in) :: s_parent, s_child 285 REAL(kind=8), intent(in) :: ds_parent, ds_child 285 286 INTEGER, intent(in) :: dir 286 287 ! … … 311 312 ELSE 312 313 ! 313 !CDIR NOLOOPCHG 314 do j = 1,coeffraf 315 !CDIR VECTOR 316 do i= 1,np 314 315 do i = 1,np 316 do j = 1,coeffraf 317 317 x(i) = x(i) + y(indchildaverage(i,dir) + j-1 ) 318 318 enddo 319 319 enddo 320 320 IF (.not.Agrif_Update_Weights) THEN … … 338 338 real, dimension(nc), intent(in) :: y 339 339 integer, intent(in) :: np, nc 340 real , intent(in) :: s_parent, s_child341 real , intent(in) :: ds_parent, ds_child342 !--------------------------------------------------------------------------------------------------- 343 REAL :: xpos, xposfin340 real(kind=8), intent(in) :: s_parent, s_child 341 real(kind=8), intent(in) :: ds_parent, ds_child 342 !--------------------------------------------------------------------------------------------------- 343 REAL(kind=8) :: xpos, xposfin 344 344 INTEGER :: i, ii, diffmod 345 345 INTEGER :: it1, it2
Note: See TracChangeset
for help on using the changeset viewer.