diff --git a/matlab/JET/gdat_jet.m b/matlab/JET/gdat_jet.m
index d522d546b81f183e28ccfc711951c545308688c7..915e0152cc33f4fc0cac71ac110a085559f94857 100644
--- a/matlab/JET/gdat_jet.m
+++ b/matlab/JET/gdat_jet.m
@@ -625,6 +625,9 @@ elseif strcmp(mapping_for_jet.method,'switchcase')
       end
       gdat_data.gdat_params.time = time;
       gdat_data.t = time;
+      if ~isfield(gdat_data.gdat_params,'equil') || isempty(gdat_data.gdat_params.equil)
+        gdat_data.gdat_params.equil = 'EFIT';
+      end
       zshift = 0.;
       if isfield(gdat_data.gdat_params,'zshift') && ~isempty(gdat_data.gdat_params.zshift)
         zshift = gdat_data.gdat_params.zshift;
@@ -642,7 +645,8 @@ elseif strcmp(mapping_for_jet.method,'switchcase')
       if ~isfield(gdat_data.gdat_params,'nz') || isempty(gdat_data.gdat_params.nz)
         gdat_data.gdat_params.nz = 129;
       end
-      [efitdata,eqd]=geteqdskJET(shot,time,gdat_data.gdat_params.nr,gdat_data.gdat_params.nz,[],zshift_eff,[],[],[],100,[],gdat_data.gdat_params.jet_user);
+      [efitdata,eqd]=geteqdskJET(shot,time,gdat_data.gdat_params.nr,gdat_data.gdat_params.nz,[],zshift_eff, ...
+                                      gdat_data.gdat_params.equil,[],[],100,[],gdat_data.gdat_params.jet_user);
       if length(time) > 1
         gdat_data.eqdsk = eqd;
         for itime=1:length(time)