Changeset 3766


Ignore:
Timestamp:
04/15/18 20:29:52 (2 years ago)
Author:
oboucher
Message:

adding zoom case

File:
1 edited

Legend:

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

    r3765 r3766  
    44 
    55#--choose the resolution 
    6 #--only LR works for now ! 
     6#--only LR and zoom works fully for now ! 
    77#lmdz='VVLR' 
    88#lmdz='VLR_L79' 
    99#lmdz='VLR' 
    10 lmdz='LR' 
     10#lmdz='LR' 
    1111#lmdz='MR' 
     12#--added for Frederique Cheruy 
     13lmdz='zoom_128x89' 
    1214 
    1315dirout='/data/'${USER}'/CMIP6/VOLC/'${lmdz}'_'${vv}'/' 
     
    2931dimz=19 
    3032output='/data/${USER}/CMIP6/VOLC/VVLR_${vv}/' 
     33latfirst=90. 
     34latinc=-180./float(dimlatlmdz-1) 
     35latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 
    3136endif 
    3237 
     
    3641dimz=39 
    3742output='/data/${USER}/CMIP6/VOLC/VLR_${vv}/' 
     43latfirst=90. 
     44latinc=-180./float(dimlatlmdz-1) 
     45latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 
    3846endif 
    3947 
     
    4351dimz=79 
    4452output='/data/${USER}/CMIP6/VOLC/VLR_L79_${vv}/' 
     53latfirst=90. 
     54latinc=-180./float(dimlatlmdz-1) 
     55latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 
    4556endif 
    4657 
     
    5061dimz=79 
    5162output='/data/${USER}/CMIP6/VOLC/LR_${vv}/' 
     63latfirst=90. 
     64latinc=-180./float(dimlatlmdz-1) 
     65latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 
    5266endif 
    5367 
     
    5771dimz=79 
    5872output='/data/${USER}/CMIP6/VOLC/MR_${vv}/' 
    59 endif 
    60  
    6173latfirst=90. 
    6274latinc=-180./float(dimlatlmdz-1) 
    63 lonfirst=-180. 
    64 loninc=360./float(dimlonlmdz-1) 
    6575latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 
    66 longitudelmdz=lonfirst+loninc*indgen(dimlonlmdz) 
     76endif 
     77 
     78if (lmdz eq 'zoom_128x89') then begin 
     79dimlonlmdz=128 
     80dimlatlmdz=89 
     81dimz=79 
     82output='/data/${USER}/CMIP6/VOLC/${lmdz}_${vv}/' 
     83latitudelmdz=[90., 86.85495, 83.8430634, 80.9792252, 78.2779922, 75.7526321, 73.4137878, 71.2680206, 69.3164978, $ 
     84              67.5542145,   65.9702301, 64.5487747, 63.2709846, 62.1167107, 61.0660896, 60.1005592, 59.2034454, $  
     85              58.360157,    57.5581551, 56.7868576, 56.0374718, 55.302887,  54.5775719, 53.8574829, 53.1399384, $ 
     86              52.4233971,   51.7071419, 50.9909439, 50.2747459, 49.5585518, 48.8423538, 48.1261597, 47.4099655, $ 
     87              46.6937675,   45.9775734, 45.2613792, 44.5451813, 43.8289604, 43.1125679, 42.395504,  41.6765518, $  
     88              40.9534645,   40.2226982, 39.4792557, 38.7165909, 37.9265213, 37.099102,  36.2225075, 35.2828789, $  
     89              34.2642632,   33.1487236, 31.9167614, 30.5481758, 29.0234509, 27.32551,   25.4415455, 23.3643932, $  
     90              21.0930786,   18.6323872, 15.9917059, 13.1835403, 10.2221146, 7.12228632, 3.89880991, 0.565921485,$ 
     91              -2.86286402, -6.37483072, -9.95808029, -13.6015415, -17.2949867, -21.029068, -24.7953758, -28.5864811, $  
     92             -32.3959961, -36.2186012, -40.0500298, -43.8870316, -47.7272758, -51.5692215, -55.4119453, -59.2549667, $ 
     93             -63.0980797, -66.9412079, -70.7843399, -74.6274719, -78.4706039, -82.313736, -86.156868, -90.] 
     94endif 
     95 
    6796NSW=6 
    6897NLW=16 
     
    141170; 
    142171;---exact altitudes of LMDZ -- L79 
    143 ;--only exists for LR at the moment 
     172;--only exists for LR and zoom_128x89 at the moment 
    144173if (dimz eq 79) then begin 
    145   filename='./zalt_zonmean_LR_l79_rev.nc' 
     174  filename='./zalt_zonmean_${lmdz}_l79_rev.nc' 
    146175  NETCDFREAD,filename,'GEOP',zz,dimzz 
    147176  NETCDFREAD,filename,'LAT',zzlat,dimzzlat 
     
    209238; 
    210239for year=0, dimtime/month_in_year-1 do begin 
    211 ;for year=141, 142 do begin ;--Pinatubo 
     240;for year=141, 141 do begin ;--Pinatubo 
    212241; 
    213242tau_550_lmdz(*,*,*)=1.e-15 
Note: See TracChangeset for help on using the changeset viewer.