#!/bin/sh ######## code REDEC (RE = retrieve; DEC = decode) # @(#)retpak July 95 - ECMWF # 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. # # Procédure à faire exécuter à Reading pour extraire les données du CEPMMT # sur la station ecgate1 A.M. LANQUETTE Juillet 1996 # Guy Cautenet, Anne-Marie Lanquette, Y.P. # (On appelle grista1.out (déjà compilé) dans $HOME) # Changer le nom d'utilisateur {tpp} en votre nom utilisateur #MODELE mis à jour le 16 aout 2000 pour analyse sur niveaux AML #################################################################### # setting the parameters for NQE # --------------------------- #QSUB -s /bin/sh # Change to Bourne Shell #QSUB -q ecgate # Exécution sur ecgate1 ou ecgate2 #%%% QSUB -a "tomorrow 6am" # Exécution demain à 6 h locale #QSUB # #------------------------------- # 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 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 données # # repres=sh harmoniques sphériques (pour données 3D) # repres=gg Grille gaussienne (pour données de surface) # # troncation initiale des données spectrales (réduit le temps calcul) # RESOL=106, # grid >= 1.5 # RESOL=213, # grid >= 0.6 # RESOL=319, # grid >= 0.4 # RESOL=639, # grid < 0.4 # RESOL=AUTO, # adapte la résolution à la maille demandée (défaut actuel) # RESOL=AV, # (or ARCHIVED VALUE) full resolution # # maille de la grille # grid=.75/.75, Résolution en latitude/longitude de la grille (3 chiffres) ################################################################# # 2 - Sélection 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 Dôme (81*43) # area=50./-80./30./-60.,grid=0.5/0.5, # aire NARE (41*41) # # (latitude négative vers le sud, positive vers le nord) # (longitude négative vers l'ouest, positive vers l'est) # # Provisoire : en raison de problème de mémoire sur opgcf4, la zone couverte # et la résolution des données doivent conduire à moins de 85 points en # longitude et 45 points en latitude (ou à moins de 75 points en # longitude et 51 points en latitude en utilisant les programmes modifiés # à cet effet et lancés avec l'option "guy" par les procédures mto_exec.sh5 # cep_traitement.sh5 ou cep_binaire.sh5) (en tout cas, le produit de ces 2 # nombres doit être inférieur ou égal à 3825). De même, le nombre de # paramètres demandés doit être inférieur ou égal à 7 pour chaque demande. ################################################################# # 3 - Type de données (parmi d'autres) # # type=init, données initialisées # type=analysis, données d'analyse non initialisées # type=fg, données "First Guess" # type=FC,step=24,time=12, données basées sur le Forecast de 12 h UTC # après 24 h (par exemple) de simulation # # NOTE : Après le 31/01/96, l'analyse non initialisée a été supprimée. # Il faut utiliser # type=analysis, # pour les dates postérieures au 31/01/96, sachant toutefois qu'il s'agit # bien d'analyse initialisée # # éventuellement conditions spéciales # (Research au lieu d'Operation en standard) # CLASS=RD,EXPVER="ABOI", ################################################################# # 4 - Données 2D ou 3D (parmi d'autres) ####### soit Champs de surface, 2D: ############################ # # param = nom des paramètres (Provisoire : 7 au plus, au choix à chaque # demande) # # SP jpar = 134 Pression de Surface (Pa) # ST jpar = 139 Température de surface (K) # SSW jpar = 140 Surface soil wetness (m d'eau) (SWL1 après 4/8/1993) ** # LSP jpar = 142 Large scale precipitation (m accumulé) # CP jpar = 143 Convective precipitation (m accumulé) # SF jpar = 144 Snowfall (m d'eau accumulé) # TCC jpar = 164 Total cloud cover (0-1) # 10U jpar = 165 U à 10 mètres (m s-1) # 10V jpar = 166 V à 10 mètres (m s-1) # 2T jpar = 167 Température à 2 mètres (K) # 2D jpar = 168 Température de rosée à 2 mètres (K) # SR jpar = 173 Surface Roughness (m) # AL jpar = 174 Albedo # SSR jpar = 176 Surface solar radiation (W/m^2 s accumulé) # 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) # + paramètres qui ne dépendent que de la résolution demandée # Z jpar = 129 Géopotentiel 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=an, # # Pour les données accumulées, => données de Firstguess après 6 h de simulation # param=LSP/CP/SF/SSR # type=FG,step=6, # # ** 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 résultat : Yp02081100sur_don ####### ou Champs atmosphériques, 3D: ######################## # Niveaux pressure level = "pl" # level=all, # # NOTE : après le 9 mars 1999, il y a 21 niveaux au lieu de 15 et les # niveaux supplémentaires 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 paramètres (Provisoire : 7 au plus, au choix à chaque # demande) # # Z jpar = 129 Géopotentiel (m^2/s^2) # T jpar = 130 Température (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 Humidité spécifique (en kg/kg) (jusqu'au 3/04/1995) # W jpar = 135 Vent: w (verticale, en Pa/s) # VO jpar = 138 Vorticité (s-1) # D jpar = 155 Divergence du vent (s-1) # R jpar = 157 Humidité relative (%) # # repres=sh,levtype=pl, # param=t/q/u/v/w/z # -------------------------------------------------------------- # NOTE : après le 04/04/1995 # le paramètre Q n'est plus dans sh (harmoniques sphériques) # mais dans gg (Grille gaussienne) # Il doit donc être demandé séparément par : # # param=q,repres=gg,levtype=pl, # # ou recalculé à partir de R (Humidité relative (%)) par cep_Creatpvor.out # # repres=sh,levtype=pl, # param=t/r/u/v/w/z # ################################################################# # 5 - Date(s) et Heure(s) # # plusieurs possibilités au choix : # date=19961104,time=00, # date=19961104,time=00/to/18/by/6, # date=19961104/to/19920906,time=00, # date=19961104/to/19920907/by/3,time=00, # date=19961104/to/19920906,time=00/to/18/by/6, # date=19961104/to/19920907/by/3,time=00/to/18/by/6, 1############################################################### # I - Extraction des données (code Mars) # ------------------------------------------------------------- # Retrieve the data from MARS in the data file "data_02081100". # Link data file to target file. mars << EOF retrieve,date=20020811/to/20020812,time=00/to/18/by/6,target="data_02081100", RESOL=AV, area=25/-25/5/-5,grid=0.5/0.5, type=an, # 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 EOF if [ $? != 0 ] then echo " The MARS request failed." exit 1 fi ################################################################ # II - Données MARS dans input_datafile # ------------------------------------------------------------- rm input_datafile ln -s data_02081100 input_datafile ################################################################ # III - Exécution du programme de décodage #---------------------------------------------------------------------- $HOME/grista1.out ################################################################ # IV - Sauvegarde des données dans le répertoire $SCRATCH # /scratch/ms/fr/{tpp} #--------------------------------------------------------------- rm data_02081100 compress fort.12 ls -al fort.12.Z # pour des données d'altitude : mv fort.12.Z $SCRATCH/Yp02081100_don.Z # pour des données de surface : # mv fort.12.Z $SCRATCH/Yp02081100sur_don.Z ################################################################# ################################################# # tidy up by deleting unwanted files echo " Fin de la procédure cep_reading.sh5 données sur niveaux" #------------------------------- # This is done by the trap command. exit 0 # # End of example script cep_reading.sh5 du 1 aout 2000 recopie dans pnivmod #################################################################