Changeset 2146


Ignore:
Timestamp:
08/08/13 14:18:41 (8 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:
2 added
1 edited

Legend:

Unmodified
Added
Removed
  • TOOLS/CREATE_INTERP_AEROSOLS/prepare_aerosol_climato_annuel.job

    r2141 r2146  
    44#les donnees necessaires a lmdz pour les concentrations d'aerosols 
    55 
    6 ##Attention les fichiers concdust_aero_IPSL-CM5A-LR_amip_r1i1p1_197901-200912.nc et climato_1944_LOI_IPCC_HISTOR.nc sont disponibles sous svn  
     6##Attention les fichiers lev_file.nc et climato_1944_LOI_IPCC_HISTOR.nc sont disponibles sous svn  
    77# au meme emplacement que ce script.  
    8 # les variables issues de  concdust_aero_IPSL-CM5A-LR_amip_r1i1p1_197901-200912.nc ne sont valables que pour la grille 96x95x39 
     8# les variables issues de  lev_file.nc ne sont valables que pour la grille 96x95x39 
    99 
    1010set -ex 
     11 
     12############################################## 
     13## repertoire de travail sur le scratchdir 
     14############################################## 
    1115 
    1216DIRWORK=$SCRATCHDIR/tmp 
     
    1519cd ${DIRWORK} 
    1620\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 ############################################## 
     21cp /ccc/scratch/cont003/dsm/p86cozic/CLIMATO/FILES/lev_file.nc exemple.nc 
     22cp  /ccc/scratch/cont003/dsm/p86cozic/CLIMATO/climato_1944_LOI_IPCC_HISTOR.nc exemple2.nc 
    2623 
    2724 
     
    3128DIREXP=$STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/ESM1989.4/CHM/Analyse/MM 
    3229 
     30mkdir -p  $STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/${EXPER}/CHM/Analyse/CLIMATO 
    3331 
     32for ANNEE in   1985 1986 1987 1988 1989 1990 1991 1992 1993 1994 1995 1996 1997 1998 1999 2000 2001 2002 2003 2004 2005 
    3433 
     34do 
    3535 
    3636    FICHIER=$DIREXP/${EXPER}_${ANNEE}_allMM_inca_avgr.nc 
    3737 
    38  
    3938## extraction des especes 
    4039############################################# 
    41     ccc_hsm get $FICHIER 
    4240    cp $FICHIER tmp.nc 
    4341    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 
     
    4947    ncpdq -O -a -lat tmp1bis.nc tmp2.nc  
    5048 
    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.nc  
    52  
     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     
    5351    ncatted -O -a units,AIBCM,o,c,"ug/m3" tmp3.nc tmp3.nc 
    5452    ncatted -O -a units,AIPOMM,o,c,"ug/m3" tmp3.nc tmp3.nc 
     
    5957    ncatted -O -a units,ASSSM,o,c,"ug/m3" tmp3.nc tmp3.nc 
    6058    ncatted -O -a units,SSSSM,o,c,"ug/m3" tmp3.nc tmp3.nc 
    61  
    6259    ncatted -O -a units,CSSO4M,o,c,"ug/m3" tmp3.nc tmp3.nc 
    6360    ncatted -O -a units,ASSO4M,o,c,"ug/m3" tmp3.nc tmp3.nc 
     61     
    6462    ncrename -O -v LOAD_ASSO4M,load_SO4  tmp3.nc tmp3.nc 
    6563    ncrename -O -v ASSO4M,SO4  tmp3.nc tmp3.nc 
     
    7573    ncrename -O -d time_counter,TIME  tmp3.nc tmp3.nc 
    7674    ncrename -O -d presnivs,lev  tmp3.nc tmp3.nc 
    77     ncks -A -v b,ap exemple.nc tmp3.nc 
    78     ncks -A -v lev,lev_bnds exemple.nc tmp3.nc 
     75     
     76    ncks -A -v lev,b,ap,lev_bnds exemple.nc tmp3.nc 
    7977    ncks -A -v p0 exemple2.nc tmp3.nc 
    80  
     78     
    8179    mv tmp3.nc $STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/historical/${EXPER}/CHM/Analyse/CLIMATO/climato_${ANNEE}_${EXPER}.nc 
    82  
    8380    \rm  tmp*.nc 
    8481 
     82done 
    8583 
    86 done 
Note: See TracChangeset for help on using the changeset viewer.