Changeset 6488

Timestamp:

11/04/10 07:53:15 (14 years ago)

Author:

Mathieu Morlighem

Message:

Simplified Radarpower

File:

• issm/trunk/src/m/model/radarpower.m (modified) (6 diffs)

issm/trunk/src/m/model/radarpower.m

@@ -1 @@
-function md=radarpower(md,hemisphere,xlim,ylim,highres)
+function md=radarpower(md,xlim,ylim,highres)
 %RADARPOWER - overlay a power radar image on an existing mesh
 %
@@ -7 +7 @@
 %
 %   Usage:
-%      md=radarpower(md,hemisphere,xlim,ylim,highres)
+%      md=radarpower(md,xlim,ylim,highres)
+%      md=radarpower(md,xlim,ylim)
+%      md=radarpower(md)
 
 global ISSM_TIER
@@ -15 +17 @@
 if isempty(MODELDATA), MODELDATA='/u/astrid-r1b/larour/ModelData'; end
 
+%Parse inputs
+if nargin<4,
+        highres=0;
+end
+if nargin<3,
+        xlim=[min(md.x) max(md.x)];
+        ylim=[min(md.y) max(md.y)];
+end
+
 %find gdal coordinates
-x0=min(xlim);
-x1=max(xlim);
-
-y0=min(ylim);
-y1=max(ylim);
+x0=min(xlim); x1=max(xlim);
+y0=min(ylim); y1=max(ylim);
 
 %the geotiff image is either 200m or 1km accuracy. 
-if strcmpi(hemisphere,'n'),
+if strcmpi(md.hemisphere,'n'),
         if ~exist([MODELDATA '/MOG/mog150_greenland_map.jpg']),
                 error(['radarpower error message: file ' MODELDATA '/MOG/mog150_greenland_map.jpg not found. Check MODELDATA variable..']);
@@ -35 +43 @@
 
         %Read and crop file
+        disp('Warning: expecting coordinates in polar stereographic (Std Latitude: 70ºN Meridian: 45º)');
         im=imread(jpgim);
         im=im(toplefty:toplefty+sizey,topleftx:topleftx+sizex);
@@ -41 +50 @@
         md.sarym=(y0:(y1-y0)/(size(md.sarpwr,1)-1):y1);
 
-else
+elseif strcmpi(md.hemisphere,'s'),
         if highres,
                 if ~exist([MODELDATA '/MosaicTiffRsat/amm125m_v2_200m.tif']),
@@ -66 +75 @@
         system('rm -rf ./temp.tif');
 
+else
+        error('field hemisphere should either be ''n'' or ''s''');
 end