Ignore:
Timestamp:
Dec 17, 2024, 2:02:51 PM (34 hours ago)
Author:
Laurent Fairhead
Message:

New version of LMDZ_Setup as rewritten by A. Barral

File:
1 edited

Legend:

Unmodified
Added
Removed
  • BOL/LMDZ_Setup/interp_aerosols.sh

    r4615 r5415  
    66
    77year=$1
    8 echo Dans interp_aerosols.sh year= $year
     8echo "Dans interp_aerosols.sh year= $year"
    99
    1010# Extraction d'un fichier de grille physique a partir de grilles_gcm.nc
     
    1212
    1313\rm -rf tmp*nc
    14 if [ 0 = 0 ] ; then
    15 imp1=`ncdump -h grilles_gcm.nc | grep lonv | head -1 | awk ' { print $3 } '`
    16 echo $imp1
    17 (( imm1 = $imp1 - 2 ))
    18 echo $imm1
    19 ncks -d lonv,0,$imm1 grilles_gcm.nc -v phis -O grille_phys.nc
    20 ncrename -v lonv,lon -v latu,lat -d lonv,lon -d latu,lat -O grille_phys.nc
    21 ncap2 -s "lon=-360.+lon" grille_phys.nc -O tmp.nc
    22 \mv -f tmp.nc grille_phys.nc
    23 #exit
     14if [[ 0 = 0 ]]; then
     15  imp1=$(ncdump -h grilles_gcm.nc | grep lonv | head -1 | awk ' { print $3 } ')
     16  echo "$imp1"
     17  (( imm1 = imp1 - 2 ))
     18  echo "$imm1"
     19  ncks -d "lonv,0,$imm1" grilles_gcm.nc -v phis -O grille_phys.nc
     20  ncrename -v lonv,lon -v latu,lat -d lonv,lon -d latu,lat -O grille_phys.nc
     21  ncap2 -s "lon=-360.+lon" grille_phys.nc -O tmp.nc
     22  \mv -f tmp.nc grille_phys.nc
     23  #exit 1
    2424fi
    2525
     
    2828
    2929orig=aerosols${year}_from_inca.nc
    30 ln -s $LMDZ_Init/$orig
     30ln -s "$LMDZ_INIT/$orig" .
    3131
    3232
    33 # Interpollation horizontale sur la grille cible
     33# Interpolation horizontale sur la grille cible
    3434# ----------------------------------------------
    3535
    3636# Un cdo direct ne marche pas. Il faut donner la liste des variables ...
    37 vars="" ; for var in `ncdump -h $orig | grep float | sed -e 's/^.*.float //' | cut -d'(' -f1 | sed -e 's/ ;//' -e /lat/d -e /lon/d` ; do  vars="$vars,$var" ; done
     37vars=""
     38for var in $(ncdump -h "$orig" | grep float | sed -e 's/^.*.float //' | cut -d'(' -f1 | sed -e 's/;//' -e /lat/d -e /lon/d); do
     39  vars="$vars,$var"
     40done
    3841
    39 echo cdo remapcon,grille_phys.nc -selvar$vars $orig tmp.nc
    40 cdo remapcon,grille_phys.nc -selvar$vars $orig tmp.nc
     42echo cdo remapcon,grille_phys.nc "-selvar$vars" "$orig" tmp.nc
     43cdo remapcon,grille_phys.nc "-selvar$vars" "$orig" tmp.nc
    4144#ncks -v ps,ap,b,presnivs $orig -A tmp.nc
    42 nccopy -k classic tmp.nc aerosols.${year}.nc
     45nccopy -k classic tmp.nc "aerosols.$year.nc"
Note: See TracChangeset for help on using the changeset viewer.