Changeset 2146

Timestamp:

08/08/13 14:18:41 (11 years ago)

Author:

acosce

Message:

• make some correction in prepare_aerosol_climato_annuel
• lev_file.nc is the netcdf file with lev / ap / b variables use in aerosol concentration file in 969539
• correction_axe.py is a script use to invert lev axis to create previous file (lev_file.nc) from cmor output

Location:

TOOLS/CREATE_INTERP_AEROSOLS

Files:

• correction_axe.py (added)
• lev_file.nc (added)
• prepare_aerosol_climato_annuel.job (modified) (6 diffs)

TOOLS/CREATE_INTERP_AEROSOLS/prepare_aerosol_climato_annuel.job

@@ -4 +4 @@
 #les donnees necessaires a lmdz pour les concentrations d'aerosols
 
-##Attention les fichiers concdust_aero_IPSL-CM5A-LR_amip_r1i1p1_197901-200912.nc et climato_1944_LOI_IPCC_HISTOR.nc sont disponibles sous svn 
+##Attention les fichiers lev_file.nc et climato_1944_LOI_IPCC_HISTOR.nc sont disponibles sous svn 
 # au meme emplacement que ce script. 
-# les variables issues de  concdust_aero_IPSL-CM5A-LR_amip_r1i1p1_197901-200912.nc ne sont valables que pour la grille 96x95x39
+# les variables issues de  lev_file.nc  ne sont valables que pour la grille 96x95x39
 
 set -ex
+
+##############################################
+## repertoire de travail sur le scratchdir
+##############################################
 
 DIRWORK=$SCRATCHDIR/tmp
@@ -15 +19 @@
 cd ${DIRWORK}
 \rm -f *
-cp /..../concdust_aero_IPSL-CM5A-LR_amip_r1i1p1_197901-200912.nc exemple.nc
-cp  /.../climato_1944_LOI_IPCC_HISTOR.nc exemple2.nc
-
-for ANNEE in 1999 2000 2001 2002 2003 2004 2005
-do
-
-##############################################
-## repertoire de travail sur le scratchdir
-##############################################
+cp /ccc/scratch/cont003/dsm/p86cozic/CLIMATO/FILES/lev_file.nc exemple.nc
+cp  /ccc/scratch/cont003/dsm/p86cozic/CLIMATO/climato_1944_LOI_IPCC_HISTOR.nc exemple2.nc
 
 
@@ -31 +28 @@
 DIREXP=$STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/ESM1989.4/CHM/Analyse/MM
 
+mkdir -p  $STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/${EXPER}/CHM/Analyse/CLIMATO
 
+for ANNEE in   1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005
 
+do
 
     FICHIER=$DIREXP/${EXPER}_${ANNEE}_allMM_inca_avgr.nc
 
-
 ## extraction des especes
 #############################################
-    ccc_hsm get $FICHIER
     cp $FICHIER tmp.nc
     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 
 
-    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.nc 
-
+    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 
+    
     ncatted -O -a units,AIBCM,o,c,"ug/m3" tmp3.nc tmp3.nc
     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
     ncatted -O -a units,SSSSM,o,c,"ug/m3" tmp3.nc tmp3.nc
-
     ncatted -O -a units,CSSO4M,o,c,"ug/m3" tmp3.nc tmp3.nc
     ncatted -O -a units,ASSO4M,o,c,"ug/m3" tmp3.nc tmp3.nc
+    
     ncrename -O -v LOAD_ASSO4M,load_SO4  tmp3.nc tmp3.nc
     ncrename -O -v ASSO4M,SO4  tmp3.nc tmp3.nc
@@ -75 +73 @@
     ncrename -O -d time_counter,TIME  tmp3.nc tmp3.nc
     ncrename -O -d presnivs,lev  tmp3.nc tmp3.nc
-    ncks -A -v b,ap exemple.nc tmp3.nc
-    ncks -A -v lev,lev_bnds exemple.nc tmp3.nc
+    
+    ncks -A -v lev,b,ap,lev_bnds exemple.nc tmp3.nc
     ncks -A -v p0 exemple2.nc tmp3.nc
-
+    
     mv tmp3.nc $STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/${EXPER}/CHM/Analyse/CLIMATO/climato_${ANNEE}_${EXPER}.nc
-
     \rm  tmp*.nc
 
+done
 
-done