Changeset 28097

Timestamp:

02/15/24 05:21:30 (13 months ago)

Author:

Mathieu Morlighem

Message:

BUG: fixing library problems on totten

File:

• issm/trunk-jpl/src/m/plot/radarpower.m (modified) (8 diffs)

issm/trunk-jpl/src/m/plot/radarpower.m

@@ -11 +11 @@
 
 %Parse inputs
-if nargin==1,
+if nargin==1
         options=pairoptions;
 else
@@ -20 +20 @@
 end
 
-highres = getfieldvalue(options,'highres',0);
 xlim    = getfieldvalue(options,'xlim',[min(md.mesh.x) max(md.mesh.x)]);
 ylim    = getfieldvalue(options,'ylim',[min(md.mesh.y) max(md.mesh.y)]);
@@ -28 +27 @@
 c = getfieldvalue(options,'overlay_adjust_c',0);
 d = getfieldvalue(options,'overlay_adjust_d',1);
+if diff(xlim)<1000e3 && diff(ylim)<1000e3
+        highres = getfieldvalue(options,'highres',1);
+else
+        highres = getfieldvalue(options,'highres',0);
+end
 
 %find GDAL coordinates
@@ -33 +37 @@
 y0=min(ylim); y1=max(ylim);
 
-if ~exist(options,'overlay_image'), % no image provided, go look into ModelData for one!{{{
+%Fix libraries if need be
+if strcmp(oshostname(),'totten')
+        setenv('LD_LIBRARY_PATH', '/usr/lib/x86_64-linux-gnu/libgeos_c.so.1')
+        setenv('PROJ_LIB', '/usr/share/proj/');
+end
+
+if ~exist(options,'overlay_image'), % no image provided, go look into ModelData for one
         if exist(options,'geotiff_name'),
                 paths = {getfieldvalue(options,'geotiff_name')};
@@ -83 +93 @@
         %Crop radar image from xylim
         filename='./temp.tif';
-        eval(['!gdal_translate -quiet -projwin ' num2str(x0) ' ' num2str(y1) ' ' num2str(x1) ' ' num2str(y0) ' ' geotiff_name ' ' filename ]);
+
+
+        system(['gdal_translate -quiet -projwin ' num2str(x0) ' ' num2str(y1) ' ' num2str(x1) ' ' num2str(y0) ' ' geotiff_name ' ' filename ]);
 
         %Read in temp.tif:
@@ -98 +110 @@
                 system('rm -rf ./temp.tif');
         end
-        %}}}
-else %user provided image {{{
+else %user provided image
 
         %user provided an image. check we also have overlay_xlim and overlay_ylim  options, to know what range of coordinates the image covers.
@@ -154 +165 @@
                 md.radaroverlay.y=(y0:(y1-y0)/(size(md.radaroverlay.pwr,1)-1):y1);
         end
-end %}}}
-
-%Was a triangulation requested for the area of the image that is not covered by the mesh? %{{{
+end
+
+%Was a triangulation requested for the area of the image that is not covered by the mesh?
 if strcmpi(getfieldvalue(options,'outertriangulation','no'),'yes'),
 
@@ -194 +205 @@
         md.radaroverlay.outery=outermd.mesh.y;
 
-end %}}}
+end