#
#######################################################################################################################################################
#
# Set namelist and call model:  

YY=${RESTARTYEAR}

while [ ${YY} -le ${EXPERIMENTEYEAR} ]; do

echo $EXPERIMENTEYEAR
    echo 'rabotaet'

    for ((MM=${RESTARTMONTH}; MM<=12; MM++)); do
       MM2=${MM}
       if [ ${MM} -le 9 ]; then
	   MM2=0${MM}
fi

       echo ${MM}!!!!!!!!!!!! ${EXPNO}_${YY}${MM2}.01_chem_m.nc

	if [ ${YY} -gt ${RESTARTYEAR} -o ${MM} -ge ${RESTARTMONTH} ]; then

	    if [ -s namelist.echam ]; then
		rm -f namelist.echam
	fi
		
#  namelist control variables and output control for grid variables
#  spectral variables are written out by default except liq. water
#  for production runs set LABORT=.FALSE.
#
	    cat > namelist.echam << EOF
&RUNCTL
  LSOCOL=${LSOCOL}
  LRESUME=$RERUN,
  OUT_DATAPATH = "$DPATH"
  OUT_EXPNAME  = "$EXPNO"
  OUT_FILETYPE = 2
  TRAC_FILETYPE = 2
  DT_START  = ${EXPERIMENTSYEAR},01,01,12,0,0
  DT_STOP   = ${EXPERIMENTEYEAR},01,01,12,0,0
  NO_CYCLES = 1
  PUTDATA   = 6, 'hours', 'first', 0
  PUTRERUN  = 1, 'months', 'last', 0
  DELTA_TIME = 900
  LAMIP=.TRUE.
  LMIDATM=.TRUE.
  LABORT=.FALSE.
  NPROCA=${NPROCA}
  NPROCB=${NPROCB}
  NPROMA=${NPROMA}
/
&SOCOLCTL
  LCHEM=${LCHEM}
  lradcoup=.true.
  lo3_coupl=.true.
  lch4_coupl=.true.
  ln2o_coupl=.true.
  lodscl_coupl=.true.
  lh2o_coupl=.true.
  lsrb_lya_heating=.true.
  nastep=20
  CYEAR=2000
  lqbonudg=.true.
  lhetchem=.true.
  lvolcano=.false.
  linteractive_drydep=.true.
  linteractive_wetdep=.true.
  lscav_gas=.true.
  lscav_cv=.true.
  lscav_ls=.true.
  lscav_aer=.true.
  iscav_easy=3
  coeff_para=2
  interactivelnox=.true.
  nlev_fbc=1
  nlev_mbc=1
  nlev_drydep=1
  ln2o5hydro = .true.  
  lgn2o5_const = .false.
  ln2o5_dry = .true.
/
&PHYSCTL
  LCOND= .true.
  LCONV= .true.
/
&DYNCTL
  VCHECK=235.
  SPDRAG=0.926E-4
/
&RADCTL
  IAERO=10
  LGADSRH=.TRUE.
/
EOF


	    # Call model:
	    date
	    #$MODEL
	    mpirun $MODEL


	    # Save restart files:
            # a) ECHAM5:
	    cp -f rerun_${EXPNO}_echam ${EXPNO}_rerun_${YY}${MM2}_echam.nc
	    gzip -f ${EXPNO}_rerun_${YY}${MM2}_echam.nc

            # b) MEZON:
	    if [ ${LSOCOL} = .TRUE. -a ${LCHEM} = .TRUE. ]; then
		cp -f rerun_${EXPNO}_tracer ${EXPNO}_rerun_${YY}${MM2}_chem1.nc
		gzip -f ${EXPNO}_rerun_${YY}${MM2}_chem1.nc
	    fi

	    RERUN=.true.

           #conversion to monthly means    
           cdo monmean ${EXPNO}_${YY}${MM2}.01_aero2.nc ${EXPNO}_${YY}${MM2}.01_aero2_m.nc
           rm ${EXPNO}_${YY}${MM2}.01_aero2.nc
           cdo monmean ${EXPNO}_${YY}${MM2}.01.nc ${EXPNO}_${YY}${MM2}.01_m.nc
           rm ${EXPNO}_${YY}${MM2}.01.nc
           cdo monmean ${EXPNO}_${YY}${MM2}.01_tracer.nc ${EXPNO}_${YY}${MM2}.01_tracer_m.nc
           rm ${EXPNO}_${YY}${MM2}.01_tracer.nc
           echo aero2 is converted to monthly

	fi

    done

    YY=`expr ${YY} + 1`
    RESTARTMONTH=1

done

ls -lt

