Changeset 3938
- Timestamp:
- 06/22/18 19:10:19 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TOOLS/CMIP6_FORCING/AER_STRAT/volc.sh
r3767 r3938 1 1 dirpwd=`pwd` 2 2 3 vv='v 4'3 vv='v3' 4 4 5 5 #--choose the resolution … … 8 8 #lmdz='VLR_L79' 9 9 #lmdz='VLR' 10 #lmdz='LR'10 lmdz='LR' 11 11 #lmdz='MR' 12 12 #--added for Frederique Cheruy 13 lmdz='zoom_128x89'13 #lmdz='zoom_128x89' 14 14 15 15 dirout='/data/'${USER}'/CMIP6/VOLC/'${lmdz}'_'${vv}'/' … … 227 227 ggg_sun_lmdz_ave=fltarr(dimlatlmdz,dimz,NSW,month_in_year) 228 228 ; 229 tau_sun_lmdz_tr=fltarr(dimlatlmdz,dimz,NSW,month_in_year) 230 ome_sun_lmdz_tr=fltarr(dimlatlmdz,dimz,NSW,month_in_year) 231 ggg_sun_lmdz_tr=fltarr(dimlatlmdz,dimz,NSW,month_in_year) 232 ; 229 233 ;--LW properties tau_abs 230 234 tau_ear_lmdz=fltarr(dimlatlmdz,dimz,NLW,month_in_year) 231 235 ; 232 236 tau_ear_lmdz_ave=fltarr(dimlatlmdz,dimz,NLW,month_in_year) 237 ; 238 tau_ear_lmdz_tr=fltarr(dimlatlmdz,dimz,NLW,month_in_year) 233 239 ; 234 240 tau_sun_lmdz_ave(*,*,*,0)=1.e-15 … … 238 244 ; 239 245 for year=0, dimtime/month_in_year-1 do begin 240 ;for year=141, 141 do begin ;--Pinatubo241 246 ; 242 247 tau_550_lmdz(*,*,*)=1.e-15 … … 402 407 ; 403 408 for 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) 408 413 endfor 409 414 ; … … 457 462 netcdfwrite,output+'taulwstrat.2D.ave.nc',opticstruct,clobber=1, $ 458 463 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 ; 468 for year=2015,2023 do begin 469 ; 470 chyr=strcompress(year,/rem) 471 ;--compute weights 472 w1=float(2024-year)/10. 473 w2=1-w1 474 ;--SW properties 475 tau_sun_lmdz_tr=w1*tau_sun_lmdz+w2*tau_sun_lmdz_ave 476 ome_sun_lmdz_tr=w1*tau_sun_lmdz*ome_sun_lmdz+w2*tau_sun_lmdz_ave*ome_sun_lmdz_ave 477 ggg_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 479 ggg_sun_lmdz_tr=ggg_sun_lmdz_tr/ome_sun_lmdz_tr 480 ome_sun_lmdz_tr=ome_sun_lmdz_tr/tau_sun_lmdz_tr 481 ; 482 ;--LW properties 483 tau_ear_tr=w1*tau_ear_lmdz+w2*tau_ear_lmdz_ave 484 ; 485 ;--prepare SW output 486 opticstruct={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 ; 492 opticstruct.lat=latitudelmdz 493 opticstruct.lev=lev 494 opticstruct.wav=(wl1_sun(0:NSW-1)+wl2_sun(0:NSW-1))/2. 495 opticstruct.time=float(indgen(month_in_year)+1) 496 opticstruct.tau_sun=tau_sun_lmdz_tr 497 opticstruct.ome_sun=ome_sun_lmdz_tr 498 opticstruct.ggg_sun=ggg_sun_lmdz_tr 499 ; 500 attributes = {units:strarr(7),long_name:strarr(7)} 501 attributes.units = ['degrees_north','level','meters','month','-','-','-'] 502 attributes.long_name = ['latitude','level','wavelength','time','tau_sun','ome_sun','g_sun'] 503 ; 504 dimensions = {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 ; 510 netcdfwrite,output+'tauswstrat.2D.'+chyr+'.nc',opticstruct,clobber=1, $ 511 attributes=attributes, dimensions=dimensions 512 ; 513 ;--prepare LW output 514 opticstruct={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 ; 518 opticstruct.lat=latitudelmdz 519 opticstruct.lev=lev 520 opticstruct.wav=(wl1_earth+wl2_earth)/2. 521 opticstruct.time=float(indgen(month_in_year)+1) 522 opticstruct.tau_ear=tau_ear_lmdz_tr 523 ; 524 attributes = {units:strarr(5),long_name:strarr(5)} 525 attributes.units = ['degrees_north','level','cm-1','month','-'] 526 attributes.long_name = ['latitude','level','wavenumber','time','tau_ear'] 527 ; 528 dimensions = {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 ; 534 netcdfwrite,output+'taulwstrat.2D.'+chyr+'.nc',opticstruct,clobber=1, $ 535 attributes=attributes, dimensions=dimensions 536 ; 537 endfor 538 ;--end loop on years 459 539 ; 460 540 end … … 479 559 EOF2 480 560 481 #qsub volc.job561 qsub volc.job
Note: See TracChangeset
for help on using the changeset viewer.