Changeset 2044 in Sophya for trunk/SophyaProg/Tests/sph2lm.cc

Timestamp:

Jun 5, 2002, 1:58:01 PM (23 years ago)

Author:

ansari

Message:

Correction fichier Makefile et adaptation sph2lm aux modifs de GLM - Reza 5/6/2002

File:

• trunk/SophyaProg/Tests/sph2lm.cc (modified) (3 diffs)

trunk/SophyaProg/Tests/sph2lm.cc

@@ -376 +376 @@
   }
   int filw = 0;
-  if (narg > 4) filw = atoi(arg[3]);
+  if (narg > 4) filw = atoi(arg[4]);
 
   cout << "\n >>>> Starting sph2lm: arg[1]= " << arg[1] 
@@ -432 +432 @@
          << " Longitude= " << longitude_map[nlm] << endl;
 
-    LocalMap<r_8> lm(nx_map[nlm], ny_map[nlm]);
+    LocalMap<r_8> lm(nx_map[nlm], ny_map[nlm],
+                     nx_map[nlm]*resolx_map[nlm]/60.,
+                     ny_map[nlm]*resoly_map[nlm]/60.,
+                     thetadeg_map[nlm], longitude_map[nlm],
+                     nx_map[nlm]/2, ny_map[nlm]/2);
     TMatrix<r_8> mtx(nx_map[nlm], ny_map[nlm]);
 
-    lm.SetOrigin(thetadeg_map[nlm], longitude_map[nlm]);
-    lm.SetSize(nx_map[nlm]*resolx_map[nlm]/60.,
-               ny_map[nlm]*resoly_map[nlm]/60.);
+    //    lm.SetOrigin(thetadeg_map[nlm], longitude_map[nlm]);
+    //    lm.SetSize(nx_map[nlm]*resolx_map[nlm]/60.,
+    //       ny_map[nlm]*resoly_map[nlm]/60.);
   
     lm.SetPixels((nlm+1)*100.);
@@ -447 +451 @@
     
     cout << " Doing lm.Project(sphck) ... " << endl;
-    lm.Project(sphck);
+    lm.ProjectionToSphere(sphck);
     cout << " Doing mtx.Project(sphck) ... " << endl;
     ProjectMatrix2Sph(mtx, sphckmtx, userot);