#!/bin/sh datey=`date +%y` datem=`date +%m` dated=`date +%d` dateH=`date +%H` min00=`echo 00` min05=`echo 05` min10=`echo 10` min15=`echo 15` min20=`echo 20` min25=`echo 25` yeard=`echo 20` endat=`echo .00` resol='crugrid_96x73' #resol='crugrid' scale='12' stashid='05216' stasholr='02205' stashu='15201' stashv='15202' latid='5' latid2='10' latid3='15' for jobid in 'aabaaa' 'exeaua' 'gcnaaa' 'gcoaaa' 'exeava' 'exeawa' 'aacaaa' do cp ${jobid}.pa1988.${stashid}.${resol}.bint.nc toto1.nc # Std cp ${jobid}.pa1988.${stasholr}.${resol}.bint.nc toto1.${stasholr}.nc # Std cp ${jobid}.pa1988.${stashu}.${resol}.bint.nc toto1.${stashu}.nc # Std cp ${jobid}.pa1988.${stashv}.${resol}.bint.nc toto1.${stashv}.nc # Std #latid='14' #latid2='16' # first experiment: latitude 15 N ferret_v581 -memsize 9 -script > timeser_waf_${latid}N.jnl << EOH set mode metafile use toto1.nc use toto1.${stasholr}.nc use toto1.${stashu}.nc set window/aspect=.6 set region/x=10W:370E/Y=${latid}N:${latid}N let rra=86400.*(PRECIP[d=toto1.nc]) let rrb=1.*(OLR[d=toto1.${stasholr}.nc]) let rrc=1.*(U[d=toto1.${stashu}.nc]) let rrd=1.*(U[d=toto1.${stashv}.nc]) let modu=(rrc*rrc + rrd*rrd)^0.5 let rainfall=rra[x=@ave,l=@sbx:10] let volr=rrb[x=@ave,l=@sbx:10] let uwind=rrc[x=@ave,z=925,l=@sbx:10] SET VARIABLE/TITLE="Precipitation"/Unit="mm/day" rainfall SET VARIABLE/TITLE="OLR"/Unit="w/m^2" volr SET VARIABLE/TITLE="Zonal wind"/Unit="m/s" uwind go left_axis_plot/nolabel/col=2/thick=3 volr 2 go right_axis_plot/nolabel/col=3/thick=3 uwind !plot/ov/col=4/nolabel/thick=3/SIZE=0.12 vu !plot/nolabel/thick=3 vprec !plot/ov/col=2/nolabel/thick=3 volr !plot/ov/col=4/nolabel/thick=3/SIZE=0.12 vu exit EOH Fprint -o fer_SM.ps -p portrait metafile.plt -l cps rm -f metafile_*.plt ferret.jnl* dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT fer_SM.x dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT timeser_waf_${latid}N.jnl dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT fer_SM.ps . mv fer_SM.ps ${jobid}_${latid}N.${resol}.1988.ps ~/bin/ps_thicken ${jobid}_${latid}N.${resol}.1988.ps 4 > file_thick.ps mv file_thick.ps ${jobid}_${latid}N.${resol}.1988.ps # 2nd experiment: latitude 10 N ferret_v581 -memsize 9 -script > timeser_waf_${latid2}N.jnl << EOH set mode metafile use toto1.nc use toto1.${stasholr}.nc use toto1.${stashu}.nc set window/aspect=.6 set region/x=10W:370E/Y=${latid2}N:${latid2}N let rra=86400.*(PRECIP[d=toto1.nc]) let rrb=1.*(OLR[d=toto1.${stasholr}.nc]) let rrc=1.*(U[d=toto1.${stashu}.nc]) let rrd=1.*(U[d=toto1.${stashv}.nc]) let modu=(rrc*rrc + rrd*rrd)^0.5 let rainfall=rra[x=@ave,l=@sbx:10] let volr=rrb[x=@ave,l=@sbx:10] let uwind=rrc[x=@ave,z=925,l=@sbx:10] SET VARIABLE/TITLE="Precipitation"/Unit="mm/day" rainfall SET VARIABLE/TITLE="OLR"/Unit="w/m^2" volr SET VARIABLE/TITLE="Zonal wind"/Unit="m/s" uwind go left_axis_plot/nolabel/col=2/thick=3 volr go right_axis_plot/nolabel/col=3/thick=3 uwind exit EOH Fprint -o fer_SM.ps -p portrait metafile.plt -l cps rm -f metafile_*.plt ferret.jnl* dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT fer_SM.x dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT timeser_waf_${latid2}N.jnl dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT fer_SM.ps . mv fer_SM.ps ${jobid}_${latid2}N.${resol}.1988.ps ~/bin/ps_thicken ${jobid}_${latid2}N.${resol}.1988.ps 4 > file_thick.ps mv file_thick.ps ${jobid}_${latid2}N.${resol}.1988.ps # 3rd experiment: latitude 10 N ferret_v581 -memsize 9 -script > timeser_waf_${latid3}N.jnl << EOH set mode metafile use toto1.nc use toto1.${stasholr}.nc use toto1.${stashu}.nc set window/aspect=.6 set region/x=10W:370E/Y=${latid3}N:${latid3}N let rra=86400.*(PRECIP[d=toto1.nc]) let rrb=1.*(OLR[d=toto1.${stasholr}.nc]) let rrc=1.*(U[d=toto1.${stashu}.nc]) let rrd=1.*(U[d=toto1.${stashv}.nc]) let modu=(rrc*rrc + rrd*rrd)^0.5 let rainfall=rra[x=@ave,l=@sbx:10] let volr=rrb[x=@ave,l=@sbx:10] let uwind=rrc[x=@ave,z=925,l=@sbx:10] SET VARIABLE/TITLE="Precipitation"/Unit="mm/day" rainfall SET VARIABLE/TITLE="OLR"/Unit="w/m^2" volr SET VARIABLE/TITLE="Zonal wind"/Unit="m/s" uwind go left_axis_plot/nolabel/col=2/thick=3 volr go right_axis_plot/nolabel/col=3/thick=3 uwind exit EOH Fprint -o fer_SM.ps -p portrait metafile.plt -l cps rm -f metafile_*.plt ferret.jnl* dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT fer_SM.x dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT timeser_waf_${latid3}N.jnl dateT=`echo $yeard$datey$dated$dateH$min00$endat` touch -t $dateT fer_SM.ps . mv fer_SM.ps ${jobid}_${latid3}N.${resol}.1988.ps ~/bin/ps_thicken ${jobid}_${latid3}N.${resol}.1988.ps 4 > file_thick.ps mv file_thick.ps ${jobid}_${latid3}N.${resol}.1988.ps done ls -lrt