Changeset 4179 for trunk/LMDZ.GENERIC/libf/phygeneric/util_linear_interpolation.F90

Timestamp:

Apr 10, 2026, 5:31:47 PM (33 hours ago)

Author:

emillour

Message:

Generic PCM:
Clean up code: put inter_line.F in a module util_linear_interpolation.F90 and
also turn soil.F and soil_settings.F into modules.
EM

File:

• trunk/LMDZ.GENERIC/libf/phygeneric/util_linear_interpolation.F90 (moved) (moved from trunk/LMDZ.GENERIC/libf/phygeneric/interp_line.F) (3 diffs)

trunk/LMDZ.GENERIC/libf/phygeneric/util_linear_interpolation.F90

@@ -1 @@
-      subroutine interp_line(x1,y1,len1,x2,y2,len2)
+module util_linear_interpolation_mod
+
+implicit none
+
+contains
+
+    subroutine interp_line(x1,y1,len1,x2,y2,len2)
       implicit none
 !-----------------------------------------------------------------------
@@ -10 +16 @@
 !  the nearby x1(:) point
 !  
-c-----------------------------------------------------------------------
+!-----------------------------------------------------------------------
 !  arguments
 !  ---------
 !  inputs:
-      real x1(len1) ! ordered list of abscissas
-      real y1(len1) ! values at x1(:)
-      integer len1  ! length of x1(:) and y1(:)
-      real x2(len2) !ordered list of abscissas at which interpolation is done
-      integer len2  ! length of x2(:) and y2(:)
+      real,intent(in) :: x1(len1) ! ordered list of abscissas
+      real,intent(in) :: y1(len1) ! values at x1(:)
+      integer,intent(in) :: len1  ! length of x1(:) and y1(:)
+      real,intent(in) :: x2(len2) ! ordered list of abscissas at which interpolation is done
+      integer,intent(in) :: len2  ! length of x2(:) and y2(:)
 !  outputs:
-      real y2(len2) ! interpolated values
+      real,intent(out) :: y2(len2) ! interpolated values
 !-----------------------------------------------------------------------
 
 ! local variables:
-      integer i,j
+      integer :: i,j
       
 
       do i=1,len2
         ! check if x2(i) lies outside of the interval covered by x1()
-        if(((x2(i).le.x1(1)).and.(x2(i).le.x1(len1))).or.
-     &     ((x2(i).ge.x1(1)).and.(x2(i).ge.x1(len1)))) then
+        if(((x2(i).le.x1(1)).and.(x2(i).le.x1(len1))).or. &
+           ((x2(i).ge.x1(1)).and.(x2(i).ge.x1(len1)))) then
           ! set y2(i) to y1(1) or y1(len1)
           if (abs(x2(i)-x1(1)).lt.abs(x2(i)-x1(len1))) then
@@ -43 +49 @@
         ! find the nearest neigbours and do a linear interpolation
          do j=1,len1-1
-          if(((x2(i).ge.x1(j)).and.(x2(i).le.x1(j+1))).or.
-     &       ((x2(i).le.x1(j)).and.(x2(i).ge.x1(j+1)))) then
-            y2(i)=((x2(i)-x1(j))/(x1(j+1)-x1(j)))*y1(j+1)+
-     &            ((x2(i)-x1(j+1))/(x1(j)-x1(j+1)))*y1(j)
+          if(((x2(i).ge.x1(j)).and.(x2(i).le.x1(j+1))).or. &
+             ((x2(i).le.x1(j)).and.(x2(i).ge.x1(j+1)))) then
+            y2(i)=((x2(i)-x1(j))/(x1(j+1)-x1(j)))*y1(j+1)+ &
+                  ((x2(i)-x1(j+1))/(x1(j)-x1(j+1)))*y1(j)
           endif
          enddo
         endif
 
-      enddo
+      enddo ! of do i=1,len2
 
-      end
+    end subroutine interp_line
+
+end module util_linear_interpolation_mod