#rename variable
ncrename -O -v LATITUDE,y              ${file} 

#rename dimension
ncrename -O -d LATITUDE,y       temp.nc       

#sums
ncap -O -s "LOAD_BCPOM[time,x,y]=LOAD_BC+LOAD_POM" ${file} ${file} 

# compute vertical integrals 
\rm dummy.nc
ncwa -O -v AERH2O3D_AER   -a z -y ttl     ${file} dummy.nc  
ncrename -O -v AERH2O3D_AER,AERH2O_AER       dummy.nc 
ncks -F -A -v AERH2O_AER dummy.nc ${file} 

#concatenate variable
ncks -F -A -v AIRE temp.nc ${file}


#transpose axis
ncecat -O $file $file
ncpdq -O -a time,record $file $file
ncwa -O -a record $file $file
ncpdq -O -a time,y,x $file $file

# turn z 
ncpdq -O -a -z $file $file

# create yearly mean file from 12 months in one file
orgfile=MM_${model}.an${an}.mALLYEAR_an.nc 
ncwa -F -O -a time ${file} ${orgfile}

# calculate yearly average fom 12 files each containing 1 month
ncra -F -O  modis_2000??.nc ${filer}


# create record dimension 
ncecat -O $file $file 
# make x unlimited
ncpdq -O -a x,record $file $file 
# delete empty record
ncwa -O -a record $file $file 

# correct axis 
file1=aire_1.nc
ncks -h -F -O -d x,1,180 $file $file1

file2=aire_2.nc
ncks -h -F -O -d x,181,360 $file $file2
ncap -O -s "x[x]=x-360" $file2 $file2

ncrcat -O $file2 $file1 $file
\rm  $file1 $file2

# turn lat to be consistent with interpolated files 
	ncpdq -O -a "-y" $file $file 

# reorder longitude
	ncecat -O $file $file 
	ncpdq -O -a x,record $file $file 
	ncwa -O -a record $file $file 

	file1=DD_${model}.an${AN}.mALLYEAR_1.nc
	ncks -F -O -d x,1,180 $file $file1

	file2=DD_${model}.an${AN}.mALLYEAR_2.nc
	ncks -F -O -d x,181,360 $file $file2
	ncap -O -s "x[x]=x-360" $file2 $file2
	
	ncrcat -O $file2 $file1 $file
	\rm  $file1 $file2	

	ncpdq -O -a time,y,x $file $file 

# identify missing values 
	ncatted -O -a missing_value,$var,o,f,0 ${file} 
	ncatted -O -a missing_value,$var,o,f,-9999 ${file} 



for AN in ${ANN[*]} ; do
# compute monthly mean files
    for IMTH in ${MTH[*]} ; do 
	file=DD_${model}.an${AN}.m${IMTH}.nc
	filemean=DD_${model}.an${AN}.m${IMTH}_m.nc 
	ncra -F -O  ${file} ${filemean}
    done

# concatenate monthly mean files
    monthfile=DD_${model}.an${AN}.mALLYEAR_m.nc    
    ncrcat -F -O DD_${model}.an${AN}.m??_m.nc ${monthfile}

#
# calculate yearly mean file
#
    orgfile=DD_${model}.an${AN}.mALLYEAR_an.nc
    ncra -F -O ${monthfile} ${orgfile}

done

#delete attributes
    ncatted -O -a ,,d,, interdummy.nc
    ncatted -h -a ,global,d,, interdummy.nc
    ncatted -h -a missing_value,,d,, interdummy.nc