Changeset 22 for trunk/src


Ignore:
Timestamp:
12/20/10 15:52:16 (13 years ago)
Author:
pinsard
Message:

remove most hard coded directories from compute_erai_daily_region_2d.sh

File:
1 edited

Legend:

Unmodified
Added
Removed

  • trunk/src/compute_erai_daily_region_2d.sh

    r11 r22  
    2020# DESCRIPTION 
    2121# =========== 
     22# 
     23# Put in ${TROPFLUX_ID} ERA-Intermin reference files 
     24# 
     25# Log file is written on :file:`${TROPFLUX_LOG}/compute_erai_daily_region_2d.log.{YYYYMMDDTHHMMSSZ}` 
    2226# 
    2327#     .. graphviz:: 
    … …  
    4650#           filein_v10 [shape=ellipse,fontname=Courier,label="/bdd/ERAI/NETCDF/GLOBAL_075/4xdaily/AN_SF/yyyy/mm/v10.yyyymm.ashei.GLOBAL_075.nc"]; 
    4751# 
    48 #           fileout_lsp [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_lsp_gridOrig.nc"]; 
    49 #           fileout_cp [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_cp_gridOrig.nc"]; 
    50 #           fileout_e [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_e_gridOrig.nc"]; 
    51 #           fileout_slhf [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_slhf_gridOrig.nc"]; 
    52 #           fileout_ssr [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_ssr_gridOrig.nc"]; 
    53 #           fileout_trd [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_trd_gridOrig.nc"]; 
    54 #           fileout_sshf [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_sshf_gridOrig.nc"]; 
    55 #           fileout_ssrd [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_ssrd_gridOrig.nc"]; 
    56 #           fileout_str [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_str_gridOrig.nc"]; 
    57 #           fileout_ewss [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_ewss_gridOrig.nc"]; 
    58 #           fileout_nsss [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_nsss_gridOrig.nc"]; 
    59 #           fileout_geopt [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_geopt_gridOrig.nc"]; 
    60 #           fileout_msl [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_msl_gridOrig.nc"]; 
    61 #           fileout_sstk [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_sstk_gridOrig.nc"]; 
    62 #           fileout_t2 [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_t2_gridOrig.nc"]; 
    63 #           fileout_d2 [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_d2_gridOrig.nc"]; 
    64 #           fileout_u10 [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_u10_gridOrig.nc"]; 
    65 #           fileout_v10 [shape=ellipse,fontname=Courier,label="/bdd/IPCC/Reanalysis/ERAI/${reg}/erai_GLOB_1d_yyyy0101_yyyy1231_v10_gridOrig.nc"]; 
     52#           fileout_lsp [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_lsp_gridOrig.nc"]; 
     53#           fileout_cp [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_cp_gridOrig.nc"]; 
     54#           fileout_e [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_e_gridOrig.nc"]; 
     55#           fileout_slhf [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_slhf_gridOrig.nc"]; 
     56#           fileout_ssr [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_ssr_gridOrig.nc"]; 
     57#           fileout_trd [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_trd_gridOrig.nc"]; 
     58#           fileout_sshf [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_sshf_gridOrig.nc"]; 
     59#           fileout_ssrd [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_ssrd_gridOrig.nc"]; 
     60#           fileout_str [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_str_gridOrig.nc"]; 
     61#           fileout_ewss [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_ewss_gridOrig.nc"]; 
     62#           fileout_nsss [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_nsss_gridOrig.nc"]; 
     63#           fileout_geopt [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_geopt_gridOrig.nc"]; 
     64#           fileout_msl [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_msl_gridOrig.nc"]; 
     65#           fileout_sstk [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_sstk_gridOrig.nc"]; 
     66#           fileout_t2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_t2_gridOrig.nc"]; 
     67#           fileout_d2 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_d2_gridOrig.nc"]; 
     68#           fileout_u10 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_u10_gridOrig.nc"]; 
     69#           fileout_v10 [shape=ellipse,fontname=Courier,label="${TROPFLUX_ID}/erai_GLOB_1d_yyyy0101_yyyy1231_v10_gridOrig.nc"]; 
    6670# 
    6771#           compute_erai_daily_region_2d [shape=box, 
    … …  
    7579#       } 
    7680# 
     81# EXAMPLES 
     82# ======== 
     83# 
     84# On climserv only if you don't have any ERA-I reference data, 
     85# you just have to run this tool :: 
     86# 
     87#  $ compute_erai_daily_region_2d.sh 
     88# 
     89# And look at log file with :: 
     90# 
     91#  $ tlogd.sh compute_erai_daily_region_2d 
     92# 
     93# and of course on files in ${TROPFLUX_ID}. 
     94 
     95# 
    7796# SEE ALSO 
    7897# ======== 
    … …  
    8099# :ref:`guide data ERA-I <data_in_erai>` 
    81100# 
     101# http://climserv.ipsl.polytechnique.fr/fr/les-donnees/era-interim-4.html 
     102# for variables names 
     103# 
     104# :ref:`tropflux_profile.sh` 
     105# 
    82106# TODO 
    83107# ==== 
    84108# 
    85 # trace 
     109# split get and process files 
    86110# 
    87111# hard coded directory - usage of ${TROPFLUX_ID} 
    88112# 
    89113# coding rules 
     114# 
     115# resolve missing long name variable for Netcdfname in http://climserv.ipsl.polytechnique.fr/fr/les-donnees/era-interim-4.html 
     116# 
     117#  - slhf 
     118#  - trd  
     119#  - sshf 
    90120# 
    91121# EVOLUTIONS 
    92122# ========== 
    93123# 
     124# - fplod 20101220T142107Z aedon.locean-ipsl.upmc.fr (Darwin) 
     125# 
     126#   * add Long name variables 
     127#   * replace /bdd/IPCC/Reanalysis/ERAI/${reg}/ by ${TROPFLUX_ID} 
     128#   * add examples 
     129#   * add log management 
     130#   * dirtemp not hard coded values 
     131#   * check dirin permission 
     132# 
    94133# - fplod 20101216T110626Z aedon.locean-ipsl.upmc.fr (Darwin) 
    95134# 
    … …  
    98137# 
    99138#- 
     139system=$(uname) 
     140case "${system}" in 
     141   AIX|IRIX64) 
     142      echo "www : no specific posix checking" 
     143   ;; 
     144   *) 
     145     set -o posix 
     146   ;; 
     147esac 
     148unset system 
     149# 
     150LANG=POSIX 
     151# 
     152command=$(basename ${0}) 
     153log_date=$(date -u +"%Y%m%dT%H%M%SZ") 
     154# 
     155usage=" Usage : ${command}" 
     156# 
     157hostname=$(hostname) 
     158# default 
     159# N.A. because no parameters 
     160# 
     161set -u 
     162# 
     163# check for ${TROPFLUX_LOG} definition 
     164if [ "${TROPFLUX_LOG}" = "" ] 
     165then 
     166   echo "${command} : eee : \${TROPFLUX_LOG} not defined" 
     167   exit 1 
     168fi 
     169# 
     170# check for ${TROPFLUX_LOG} existence 
     171if [ ! -d ${TROPFLUX_LOG} ] 
     172then 
     173   echo "${command} : eee : ${TROPFLUX_LOG} not found" 
     174   exit 1 
     175fi 
     176# 
     177# check for permission access on TROPFLUX_LOG 
     178if [ ! -x ${TROPFLUX_LOG} ] 
     179then 
     180   echo "${command} : eee : ${TROPFLUX_LOG} not reachable" 
     181   exit 1 
     182fi 
     183# 
     184# check for write permission on TROPFLUX_LOG 
     185if [ ! -w ${TROPFLUX_LOG} ] 
     186then 
     187   echo "${command} : eee : ${TROPFLUX_LOG} not writable" 
     188   exit 1 
     189fi 
     190# 
     191dirtemp=${TROPFLUX_LOG}/$(basename ${0} .sh).${log_date}.temp/ 
     192mkdir -p ${dirtemp} 
     193# 
     194log=${TROPFLUX_LOG}/$(basename ${0} .sh).log.${log_date} 
     195echo "[Context]" 1>> ${log} 
     196echo "command=$(basename ${0})" 1>>${log} 
     197echo "hostname=${hostname}" 1>> ${log} 
     198echo "runtime=${log_date}" 1>> ${log} 
     199unset log_date 
     200echo "" 1>> ${log} 
    100201# 
    101202minlat=-90. 
    … …  
    104205maxlon=360. 
    105206reg=GLOB 
    106 dirtemp=/home/mlengaig/ERAI/1d_2d_${reg}/ 
    107207dirin=/bdd/ERAI/NETCDF/GLOBAL_075/4xdaily/ 
    108 dirout=/bdd/IPCC/Reanalysis/ERAI/${reg}/ 
    109 mkdir ${dirout} 
    110 mkdir ${dirtemp} 
    111 \rm ${dirtemp}/* 
     208# check for dirin existence 
     209if [ ! -d ${dirin} ] 
     210then 
     211   echo "${command} : eee : ${dirin} not found" 
     212   exit 1 
     213fi 
     214# 
     215# check for permission access on dirin 
     216if [ ! -x ${dirin} ] 
     217then 
     218   echo "${command} : eee : ${dirin} not reachable" 
     219   exit 1 
     220fi 
     221# 
     222# lsp = Large scale precipitation 
     223# cp = Convective precipitation 
    112224for var in lsp cp; do 
    113225for 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 
    … …  
    125237\rm ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc 
    126238done 
    127 ncrcat ${dirtemp}/temp_int_${var}_${date}* /bdd/IPCC/Reanalysis/ERAI/${reg}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc 
     239ncrcat ${dirtemp}/temp_int_${var}_${date}* ${TROPFLUX_ID}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc 
    128240\rm ${dirtemp}/temp_int_${var}_${date}* 
    129241done 
    130242done 
    131243 
     244# e = Evaporation 
     245# slhf = ? 
     246# ssr = Surface solar radiation 
     247# trd = ? 
     248# sshf = ? 
     249# ssrd = Downward surface solar radiation 
     250# str = Surface thermal radiation 
     251# ewss = East /West surface stress 
     252# nsss = North/South surface stress 
    132253for var in e slhf ssr trd sshf ssrd str ewss nsss ; do 
    133254for 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 
    … …  
    145266\rm ${dirtemp}/${var}.${date}${month}.fshei.GLOBAL_075.nc 
    146267done 
    147 ncrcat ${dirtemp}/temp_int_${var}_${date}* /bdd/IPCC/Reanalysis/ERAI/${reg}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc 
     268ncrcat ${dirtemp}/temp_int_${var}_${date}* ${TROPFLUX_ID}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc 
    148269\rm ${dirtemp}/temp_int_${var}_${date}* 
    149270done 
    … …  
    151272 
    152273 
     274# geopt = Geopotentiel 
     275# msl = Mean sea level pressure 
     276# sstk = Sea surface temperature 
     277# t2 = Temperature at 2 meters 
     278# d2 = Dew point at 2 meters 
     279# u10 = 10m U wind component 
     280# v10 = 10m V wind component 
    153281for var in geopt msl sstk t2 d2 u10 v10 ; do 
    154282for 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 
    … …  
    164292\rm ${dirtemp}/${var}.${date}${month}.ashei.GLOBAL_075.nc 
    165293done 
    166 ncrcat ${dirtemp}/temp_int_${date}* /bdd/IPCC/Reanalysis/ERAI/${reg}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc 
     294ncrcat ${dirtemp}/temp_int_${date}* ${TROPFLUX_ID}/erai_${reg}_1d_${date}0101_${date}1231_${var}_gridOrig.nc 
    167295 
    168296\rm ${dirtemp}/temp_int_${date}* 
    169297done 
    170298done 
     299# 
     300set 
     301# end 
     302exit 0 
Note: See TracChangeset for help on using the changeset viewer.