varamma_profile.sh - define VARAMMA environnement
Online usage
$ . ./varamma_profile.sh -d directory -i indir -o outdir -t tempdir
In ${HOME}/.profile, add the following line
. varamma_profile.sh -d directory -i indir -o outdir -t tempdir
define VARAMMA environnement
${VARAMMA}
${VARAMMA_LOG} ++
${VARAMMA_ID}
${VARAMMA_OD}
MANPATH++
For fplod, on aedon.locean-ipsl.upmc.fr:
$ cd /usr/home/fplod/incas/varamma/varamma_ws/src/ $ . ./varamma_profile.sh \ -d $(pwd) \ -i /usr/temp/${LOGNAME}/varamma_d/ \ -o /usr/temp/${LOGNAME}/varamma_d/ \ -t /usr/temp/${LOGNAME}/log/
For fplod on zeus.locean-ipsl.umpc.fr:
$ cd ${HOME}/incas/varamma/varamma_ws/src/ $ . ./varamma_profile.sh \ -d $(pwd) \ -i /usr/temp/${LOGNAME}/varamma_d/ \ -o /usr/temp/${LOGNAME}/varamma_d/ \ -t /usr/temp/${LOGNAME}/log/
/usr/home/fplod/incas/varamma/varamma_ws/src/varamma_profile.sh sur aedon.locean-ipsl.upmc.fr
++ option bavarde
++ machine dependant
++ besoin de posix
++ pas de MANPATH defini par défaut sur zeus
++ climserv
$Id: varamma_profile.sh 2 2008-12-18 16:56:52Z pinsard $
- move all *.sh and *.pro files to src/
dlogd.sh - prompt for destruction of the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ
$ dlogd.sh action
dlogd.sh prompt for destruction of the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ.
/usr/home/fplod/incas/varamma/varamma_ws/src/dlogd.sh sur aedon.locean-ipsl.upmc.fr
$Id: dlogd.sh 2 2008-12-18 16:56:52Z pinsard $
elogd.sh - edit the most recent ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ
$ elogd.sh action
elogd.sh launch ${EDITOR} on the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ.
/usr/home/fplod/incas/varamma/varamma_ws/src/elogd.sh sur aedon.locean-ipsl.upmc.fr
$Id: elogd.sh 2 2008-12-18 16:56:52Z pinsard $
plogd.sh - prints the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ
$ plogd.sh action
plogd.sh prints the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ.
It uses a2ps application.
Print is one on the default printer.
/usr/home/fplod/incas/varamma/varamma_ws/src/plogd.sh sur aedon.locean-ipsl.upmc.fr
$Id: plogd.sh 2 2008-12-18 16:56:52Z pinsard $
tlogd.sh - shows the most recent release of ${VARAMMA_LOG}/<action>.log.YYYY-MM-DDTHH:MM:SSZ
$ tlogd.sh action
tlogd.sh shows the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ.
/usr/home/fplod/incas/varamma/varamma_ws/src/tlogd.sh sur aedon.locean-ipsl.upmc.fr
$Id: tlogd.sh 2 2008-12-18 16:56:52Z pinsard $
cat_clim_ecmwf.sh - build ECMWF climatology file
cat_clim_ecmwf.sh
Build ECMWF climatology file of relative humidity ++using mean daily file by year and reduce geographical domain ++to [60 ° W, 50 ° E] [30 ° S, 45 ° N].
To build ECMWF climatology :
$ cat_clim_ecmwf.sh
/usr/home/fplod/incas/varamma/varamma_ws/src/cat_clim_ecmwf.sh sur aedon.locean-ipsl.upmc.fr
If yearmin and yearmax are the same one can compare input and output files like this:
For yearmin=yearmax=2000
$ cdo diffv ${VARAMMA_ID}/++_ng.nc \ ${VARAMMA_OD}/ecmwf.day.mean_climato_ng.nc
There must be a difference because 2000 is a leap year.
For yearmin=yearmax=2001
$ cdo diffv ${VARAMMA_ID}/ecmwf.day.mean_2001_ng.nc \ ${VARAMMA_OD}/ecmwf.day.mean_climato_ng.nc
There must be no difference on values of ecmwf and info fields ++ à vérifier
++ unset
++ option debug/verbose
++ gestion log
++ vérification accès filein
++ pb levels operationel vs reanalyse
++ 29 février
$Id: cat_clim_ecmwf.sh 2 2008-12-18 16:56:52Z pinsard $
cat_clim_olr.sh - build OLR climatology file
cat_clim_olr.sh
Build OLR climatology file using mean daily file by year and reduce geographical domain to [60 ° W, 50 ° E] [30 ° S, 45 ° N].
To build OLR climatology :
$ cat_clim_olr.sh
/usr/home/fplod/incas/varamma/varamma_ws/src/cat_clim_olr.sh sur aedon.locean-ipsl.upmc.fr
If yearmin and yearmax are the same one can compare input and output files like this:
For yearmin=yearmax=2000
$ cdo diffv ${VARAMMA_ID}/olr.day.mean_2000_ng.nc \ ${VARAMMA_OD}/olr.day.mean_climato_ng.nc
There must be a difference because 2000 is a leap year.
For yearmin=yearmax=2001
$ cdo diffv ${VARAMMA_ID}/olr.day.mean_2001_ng.nc \ ${VARAMMA_OD}/olr.day.mean_climato_ng.nc
There must be no difference on values of olr and info fields ++ à vérifier
++ unset
++ option debug/verbose
++ gestion log
++ vérification accès filein
!! pb zeus et dedale signalé à reseau le 200807
$Id: cat_clim_olr.sh 2 2008-12-18 16:56:52Z pinsard $
is_leap_year.sh - tell if year is leap
is_leap_year.sh [yyyy]
Tell if year is leap by exit status.
If no argument is given, current year will be used.
Is 2000 leap ?
$ is_leap_year.sh 2000 $ echo ${?} 0
Is 2001 leap ?
$ is_leap_year.sh 2001 $ echo ${?} 1
/usr/home/fplod/incas/varamma/varamma_ws/is_leap_year.sh sur aedon.locean-ipsl.upmc.fr
$Id: is_leap_year.sh 2 2008-12-18 16:56:52Z pinsard $
olr_split.sh - split OLR mean daily file
olr_split.sh [-g] [-r] [-f] [-y year]
Split OLR mean daily file by year and reduce geographical domain to [60 ° W, 50 ° E] [30 ° S, 45 ° N].
If -g option is set, there will be no geographical domain reduction.
If -f option is set, previously output files will be overwritten.
If -r option is set, output files will be concatenated in one (might be usefull, combinated with -g option to check).
If -y option is set with a year, only this year will be splitted.
To split ${VARAMMA_ID}/olr.day.mean.nc between [2000,2007] and between [60 ° W, 50 ° E] [30 ° S, 45 ° N]
$ olr_split.sh
To split ${VARAMMA_ID}/olr.day.mean.nc in 2006 between [60 ° W, 50 ° E] [30 ° S, 45 ° N]
$ olr_split.sh -y 2006
To split globaly and build a recomposed file
$ olr_split.sh -g -r
one might compare original and rebuild file
$ cdo diffv ${VARAMMA_ID}/olr.day.mean.nc ${VARAMMA_OD}/olr.day.mean_rebuild.nc
diff should only show differences in global attributes
Does'nt work on zeus because of UDUnits library problems
++++ GROS DOUTE SUR LES LONGITUDES DANS CES FICHIERS +++++++
/usr/home/fplod/incas/varamma/varamma_ws/olr_split.sh sur aedon.locean-ipsl.upmc.fr
++ unset
++ option debug/verbose
++ gestion log
++ vérification accès filein
++ found none interactive actionof ncks if output file exists
$Id: olr_split.sh 2 2008-12-18 16:56:52Z pinsard $
paper01_pre.sh - data prerequisite for paper #1
paper01_pre.sh
Put in ${VARAMMA_ID} all the data files necessary to plot all the figures for paper1.
++ les fichiers de climato sont produits à partir des fichiers de référence donc il faudrait plutôt lancer clim++ avant ou dedans
++ quels sont les fichiers de références AMSU
++ quels sont les fichiers de références ECMWF
++ comment ça marche sur climserv
Does'nt work on zeus because of UDUnits library problems
/usr/home/fplod/incas/varamma/varamma_ws/paper01_pre.sh sur aedon.locean-ipsl.upmc.fr
returns a file name writable in ${VARAMMA_OD}
Categories: | outputs |
---|---|
Uses: | |
Params: |
|
Keywords: |
|
Returns: | filename (dirname and basename) or -1 if error |
Restrictions: |
|
Pre: | |
Post: | |
Todo: | |
Examples: |
To prepare the production of ${VARAMMA_OD}/figure08.png:
IDL> fullfilename_out=forout('figure08.png',/OVERWRITE) IDL> print, fullfilename_out /usr/temp/fplod/varamma_d/figure08.png
If asked an other time without /OVERWRITE, the file is not anymore writable
IDL> fullfilename_out=forout('figure08.png') IDL> print, fullfilename_out -1
History: |
---|
plot hovmuller
Categories: | PNG hoevmuller AMSU OLR |
---|---|
Params: |
|
Keywords: |
|
Restrictions: |
|
Returns: | |
Uses: | |
Pre: |
for AMSU dataset be sure to have cx.anyyyy.nc in the directory defined in ${VARAMMA_ID}/ for OLR dataset be sure to have olr.day.mean_yyyy_ng.nc in the directory defined in ${VARAMMA_ID}/ |
Post: | |
Todo: | ++ parametrisation date début date fin ++ exploitation du champ info des fichiers OLR ++ gros souci avec les longitude OLR !!!! ++ split read and plot ++ saveimage or openps (for image to provide to publisher) ++ filename si small idiot car peu contenir plusieurs figures donc pas canal+year... |
Examples: |
To plot AMSU a4 hovmuller latitude plot for 2006 and produce a PNG file:
IDL> hovmuller, 'AMSU', 'a4', 2006L, AXE='latitude'
To plot AMSU a4 and a5 hovmuller latitude plots for 2006 on the same figure:
IDL> hovmuller, 'AMSU', 'a4', 2006L, AXE='latitude', /nosave,small=[2,1,1] IDL> hovmuller, 'AMSU', 'a5', 2006L, AXE='latitude', /noerase,/nocolorbar,small=[2,1,2]
++ les deux dessins n'occupent pas la même taille.
To plot OLR hovmuller latitude plot for 2006 and produce a PNG file:
IDL> hovmuller, 'OLR', 'n.a.', 2006L, AXE='latitude'
To plot OLR hovmuller longitude plot for 2006 and produce a PNG file:
IDL> hovmuller, 'OLR', 'n.a', 2006L, AXE='longitude'
To plot AMSU a4 hovmuller latitude plot climatology and produce a PNG file:
IDL> hovmuller, 'AMSU', 'a4', 0000L, DATASETTYPE='varamma_t3',AXE='latitude'
History: |
---|
hovmuller latitude canal='a7' mois=['Jan','Feb','Mar','Apr','May','Jun','Jul','Aug','Sep','Oct','Nov','Dec'] mo_lon=[31,28,31,30,31,30,31,31,30,31,30,31] imo=4 &im1=10 ; no du mois file='stagiaire/'+canal+'.an2006.nc' initncdf, file, xaxisname = 'xlon', yaxisname = 'ylat' domdef,-10,5,-30,45 ;domdef,0,5,-30,45 j1=total(mo_lon(0:imo-2)) & J2=total(mo_lon(0:im1-1)) data=read_ncdf('moyenne_tb', j1,j2,/timestep,timevar = 'jours', file = file) time=julday(imo,1,2006)+lindgen(jpt) ; redefinition axe temps pltt,window=0, data, 'yt',title='latitude - time '+canal;,min=250,max=295 saveimage, 'sorties/'+canal+'20060410-hov-10-5.png',/png ;capture d'ecran
Version: | $Id: hovmuller.pro 2 2008-12-18 16:56:52Z pinsard $ |
---|
plot all figures for paper #1
Categories: | |
---|---|
Uses: | |
Restrictions: |
|
Pre: | |
Post: | |
Todo: | |
Examples: |
To produce PNG files :
IDL> paper01
History: |
---|
plot time series
Categories: | PNG Climatology AMSU OLR2 |
---|---|
Params: |
|
Keywords: | |
Restrictions: |
|
Returns: | |
Uses: | |
Pre: |
for AMSU dataset be sure to have cx.an*yyyy*.climato.nc ++ in the directory defined in ${VARAMMA_ID}/ for OLR dataset be sure to have olr.day.mean_climato_ng.nc in the directory defined in ${VARAMMA_ID}/ |
Post: |
PNG file is now present in ${VARAMMA_OD} |
Todo: | ++ exploitation du champ info des fichiers OLR ++ split read and plot ++ saveimage or openps (for image to provide to publisher) ++ add _EXTRA keyword |
Examples: |
To plot AMSU a4 climatology
IDL> timeserie, 'AMSU', 'a4'
To plot OLR climatology
IDL> timeserie, 'OLR', 'n.a.'
History: |
|
---|---|
Version: | $Id: timeserie.pro 2 2008-12-18 16:56:52Z pinsard $ |
plot figure 8 of paper1
Categories: | PNG hoevmuller VARAMMA |
---|---|
Uses: | |
Restrictions: |
|
Pre: |
for AMSU dataset be sure to have cx.anyyyy.nc in the directory defined in ${VARAMMA_ID}/ |
Post: | |
Todo: |
++ different min et max pour chaque plot !
++ write a) b) c) d) near each plot
Examples: |
---|
To produce a PNG file :
IDL> figure08
History: |
---|
- change hovmuller_latitude to hovmuller and the colling sequence
plot figure 9 of paper1
Categories: | PNG hoevmuller VARAMMA |
---|---|
Uses: | |
Restrictions: |
|
Pre: |
for AMSU dataset be sure to have cx.anyyyy.nc in the directory defined in ${VARAMMA_ID}/ |
++ OLR ++ ECMWF
Post: | varamma_profile.sh |
---|---|
Todo: | |
Examples: |
To produce a PNG file :
IDL> figure09
History: |
---|
plot figure 11 of paper1
Uses: | |
---|---|
Restrictions: |
|
Pre: |
++ OLR ++ ECMWF
Post: |
++ jpg is now present in ++ |
---|---|
Todo: | |
Examples: |
To produce a PNG file :
IDL> figure11
History: |
---|
plot figure 12 of paper1
Categories: | AMSU |
---|---|
Uses: | |
Restrictions: |
|
Pre: |
for AMSU dataset be sure to have cx.anyyyyclimato_ng.nc in the directory defined in ${VARAMMA_ID}/ |
Post: | |
Todo: | |
Examples: |
To produce a PNG file :
IDL> figure12
History: |
---|
plot figure 13 of paper1
Categories: | AMSU |
---|---|
Uses: | |
Restrictions: |
|
Pre: |
for AMSU dataset be sure to have cx.anyyyy_climato_ng.nc in the directory defined in ${VARAMMA_ID}/ |
Post: |
++ jpg is now present in ++ |
Todo: | |
Examples: |
To produce a PNG file :
IDL> figure13
History: |
---|
plot figure 14 of paper1
Categories: | |
---|---|
Uses: | |
Restrictions: |
|
Pre: |
for AMSU dataset be sure to have cx.anyyyy_climato.nc in the directory defined in ${VARAMMA_ID}/ |
Post: | |
Todo: | |
Examples: |
To produce a PNG file :
IDL> figure14
History: |
---|