Observation data download and process

1. https://coastwatch.glerl.noaa.gov/
2. The average surface temperature
   
   https://coastwatch.glerl.noaa.gov/ftp/glsea/avgtemps/1995/glsea-temps1995_1024.dat
   
   The data can be downloaded by the following command:

wget https://coastwatch.glerl.noaa.gov/ftp/glsea/avgtemps/1995/glsea-temps1995_1024.dat

or we can use loop to download all the data of different years:

for yr in{1995..2022}
do
echo $yr
wget https://coastwatch.glerl.noaa.gov/ftp/glsea/avgtemps/${yr}/glsea-temps${yr}_1024.dat
done

2. The ice coverage
   https://coastwatch.glerl.noaa.gov/statistic/ice/dat/g2008_2009_ice.dat
   
3. The spatial surface average temperature
   
   https://coastwatch.glerl.noaa.gov/erddap/index.html
   
   https://coastwatch.glerl.noaa.gov/erddap/search/index.html?page=1&itemsPerPage=1000&searchFor=glsea

https://coastwatch.glerl.noaa.gov/erddap/files/GLSEA_GCS/

4. The spatial ice coverage
https://coastwatch.glerl.noaa.gov/erddap/search/index.html?page=1&itemsPerPage=1000&searchFor=glsea

The data of step 3 & 4 can also be downloaded by commands:

yr=2022
nday=365
for ((i=1;i<=${nday};i++));
do
 printf -v jd “%03d” ${i}
 mm=$(date -d “${yr}-01-01 +$(( ${i} - 1 ))day” “+%m”)
 dd=$(date -d “${yr}-01-01 +$(( ${i} - 1 ))day” “+%d”)
xtmp=“wget -N https://coastwatch.glerl.noaa.gov/erddap/files/GLSEA_GCS/${yr}/${mm}/${yr}_${jd}_glsea_sst.nc”
echo ${xtmp}
eval ${xtmp}
#
xtmp=“wget -N https://coastwatch.glerl.noaa.gov/erddap/files/GL_Ice_Concentration_GCS/${yr}/${mm}/${yr}_${jd}_glsea_ice.nc”
echo ${xtmp}
eval ${xtmp}
done

5. The process of the observation data
   The main thought is to interpolate the origin observation data into fvcom grids.

clear all;close all;clc;

map=load('/Users/chuyanzhao/Desktop/observation_data/process_code/boundary_superior.dat');

%input yy/mm/dd
yr=2019;
mm=1;
dd=1;

%
bst0=datenum(yr,1,1);
bst2=datenum(yr,mm,dd);
xday=bst2-bst0+1;

ncload(['/Users/chuyanzhao/Desktop/observation_data/obs_glsea_SST/2019_glsea_sst/' num2str(yr) '_' num2str(xday,'%.3i') '_glsea_sst.nc'],'lon','lat','sst');
[xlon0,xlat0]=meshgrid(lon,lat);
[ny0,nx0]=size(xlon0);
sst_glsea=sst;

ncload(['/Users/chuyanzhao/Desktop/output/superior_0001.nc'],'lon','lat','temp','time','zeta');
xlon=lon-360;
xlat=lat;

% interpolate the origin data to fvcom grids
clear ind
ind=find(sst_glsea>=0);% export the location of lake
clear tmp
tmp=interp_scatter_quick(xlon0(ind),xlat0(ind),sst_glsea(ind),xlon,xlat,0.5,1e-5,2,2);
temp=tmp';

% figure with original data
figure(1);clf;hold on;
scattercontourf(xlon0(ind),xlat0(ind),sst_glsea(ind),[0:1:10]);
plot_fvcom_obc(map,[210 180 140]/255);
axis([-92.2 -84.4 46.2 49.2]);

% figure with interpolated data
figure(2);clf;hold on;
scattercontourf(xlon,xlat,temp',[0:1:10]);
plot_fvcom_obc(map,[210 180 140]/255);    
axis([-92.2 -84.4 46.2 49.2]);