Changeset 105

Timestamp:

06/11/13 14:28:57 (11 years ago)

Author:

leila_ocean

Message:

weights

File:

• altifloat/src/lapfilu.h (modified) (2 diffs)

altifloat/src/lapfilu.h

@@ -13 +13 @@
       //Example: a weighted mean
       YS1=c1*fij + c2*fim1j + c3*fip1j + c4*fijm1 + c5*fijp1;
+   }
+}
 
-    }
-}
 
 backward (YREAL xu_d, YREAL fij, YREAL fim1j, YREAL fip1j, YREAL fijm1, YREAL fijp1)
@@ -38 +38 @@
 
 }
+
+
+
+//LLL START
+forward (YREAL xu_d, YREAL fij, YREAL fim1j, YREAL fip1j, YREAL fijm1, YREAL fijp1, YREAL wij, YREAL wim1j, YREAL wip1j, YREAL wijm1, YREAL wijp1)
+{
+ 
+  if (Yt==1) { //first step time
+    YS1=xu_d;
+  }
+  else
+    {
+      /* Boundary condition ? 
+         if(Yi==0 || Yi==nlon-1 || Yj==0 || Yj==nlat-1)
+      */
+
+//LLL The filter is actually L = I + kp Dt D, where D is the Laplacian Operator, I is the identity matrix, kp is the diffusion coefficient, Dt is the time step. kp*Dt is included in the weight matrix w among possibly other things to be determined soon
+/*      if(Yi==Yj) {
+        YS1 = 1.0 + c1*fij*wij + c2*fim1j*wim1j + c3*fip1j*wip1j + c4*fijm1*wijm1 + c5*fijp1*wijp1;
+    }
+    else
+    {
+        YS1 = c1*fij*wij + c2*fim1j*wim1j + c3*fip1j*wip1j + c4*fijm1*wijm1 + c5*fijp1*wijp1;   
+    }
+   }
+}
+
+
+backward (YREAL xu_d, YREAL fij, YREAL fim1j, YREAL fip1j, YREAL fijm1, YREAL fijp1, YREAL wij, YREAL wim1j, YREAL wip1j, YREAL wijm1, YREAL wijp1)
+{
+  
+  if (Yi==9 && Yj==16)
+    {
+      printf("YG=%g\n",YG1_u_d(Yi,Yj));
+    }
+  
+if (Yt==1) { //first step time
+    YJ1I1=1;
+  } 
+  else
+    {
+      YJ1I2=c1*wij;
+      YJ1I3=c2*wim1j;
+      YJ1I4=c3*wip1j;
+      YJ1I5=c4*wijm1;
+      YJ1I6=c5*wijp1;
+      YJ1I7=c1*fij;
+      YJ1I8=c2*fim1j;
+      YJ1I9=c3*fip1j;
+      YJ1I10=c4*fijm1;
+      YJ1I11=c5*fijp1;
+    }
+}
+*/
+//LLL end