Some correction in cse of iteration of hindcast

Some correction for plotting

Forecast_tools/iteration.m

@@ -3,7 +3,9 @@
 start
 crocotools_param
 
-time=(floor(now)-datenum(Yorig,1,1)-(hdays-1)-timezone/24)*86400;
+%time=(floor(now)-datenum(Yorig,1,1)-(hdays-1)-timezone/24)*86400;
+time=(floor(now)-datenum(Yorig,1,1)-(hdays)-timezone/24)*86400;
+time
 
 % shift back restart time after first iteration
 eval(['!cp ','croco_rst.nc',' ','croco_ini.nc'])
@@ -12,9 +14,10 @@ nc{'scrum_time'}(1)=time;
 close(nc)
 
 for i=1:3
+%for i=1:1
   disp(['I=',num2str(i+1)])
-  %eval(['!export OMP_NUM_THREADS=4 ; ./croco ','croco_hindcast.in ',' > ','croco_hindcast_',num2str(date),'.out']);
-  eval(['!mpirun -np 4 ./croco ','croco_hindcast.in ',' > ','croco_hindcast_',num2str(date),'.out']);
+  %eval(['!export OMP_NUM_THREADS=4 ; ./croco croco_hindcast.in > croco_hindcast_',datestr(date,'yyyy-mm-dd'),'.out']);
+  eval(['!mpirun -np 4 ./croco croco_hindcast.in > croco_hindcast_',datestr(date,'yyyy-mm-dd'),'.out']);
 
   % shift back restart time for next iteration
   eval(['!cp ','croco_rst.nc',' ','croco_ini.nc'])

Forecast_tools/plot_forecast_croco.m

@@ -33,7 +33,7 @@ mean_currents = 0;  % 1: plot 0-500 mean currents
                     % 0: surface currents
 
 plot_fig_1 = 0; % wind
-plot_fig_2 = 0; % surface currents
+plot_fig_2 = 1; % surface currents
 plot_fig_3 = 1; % SST + surface currents
 plot_fig_4 = 1; % SSH
 plot_fig_5 = 0; % HBL
@@ -160,7 +160,7 @@ if plot_fig_1,
 
   % export_fig -transparent file.pdf
   % eval(['! mv -f file.pdf ',dirout,outname,'.pdf']);
-  eval(['print -dpdf ',dirout,outname,'.pdf']);
+  eval(['print -bestfit -dpdf',dirout,outname,'.pdf']);
  end
 close(nc)
 end % plot_fig_1
@@ -249,7 +249,7 @@ for tndx=hdays+1:nx+hdays
 
     % export_fig -transparent file.pdf
     % eval(['! mv -f file.pdf ',dirout,outname,'.pdf']);
-    eval(['print -dpdf ',dirout,outname,'.pdf']);
+    eval(['print -bestfit -dpdf ',dirout,outname,'.pdf']);
   end
 end
 end % plot_fig_2
@@ -327,7 +327,7 @@ for tndx=hdays+1:nx+hdays
 
     % export_fig -transparent file.pdf
     % eval(['! mv -f file.pdf ',dirout,outname,'.pdf']);
-    eval(['print -dpdf ',dirout,outname,'.pdf']);
+    eval(['print -bestfit -dpdf ',dirout,outname,'.pdf']);
   end
 end
 end % plot_fig_3
@@ -383,7 +383,7 @@ for tndx=hdays+1:nx+hdays
 
     % export_fig -transparent file.pdf
     % eval(['! mv -f file.pdf ',dirout,outname,'.pdf']);
-    eval(['print -dpdf ',dirout,outname,'.pdf']);
+    eval(['print -bestfit -dpdf ',dirout,outname,'.pdf']);
   end
 end
 end % plot_fig_4
@@ -437,7 +437,7 @@ for tndx=hdays+1:nx+hdays
 
     % export_fig -transparent file.pdf
     % eval(['! mv -f file.pdf ',dirout,outname,'.pdf']);
-    eval(['print -dpdf ',dirout,outname,'.pdf']);
+    eval(['print  -bestfit -dpdf ',dirout,outname,'.pdf']);
   end
 end
 end % plot_fig_5