Changeset 18324


Ignore:
Timestamp:
08/03/14 19:08:20 (11 years ago)
Author:
Eric.Larour
Message:

CHG: minor changes for Greenland.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • issm/branches/trunk-jpl-ad2/src/m/plot/prepare_gl.m

    r18323 r18324  
    1313%database: {{{
    1414if strcmpi(md.mesh.hemisphere','s'),
    15         geotiff_name_earth=['/u/astrid-r1b/ModelData/MosaicTiffRsat/amm125m_v2_200m.wgs84.tif'];
    1615        radar_image_earth='/Users/larour/ModelData/TrueMarble_GeoTIFF/radar_no_antarctica.png';
    17         triangulation_earth='/Users/larour/ModelData/DEM_geotiff/processed/triangulation.mat';
     16        triangulation_earth='/Users/larour/ModelData/DEM_geotiff/processed/triangulation_antarctica.mat';
    1817        triangulation='/Users/larour/ModelData/BamberDEMAntarctica1km/triangulated/triangulation_15.mat';
    1918        radar_image='/Users/larour/ModelData/BamberDEMAntarctica1km/triangulated/radar.latlong.png';
    2019
    2120elseif strcmpi(md.mesh.hemisphere','n'),
    22         geotiff_name_earth=['/u/astrid-r1b/ModelData/MosaicTiffRsat/amm125m_v2_200m.wgs84.tif'];
    2321        radar_image_earth='/Users/larour/ModelData/TrueMarble_GeoTIFF/radar_no_greenland.png';
     22        triangulation_earth='/Users/larour/ModelData/DEM_geotiff/processed/triangulation_greenland.mat';
    2423        triangulation='/u/astrid-r1b/ModelData/HowatDEMGreenland2012/triangulated/triangulation_15.mat';
    25         radar_image='/u/astrid-r1b/ModelData/HowatDEMGreenland2012/triangulated/radar.png';
    26         triangulation_hd='/u/astrid-r1b/ModelData/HowatDEMGreenland2012/triangulated/triangulation_2.mat';
     24        radar_image='/u/astrid-r1b/ModelData/HowatDEMGreenland2012/triangulated/radar.latlong.png';
    2725else
    2826        error('prepare_gl error usage: ''hemisphere'' should be set to ''s'' or ''n''');
     
    147145R=6371000*ones(length(x),1)+scaling_factor*surface;
    148146
    149 [lat,long]=xy2ll(x,y,-1);
     147if strcmpi(md.mesh.hemisphere','s'),
     148        [lat,long]=xy2ll(x,y,-1);
     149elseif strcmpi(md.mesh.hemisphere','n'),
     150        [lat,long]=xy2ll(x,y,1);
     151end
    150152
    151153x = R .* cosd(lat) .* cosd(long);
Note: See TracChangeset for help on using the changeset viewer.