Changeset 3938


Ignore:
Timestamp:
06/22/18 19:10:19 (2 years ago)
Author:
oboucher
Message:

Adding transition period 2015-2023 for strat aerosols

File:
1 edited

Legend:

Unmodified
Added
Removed
  • TOOLS/CMIP6_FORCING/AER_STRAT/volc.sh

    r3767 r3938  
    11dirpwd=`pwd` 
    22 
    3 vv='v4' 
     3vv='v3' 
    44 
    55#--choose the resolution 
     
    88#lmdz='VLR_L79' 
    99#lmdz='VLR' 
    10 #lmdz='LR' 
     10lmdz='LR' 
    1111#lmdz='MR' 
    1212#--added for Frederique Cheruy 
    13 lmdz='zoom_128x89' 
     13#lmdz='zoom_128x89' 
    1414 
    1515dirout='/data/'${USER}'/CMIP6/VOLC/'${lmdz}'_'${vv}'/' 
     
    227227ggg_sun_lmdz_ave=fltarr(dimlatlmdz,dimz,NSW,month_in_year) 
    228228; 
     229tau_sun_lmdz_tr=fltarr(dimlatlmdz,dimz,NSW,month_in_year) 
     230ome_sun_lmdz_tr=fltarr(dimlatlmdz,dimz,NSW,month_in_year) 
     231ggg_sun_lmdz_tr=fltarr(dimlatlmdz,dimz,NSW,month_in_year) 
     232; 
    229233;--LW properties tau_abs 
    230234tau_ear_lmdz=fltarr(dimlatlmdz,dimz,NLW,month_in_year) 
    231235; 
    232236tau_ear_lmdz_ave=fltarr(dimlatlmdz,dimz,NLW,month_in_year) 
     237; 
     238tau_ear_lmdz_tr=fltarr(dimlatlmdz,dimz,NLW,month_in_year) 
    233239; 
    234240tau_sun_lmdz_ave(*,*,*,0)=1.e-15 
     
    238244; 
    239245for year=0, dimtime/month_in_year-1 do begin 
    240 ;for year=141, 141 do begin ;--Pinatubo 
    241246; 
    242247tau_550_lmdz(*,*,*)=1.e-15 
     
    402407; 
    403408for mth=1, month_in_year-1 do begin 
    404 ggg_sun_lmdz_ave(*,*,*,mth)=ggg_sun_lmdz_ave(*,*,*,0) 
    405 ome_sun_lmdz_ave(*,*,*,mth)=ome_sun_lmdz_ave(*,*,*,0) 
    406 tau_sun_lmdz_ave(*,*,*,mth)=tau_sun_lmdz_ave(*,*,*,0) 
    407 tau_ear_lmdz_ave(*,*,*,mth)=tau_ear_lmdz_ave(*,*,*,0) 
     409  ggg_sun_lmdz_ave(*,*,*,mth)=ggg_sun_lmdz_ave(*,*,*,0) 
     410  ome_sun_lmdz_ave(*,*,*,mth)=ome_sun_lmdz_ave(*,*,*,0) 
     411  tau_sun_lmdz_ave(*,*,*,mth)=tau_sun_lmdz_ave(*,*,*,0) 
     412  tau_ear_lmdz_ave(*,*,*,mth)=tau_ear_lmdz_ave(*,*,*,0) 
    408413endfor 
    409414; 
     
    457462netcdfwrite,output+'taulwstrat.2D.ave.nc',opticstruct,clobber=1, $ 
    458463            attributes=attributes, dimensions=dimensions 
     464; 
     465;--finally prepare the transition period 2015-2023 
     466;--we mix the last year of the period with the climatology 
     467; 
     468for year=2015,2023 do begin 
     469; 
     470chyr=strcompress(year,/rem) 
     471;--compute weights 
     472w1=float(2024-year)/10. 
     473w2=1-w1 
     474;--SW properties 
     475tau_sun_lmdz_tr=w1*tau_sun_lmdz+w2*tau_sun_lmdz_ave 
     476ome_sun_lmdz_tr=w1*tau_sun_lmdz*ome_sun_lmdz+w2*tau_sun_lmdz_ave*ome_sun_lmdz_ave 
     477ggg_sun_lmdz_tr=w1*tau_sun_lmdz*ome_sun_lmdz*ggg_sun_lmdz+ $  
     478                w2*tau_sun_lmdz_ave*ome_sun_lmdz_ave*ggg_sun_lmdz_ave 
     479ggg_sun_lmdz_tr=ggg_sun_lmdz_tr/ome_sun_lmdz_tr 
     480ome_sun_lmdz_tr=ome_sun_lmdz_tr/tau_sun_lmdz_tr 
     481; 
     482;--LW properties 
     483tau_ear_tr=w1*tau_ear_lmdz+w2*tau_ear_lmdz_ave 
     484; 
     485;--prepare SW output 
     486opticstruct={lat:fltarr(dimlatlmdz),lev:fltarr(dimz),             $  
     487             wav:fltarr(NSW),time:fltarr(month_in_year),          $  
     488             tau_sun:fltarr(dimlatlmdz,dimz,NSW,month_in_year),   $ 
     489             ome_sun:fltarr(dimlatlmdz,dimz,NSW,month_in_year),   $ 
     490             ggg_sun:fltarr(dimlatlmdz,dimz,NSW,month_in_year) } 
     491; 
     492opticstruct.lat=latitudelmdz 
     493opticstruct.lev=lev 
     494opticstruct.wav=(wl1_sun(0:NSW-1)+wl2_sun(0:NSW-1))/2. 
     495opticstruct.time=float(indgen(month_in_year)+1) 
     496opticstruct.tau_sun=tau_sun_lmdz_tr 
     497opticstruct.ome_sun=ome_sun_lmdz_tr 
     498opticstruct.ggg_sun=ggg_sun_lmdz_tr 
     499; 
     500attributes = {units:strarr(7),long_name:strarr(7)} 
     501attributes.units =     ['degrees_north','level','meters','month','-','-','-'] 
     502attributes.long_name = ['latitude','level','wavelength','time','tau_sun','ome_sun','g_sun'] 
     503; 
     504dimensions = {isdim:intarr(7), links:intarr(4,7)} 
     505       dimensions.isdim =  [1,1,1,1,0,0,0]  ; (1=dimension, 0=variable) 
     506       dimensions.links = [ [-1,-1,-1,-1],[-1,-1,-1,-1],   $ 
     507                            [-1,-1,-1,-1],[-1,-1,-1,-1],   $ 
     508                            [0,1,2,3],[0,1,2,3],[0,1,2,3]  ] 
     509; 
     510netcdfwrite,output+'tauswstrat.2D.'+chyr+'.nc',opticstruct,clobber=1, $ 
     511            attributes=attributes, dimensions=dimensions 
     512; 
     513;--prepare LW output 
     514opticstruct={lat:fltarr(dimlatlmdz),lev:fltarr(dimz),             $  
     515             wav:fltarr(NLW),time:fltarr(month_in_year),          $ 
     516             tau_ear:fltarr(dimlatlmdz,dimz,NLW,month_in_year) } 
     517; 
     518opticstruct.lat=latitudelmdz 
     519opticstruct.lev=lev 
     520opticstruct.wav=(wl1_earth+wl2_earth)/2. 
     521opticstruct.time=float(indgen(month_in_year)+1) 
     522opticstruct.tau_ear=tau_ear_lmdz_tr 
     523; 
     524attributes = {units:strarr(5),long_name:strarr(5)} 
     525attributes.units = ['degrees_north','level','cm-1','month','-'] 
     526attributes.long_name = ['latitude','level','wavenumber','time','tau_ear'] 
     527; 
     528dimensions = {isdim:intarr(5), links:intarr(4,5)} 
     529       dimensions.isdim =  [1,1,1,1,0]  ; (1=dimension, 0=variable) 
     530       dimensions.links = [ [-1,-1,-1,-1],[-1,-1,-1,-1],   $ 
     531                            [-1,-1,-1,-1],[-1,-1,-1,-1],   $ 
     532                            [0,1,2,3]  ] 
     533; 
     534netcdfwrite,output+'taulwstrat.2D.'+chyr+'.nc',opticstruct,clobber=1, $ 
     535            attributes=attributes, dimensions=dimensions 
     536; 
     537endfor 
     538;--end loop on years 
    459539; 
    460540end 
     
    479559EOF2 
    480560 
    481 #qsub volc.job 
     561qsub volc.job 
Note: See TracChangeset for help on using the changeset viewer.