source: CONFIG/LMDZOR/branches/LMDZOR_v4/CREATE/SCRIPT/era2gcm.x @ 1190

Last change on this file since 1190 was 1190, checked in by jgipsl, 14 years ago

Modifications in config LMDZOR_v4 (config under developement) :

  • added resolution 280x280x39 with zoom Afric in LMDZOR experience
  • added guide.def in in LMDZOR experience
  • nudging files have to be added manualy in COMP/lmdz.card
  • Property svn:executable set to *
File size: 3.5 KB
Line 
1#!/bin/csh
2
3set mois=$1
4set an=$2
5set tmin=$3
6set tmax=$4
7set GRILLE_FI=$5
8set outsuf=_ERAI_$an-$mois.nc
9set climserv=local
10
11#####################################################################
12# Grille
13#####################################################################
14if ( ! -f $GRILLE_FI ) then
15    echo fichier start.nc must exist
16    exit
17endif
18
19set iip1=`ncdump -h $GRILLE_FI | grep rlonu | head -1 | awk ' { print $3 } '`
20set jjm=`ncdump -h $GRILLE_FI | grep rlatv | head -1 | awk ' { print $3 } '`
21set jjp1=`expr $jjm + 1`
22set iim=`expr $iip1 - 1`
23
24echo iim=$iim jjm=$jjm
25
26###################################################################
27# Choix de la periode temporelle
28###################################################################
29set t0="l=$tmin"
30set t1tn="l=${tmin}:${tmax}"
31
32###################################################################
33# Fichiers donnees
34###################################################################
35if ( 1 == 0 ) then #ncep
36
37    set varu=ua
38    set varv=va
39    set vart=ta
40    set varq=hur
41    set varp=ps
42   
43    set fileu="$NCEP_DIR/ua.nc"
44    set filev="$NCEP_DIR/va.nc"
45    set filet="$NCEP_DIR/ta.nc"
46    set fileq="$NCEP_DIR/hur.nc"
47    set filep="$NCEP_DIR/ps.nc"
48
49else # ERAI
50
51    set suf=ei.GLOBAL_075
52
53    if ( $climserv == poly ) then
54        set ANA_DIR='http://dods-ipsl.ipsl.polytechnique.fr/ipsl-bin/nph-dods/ERAI/NETCDF/GLOBAL_075/4xdaily'
55    else
56        set ANA_DIR='http://dodsp.idris.fr:81/cgi-bin/nph-dods/rpsl376/ERAI/NETCDF/GLOBAL_075/4xdaily'
57# at ulam: set ANA_DIR='/u/rech/psl/rpsl376/ERAI/NETCDF/GLOBAL_075/4xdaily'
58    endif
59
60    set vart=ta
61    set varu=u
62    set varv=v
63    set varq=r
64    set varp=msl
65    set filet="$ANA_DIR/AN_PL/$an/$vart.$an$mois.aph$suf.nc"
66    set fileu="$ANA_DIR/AN_PL/$an/$varu.$an$mois.aph$suf.nc"
67    set filev="$ANA_DIR/AN_PL/$an/$varv.$an$mois.aph$suf.nc"
68    set fileq="$ANA_DIR/AN_PL/$an/$varq.$an$mois.aph$suf.nc"
69    set filep="$ANA_DIR/AN_SF/$an/$varp.$an$mois.ash$suf.nc"
70
71endif
72
73###################################################################
74# script ferret pour interpolation
75###################################################################
76
77cat <<eod> ! tmp.jnl
78set memory/size=50
79use "$GRILLE_FI"
80use "$fileu"
81use "$filev"
82use "$filet"
83use "$fileq"
84use "$filep"
85
86let uwnd=$varu
87let vwnd=$varv
88let tempair='T'
89let air=$vart
90let rh=$varq
91let sp=$varp
92let year=$an
93
94define axis/t=1-jan-${an}:31-dec-${an}:6/units=hours thour
95define grid/like=uwnd[d=2]/x=cu[d=1]/y=cu[d=1] grille_u
96define grid/like=vwnd[d=3]/x=cv[d=1]/y=cv[d=1] grille_v
97define grid/like=air[d=4]/x=cv[d=1]/y=cu[d=1] grille_T
98
99!save/clobber/file=T$outsuf air[d=4,g=grille_T,i=1:$iip1,j=1:$jjp1,$t0,gt=thour@asn]
100!repeat/$t1tn save/file="T$outsuf"/append air[d=4,g=grille_T,i=1:$iip1,j=1:$jjp1,gt=thour@asn]
101
102save/clobber/file=u$outsuf uwnd[d=2,g=grille_u,i=1:$iip1,j=1:$jjp1,$t0,gt=thour@asn]
103repeat/$t1tn save/file="u$outsuf"/append uwnd[d=2,g=grille_u,i=1:$iip1,j=1:$jjp1,gt=thour@asn]
104
105save/clobber/file=v$outsuf vwnd[d=3,g=grille_v,i=1:$iip1,j=1:$jjm,$t0,gt=thour@asn]
106repeat/$t1tn save/file="v$outsuf"/append vwnd[d=3,g=grille_v,i=1:$iip1,j=1:$jjm,gt=thour@asn]
107
108!save/clobber/file=hur$outsuf rh[d=5,g=grille_T,i=1:$iip1,j=1:$jjp1,$t0,gt=thour@asn]
109!repeat/$t1tn save/file="hur$outsuf"/append rh[d=5,g=grille_T,i=1:$iip1,j=1:$jjp1,gt=thour@asn]
110
111!save/clobber/file=ps$outsuf sp[d=6,g=grille_T,i=1:$iip1,j=1:$jjp1,$t0,gt=thour@asn]
112!repeat/$t1tn save/file="ps$outsuf"/append sp[d=6,g=grille_T,i=1:$iip1,j=1:$jjp1,gt=thour@asn]
113
114eod
115
116ferret <<eod > & /dev/null
117go tmp.jnl
118quit
119eod
Note: See TracBrowser for help on using the repository browser.