Changeset 23 for trunk/src

Timestamp:

12/20/10 18:21:39 (13 years ago)

Author:

pinsard

Message:

while vs for in compute_erai_daily_region_2d.sh

File:

• trunk/src/compute_erai_daily_region_2d.sh (modified) (6 diffs)

trunk/src/compute_erai_daily_region_2d.sh

@@ -107 +107 @@
 # ====
 #
+# make it work 
+#
+# avoid ncrcat interactive question
+#
 # split get and process files
 #
@@ -121 +125 @@
 # EVOLUTIONS
 # ==========
+#
+# - fplod 20101220T150356Z aedon.locean-ipsl.upmc.fr#
+#
+#   * indentation
+#   * replace date by yyyy (date is a shell command)
+#   * replace for yyyy by while yyyy
+#   * replace for month by while month
+#   * replace for day by while day
+#   * work only on existing dates
 #
 # - fplod 20101220T142107Z aedon.locean-ipsl.upmc.fr (Darwin)
@@ -205 +218 @@
 maxlon=360.
 reg=GLOB
+yearmin=1989
+yearmax=2010
+monthmin=1
+monthmax=12
+cmonthmin=$(printf "%2.2d" ${monthmin})
+cmonthmax=$(printf "%2.2d" ${monthmax})
 dirin=/bdd/ERAI/NETCDF/GLOBAL_075/4xdaily/
 # check for dirin existence
@@ -223 +242 @@
 # cp = Convective precipitation
 for var in lsp cp; do
-for date in 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010; do
-for month in 01 02 03 04 05 06 07 08 09 10 11 12; do
-for day in 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; do
-ncpdq -U ${dirin}/FC_SF/${date}/${month}/${var}.${date}${month}${day}* ${dirtemp}/${var}.${date}${month}${day}.temp.GLOBAL_075.nc
+   yyyy=${yearmin}
+   while [ ${yyyy} -le ${yearmax} ]
+   do
+      month=${monthmin}
+      while [ ${month} -le ${monthmax} ]
+      do
+         cmonth=$(printf "%2.2d" ${month})
+         daymin=1
+         daymax=$(cal ${month} ${yyyy} | grep . | fmt -1 | tail -1)
+         day=${daymin}
+         while [ ${day} -le ${daymax} ]
+         do
+            cday=$(printf "%2.2d" ${day})
+            if [ -f ${dirin}/FC_SF/${yyyy}/${cmonth}/${var}.${yyyy}${cmonth}${cday}* ] 
+            then
+               ncpdq -U ${dirin}/FC_SF/${yyyy}/${cmonth}/${var}.${yyyy}${cmonth}${cday}* ${dirtemp}/${var}.${yyyy}${cmonth}${cday}.temp.GLOBAL_075.nc
+            else
+               echo "${command} : iii : no files ${dirin}/FC_SF/${yyyy}/${cmonth}/${var}.${yyyy}${cmonth}${cday}*" >> ${log}
+            fi
+            day=$(( ${day} + 1 ))
+         done
+         exist_temp_files=$(find ${dirtemp} -name "${var}.${yyyy}${cmonth}*temp*")
+         if [ "${exist_temp_files}" != "" ]
+         then 
+            ncrcat ${dirtemp}/${var}.${yyyy}${cmonth}*temp*  ${dirtemp}/${var}.${yyyy}${cmonth}.fshei.GLOBAL_075.nc
+            rm ${dirtemp}/${var}.${yyyy}${cmonth}*temp* 2> /dev/null
+            ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,0,,5 ${dirtemp}/${var}.${yyyy}${cmonth}.fshei.GLOBAL_075.nc ${dirtemp}/temp_time0.nc
+            ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,4,,5 ${dirtemp}/${var}.${yyyy}${cmonth}.fshei.GLOBAL_075.nc ${dirtemp}/temp_time4.nc
+            ncflint -w -1,1 ${dirtemp}/temp_time0.nc ${dirtemp}/temp_time4.nc ${dirtemp}/temp_time.nc
+            ncflint -w 1000.,0.0 ${dirtemp}/temp_time.nc ${dirtemp}/temp_time.nc ${dirtemp}/temp_int_${var}_${yyyy}${cmonth}.nc
+            rm ${dirtemp}/temp_time* 2> /dev/null
+            rm ${dirtemp}/${var}.${yyyy}${cmonth}.fshei.GLOBAL_075.nc 2> /dev/null
+         else
+            echo "${command} : iii : no files ${dirtemp}/${var}.${yyyy}${cmonth}*temp*" >> ${log}
+         fi
+         unset exist_temp_files
+         month=$(( ${month} + 1 ))
+      done
+      ncrcat ${dirtemp}/temp_int_${var}_${yyyy}* ${TROPFLUX_ID}/erai_${reg}_1d_${yyyy}${cmonthmin}01_${yyyy}${cmonthmax}31_${var}_gridOrig.nc
+      rm ${dirtemp}/temp_int_${var}_${yyyy}* 2> /dev/null
+      yyyy=$(( ${yyyy} + 1 ))
+   done
 done
-ncrcat ${dirtemp}/${var}.${date}${month}*temp*  ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc
-\rm ${dirtemp}/${var}.${date}${month}*temp*
-ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,0,,5 ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc ${dirtemp}/temp_time0.nc
-ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,4,,5 ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc ${dirtemp}/temp_time4.nc
-ncflint -w -1,1 ${dirtemp}/temp_time0.nc ${dirtemp}/temp_time4.nc ${dirtemp}/temp_time.nc
-ncflint -w 1000.,0.0 ${dirtemp}/temp_time.nc ${dirtemp}/temp_time.nc ${dirtemp}/temp_int_${var}_${date}${month}.nc
-\rm ${dirtemp}/temp_time*
-\rm ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc
-done
-ncrcat ${dirtemp}/temp_int_${var}_${date}* ${TROPFLUX_ID}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc
-\rm ${dirtemp}/temp_int_${var}_${date}*
-done
-done
-
+#
 # e = Evaporation
 # slhf = ?
@@ -252 +296 @@
 # nsss = North/South surface stress
 for var in e slhf ssr trd sshf ssrd str ewss nsss ; do
-for date in 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010 ; do
-for month in 01 02 03 04 05 06 07 08 09 10 11 12; do
-for day in 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; do
-ncpdq -U ${dirin}/FC_SF/${date}/${month}/${var}.${date}${month}${day}* ${dirtemp}/${var}.${date}${month}${day}.temp.GLOBAL_075.nc
+   yyyy=${yearmin}
+   while [ ${yyyy} -le ${yearmax} ]
+   do
+      month=${monthmin}
+      while [ ${month} -le ${monthmax} ]
+      do
+         cmonth=$(printf "%2.2d" ${month})
+         daymin=1
+         daymax=$(cal ${month} ${yyyy} | grep . | fmt -1 | tail -1)
+         day=${daymin}
+         while [ ${day} -le ${daymax} ]
+         do
+            cday=$(printf "%2.2d" ${day})
+            if [ -f ${dirin}/FC_SF/${yyyy}/${cmonth}/${var}.${yyyy}${cmonth}${cday}* ] 
+            then
+               ncpdq -U ${dirin}/FC_SF/${yyyy}/${cmonth}/${var}.${yyyy}${cmonth}${cday}* ${dirtemp}/${var}.${yyyy}${cmonth}${cday}.temp.GLOBAL_075.nc
+            else
+               echo "${command} : iii : no files ${dirin}/FC_SF/${yyyy}/${cmonth}/${var}.${yyyy}${cmonth}${cday}*" >> ${log}
+            fi
+            day=$(( ${day} + 1 ))
+         done
+         exist_temp_files=$(find  ${dirtemp} -name "${var}.${yyyy}${cmonth}*temp*")
+         if [ "${exist_temp_files}" != "" ]
+         then
+            ncrcat ${dirtemp}/${var}.${yyyy}${cmonth}*temp*  ${dirtemp}/${var}.${yyyy}${cmonth}.fshei.GLOBAL_075.nc
+            rm ${dirtemp}/${var}.${yyyy}${cmonth}*temp* 2> /dev/null
+            ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,0,,5 ${dirtemp}/${var}.${yyyy}${cmonth}.fshei.GLOBAL_075.nc ${dirtemp}/temp_time0.nc
+            ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,4,,5 ${dirtemp}/${var}.${yyyy}${cmonth}.fshei.GLOBAL_075.nc ${dirtemp}/temp_time4.nc
+            ncflint -w 1.1574074,-1.1574074 ${dirtemp}/temp_time0.nc ${dirtemp}/temp_time4.nc ${dirtemp}/temp_time.nc
+            ncflint -w 0.00001,0.0 ${dirtemp}/temp_time.nc ${dirtemp}/temp_time.nc ${dirtemp}/temp_int_${var}_${yyyy}${cmonth}.nc
+            rm ${dirtemp}/temp_time* 2> /dev/null
+            rm ${dirtemp}/${var}.${yyyy}${cmonth}.fshei.GLOBAL_075.nc 2> /dev/null
+         else
+            echo "${command} : iii : no files ${dirtemp}/${var}.${yyyy}${cmonth}*temp*" >> ${log}
+         fi
+         unset exist_temp_files
+         month=$(( ${month} + 1 ))
+      done
+      ncrcat ${dirtemp}/temp_int_${var}_${yyyy}* ${TROPFLUX_ID}/erai_${reg}_1d_${yyyy}${cmonthmin}01_${yyyy}${cmonthmax}31_${var}_gridOrig.nc
+      rm ${dirtemp}/temp_int_${var}_${yyyy}* 2> /dev/null
+      yyyy=$(( ${yyyy} + 1 ))
+   done
 done
-ncrcat ${dirtemp}/${var}.${date}${month}*temp*  ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc
-\rm ${dirtemp}/${var}.${date}${month}*temp*
-ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,0,,5 ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc ${dirtemp}/temp_time0.nc
-ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,4,,5 ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc ${dirtemp}/temp_time4.nc
-ncflint -w 1.1574074,-1.1574074 ${dirtemp}/temp_time0.nc ${dirtemp}/temp_time4.nc ${dirtemp}/temp_time.nc
-ncflint -w 0.00001,0.0 ${dirtemp}/temp_time.nc ${dirtemp}/temp_time.nc ${dirtemp}/temp_int_${var}_${date}${month}.nc
-\rm ${dirtemp}/temp_time*
-\rm ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc
-done
-ncrcat ${dirtemp}/temp_int_${var}_${date}* ${TROPFLUX_ID}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc
-\rm ${dirtemp}/temp_int_${var}_${date}*
-done
-done
-
-
+#
 # geopt = Geopotentiel
 # msl = Mean sea level pressure
@@ -280 +348 @@
 # v10 = 10m V wind component
 for var in geopt msl sstk t2 d2 u10 v10 ; do
-for date in 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 2007 2008 2009 2010; do
-for month in 01 02 03 04 05 06 07 08 09 10 11 12; do
-
-ncpdq -U ${dirin}/AN_SF/${date}/${var}.${date}${month}.ashei.GLOBAL_075.nc ${dirtemp}/${var}.${date}${month}.ashei.GLOBAL_075.nc
-ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,0,,4 ${dirtemp}/${var}.${date}${month}.ashei.GLOBAL_075.nc ${dirtemp}/temp_time1.nc
-ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,1,,4 ${dirtemp}/${var}.${date}${month}.ashei.GLOBAL_075.nc ${dirtemp}/temp_time2.nc
-ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,2,,4 ${dirtemp}/${var}.${date}${month}.ashei.GLOBAL_075.nc ${dirtemp}/temp_time3.nc
-ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,3,,4 ${dirtemp}/${var}.${date}${month}.ashei.GLOBAL_075.nc ${dirtemp}/temp_time4.nc
-ncea ${dirtemp}/temp_time* ${dirtemp}/temp_int_${date}${month}.nc
-\rm ${dirtemp}/temp_time*
-\rm ${dirtemp}/${var}.${date}${month}.ashei.GLOBAL_075.nc
+   yyyy=${yearmin}
+   while [ ${yyyy} -le ${yearmax} ]
+   do
+      month=${monthmin}
+      while [ ${month} -le ${monthmax} ]
+      do
+         cmonth=$(printf "%2.2d" ${month})
+         ncpdq -U ${dirin}/AN_SF/${yyyy}/${var}.${yyyy}${cmonth}.ashei.GLOBAL_075.nc ${dirtemp}/${var}.${yyyy}${cmonth}.ashei.GLOBAL_075.nc
+         ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,0,,4 ${dirtemp}/${var}.${yyyy}${cmonth}.ashei.GLOBAL_075.nc ${dirtemp}/temp_time1.nc
+         ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,1,,4 ${dirtemp}/${var}.${yyyy}${cmonth}.ashei.GLOBAL_075.nc ${dirtemp}/temp_time2.nc
+         ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,2,,4 ${dirtemp}/${var}.${yyyy}${cmonth}.ashei.GLOBAL_075.nc ${dirtemp}/temp_time3.nc
+         ncks -d lat,${minlat},${maxlat} -d lon,${minlon},${maxlon} -d time,3,,4 ${dirtemp}/${var}.${yyyy}${cmonth}.ashei.GLOBAL_075.nc ${dirtemp}/temp_time4.nc
+         ncea ${dirtemp}/temp_time* ${dirtemp}/temp_int_${yyyy}${cmonth}.nc
+         rm ${dirtemp}/temp_time* 2> /dev/null
+         rm ${dirtemp}/${var}.${yyyy}${cmonth}.ashei.GLOBAL_075.nc 2> /dev/null
+         month=$(( ${month} + 1 ))
+      done
+      ncrcat ${dirtemp}/temp_int_${yyyy}* ${TROPFLUX_ID}/erai_${reg}_1d_${yyyy}${cmonthmin}01_${yyyy}${cmonthmax}31_${var}_gridOrig.nc
+      rm ${dirtemp}/temp_int_${yyyy}* 2> /dev/null
+      yyyy=$(( ${yyyy} + 1 ))
+   done
 done
-ncrcat ${dirtemp}/temp_int_${date}* ${TROPFLUX_ID}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc
-
-\rm ${dirtemp}/temp_int_${date}*
-done
-done
-#
-set
+#
 # end
 exit 0