#!/bin/ksh ######## code REDEC (RE = retrieve; DEC = decode) # This shell-script runs the following commands on a workstation : # Retrieve and decode mars packed data. # All the job output is directed to the standard output. # # Procedure a faire executer a Reading pour extraire les donnees du CEPMMT # sur la station ecgate1 A.M. LANQUETTE Juillet 1996 # Guy Cautenet, Anne-Marie Lanquette, Y.P. # (On appelle grista1.out (deja compile) dans $HOME) # Changer le nom d'utilisateur {tpp} en votre nom utilisateur # #################################################################### # setting the parameters for LoadLeveler # --------------------------- #@ shell = /usr/bin/ksh #@ class = normal #@ queue # #------------------------------- # commands to be executed #------------------------------- set -x PATH=$PATH:. # allows you to run any of your programs or # scripts held in the current directory export PATH # modifier en fonction du nom {auc} utilisateur SCRATCH="/ws/scratch/ms/fr/auc/scratchdir" cd $SCRATCH # # Tidy up at the end of the script or if the scripts end abnormally. trap 'eval cd $CURRENT_DIR;eval /bin/rm -r $SCRATCH/input_datafile ' 0 1 2 3 6 15 # ################################################################# # MARS request ################################################################# # 1 - Grille des donnees # # repres=sh harmoniques spheriques (pour donnees 3D) # repres=gg Grille gaussienne (pour donnees de surface) # # troncation initiale des donnees spectrales (reduit le temps calcul) # RESOL=106, # grid >= 1.5 - avant sep. 1991 # RESOL=213, # grid >= 0.6 - apres sep. 1991 # RESOL=319, # grid >= 0.4 - apres avril 1998 # RESOL=511, # grid < 0.3 - apres nov. 2000 # RESOL=AUTO, # adapte la resolution a la maille demandee (defaut actuel) # RESOL=AV, # (or ARCHIVED VALUE) full resolution # # maille de la grille # grid=.75/.75, Resolution en latitude/longitude de la grille (3 chiffres) ################################################################# # 2 - Selection de la zone # # area=60/-40.5/30/20.25, =latitude Nord/longitude Ouest/latitude # Sud/longitude Est # # area=75/-81/30/30,grid=1.5/1.5, # aire atlantique (75*31) # area=60/-40.5/30/20.25,grid=.75/.75, # FRONT-PYREX etc (82*41) # area=60/-30./30/30.,grid=.75/.75, # ESTIME 2 et 3 (81*41) # area=56./-15./35./25.,grid=0.5/0.5, # aire Puy de Dome (81*43) # area=50./-80./30./-60.,grid=0.5/0.5, # aire NARE (41*41) # area=50./120./20./150.,grid=0.5/0.5, # aire Japon # area=45/25/10/70,grid=.5/.5, # aire Kowet # # (latitude negative vers le sud, positive vers le nord) # (longitude negative vers l'ouest, positive vers l'est) # # Note : Par defaut, la zone couverte et la resolution des donnees # doivent conduire a moins de 180 points en longitude et 80 points en # latitude. De meme, le nombre de parametres demandes doit etre # inferieur ou egal a 7 pour chaque demande. ################################################################# # 3 - Type de donnees (parmi d'autres) # # type=init, donnees initialisees # type=analysis, donnees d'analyse non initialisees # type=FG, donnees "First Guess" (supprimee le 12/09/2000) # type=FC,step=24,time=12, donnees basees sur le Forecast de 12 h UTC # apres 24 h (par exemple) de simulation # # NOTE : Apres le 31/01/96, l'analyse non initialisee a ete supprimee. # Il faut utiliser # type=analysis, # pour les dates posterieures au 31/01/96, sachant toutefois qu'il s'agit # bien d'analyse initialisee # # NOTE ** Depuis le 12/09/2000: # Users of the archive may wish to take note of the removal of the TYPE=FG # 'first-guess fields' from the archive. To compute the analysis increments # one should subtract the most recent forecast (base 00 or 12 UTC) from the # analysis. Users of the TYPE=FG data should retrieve the corresponding # forecast from the previous cycles instead: # # DATE=0,TYPE=FG,TIME=00,STEP=06 -> DATE=-1,TYPE=FC,TIME=12,STEP=12 # DATE=0,TYPE=FG,TIME=06,STEP=06 -> DATE= 0,TYPE=FC,TIME=00,STEP=06 # DATE=0,TYPE=FG,TIME=12,STEP=06 -> DATE= 0,TYPE=FC,TIME=00,STEP=12 # DATE=0,TYPE=FG,TIME=18,STEP=06 -> DATE= 0,TYPE=FC,TIME=12,STEP=06 # # eventuellement conditions speciales # (Research au lieu d'Operation en standard) # CLASS=RD,EXPVER="ABOI", ################################################################# # 4 - Donnees 2D ou 3D (parmi d'autres) ####### soit Champs de surface, 2D: ############################ # # http://wms.ecmwf.int/publications/manuals/libraries/tables/localTable128.html # # param = nom des parametres (Provisoire : 7 au plus, au choix a chaque # demande) # # SP jpar = 134 Pression de Surface (Pa) # ST jpar = 139 Temperature de surface (K) # SSW jpar = 140 Surface soil wetness (m d'eau) (SWL1 apres 4/8/1993) ** # n'existe plus depuis le 27/6/2000 remplace par # 39 Volumetric soil water layer 1 SWVL1 (m^3 m^-3) # LSP jpar = 142 Large scale precipitation (m accumule) # CP jpar = 143 Convective precipitation (m accumule) # SF jpar = 144 Snowfall (m d'eau accumule) issu de la prevision forecast # peut etre remplace par 141 Snow depth SD (m of water equivalent) # TCC jpar = 164 Total cloud cover (0-1) # 10U jpar = 165 U a 10 metres (m s-1) # 10V jpar = 166 V a 10 metres (m s-1) # 2T jpar = 167 Temperature a 2 metres (K) # 2D jpar = 168 Temperature de rosee a 2 metres (K) # SR jpar = 173 Surface Roughness (m) # AL jpar = 174 Albedo # SSR jpar = 176 Surface solar radiation (W/m^2 s accumule) # CCC jpar = 185 Convective cloud cover (0-1) ** # LCC jpar = 186 Low cloud cover (0-1) # MCC jpar = 187 Medium cloud cover (0-1) # HCC jpar = 188 High cloud cover (0-1) # SKT jpar = 235 Skin Temperature (K) # STL2 jpar = 170 Soil temperature level 2 (K) # SWL2 jpar = 171 Soil wetness level 2 (m d'eau) # STL3 jpar = 183 Soil temperature level 3 (K) # SWL3 jpar = 184 Soil wetness level 3 (m d'eau) # STL4 jpar = 236 Soil temperature level 4 (K) # SWL4 jpar = 237 Soil wetness level 4 (m d'eau) # SWVL1 jpar = 39 Volumetric soil water layer 1 (m^3 m^-3) # SWVL2 jpar = 40 Volumetric soil water layer 2 (m^3 m^-3) # SWVL3 jpar = 41 Volumetric soil water layer 3 (m^3 m^-3) # SWVL4 jpar = 42 Volumetric soil water layer 4 (m^3 m^-3) # + parametres qui ne dependent que de la resolution demandee # Z jpar = 129 Geopotentiel du sol (m^2/s^2) # LSM jpar = 172 Land Sea Mask (0 = mer, 1 = terre) # # repres=gg,levtype=sfc,level=off, # param=SP/ST # param=SSW/TCC/AL/SR # param=10U/10V/2T/2D # param=Z/LSM # type=init, # # Pour les donnees accumulees, => donnees de Firstguess apres 6 h de simulation # param=LSP/CP/SF/SSR # type=FG,step=6, (avant le 12/09/2000) # type=FC,step=12, (apres le 12/09/2000) Le Forecast n'existe que pour # 0 h et 12 h, pas 6 h ni 18 h ! # # ** An upgraded version of the IFS model (Cycle 22r3, atmospheric model # identification number 196 in the GRIB headers) was implemented on # Tuesday, 27 June 2000. This version includes important changes to the # parametrization schemes for the land surface, lying snow and sea-ice. # Only the replacement parameters (no. 39, 40, 41, 42) will be available in # MARS from now on. The old products SWL1 to SWL4 (no. 140, 171, 184, 237) # will no longer be archived in MARS. In addition, the parameters CCC, VEG, # and ASQ (no. 185, 199, 233) will no longer be available. # # attention fichier de resultat : Yp01080600sur_don ####### ou Champs atmospheriques, 3D: ######################## # Niveaux pressure level = "pl" # level=all, # # NOTE : apres le 9 mars 1999, il y a 21 niveaux au lieu de 15 et les # niveaux supplementaires sont au-dessus de 30 hPa. Il convient alors # de demander (identique aux niveaux depuis 1989) : # # level=10/30/50/70/100/150/200/250/300/400/500/700/850/925/1000, # # param = nom des parametres (Provisoire : 7 au plus, au choix a chaque # demande) # # Z jpar = 129 Geopotentiel (m^2/s^2) # T jpar = 130 Temperature (en K) # U jpar = 131 Vent: u composante est-ouest en m/s) # V jpar = 132 : v composante sud-nord, en m/s) # Q jpar = 133 Humidite specifique (en kg/kg) (jusqu'au 3/04/1995) # W jpar = 135 Vent: w (verticale, en Pa/s) # VO jpar = 138 Vorticite (s-1) # D jpar = 155 Divergence du vent (s-1) # R jpar = 157 Humidite relative (%) # # repres=sh,levtype=pl, # param=t/q/u/v/w/z # -------------------------------------------------------------- # NOTE : apres le 04/04/1995 # le parametre Q n'est plus dans sh (harmoniques spheriques) # mais dans gg (Grille gaussienne) # Il doit donc etre demande separement par : # # param=q,repres=gg,levtype=pl, # # ou recalcule a partir de R (Humidite relative (%)) par cep_Creatpvor.out # # repres=sh,levtype=pl, # param=t/r/u/v/w/z # ################################################################# # 5 - Date(s) et Heure(s) # # plusieurs possibilites au choix : # date=20010806,time=00, # date=20010806,time=00/to/18/by/6, # date=20010806/to/20001107,time=00, # date=20010806/to/20001107/by/3,time=00, # date=20010806/to/20001107,time=00/to/18/by/6, # date=20010806/to/20001107/by/3,time=00/to/18/by/6, # # N.B. : pour des moyennes mensuelles, preciser : # stream=monthly, # date=19950100,time=00/06/12/18, jour a 00 ! et modifier # manuellement Yp{anmmjjhh{_don pour mettre le jour a 0 ################################################################ # I - Extraction des donnees (code Mars) # ------------------------------------------------------------- # Retrieve the data from MARS in the data file "data_01080600". # Link data file to target file. mars << *eor retrieve,date=20010806/to/20010810,time=00/to/18/by/6,target="data_01080600", RESOL=AV, area=40/-45/-5/45,grid=0.5/0.5, type=analysis, repres=sh,levtype=pl, level=10/30/50/70/100/150/200/250/300/400/500/700/850/925/1000, param=t/r/u/v/w/z end *eor if [ $? != 0 ] then echo " The MARS request failed." exit 1 fi ################################################################ # II - Donnees MARS dans input_datafile # ------------------------------------------------------------- ln -fs data_01080600 input_datafile ################################################################ # III - Execution du programme de decodage #---------------------------------------------------------------------- $HOME/grista1.out ################################################################ # IV - Sauvegarde des donnees dans le repertoire $SCRATCH #--------------------------------------------------------------- rm data_01080600 compress fort.12 ls -al fort.12.Z # pour des donnees d'altitude : mv fort.12.Z Yp01080600_don.Z # pour des donnees de surface : # mv fort.12.Z $SCRATCH/Yp01080600sur_don.Z ################################################################# #------------------------------- # tidy up by deleting unwanted files echo " Fin de la procedure cep_reading.sh " #------------------------------- # This is done by the trap command. exit 0 # # End of example script retpak modifie le 2 juillet 1996 #################################################################