#! /bin/bash
# Global 8-day LAI data, 1-km resolution; conservative remapping to 0.5 degree; from Dai group
# huge, amounts to 1TB.
export modis_ROOT=/home/gwang/Data/modis
cd $modis_ROOT
mkdir LAI_monthly

for year in 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009
do

cd $modis_ROOT

tar -xvf global_30s_${year}.tar
bunzip2 *${year}*.bz2
mkdir global_30s_${year}
mv *${year}*.nc $modis_ROOT/global_30s_${year}
cp $modis_ROOT/weights.nc $modis_ROOT/global_30s_${year}
cd $modis_ROOT/global_30s_${year}

# mygrid
cat > mygrid << EOF
gridtype = lonlat
xsize    = 140
ysize    = 160
xfirst   = -19.995825
xinc     = 0.5
yfirst   = -40.0041693
yinc     = 0.5
EOF

# January
for dates in 001 009 017 025
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
# cdo gencon,mygrid africa_30s_${year}_${dates}.nc weights.nc
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-01-15 temp2.nc africa_${year}_jan.nc
cp africa_${year}_jan.nc $modis_ROOT/LAI_monthly
rm temp* africa*

# Febuary
for dates in 033 041 049 057
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-02-15 temp2.nc africa_${year}_feb.nc
cp africa_${year}_feb.nc $modis_ROOT/LAI_monthly
rm temp* africa*

# March
for dates in 065 073 081 089
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-03-15 temp2.nc africa_${year}_mar.nc
cp africa_${year}_mar.nc $modis_ROOT/LAI_monthly
rm temp* africa*

# April
for dates in 097 105 113 121
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-04-15 temp2.nc africa_${year}_apr.nc
cp africa_${year}_apr.nc $modis_ROOT/LAI_monthly
rm temp* africa*  

# May
for dates in 129 137 145
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-05-15 temp2.nc africa_${year}_may.nc
cp africa_${year}_may.nc $modis_ROOT/LAI_monthly
rm temp* africa*  


# June
for dates in 153 161 169 177
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-06-15 temp2.nc africa_${year}_jun.nc
cp africa_${year}_jun.nc $modis_ROOT/LAI_monthly
rm temp* africa*  


# July
for dates in 185 193 201 209
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-07-15 temp2.nc africa_${year}_jul.nc
cp africa_${year}_jul.nc $modis_ROOT/LAI_monthly
rm temp* africa*  


# August
for dates in 217 225 233 241
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-08-15 temp2.nc africa_${year}_aug.nc
cp africa_${year}_aug.nc $modis_ROOT/LAI_monthly
rm temp* africa*  


# September
for dates in 249 257 265 273
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-09-15 temp2.nc africa_${year}_sept.nc
cp africa_${year}_sept.nc $modis_ROOT/LAI_monthly
rm temp* africa*


# October
for dates in 281 289 297
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-10-15 temp2.nc africa_${year}_oct.nc
cp africa_${year}_oct.nc $modis_ROOT/LAI_monthly
rm temp* africa*


# November
for dates in 305 313 321 329
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-11-15 temp2.nc africa_${year}_nov.nc
cp africa_${year}_nov.nc $modis_ROOT/LAI_monthly
rm temp* africa* 


# December
for dates in 337 345 353 361
do
cdo sellonlatbox,-20,50,-40,40 global_30s_${year}_${dates}.nc africa_30s_${year}_${dates}.nc # Africa domain
cdo remap,mygrid,weights.nc africa_30s_${year}_${dates}.nc africa_30s_${year}_${dates}_mygrid.nc
done

cdo mergetime africa_30s_*_mygrid.nc temp1.nc
cdo timmean temp1.nc temp2.nc
cdo setdate,${year}-12-15 temp2.nc africa_${year}_dec.nc
cp africa_${year}_dec.nc $modis_ROOT/LAI_monthly
rm temp* africa*

rm global_30s_${year}_*

done

rmdir global*


转载本文请联系原作者获取授权，同时请注明本文来自王国杰科学网博客。
链接地址：https://blog.sciencenet.cn/blog-569118-541704.html

上一篇：外面鞭炮很响
下一篇：install RegCM4 on HUYGENS