MOM-TOPAZ 대기 강제력 입력을 위한 ERA5 데이터 변환 스크립트

[HCJung-jbnu]

#!/bin/ksh

YEAR=1993

while [ ${YEAR} -le 2023 ]
do

  for VAR in u10 v10 qair t2m msl msr mtpr ssrd strd
  do

    ls /data/kong/knusml5/ERA5_day/ERA5_y${YEAR}m??d??.nc > temp.list

    mkdir -p ${VAR}

    while read fname
    do

      fname2=`basename ${fname}`

      echo "${fname2} -> ${VAR}"

      cdo select,name=${VAR} ${fname} ./${VAR}/ERA5_${VAR}_${fname2:5:11}.nc

    done < temp.list

    rm -f temp.list

    cdo -f nc2 -b F64 mergetime ./${VAR}/ERA5_${VAR}_y${YEAR}m??d??.nc ERA5_${VAR}_y${YEAR}_tmp.nc

    ncatted -a modulo,time,c,c," " ERA5_${VAR}_y${YEAR}_tmp.nc

    cdo invertlat ERA5_${VAR}_y${YEAR}_tmp.nc ERA5_${VAR}_y${YEAR}.nc

    rm -f ERA5_${VAR}_y${YEAR}_tmp.nc

  done

  ((YEAR++))

done

[HCJung-jbnu]

#!/bin/ksh

YEAR=1994

while [ ${YEAR} -le 2022 ]
do

  for VAR in u10 v10 qair t2m msl msr mtpr ssrd strd
  do

    ((YEAR2 = YEAR + 1))

    cdo invertlat ./${VAR}/ERA5_${VAR}_y${YEAR2}m01d01.nc ./invert.nc

    cdo -f nc2 -b F64 mergetime ./ERA5_${VAR}_y${YEAR}.nc ./invert.nc ./ERA5_${VAR}_y${YEAR}_final.nc

    rm -f ./invert.nc

#    ncatted -a modulo,time,c,c," " ./ERA5_${VAR}_y${YEAR}_final.nc

  done

  ((YEAR++))

done

[HCJung-jbnu]

#!/bin/ksh

YEAR=2011

while [ ${YEAR} -le 2022 ]
do

  for VAR in u10 v10 qair t2m msl msr mtpr ssrd strd
  do

     cdo daymean -cat "ERA5_${VAR}_y${YEAR}_final_sm.nc" temp1.nc

     cdo seltimestep,1 ERA5_${VAR}_y${YEAR}_final_sm.nc temp2.nc

     cdo mergetime temp2.nc temp1.nc ERA5_${VAR}_y${YEAR}_final_sm_daily.nc

     rm -rf ./temp?.nc

  done

  ((YEAR++))

done