VARAMMA manuals

Shell scripts

varamma_profile.sh

dlogd.sh

elogd.sh

plogd.sh

tlogd.sh

cat_clim_ecmwf.sh

cat_clim_olr.sh

is_leap_year.sh

olr_split.sh

paper01_pre.sh

IDL scripts

forout.pro

hovmuller.pro

paper01.pro

timeserie.pro

figure08.pro

figure09.pro

figure11.pro

figure12.pro

figure13.pro

figure14.pro

varamma_profile.sh

NAME

varamma_profile.sh - define VARAMMA environnement

SYNOPSYS

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

DESCRIPTION

define VARAMMA environnement

${VARAMMA}

${VARAMMA_LOG} ++

${VARAMMA_ID}

${VARAMMA_OD}

MANPATH++

EXAMPLES

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/

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/src/varamma_profile.sh sur aedon.locean-ipsl.upmc.fr

EVOLUTIONS

++ 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 $

  • fplod 2008-08-14T12:28:16Z aedon.locean-ipsl.upmc.fr (Darwin)
  • move all *.sh and *.pro files to src/
  • fplod 2008-07-21T10:58:47Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation

dlogd.sh

NAME

dlogd.sh - prompt for destruction of the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ

SYNOPSIS

$ dlogd.sh action

DESCRIPTION

dlogd.sh prompt for destruction of the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ.

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/src/dlogd.sh sur aedon.locean-ipsl.upmc.fr

EVOLUTIONS

$Id: dlogd.sh 2 2008-12-18 16:56:52Z pinsard $

  • fplod 2008-11-26T11:26:40Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation

elogd.sh

NAME

elogd.sh - edit the most recent ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ

SYNOPSIS

$ elogd.sh action

DESCRIPTION

elogd.sh launch ${EDITOR} on the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ.

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/src/elogd.sh sur aedon.locean-ipsl.upmc.fr

EVOLUTIONS

$Id: elogd.sh 2 2008-12-18 16:56:52Z pinsard $

  • fplod 2008-11-26T11:26:40Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation

plogd.sh

NAME

plogd.sh - prints the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ

SYNOPSIS

$ plogd.sh action

DESCRIPTION

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.

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/src/plogd.sh sur aedon.locean-ipsl.upmc.fr

EVOLUTIONS

$Id: plogd.sh 2 2008-12-18 16:56:52Z pinsard $

  • fplod 2008-11-26T11:26:40Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation

tlogd.sh

NAME

tlogd.sh - shows the most recent release of ${VARAMMA_LOG}/<action>.log.YYYY-MM-DDTHH:MM:SSZ

SYNOPSIS

$ tlogd.sh action

DESCRIPTION

tlogd.sh shows the most recent release of ${VARAMMA_LOG}/action.log.YYYY-MM-DDTHH:MM:SSZ.

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/src/tlogd.sh sur aedon.locean-ipsl.upmc.fr

EVOLUTIONS

$Id: tlogd.sh 2 2008-12-18 16:56:52Z pinsard $

  • fplod 2008-11-26T11:26:40Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation

cat_clim_ecmwf.sh

NAME

cat_clim_ecmwf.sh - build ECMWF climatology file

SYNOPSYS

cat_clim_ecmwf.sh

DESCRIPTION

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].

EXAMPLES

To build ECMWF climatology :

$ cat_clim_ecmwf.sh

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/src/cat_clim_ecmwf.sh sur aedon.locean-ipsl.upmc.fr

COMMENTS

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

SEE ALSO

EVOLUTIONS

++ 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 $

  • fplod 2008-11-20T09:36:58Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation from /homedata/eymard/scripts/script_cat_clim_ECMWF3d

cat_clim_olr.sh

NAME

cat_clim_olr.sh - build OLR climatology file

SYNOPSYS

cat_clim_olr.sh

DESCRIPTION

Build OLR climatology file using mean daily file by year and reduce geographical domain to [60 ° W, 50 ° E] [30 ° S, 45 ° N].

EXAMPLES

To build OLR climatology :

$ cat_clim_olr.sh

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/src/cat_clim_olr.sh sur aedon.locean-ipsl.upmc.fr

COMMENTS

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

SEE ALSO

EVOLUTIONS

++ 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 $

  • fplod 2008-08-05T09:37:14Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation from /homedata/eymard/scripts/script_cat_clim_AMSU

is_leap_year.sh

NAME

is_leap_year.sh - tell if year is leap

SYNOPSYS

is_leap_year.sh [yyyy]

DESCRIPTION

Tell if year is leap by exit status.

If no argument is given, current year will be used.

EXAMPLES

Is 2000 leap ?

$ is_leap_year.sh 2000
$ echo ${?}
0

Is 2001 leap ?

$ is_leap_year.sh 2001
$ echo ${?}
1

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/is_leap_year.sh sur aedon.locean-ipsl.upmc.fr

EVOLUTIONS

$Id: is_leap_year.sh 2 2008-12-18 16:56:52Z pinsard $

olr_split.sh

NAME

olr_split.sh - split OLR mean daily file

SYNOPSYS

olr_split.sh [-g] [-r] [-f] [-y year]

DESCRIPTION

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.

EXAMPLES

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

CAUTIONS

Does'nt work on zeus because of UDUnits library problems

++++ GROS DOUTE SUR LES LONGITUDES DANS CES FICHIERS +++++++

FILES

original location

/usr/home/fplod/incas/varamma/varamma_ws/olr_split.sh sur aedon.locean-ipsl.upmc.fr

TODO

++ unset

++ option debug/verbose

++ gestion log

++ vérification accès filein

++ found none interactive actionof ncks if output file exists

EVOLUTIONS

$Id: olr_split.sh 2 2008-12-18 16:56:52Z pinsard $

  • fplod 2008-08-14T14:53:04Z aedon.locean-ipsl.upmc.fr (Darwin)
  • fplod 2008-08-13T14:26:41Z aedon.locean-ipsl.upmc.fr (Darwin)
    • add -f option
    • add -y option
    • add 180 degre to longitude centered on meridien 0 (<0 =W , >0=E)
  • fplod 2008-07-24T15:02:12Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation

paper01_pre.sh

NAME

paper01_pre.sh - data prerequisite for paper #1

SYNOPSIS

paper01_pre.sh

DESCRIPTION

Put in ${VARAMMA_ID} all the data files necessary to plot all the figures for paper1.

TODO

++ 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

CAUTIONS

Does'nt work on zeus because of UDUnits library problems

FILES

/usr/home/fplod/incas/varamma/varamma_ws/paper01_pre.sh sur aedon.locean-ipsl.upmc.fr

EVOLUTIONS

  • fplod 2008-08-13T12:27:41Z aedon.locean-ipsl.upmc.fr (Darwin)
    • add OLR reference file and yearly-afratl files production
  • fplod 2008-08-12T13:55:32Z aedon.locean-ipsl.upmc.fr (Darwin)
    • creation. not finish

forout.pro

returns a file name writable in ${VARAMMA_OD}

Categories:

outputs

Uses:
Params:
name : in, required, type=string

basename of the file to create in ${VARAMMA_OD}/

Keywords:
OVERWRITE : default=false

to overwrite the output image file if exists.

Returns:

filename (dirname and basename) or -1 if error

Restrictions:
  • Requires SAXO
Pre:

varamma_profile.sh

Post:

varamma_profile.sh

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:

hovmuller.pro

plot hovmuller

Categories:

PNG hoevmuller AMSU OLR

Params:
dataset : in, required, type='string'

dataset to plot

canal : in, optional, type='string'

canal of AMSU

only used if AMSU dataset

iyear : in, required, type='integer'

year

Keywords:
DATASETTYPE : default='for AMSU ++, for OLR ++'

dataset type name. For each dataset, there are several file type (year vs total, global vs regional). ++

AXE : default='latitude'

geographical axe to be use for hovmuller can be either 'latitude' or 'longitude'

OVERWRITE : default=false

to overwrite the output image file if exists.

Useless if NOSAVE is set.

NOSAVE : default=false

to NOT save the output image in a file

might be usefull if several plots are needed on the same figure.

file if exist might be overwritten if OVERWRITE is set to 1.

_EXTRA :

Used to declare that this routine accepts inherited keywords

see example with SMALL keyword

Restrictions:
  • Requires SAXO
Returns:
Uses:

pltt_

cm_4cal_

forout.pro

Pre:

varamma_profile.sh

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:

varamma_profile.sh

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 $

paper01.pro

plot all figures for paper #1

Categories:
Uses:

figure08.pro figure09.pro figure14.pro

Restrictions:
  • Requires SAXO
Pre:

varamma_profile.sh

paper01_pre.sh

Post:

varamma_profile.sh

Todo:
Examples:

To produce PNG files :

IDL> paper01
History:

timeserie.pro

plot time series

Categories:

PNG Climatology AMSU OLR2

Params:
dataset : in, required, type='string: AMSU OLR'

dataset to plot

canal : in, optional, type='string: ax bx'

canal of AMSU

only used if AMSU dataset

Keywords:
Restrictions:
  • Requires SAXO
Returns:
Uses:

pltt_

cm_4cal_

Pre:

varamma_profile.sh

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:

varamma_profile.sh

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:
  • fplod 2008-08-05T14:03:03Z zeus.locean-ipsl.upmc.fr (Linux)

    • creation from hovmuller_latitude.pro

    • extract from ananewvaramma3.pro 20080718 :

      ; series temporelles
      

      System Message: ERROR/3 (varamma_manuals.rst, line 1253)

      Inconsistent literal block quoting.

      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] canal='a5' file='stagiaire/'+canal+'.climato.nc' initncdf, file, xaxisname = 'xlon', yaxisname = 'ylat'

      domdef,-10,5,15,25 imo=4 ; no du mois im1=9 j1=total(mo_lon(0:imo-2)) & J2=total(mo_lon(0:im1-2))+mo_lon(im1-1)-1 ;imo=1 ; no du mois ;im1=12 ;j1=0 & j2=364 ;J2=total(mo_lon(0:im1-2))+mo_lon(im1-1)-1 print,j1,j2 data=read_ncdf('moyenne_tb', j1,j2,/timestep,timevar = 'jours', file = file) time=julday(imo,1,2000)+lindgen(jpt) seriea5=(data.arr-min(data.arr))/(max(data.arr)-min(data.arr)) pltt,window=1,seriea5, 't',title='climatology time series',min=0.1,max=0.9 xyouts,max(time)+20,0.1,'a5',color=0,charsize=2

Version:

$Id: timeserie.pro 2 2008-12-18 16:56:52Z pinsard $

figure08.pro

plot figure 8 of paper1

Categories:

PNG hoevmuller VARAMMA

Uses:

hovmuller.pro forout.pro

Restrictions:
  • Requires SAXO
Pre:

varamma_profile.sh

for AMSU dataset be sure to have cx.anyyyy.nc in the directory defined in ${VARAMMA_ID}/

Post:

varamma_profile.sh

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

figure09.pro

plot figure 9 of paper1

Categories:

PNG hoevmuller VARAMMA

Uses:

hovmuller.pro

Restrictions:
  • Requires SAXO
Pre:

varamma_profile.sh

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:

figure11.pro

plot figure 11 of paper1

Uses:

hovmuller.pro

Restrictions:
  • Requires SAXO
Pre:

varamma_profile.sh

++ OLR ++ ECMWF

Post:

varamma_profile.sh

++ jpg is now present in ++

Todo:
Examples:

To produce a PNG file :

IDL> figure11
History:

figure12.pro

plot figure 12 of paper1

Categories:

AMSU

Uses:
Restrictions:
  • Requires SAXO
Pre:

varamma_profile.sh

for AMSU dataset be sure to have cx.anyyyyclimato_ng.nc in the directory defined in ${VARAMMA_ID}/

Post:

varamma_profile.sh

Todo:
Examples:

To produce a PNG file :

IDL> figure12
History:

figure13.pro

plot figure 13 of paper1

Categories:

AMSU

Uses:
Restrictions:
  • Requires SAXO
Pre:

varamma_profile.sh

for AMSU dataset be sure to have cx.anyyyy_climato_ng.nc in the directory defined in ${VARAMMA_ID}/

Post:

varamma_profile.sh

++ jpg is now present in ++

Todo:
Examples:

To produce a PNG file :

IDL> figure13
History:

figure14.pro

plot figure 14 of paper1

Categories:
Uses:

hovmuller.pro forout.pro

Restrictions:
  • Requires SAXO
Pre:

varamma_profile.sh

for AMSU dataset be sure to have cx.anyyyy_climato.nc in the directory defined in ${VARAMMA_ID}/

Post:

varamma_profile.sh

Todo:
Examples:

To produce a PNG file :

IDL> figure14
History:

Docutils System Messages

System Message: ERROR/3 (varamma_manuals.rst, line 988); backlink

Unknown target name: "pltt".

System Message: ERROR/3 (varamma_manuals.rst, line 991); backlink

Unknown target name: "cm_4cal".

System Message: ERROR/3 (varamma_manuals.rst, line 1200); backlink

Unknown target name: "pltt".

System Message: ERROR/3 (varamma_manuals.rst, line 1203); backlink

Unknown target name: "cm_4cal".