source: BOL/Multi_atlas/Utils/multi.sh @ 4479

Last change on this file since 4479 was 4439, checked in by musat, 18 months ago

Deux corrections:
1/ pour les cas où l'on renomme une simulation (i.e. le nom en 1ere colonne
du def.txt n'est pas le même que celui de la simulation)
2/ de la période, pour les cas où l'on utilise des fichiers MO plutôt que
des SE
IonelaMusat?
  • Property svn:executable set to *
File size: 21.7 KB
Line 
1#set -vx
2
3source ~/env_Multi_atlas.sh
4
5local=`pwd -P`
6cd $local
7echo REPERTOIRE LOCAL :  $local
8
9echo '############################################################'
10echo  0/ gestion des options
11echo '############################################################'
12pcmdi_metrics=0
13tuning_metrics=1
14und=0
15zon=1
16html=1
17atlas=1
18#newatlas=0
19atlas_cesmep=0
20get=0
21cmor=1
22regdyn=0
23lma=0
24trmm=0
25axe2=0
26axe4=0
27esmval2=0
28reflcld_histo=0
29vertcld_calipso=0
30stations_day=0
31stations_mth=0
32qsub=0
33mjo=0
34river=0
35dataref=1
36
37ref=OBS
38echo ref=$ref
39
40lonmin=-180
41lonmax=180
42latmin=-90
43latmax=90
44pmin=1
45
46# ./job_multi.sh -1D 0 -get 0 -atlas 0 -pcmdi 0 -tuning 0 -lma 0 -trmm 0
47runstxt=$local/runs.txt
48
49while test -n "${1}"; do
50   case $1 in
51      -1D) und=$2 ; shift ;;
52      -get) get=$2 ; shift ;;
53      -cmor) cmor=$2 ; shift ;;
54      -zon) zon=$2 ; shift ;;
55      -atlas) atlas=$2 ; shift ;;
56      -atlas_cesmep) atlas_cesmep=$2 ; shift ;;
57      -pcmdi) pcmdi_metrics=$2 ; shift ;;
58      -tuning) tuning_metrics=$2 ; shift ;;
59      -runstxt) runstxt=$2 ; shift ;;
60      -lma) lma=$2 ; shift ;;
61      -trmm) trmm=$2 ; shift ;;
62      -axe4) axe4=$2 ; shift ;;
63      -axe2) axe2=$2 ; shift ;;
64      -regdyn) regdyn=$2 ; shift ;;
65      -esmval2) esmval2=$2 ; shift ;; 
66      -stations_day) stations_day=$2 ; shift ;; 
67      -stations_mth) stations_mth=$2 ; shift ;; 
68      -reflcld_histo) refld_histo=$2 ; shift ;;
69      -vertcld_calipso) vertcld_calipso=$2 ; shift ;;
70      -mjo) mjo=$2 ; shift ;; 
71      -river) river=$2 ; shift ;;
72      -dataref) dataref=$2 ; shift ;;
73      -lonmin) lonmin=$2 ; shift ;;
74      -lonmax) lonmax=$2 ; shift ;;
75      -latmin) latmin=$2 ; shift ;;
76      -latmax) latmax=$2 ; shift ;;
77      -pmin) pmin=$2 ; shift ;;
78      -qsub) qsub=$2 ; shift ;;
79      *) comp=$1 ;;
80   esac
81   shift
82done
83
84echo multi=$comp region: lonmin=$lonmin lonmax=$lonmax latmin=$latmin latmax=$latmax
85echo multi=$comp pmin: pmin=$pmin
86rm -f out.com$$
87runstxt=$local/runs.txt
88echo runstxt= $runstxt
89
90echo qsub is $qsub
91#exit
92
93echo DODSDIR : $DODSDIR
94MAIN_SE=$DODSDIR/$login/lmdz/SE
95TS_DA=$DODSDIR/$login/lmdz/TS_DA
96TS_MO=$DODSDIR/$login/lmdz/TS_MO
97MULTIDIR=$DODSDIR/$login/lmdz/MultiSimu
98
99if [ ! -d $MULTIDIR ] ; then mkdir -p $MULTIDIR ; fi
100if [ ! -d $MAIN_SE ] ; then mkdir -p $MAIN_SE ; fi
101if [ ! -d $MAIN_MO ] ; then mkdir -p $MAIN_MO ; fi
102if [ ! -d $MAIN_DA ] ; then mkdir -p $MAIN_DA ; fi
103
104COMP_D=$MULTIDIR/$comp
105DEF_FILE=$COMP_D/def.txt ; if [ ! -f $DEF_FILE ] ; then echo renseigner $DEF_FILE ; exit ; fi
106LOGD=$COMP_D/log ; mkdir -p $LOGD
107outfile=$LOGD/out$$ ; if [ -f ${outfile} ] ; then \rm -f $outfile ; fi
108echo Ecriture des logs dans $outfile
109
110###
111#~/LMDZ/check/make_warning.sh $COMP_D
112### nom2
113listsims=""  ; for s  in `awk ' {print $1} ' $DEF_FILE | uniq`          ; do listsims="$listsims $s"    ; done
114#echo listsims $listsims
115
116listsims=""  ; for s  in `awk ' {print $1} ' $DEF_FILE | uniq`          ; do listsims="$listsims $s"    ; done
117echo listsims $listsims
118### nom1, simdir
119for s in $listsims ; do
120### simdir
121simdir=`grep -w "^$s "  $runstxt | awk ' { print $2 } '`
122### nom1
123ss=`basename $simdir `
124for year in `grep -w "^$s" $DEF_FILE | awk ' { print $2 } '` ; do
125file=${MAIN_SE}/ORIG/${s}_SE_${year}_1M_histmth.nc
126###file=${MAIN_SE}/ORIG/${ss}_SE_${year}_1M_histmth.nc
127
128echo $file $year
129
130if [ ! -f ${file} ]; then
131 if [ -f ${simdir}/ATM/Analyse/SE/${ss}_SE_${year}_1M_histmth.nc ]; then
132  ln -s ${simdir}/ATM/Analyse/SE/${ss}_SE_${year}_1M_histmth.nc $file
133 elif [ "${ss}" == "UQLMDZ6A1000" ]; then
134  echo NO need to ln -s; the *1M_histmth.nc does not exist anymore
135 else
136  echo $year>>year.txt
137  sed -i -e 's/_/ /' year.txt
138  year1=`less year.txt | awk '{print $1}'`
139  year2=`less year.txt | awk '{print $2}'`
140  ln -s ${simdir}/ATM/Output/MO/${ss}_${year1}0101_${year2}1230_1M_histmth.nc $file
141 fi
142fi
143
144\rm ${MAIN_SE}/ORIG/UQLMDZ6AW1000_SE_1979_1988_1M_histmth.nc
145done
146done
147#${local}/../check/mk_check-avail-data.sh ${comp}
148
149echo ### On supprime dans def.txt les lignes vides ########
150sed -i "/^[ \t]*$/d" $DEF_FILE
151
152echo ### On teste d abord la coherence des noms des simulations dans def.txt et $runstxt ################
153echo MULTI-ATLAS "$comp"
154for s in $listsims ; do
155srun=`grep -w "^$s "  $runstxt | awk ' { print $1 } '`
156echo 's='"$s" '    srun='  "$srun"
157if [ "$s" != "$srun" ] ; then
158   #echo Meme nom de la simulation dans def.txt et $runstxt 's=' "$s" '    srun='  "$srun"
159#else
160   echo ATTENTION : la simulation "$s" de votre def.txt n est pas repertoriee dans $runstxt >${LOGD}/error_${comp}
161   echo             ou le meme nom est attribue a plusieures simulations >> ${LOGD}/error_${comp}
162   echo Il faut la rajouter dans $runstxt >> ${LOGD}/error_${comp}
163   echo avec le chemin vers les SE et TS dans ciclad : /thredds/ipsl/fabric/IGCM_OUT >> ${LOGD}/error_${comp}
164   echo             ou renomer votre simulation >> ${LOGD}/error_${comp}
165   more ${LOGD}/error_${comp}
166exit
167fi
168echo Checking if $s simulation is referenced many times in runs.txt
169#nrun=`./how_many_times_run.sh $s`
170#nrun=`./how_many_times_run.sh $srun`
171nrun=1
172
173if [ "$nrun" != "1" ]; then
174 echo Simulation $s est repertoriee ${nrun} fois dans $runstxt >> ${LOGD}/error_${comp}
175 echo STOPping here >> ${LOGD}/error_${comp}
176 more ${LOGD}/error_${comp}
177 exit
178fi
179done
180#exit
181listruns=""
182listrunss=""
183listrunse=""
184if [ ! -d $MAIN_SE/ORIG ] ; then mkdir -p $MAIN_SE/ORIG ; fi
185for s in $listsims ; do 
186simdir=`grep -w "^$s "  $runstxt | awk ' { print $2 } '`
187echo $simdir
188( cd $MAIN_SE/ORIG ; ln -sf $simdir/ATM/Analyse/SE/*mth*nc . ) > /dev/null 2>&1
189echo $simdir
190ss=`basename $simdir `
191listruns="$listruns $ss"
192for year in `grep -w "^$s" $DEF_FILE | awk ' { print $2 } '` ; do
193listrunss="$listrunss ${ss}_$year"
194listrune="$listrune ${ss}_SE_$year"
195done
196done
197#echo listruns = $listruns
198#echo listrunss= $listrunss
199#echo listrune= $listrune
200#ls -lrt $MAIN_SE/ORIG # Eviter de laisser des ls qui rendent la sortire illisible FH
201#listsimss="" ; for ss in `awk ' {print $1"_"$2 } ' $DEF_FILE`    ; do listsimss="$listsimss $ss" ; done
202
203ns=0
204listsimss="" ; 
205for ss in `awk ' {print $1"_"$2 } ' $DEF_FILE`    ; do 
206if [ $dataref = 0 ] ; then
207   if [ ${ns} = 0 ] ; then
208      ref=$ss
209      echo ref=$ref
210   fi 
211fi
212let "ns = ${ns} + 1"
213listsimss="$listsimss $ss" ; 
214done
215echo ref=$ref
216
217listsimse="" ; for se in `awk ' {print $1"_SE_"$2 } ' $DEF_FILE` ; do listsimse="$listsimse $se" ; done
218echo listsims $listsims
219echo listsimss $listsimss
220echo listsimse $listsimse
221
222echo SIMULATIONS: $listsims
223
224RUNS=""
225RUNS=`awk ' {print $1"_"$2 } ' $DEF_FILE`
226echo listrunss=$RUNS
227
228echo '########################################################################'
229echo  1/ Rapatriement des simulations manquantes
230echo '########################################################################'
231#### AI avril 2020 ################ Avant l acces directe aux simulations sur le dods du CCRT
232# Adresses web ne sont plus a jours
233# A commenter les adresses web CCRT
234# A mettre a jours les adress web IDRIS
235# Rajouter d'autres machines avec rsynk, wget ou ssh ...
236
237if [ $get = 1 ] ; then
238
239cd $MAIN_SE
240for sim in $listsimss ; do echo ./getun.sh ${sim} >> $outfile 2>&1 ; ./getun.sh ${sim} ; done
241cd $TS_DA
242for sim in $listsimss ; do echo ./getun.sh ${sim} >> $outfile 2>&1 ; ./getun.sh ${sim} ; done
243cd $TS_MO
244for sim in $listsimss ; do echo ./getun.sh ${sim} >> $outfile 2>&1 ; ./getun.sh ${sim} ; done
245
246fi
247
248echo '########################################################################'
249echo  1.1/ Check des fichiers input et des variables demandees pour l atlas avec
250echo       des exit si un des fichiers ou des variables ne sont pas accessibles
251echo '########################################################################'
252#IM TBC ${local}/../check/make_check_fin.sh $comp ${LOGD}/check_finput$$ 2>&1
253
254if [ $cmor = 1 ] ; then
255 cp cmor.sh $MAIN_SE/.
256 cp season.sh $MAIN_SE/.
257 cd $MAIN_SE
258 chmod +x cmor.sh
259 chmod +x season.sh
260#########################################################################
261# 1.2 "Cmorisation des fichiers mensuels : traduction au format CMIP
262#########################################################################
263# Deux lignes suivantes modif FH 2016/11/09
264# Files are taken from /thredds/ipsl/fabric/lmdz/SE/ORIG and
265# results are in /thredds/ipsl/fabric/lmdz/SE/CMOR
266#########################################################################
267
268echo cmor.sh : output dans ${LOGD}/cmor$$
269./cmor.sh -runsfile $runstxt $comp > ${LOGD}/cmor$$ 2>&1
270if [ $? != 0 ] ; then exit ; fi # Probleme de histmth non present. On arrete tout
271
272#########################################################################
273# 1.3 Interpolate fields on a common horizontal grid (VLR) and
274#########################################################################
275#     average in time "YEAR", "JJA" ...
276#     Results are put on
277#     /thredds/ipsl/fabric/lmdz/VLR/YEAR/[SIMU_NAME]_[period]/NC
278#     If wanting to rerun season.sh for observations :
279#     ./season.sh -s OBS
280#########################################################################
281
282cd $MAIN_SE
283echo season.sh : output dans ${LOGD}/season$$
284./season.sh -comp $comp > ${LOGD}/season$$ 2>&1
285cd $local
286
287fi
288
289echo '########################################################################'
290echo  2/ Lancement des atlas manquant
291echo '########################################################################'
292natlas=`wc -l $COMP_D/atlas-def | awk ' { print $1 } '`
293count_atlas=0
294
295if [ $atlas = 1 ] ; then
296#if [ $newatlas = 1 ]; then
297#sed -e 's/COMPAR/'${comp}'/' ${local}/atlas/make_atlas.sh > ${local}/atlas/make_atlas_${comp}.sh
298#$submit -v COMP=${comp},DATAREF=${dataref} ${local}/atlas/make_atlas_${comp}.sh
299
300#./make_timecount.sh ${comp}
301
302#else
303
304for sim in $listsimss ; do
305   il=1
306   if [ ! -f $COMP_D/atlas-def ] ; then
307      echo Definir les atlas a construire dans $COMP_D/atlas-def
308      echo Example :
309      cat $COMP_D/../CMIP5to6/atlas-def
310   fi
311   while [ $il -le $natlas ] ; do
312      line=`sed -n -e ${il}p $COMP_D/atlas-def`
313      seas=`echo $line | awk ' { print $1 } '`
314      proj=`echo $line | awk ' { print $2 } '`
315      vars=`echo $line | awk ' { print $3 } '`
316      atlas=${sim}+${seas}+${proj}+$vars
317      cd ${local}/atlas
318      varg=`echo $vars | sed -e 's/,/+/g'`
319      echo Atlas $atlas a creer
320      if [ $qsub = 0 ] ; then
321         echo bash ${local}/atlas/atlas.sh RUN=$sim,SEASON=$seas,PROJ=$proj,OUTDIR=$COMP_D/ATLAS,REF=$ref,VARS="$varg",RUNSTXT=$runstxt,lonmin=$lonmin,lonmax=$lonmax,latmin=$latmin,latmax=$latmax,pmin=$pmin
322         bash ${local}/atlas/atlas.sh RUN=$sim,SEASON=$seas,PROJ=$proj,OUTDIR=$COMP_D/ATLAS,REF=$ref,VARS="$varg",RUNSTXT=$runstxt,lonmin=$lonmin,lonmax=$lonmax,latmin=$latmin,latmax=$latmax,pmin=$pmin >> $outfile
323         if [ "$ref" != "OBS" ]; then
324            echo bash ${local}/atlas/atlas.sh RUN=$sim,SEASON=$seas,PROJ=$proj,OUTDIR=$COMP_D/ATLAS,REF=OBS,VARS="$varg",RUNSTXT=$runstxt,lonmin=$lonmin,lonmax=$lonmax,latmin=$latmin,latmax=$latmax,pmin=$pmin
325            bash ${local}/atlas/atlas.sh RUN=$sim,SEASON=$seas,PROJ=$proj,OUTDIR=$COMP_D/ATLAS,REF=OBS,VARS="$varg",RUNSTXT=$runstxt,lonmin=$lonmin,lonmax=$lonmax,latmin=$latmin,latmax=$latmax,pmin=$pmin >> $outfile
326         fi
327      else
328         echo $submit -v RUN=$sim,SEASON=$seas,PROJ=$proj,OUTDIR=$COMP_D/ATLAS,REF=$ref,VARS="$varg",RUNSTXT=$runstxt,lonmin=$lonmin,lonmax=$lonmax,latmin=$latmin,latmax=$latmax,pmin=$pmin atlas.sh >> ${LOGD}/ATLAS$$ 2>&1
329         $submit -v RUN=$sim,SEASON=$seas,PROJ=$proj,OUTDIR=$COMP_D/ATLAS,REF=$ref,VARS="$varg",RUNSTXT=$runstxt,lonmin=$lonmin,lonmax=$lonmax,latmin=$latmin,latmax=$latmax,pmin=$pmin atlas.sh
330         if [ "$ref" != "OBS" ]; then
331            echo $submit -v RUN=$sim,SEASON=$seas,PROJ=$proj,OUTDIR=$COMP_D/ATLAS,REF=OBS,VARS="$varg",RUNSTXT=$runstxt,lonmin=$lonmin,lonmax=$lonmax,latmin=$latmin,latmax=$latmax,pmin=$pmin atlas.sh >> ${LOGD}/ATLAS$$ 2>&1
332            $submit -v RUN=$sim,SEASON=$seas,PROJ=$proj,OUTDIR=$COMP_D/ATLAS,REF=OBS,VARS="$varg",RUNSTXT=$runstxt,lonmin=$lonmin,lonmax=$lonmax,latmin=$latmin,latmax=$latmax,pmin=$pmin atlas.sh
333         fi
334      fi
335      (( count_atlas = $count_atlas + 1 ))
336      (( il = $il + 1 ))
337   done
338done
339if [ $count_atlas != 0 ] ; then
340    echo '!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
341    echo Les atlas manquant ont ete lances. Attendre leur
342    echo creation pour relancer job_multi.sh
343    echo '!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!'
344#   exit
345fi
346
347#fi #newatlas
348fi #atlas
349
350# atlas_cesmep : cles pour activer les nouveau atlas (C-ESM-EP)
351# AI avril 2020 : Passage a la version v2 de C-ESM-EP
352# FH Juin 2022 : Passage a la version v3 de C-ESM-EP (changement de prodigfs en thredds)
353if [ $atlas_cesmep = 1 ] ; then
354   cd ~/LMDZ/C-ESM-EP
355   ./AtlasLMDZ.sh -runsfile $runstxt $comp >> ${LOGD}/CESMEP$$ 2>&1
356fi # atlas_cesmep
357
358echo '########################################################################'
359echo  3/ Calcul des moyennes zonales nco
360echo '#####################################################################'
361
362
363#trace moy.zonale OBS + multi-simulations
364if [ $zon = 1 ] ; then
365   cd $local
366   mkdir -p $COMP_D/ZON
367   for sim in $listsimse OBS ; do
368      echo make_MOY for simulation $sim, log in ${LOGD}/MOY$$
369      years=`echo $sim | sed -e 's/_SE_/ /' | awk ' { print $2 } '`
370      run=`echo $sim | sed -e 's/_SE_/ /' | awk ' { print $1 } '`
371      echo ./make_MOY.sh "$years" YEAR "$run" GLOB
372      ./make_MOY.sh "$years" YEAR "$run" GLOB      >> ${LOGD}/MOY$$ 2>&1
373      ./make_MOY.sh "$years" YEAR "$run" ZON       >> ${LOGD}/MOY$$ 2>&1
374      ./make_MOY.sh "$years" JJA "$run"  ZON       >> ${LOGD}/MOY$$ 2>&1
375      ./make_MOY.sh "$years" DJF "$run"  ZON       >> ${LOGD}/MOY$$ 2>&1
376      ./make_MOY.sh "$years" JJAS "$run" AMMACROSS >> ${LOGD}/MOY$$ 2>&1
377      ./make_MOY.sh "$years" MAM "$run"  ZON       >> ${LOGD}/MOY$$ 2>&1
378      ./make_MOY.sh "$years" SON "$run"  ZON       >> ${LOGD}/MOY$$ 2>&1
379   done
380   cd $local
381### ajout nuages LMDZ dans la compa et doit faire des liens symbol. sur
382#### les OBS calipso toujurs pour le moment
383   echo mk_clouds_NOcalipso.sh, log in  ${LOGD}/calipso$$
384   ./mk_clouds_NOcalipso.sh YEAR ZON       >> ${LOGD}/calipso$$ 2>&1
385   ./mk_clouds_NOcalipso.sh DJF ZON        >> ${LOGD}/calipso$$ 2>&1
386   ./mk_clouds_NOcalipso.sh JJA ZON        >> ${LOGD}/calipso$$ 2>&1
387   ./mk_clouds_NOcalipso.sh JJAS AMMACROSS >> ${LOGD}/calipso$$ 2>&1
388   ./mk_clouds_NOcalipso.sh MAM ZON        >> ${LOGD}/calipso$$ 2>&1
389   ./mk_clouds_NOcalipso.sh SON ZON        >> ${LOGD}/calipso$$ 2>&1
390   echo multi_ZONE.sh, log in ${LOGD}/Z$$
391   ./multi_ZONE.sh YEAR $comp ZON -runsfile $runstxt >> ${LOGD}/Z$$ 2>&1
392   ./multi_ZONE.sh JJA $comp ZON -runsfile $runstxt  >> ${LOGD}/Z$$ 2>&1
393   ./multi_ZONE.sh DJF $comp ZON -runsfile $runstxt  >> ${LOGD}/Z$$ 2>&1
394   ./multi_ZONE.sh JJAS $comp AMMACROSS -runsfile $runstxt >> ${LOGD}/Z$$ 2>&1
395   ./multi_ZONE.sh MAM $comp ZON -runsfile $runstxt  >> ${LOGD}/Z$$ 2>&1
396   ./multi_ZONE.sh SON $comp ZON -runsfile $runstxt  >> ${LOGD}/Z$$ 2>&1
397fi
398
399if [ ! -f $COMP_D/1D/1D.html -a $und = 1  ] ; then
400echo '#####################################################################'
401echo 4/ Lancement du 1D
402echo '#####################################################################'
403    cd $COMP_D
404    echo log in ${LOGD}/1D$$
405    ~/LMDZ/1D/serie1d.sh $comp  > ${LOGD}/1D$$ 2>&1
406    cd $local
407fi
408
409if [ $tuning_metrics = 1 -a ! -d $COMP_D/METRICS ] ; then
410echo '#####################################################################'
411echo 5.1/ Metriques LMD : $tuning_metrics
412echo '#####################################################################'
413   echo metrics.sh, log in ${LOGD}/METRICS$$
414   cd $local/METRICS
415   #./metrics.sh -runsfile $runstxt $comp
416   echo ./metrics.sh $comp > ${LOGD}/METRICS$$ 2>&1
417   ./metrics.sh $comp > ${LOGD}/METRICS$$ 2>&1
418fi
419
420if [ $pcmdi_metrics = 1 ] ; then
421echo '#####################################################################'
422echo 5.2/ Metriques PCMDI : $pcmdi_metrics
423echo '#####################################################################'
424cd $local
425   ./make_pcmdi.sh $comp
426fi
427
428if [ $regdyn = 1 ]; then
429echo '#####################################################################'
430echo 6/ Sorting by dynamical regimes : $regdyn
431echo '#####################################################################'
432   cd ~/users/musat/regdyn/clean_regdyn
433   ./mk_all_regdyn.sh $comp 0 > ${LOGD}/REGDYN$$ 2>&1
434   cd $local
435fi
436
437if [ $esmval2 = 1 ]; then
438echo '#####################################################################'
439echo 7/ ESMvalTool : $esmval2
440echo '#####################################################################'
441   cd $local
442   ./cvdp4atlas.sh $comp
443   ./perfm4atlas.sh $comp
444fi   
445
446if [ $mjo = 1 ]; then
447echo '#####################################################################'
448echo 8.1/ Rainfall variability / MJO diagnostics : $mjo
449echo '#####################################################################'
450   cd ~/users/idelkadi/MJO
451   ./make_diagnosMJO.sh $comp
452fi
453
454if [ $river = 1 ]; then 
455   cd ~/users/idelkadi/AXE1/River
456   ./make_diagnoRiver.sh $comp
457fi 
458
459if [ $lma = 1 ]; then
460echo '#####################################################################'
461echo 8.2/ Rainfall variability / LMA analysis of the MJO : $lma
462echo '#####################################################################'
463   #IM150722 cd $local
464   #IM120722 ./make_lma.sh $comp
465   cd ~/users/musat/d_AXE2_LMA/
466   ./make_LMA.sh $comp >  ${LOGD}/LMA$$ 2>&1
467   cd ${local}
468fi
469
470if [ $axe2 = 1 ]; then
471echo '#####################################################################'
472echo 8.3/ Rainfall variability / metrics
473echo '#####################################################################'
474   cd $local
475   bash  make_axe2.sh -runstxt $runstxt $comp    >  ${LOGD}/axe2_$$ 2>&1
476   echo  make_axe2.sh -runstxt $runstxt $comp '-->' ${LOGD}/axe2_$$
477fi
478
479if [ $trmm = 1 ]; then
480echo '#####################################################################'
481echo 8.4/ Rainfall variability / Marine Bonazzola : $trmm
482echo '#####################################################################'
483   cd ~/users/musat/atlas/d_AXE2
484   ./Axe2_Marine.sh $comp
485   cd $local
486fi
487
488if [ $river = 1 ]; then 
489echo '#####################################################################'
490echo 9/ River discharges
491echo '#####################################################################'
492   cd ~/users/idelkadi/AXE1/River
493   ./make_diagnoRiver.sh $comp
494fi 
495
496if [ $axe4 = 1 ]; then
497echo '#####################################################################'
498echo 10/ Poles
499echo '#####################################################################'
500  cd $local
501   ./make_axe4.sh $comp
502fi
503
504if [ "$reflcld_histo" = "1" ] ; then
505echo '#####################################################################'
506echo 11/ reflcld_histo : $reflcld_histo
507echo '#####################################################################'
508   cd ~/users/idelkadi/AXE3/CldReflCldFrac
509   ./make_scatReflCld.sh $comp
510   cd $local
511fi
512
513if [ $vertcld_calipso = 1 ] ; then
514echo '#####################################################################'
515echo 12/ Calipso : $vertcld_calipso
516echo '#####################################################################'
517   cd ~/users/idelkadi/AXE3/Calipso/CoupeVerticale/DJF
518   ./make_VerticalMzCldCalipso.sh $comp
519   cd ~/users/idelkadi/AXE3/Calipso/CoupeVerticale/JJA
520   ./make_VerticalMzCldCalipso.sh $comp
521   cd ~/users/idelkadi/AXE3/Calipso/CoupeVerticale/MAM
522   ./make_VerticalMzCldCalipso.sh $comp
523   cd ~/users/idelkadi/AXE3/Calipso/CoupeVerticale/SON
524   ./make_VerticalMzCldCalipso.sh $comp
525   cd ~/users/idelkadi/AXE3/Calipso/Stratos
526   ./make_StratosCalipso.sh $comp
527fi
528
529
530if [ $stations_day = 1 ] ; then
531echo '#####################################################################'
532echo 13/ Stations
533echo '#####################################################################'
534   ~/users/Binta/bibi.sh $comp
535   ~/users/cheruy/bibi.sh $comp
536fi
537
538if [ $stations_mth = 1 ] ; then
539   ~/users/traore/SE_buoy.sh $comp
540fi
541
542if [ $html = 1 ] ; then
543echo '#####################################################################'
544echo 14/ generation des pages html
545echo '#####################################################################'
546   cd ${local}
547   echo ./make_html_libigcm.sh --comp $comp --runstxt $runstxt
548   bash ./make_html_libigcm.sh --comp $comp --runstxt $runstxt > ${LOGD}/html$$ 2>&1
549   echo log in ${LOGD}/html$$ 
550   bash ./html_entete.sh $comp ${ref} >> ${LOGD}/html$$ 2>&1
551   cd ${local}/html
552   bash ./make_index_html.sh ${comp} ${dataref} >> ${LOGD}/html$$ 2>&1
553   cd $local
554   bash ./html_tuning.sh  $comp           >> ${LOGD}/html$$ 2>&1
555   bash ./html_tuning_min.sh  $comp       >> ${LOGD}/html$$ 2>&1
556   bash ./html_accueil.sh  $comp          >> ${LOGD}/html$$ 2>&1
557   bash ./html_1D.sh  $comp               >> ${LOGD}/html$$ 2>&1
558   bash ./html_lma.sh $comp               >> ${LOGD}/html_LMA$$ 2>&1
559   bash ./html_pcmdi.sh $comp             >> ${LOGD}/html$$ 2>&1
560   bash ./html_all.sh >| $MULTIDIR/all.html
561
562   # html pour les métriques
563#   cat $COMP_D/entete.html doc_pcmdi.html >| $COMP_D/PCMDI.html
564#   ref=`echo $listsimss | awk ' { print $1 } '`
565#   echo "<h1> Metrics with respect to $ref </h1>" >>  $COMP_D/PCMDI.html
566#   cat $COMP_D/pcmdi_metrics.html >> $COMP_D/PCMDI.html
567#   echo "<h1> Metrics with respect to AR4.00 </h1>" >>  $COMP_D/PCMDI.html
568#   cat $COMP_D/metrics-AR4-as-ref.html >> $COMP_D/PCMDI.html
569#   echo "<h1> Metrics with respect to CMIP5/AMIP multi model </h1>" >>  $COMP_D/PCMDI.html
570#   cat $COMP_D/metrics-CMIP5-AMIP-as-ref.html >> $COMP_D/PCMDI.html
571   for i in 1 2 3 4 5 6 7 ; do
572        ./html_axe$i.sh  $comp >> ${LOGD}/html$$ 2>&1
573   done
574
575   echo Multi atlas, log in ${LOGD}/multiatlas$$ 
576#  for proj in GLOB SH NH ; do
577#  for proj in ATLTN GLOB SH NH ; do
578   for proj in GLOB SH NH ; do
579   for bias in BIAS MAPS ; do
580       for season in YEAR DJF JJA MAM SON ; do
581        ./multi_atlas.sh $comp $bias $season $proj $ref $RUNS  >> ${LOGD}/multiatlas$$ 2>&1
582       done
583   done
584   done
585#exit
586#   cd $COMP_D
587#   ln -sf BIASYEAR.html  $comp.html
588#   ln -sf ACCUEIL.html  $comp.html
589cd $local
590./webpage_multi.sh $comp
591fi
592
593echo UN NOVEL MULTI ATLAS EST DISPONIBLE sur $COMP_D/ACCUEIL.html
Note: See TracBrowser for help on using the repository browser.