Changeset 3766 for TOOLS/CMIP6_FORCING/AER_STRAT
- Timestamp:
- 04/15/18 20:29:52 (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
TOOLS/CMIP6_FORCING/AER_STRAT/volc.sh
r3765 r3766 4 4 5 5 #--choose the resolution 6 #--only LR worksfor now !6 #--only LR and zoom works fully for now ! 7 7 #lmdz='VVLR' 8 8 #lmdz='VLR_L79' 9 9 #lmdz='VLR' 10 lmdz='LR'10 #lmdz='LR' 11 11 #lmdz='MR' 12 #--added for Frederique Cheruy 13 lmdz='zoom_128x89' 12 14 13 15 dirout='/data/'${USER}'/CMIP6/VOLC/'${lmdz}'_'${vv}'/' … … 29 31 dimz=19 30 32 output='/data/${USER}/CMIP6/VOLC/VVLR_${vv}/' 33 latfirst=90. 34 latinc=-180./float(dimlatlmdz-1) 35 latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 31 36 endif 32 37 … … 36 41 dimz=39 37 42 output='/data/${USER}/CMIP6/VOLC/VLR_${vv}/' 43 latfirst=90. 44 latinc=-180./float(dimlatlmdz-1) 45 latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 38 46 endif 39 47 … … 43 51 dimz=79 44 52 output='/data/${USER}/CMIP6/VOLC/VLR_L79_${vv}/' 53 latfirst=90. 54 latinc=-180./float(dimlatlmdz-1) 55 latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 45 56 endif 46 57 … … 50 61 dimz=79 51 62 output='/data/${USER}/CMIP6/VOLC/LR_${vv}/' 63 latfirst=90. 64 latinc=-180./float(dimlatlmdz-1) 65 latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 52 66 endif 53 67 … … 57 71 dimz=79 58 72 output='/data/${USER}/CMIP6/VOLC/MR_${vv}/' 59 endif60 61 73 latfirst=90. 62 74 latinc=-180./float(dimlatlmdz-1) 63 lonfirst=-180.64 loninc=360./float(dimlonlmdz-1)65 75 latitudelmdz=latfirst+latinc*indgen(dimlatlmdz) 66 longitudelmdz=lonfirst+loninc*indgen(dimlonlmdz) 76 endif 77 78 if (lmdz eq 'zoom_128x89') then begin 79 dimlonlmdz=128 80 dimlatlmdz=89 81 dimz=79 82 output='/data/${USER}/CMIP6/VOLC/${lmdz}_${vv}/' 83 latitudelmdz=[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.] 94 endif 95 67 96 NSW=6 68 97 NLW=16 … … 141 170 ; 142 171 ;---exact altitudes of LMDZ -- L79 143 ;--only exists for LR a t the moment172 ;--only exists for LR and zoom_128x89 at the moment 144 173 if (dimz eq 79) then begin 145 filename='./zalt_zonmean_ LR_l79_rev.nc'174 filename='./zalt_zonmean_${lmdz}_l79_rev.nc' 146 175 NETCDFREAD,filename,'GEOP',zz,dimzz 147 176 NETCDFREAD,filename,'LAT',zzlat,dimzzlat … … 209 238 ; 210 239 for year=0, dimtime/month_in_year-1 do begin 211 ;for year=141, 14 2do begin ;--Pinatubo240 ;for year=141, 141 do begin ;--Pinatubo 212 241 ; 213 242 tau_550_lmdz(*,*,*)=1.e-15
Note: See TracChangeset
for help on using the changeset viewer.