Changeset 7 for Roms_tools/Aforc_NCEP/get_NCEP_grid.m

Timestamp:

05/04/11 17:34:07 (13 years ago)

Author:

cholod

Message:

patch_romstools_16_02_2011.tar : Bug in make_NCEP.m, get_NCEP_grid.m

File:

• Roms_tools/Aforc_NCEP/get_NCEP_grid.m (modified) (4 diffs)

Roms_tools/Aforc_NCEP/get_NCEP_grid.m

@@ -79 +79 @@
 lat=lat(j);
 %
+% If we are in OpenDap (Get_My_Data=0) we need
+% a shift of decal=1, because the indexes begin at 0 in OpenDap.
+% If we use local data, Get_My_Data=1, there is no shift needed, 
+% then decal=0;
+
+if Get_My_Data==1
+    decal=0;
+else
+    decal=1;
+end
+%
 if ~isempty(i1)
-  i1min=min(i1)-1;
-  i1max=max(i1)-1;
+  i1min=min(i1)-decal;
+  i1max=max(i1)-decal;
 else
   i1min=[];
@@ -87 +98 @@
 end
 if ~isempty(i2)
-  i2min=min(i2)-1;
-  i2max=max(i2)-1;
+  i2min=min(i2)-decal;
+  i2max=max(i2)-decal;
 else
   i2min=[];
@@ -94 +105 @@
 end
 if ~isempty(i3)
-  i3min=min(i3)-1;
-  i3max=max(i3)-1;
+  i3min=min(i3)-decal;
+  i3max=max(i3)-decal;
 else
   i3min=[];
@@ -101 +112 @@
 end
 %
-jmin=min(j)-1;
-jmax=max(j)-1;
+jmin=min(j)-decal;
+jmax=max(j)-decal;
 jrange=['[',num2str(jmin),':',num2str(jmax),']'];