- Timestamp:
- 08/08/13 14:18:41 (11 years ago)
- Location:
- TOOLS/CREATE_INTERP_AEROSOLS
- Files:
-
- 2 added
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TOOLS/CREATE_INTERP_AEROSOLS/prepare_aerosol_climato_annuel.job
r2141 r2146 4 4 #les donnees necessaires a lmdz pour les concentrations d'aerosols 5 5 6 ##Attention les fichiers concdust_aero_IPSL-CM5A-LR_amip_r1i1p1_197901-200912.nc et climato_1944_LOI_IPCC_HISTOR.nc sont disponibles sous svn6 ##Attention les fichiers lev_file.nc et climato_1944_LOI_IPCC_HISTOR.nc sont disponibles sous svn 7 7 # au meme emplacement que ce script. 8 # les variables issues de concdust_aero_IPSL-CM5A-LR_amip_r1i1p1_197901-200912.ncne sont valables que pour la grille 96x95x398 # les variables issues de lev_file.nc ne sont valables que pour la grille 96x95x39 9 9 10 10 set -ex 11 12 ############################################## 13 ## repertoire de travail sur le scratchdir 14 ############################################## 11 15 12 16 DIRWORK=$SCRATCHDIR/tmp … … 15 19 cd ${DIRWORK} 16 20 \rm -f * 17 cp /..../concdust_aero_IPSL-CM5A-LR_amip_r1i1p1_197901-200912.nc exemple.nc 18 cp /.../climato_1944_LOI_IPCC_HISTOR.nc exemple2.nc 19 20 for ANNEE in 1999 2000 2001 2002 2003 2004 2005 21 do 22 23 ############################################## 24 ## repertoire de travail sur le scratchdir 25 ############################################## 21 cp /ccc/scratch/cont003/dsm/p86cozic/CLIMATO/FILES/lev_file.nc exemple.nc 22 cp /ccc/scratch/cont003/dsm/p86cozic/CLIMATO/climato_1944_LOI_IPCC_HISTOR.nc exemple2.nc 26 23 27 24 … … 31 28 DIREXP=$STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/ESM1989.4/CHM/Analyse/MM 32 29 30 mkdir -p $STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/${EXPER}/CHM/Analyse/CLIMATO 33 31 32 for ANNEE in 1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 34 33 34 do 35 35 36 36 FICHIER=$DIREXP/${EXPER}_${ANNEE}_allMM_inca_avgr.nc 37 37 38 39 38 ## extraction des especes 40 39 ############################################# 41 ccc_hsm get $FICHIER42 40 cp $FICHIER tmp.nc 43 41 ncks -v lat,lon,ps,pmid,temp,AIPOMM,ASBCM,AIBCM,ASPOMM,CSSSM,ASSSM,CIDUSTM,ASSO4M,LOAD_ASSO4M,CSSO4M,LOAD_CSSO4M,SSSSM,LOAD_SSSSM,LOAD_AIPOMM,LOAD_ASBCM,LOAD_AIBCM,LOAD_ASPOMM,LOAD_CSSSM,LOAD_ASSSM,LOAD_CIDUSTM tmp.nc tmp1.nc … … 49 47 ncpdq -O -a -lat tmp1bis.nc tmp2.nc 50 48 51 ncap -s "AIPOMM=float(AIPOMM *0.7913);ASBCM=float(ASBCM);AIBCM=float(AIBCM);ASPOMM=float(ASPOMM*0.7913);CSSSM=float(CSSSM);ASSSM=float(ASSSM);CIDUSTM=float(CIDUSTM);ASSO4M=float(ASSO4M);SSSSM=float(SSSSM);CSSO4M=float(CSSO4M);LOAD_AIPOMM=LOAD_AIPOMM*0.7913;LOAD_ASPOMM=LOAD_ASPOMM*0.7913;lat=double(lat);lon=double(lon)" tmp2.nc tmp3.nc52 49 ncap -s "AIPOMM=float(AIPOMM);ASBCM=float(ASBCM);AIBCM=float(AIBCM);ASPOMM=float(ASPOMM);CSSSM=float(CSSSM);ASSSM=float(ASSSM);CIDUSTM=float(CIDUSTM);ASSO4M=float(ASSO4M);SSSSM=float(SSSSM);CSSO4M=float(CSSO4M);lat=double(lat);lon=double(lon)" tmp2.nc tmp3.nc 50 53 51 ncatted -O -a units,AIBCM,o,c,"ug/m3" tmp3.nc tmp3.nc 54 52 ncatted -O -a units,AIPOMM,o,c,"ug/m3" tmp3.nc tmp3.nc … … 59 57 ncatted -O -a units,ASSSM,o,c,"ug/m3" tmp3.nc tmp3.nc 60 58 ncatted -O -a units,SSSSM,o,c,"ug/m3" tmp3.nc tmp3.nc 61 62 59 ncatted -O -a units,CSSO4M,o,c,"ug/m3" tmp3.nc tmp3.nc 63 60 ncatted -O -a units,ASSO4M,o,c,"ug/m3" tmp3.nc tmp3.nc 61 64 62 ncrename -O -v LOAD_ASSO4M,load_SO4 tmp3.nc tmp3.nc 65 63 ncrename -O -v ASSO4M,SO4 tmp3.nc tmp3.nc … … 75 73 ncrename -O -d time_counter,TIME tmp3.nc tmp3.nc 76 74 ncrename -O -d presnivs,lev tmp3.nc tmp3.nc 77 ncks -A -v b,ap exemple.nc tmp3.nc78 ncks -A -v lev, lev_bnds exemple.nc tmp3.nc75 76 ncks -A -v lev,b,ap,lev_bnds exemple.nc tmp3.nc 79 77 ncks -A -v p0 exemple2.nc tmp3.nc 80 78 81 79 mv tmp3.nc $STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/${EXPER}/CHM/Analyse/CLIMATO/climato_${ANNEE}_${EXPER}.nc 82 83 80 \rm tmp*.nc 84 81 82 done 85 83 86 done
Note: See TracChangeset
for help on using the changeset viewer.