Changeset 22478

Timestamp:

02/26/18 22:35:18 (7 years ago)

Author:

adhikari

Message:

CHG: ordering issue regardg even yi resolved

File:

• issm/trunk-jpl/src/c/modules/FourierLoveCorex/lovenb_sub.f90 (modified) (2 diffs)

issm/trunk-jpl/src/c/modules/FourierLoveCorex/lovenb_sub.f90

@@ -378 +378 @@
    do j = 1,nbc/6
       !!! y4, y2, y6 are not scaled yet =>>> ask lambert ==>>> 
-      love_kernels((j-1)*6+1,1) = sc((j-1)*6+1,1)               ! y4 NOT SCALED YET 
-      love_kernels((j-1)*6+2,1) = sc((j-1)*6+2,1)               ! y2 NOT SCALED YET 
-      love_kernels((j-1)*6+3,1) = sc((j-1)*6+3,1)               ! y6 NOT SCALED YET 
+      !love_kernels((j-1)*6+1,1) = sc((j-1)*6+1,1)               ! y4 NOT SCALED YET 
+      !love_kernels((j-1)*6+2,1) = sc((j-1)*6+2,1)               ! y2 NOT SCALED YET 
+      !love_kernels((j-1)*6+3,1) = sc((j-1)*6+3,1)               ! y6 NOT SCALED YET 
       love_kernels((j-1)*6+4,1) = sc((j-1)*6+4,1)*ra*go_surf    ! y1 => love_h 
       love_kernels((j-1)*6+5,1) = sc((j-1)*6+5,1)*ra*go_surf    ! y3 => love_l 
@@ -388 +388 @@
       !end if 
    end do
-
+   ! some reshuffling is needed for y4, y2, y6...
+   do j = 2,nbc/6
+      !!! y4, y2, y6 are not scaled yet =>>> ask lambert ==>>> 
+      love_kernels((j-2)*6+1,1) = sc((j-1)*6+1,1)               ! y4 NOT SCALED YET 
+      love_kernels((j-2)*6+2,1) = sc((j-1)*6+2,1)               ! y2 NOT SCALED YET 
+      love_kernels((j-2)*6+3,1) = sc((j-1)*6+3,1)               ! y6 NOT SCALED YET 
+   end do
+   !!! these are surface values... 
+   love_kernels(nbc-5,1) = sc(1,1)               ! y4 NOT SCALED YET 
+   love_kernels(nbc-4,1) = sc(2,1)               ! y4 NOT SCALED YET 
+   love_kernels(nbc-3,1) = sc(3,1)               ! y4 NOT SCALED YET 
+ 
+   !print *, ' ************* ' 
    !if (deg==2) then
    !  do ibc=1,nbc/6
-   !      print *, real(sc( (ibc-1)*6+1:(ibc-1)*6+6 , 1))
-   !   end do 
-   !   print *, ' ************* ' 
-   !   do ibc=1,nbc/6
+         !print *, real(sc( (ibc-1)*6+1:(ibc-1)*6+6 , 1))
    !      print *, real(love_kernels( (ibc-1)*6+1:(ibc-1)*6+6 , 1))
-   !   end do 
-   !   print *, ' ************* ' 
-   !   print*, ra*go_surf, ra*go0, mu0
+    !  end do 
+      !print *, ' ************* ' 
+      !do ibc=1,nbc/6
+      !   print *, real(love_kernels( (ibc-1)*6+1:(ibc-1)*6+6 , 1))
+      !end do 
+      !print *, ' ************* ' 
+      !print*, ra*go_surf, ra*go0, mu0
    !end if 
    !!! >>> SA