Script zum Medokads Datenprozessing (auf wizard)

#! /bin/ksh #set -x

# # create surface monthly data files from origional location : step1 # # $1 : longitute # $2 : latutude # # H.Luthardt 28.11-2013 #

#echo " " ; echo " running j1.sh " ; echo " "

typeset -i SIZE

SW_CDO=/sw/rhel55-x64/cdo-1.5.6/bin/cdo

PROJ=medokads DIR_PROJ="/cera_in/projects/"$PROJ"/"

DIR_DATA=$DIR_PROJ"PROZ/COMPLETE/"

DIR_PROC="${DIR_PROJ}PROZ/" DIR_WRK="${DIR_PROC}WWORK/" DIR_LS="${DIR_PROC}LARGESMALL/"

DIR_CONT_LIST=$DIR_PROJ"container_lists/"

LOG_DIR="${DIR_PROC}LOGS/"

YEARS="89 90 91 92 93 94 95 96 97 98 99 00 01 02 03 04 05 06 07 08 09" #VARS="1 2 3 4 5 6 7 8 9 a i n p s t" VARS="a i n p s t" #VARS="8 9 a i n p s t"

MONTHS="01 02 03 04 05 06 07 08 09 10 11 12" DAYS="01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 18 20 21 22 23 24 25 26 27 28 29 30 31" DAYS0='01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31' DAYS1='01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30' DAYS2='01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29' DAYS3='01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28'

#YEARS="89" #VARS="4" #MONTHS="02"

crone() {

#/sw/rhel55-x64/cdo-1.6.0/bin/cdo selday,$3 allmiss_$6.nc $5$4$3_$1e_$2n_$6.nc

#/sw/rhel55-x64/cdo-1.6.0/bin/cdo,${YEAR}-${mm}-01 settaxis allmiss.nc XXXX_${TILE}_$4.nc

# # gridID 2 # gridtype = lonlat gridsize = 200000 xname = longitude xlongname = longitude xunits = degrees_east yname = latitude ylongname = latitude yunits = degrees_north xsize = 500 ysize = 400 xfirst = $1 xinc = 0.01 yfirst = $2 yinc = -0.00999832 EOF

# east # ${SW_CDO} setgrid,${DIR_PROC}grid ${DIR_PROC}temp_$1$2.nc ${DIR_WRK}$5$4$3_$1e_$2n_$6.nc

# west

}

proc_data() {

# TILE="$1"

# for VV in $VARS ; do

# for YY in $YEARS ; do # echo ${VV} ${YEARS} # done

# done

# for VV in $VARS ; do # done

# exit

# ddir_name=${TILE}"_"${VV} # ds_name="MEDOKADS_"$ddir_name # lst_file_name=${DIR_CONT_LIST}${ds_name} # lst_file=${lst_file_name}".lst"

# echo "Filename: "${lst_file}

# if [ ! -f "${lst_file}" ]; then

# else # echo "File already exists: "${lst_file} # fi

}

cleanup_before() {

# nohup rm ${DIR_WRK}??????_${TILE}_${VV}.nc & # nohup rm WWORK/??????_$1e_$2n_2.nc & # nohup rm WWORK/??????_$1e_$2n_3.nc & # nohup rm WWORK/??????_$1e_$2n_4.nc & # nohup rm WWORK/??????_$1e_$2n_5.nc & # nohup rm WWORK/??????_$1e_$2n_7.nc & # nohup rm WWORK/??????_$1e_$2n_6.nc & # nohup rm WWORK/??????_$1e_$2n_8.nc & # nohup rm WWORK/??????_$1e_$2n_9.nc & # nohup rm WWORK/??????_$1e_$2n_a.nc & # nohup rm WWORK/??????_$1e_$2n_i.nc & # nohup rm WWORK/??????_$1e_$2n_n.nc & # nohup rm WWORK/??????_$1e_$2n_p.nc & # nohup rm WWORK/??????_$1e_$2n_s.nc & # nohup rm WWORK/??????_$1e_$2n_t.nc &

}

cleanup_tile() {

# nohup rm ${DIR_WRK}??????_${TILE}_${VV}.nc &

}

proc_east_part() {

# for n in 35 39 43 47 51 55 ; do # for n in 31 ; do

}

proc_west_part() { # for e_w in 10 ; do

# for n in 31 35 39 43 47 51 55 ; do

}

cleanup_tiles() {

}

echo "Start..."

#proc_east_part proc_west_part #cleanup_tiles

echo "Done..."

#echo " " ; echo "end j1.sh " ; echo " "

exit # =============================================================================

#DIR1=/cera_in/projects/medokads/PROZ #mkdir ${DIR1}/COMPLETE/${TILE}_${VV} #mv ${DIR1}/MN/????_${TILE}_${VV}.nc ${DIR1}/COMPLETE/${TILE}_${VV} #echo "---" #done # VV

MadWiki: Datenprocessing (last edited 2014-04-16 12:34:53 by HansRamthun)