In Preprocessing_tools, during geostrophic calculation,

add the Ekman transport only of forcing file (containing wind strss) exists

Preprocessing_tools/geost_currents.m

@@ -162,8 +162,8 @@ for l=1:tlen
 %
 % Ekman transport
 %
-  if ~isempty(frcname)
-%    disp('Add the Ekman transport')
+  if exist(frcname)
+    disp('Add the Ekman transport')
     sustr(:,:)=nfrc{'sustr'}(l,:,:);
     svstr(:,:)=nfrc{'svstr'}(l,:,:);
 %    rhoA=1.3; Cd=1.4e-3;

Preprocessing_tools/geost_currents_bry.m

@@ -193,8 +193,8 @@ for l=1:tlen
 %
 % Ekman transport
 %
-  if ~isempty(frcname)
-%    disp('Add the Ekman transport')
+   if exist(frcname)
+      disp('Add the Ekman transport')
     if obcndx==1
       tmp=squeeze(nfrc{'sustr'}(l,1,:));
       sustr=0*h;