function plot_roms_tidal_forcing_latte(gfile,tfile)








lon=nc_varget(gfile,'lon_rho');
lat=nc_varget(gfile,'lat_rho');
h=nc_varget(gfile,'h');

eta=nc_varget(tfile,'tide_Eamp');
eph=nc_varget(tfile,'tide_Ephase');
cph=nc_varget(tfile,'tide_Cphase');
cag=nc_varget(tfile,'tide_Cangle');
cmn=nc_varget(tfile,'tide_Cmin');
cmx=nc_varget(tfile,'tide_Cmax');
tnames=nc_varget(tfile,'tidal_constituents');
dims=size(eta);
size(tnames)
for i=1:dims(1)
  cons=tnames(i,:)
  
  figure(1)
  clf
 % subplot(2,3,1),
 
  set(gcf,'renderer','zbuffer')
  latte_roms_map
  hold on
  
  pcolor(lon,lat,squeeze(eta(i,:,:)))
  caxis([0 0.8])
  contour(lon,lat,h,[10 20 30 40 50 60],'k-')
  shading flat
  hold off
set(gca,'xlim',[  -75.1437  -72.2179]);
  set(gca,'ylim',[38.3815   40.7964]);
  
  colorbar
  title(['EAMP - ' cons])
  print(gcf,'-djpeg99',['OTPS_EAMP_' cons '.jpg'])
  
  clf
 % subplot(2,3,2),
  set(gcf,'renderer','zbuffer')
  latte_roms_map
  hold on
  pcolor(lon,lat,squeeze(eph(i,:,:)))
  caxis([100 140])
  shading flat
  contour(lon,lat,h,[10 20 30 40 50 60],'k')
  hold off
set(gca,'xlim',[  -75.1437  -72.2179]);
  set(gca,'ylim',[38.3815   40.7964])
  title(['EPHASE - ' cons])
  colorbar
  print(gcf,'-djpeg99',['OTPS_EPHASE_' cons '.jpg'])
  
  clf
  %subplot(2,3,3),
  set(gcf,'renderer','zbuffer')
  latte_roms_map
  hold on
  pcolor(lon,lat,squeeze(cph(i,:,:)))
  caxis([0 360])
  shading flat
  contour(lon,lat,h,[10 20 30 40 50 60],'k')
  hold off
set(gca,'xlim',[  -75.1437  -72.2179]);
  set(gca,'ylim',[38.3815   40.7964])
  title(['CPHASE - ' cons])
  colorbar
  print(gcf,'-djpeg99',['OTPS_CPHASE_' cons '.jpg'])

  clf
  %subplot(2,3,4),
  set(gcf,'renderer','zbuffer')
  latte_roms_map
  hold on
  pcolor(lon,lat,squeeze(cag(i,:,:)))
  caxis([0 360])
  shading flat
  contour(lon,lat,h,[10 20 30 40 50 60],'k')
  hold off
set(gca,'xlim',[  -75.1437  -72.2179]);
  set(gca,'ylim',[38.3815   40.7964])
  title(['CANGLE - ' cons])
  colorbar
  print(gcf,'-djpeg99',['OTPS_CANGLE_' cons '.jpg'])
  clf
 % subplot(2,3,5),
  set(gcf,'renderer','zbuffer')
  latte_roms_map
  hold on
  pcolor(lon,lat,squeeze(cmx(i,:,:)))
  caxis([0 0.2])
  shading flat
  contour(lon,lat,h,[10 20 30 40 50 60],'k')
  hold off
set(gca,'xlim',[  -75.1437  -72.2179]);
  set(gca,'ylim',[38.3815   40.7964])
  title(['CMAX - ' cons])
  colorbar
  print(gcf,'-djpeg99',['OTPS_CMAX_' cons '.jpg'])
  
  clf
  %subplot(2,3,6),
  set(gcf,'renderer','zbuffer')
  latte_roms_map
  hold on
  pcolor(lon,lat,squeeze(cmn(i,:,:)))
  caxis([-0.2 0.2])
  shading flat
  contour(lon,lat,h,[10 20 30 40 50 60],'k')
  hold off
set(gca,'xlim',[  -75.1437  -72.2179]);
  set(gca,'ylim',[38.3815   40.7964])
  title(['CMIN - ' cons])
  colorbar
  print(gcf,'-djpeg99',['OTPS_CMIN_' cons '.jpg'])
  
  
  
  
  
  
  
  
end

