Changeset 3766

Timestamp:

04/15/18 20:29:52 (6 years ago)

Author:

oboucher

Message:

adding zoom case

File:

• TOOLS/CMIP6_FORCING/AER_STRAT/volc.sh (modified) (8 diffs)

TOOLS/CMIP6_FORCING/AER_STRAT/volc.sh

@@ -4 +4 @@
 
 #--choose the resolution
-#--only LR works for now !
+#--only LR and zoom works fully for now !
 #lmdz='VVLR'
 #lmdz='VLR_L79'
 #lmdz='VLR'
-lmdz='LR'
+#lmdz='LR'
 #lmdz='MR'
+#--added for Frederique Cheruy
+lmdz='zoom_128x89'
 
 dirout='/data/'${USER}'/CMIP6/VOLC/'${lmdz}'_'${vv}'/'
@@ -29 +31 @@
 dimz=19
 output='/data/${USER}/CMIP6/VOLC/VVLR_${vv}/'
+latfirst=90.
+latinc=-180./float(dimlatlmdz-1)
+latitudelmdz=latfirst+latinc*indgen(dimlatlmdz)
 endif
 
@@ -36 +41 @@
 dimz=39
 output='/data/${USER}/CMIP6/VOLC/VLR_${vv}/'
+latfirst=90.
+latinc=-180./float(dimlatlmdz-1)
+latitudelmdz=latfirst+latinc*indgen(dimlatlmdz)
 endif
 
@@ -43 +51 @@
 dimz=79
 output='/data/${USER}/CMIP6/VOLC/VLR_L79_${vv}/'
+latfirst=90.
+latinc=-180./float(dimlatlmdz-1)
+latitudelmdz=latfirst+latinc*indgen(dimlatlmdz)
 endif
 
@@ -50 +61 @@
 dimz=79
 output='/data/${USER}/CMIP6/VOLC/LR_${vv}/'
+latfirst=90.
+latinc=-180./float(dimlatlmdz-1)
+latitudelmdz=latfirst+latinc*indgen(dimlatlmdz)
 endif
 
@@ -57 +71 @@
 dimz=79
 output='/data/${USER}/CMIP6/VOLC/MR_${vv}/'
-endif
-
 latfirst=90.
 latinc=-180./float(dimlatlmdz-1)
-lonfirst=-180.
-loninc=360./float(dimlonlmdz-1)
 latitudelmdz=latfirst+latinc*indgen(dimlatlmdz)
-longitudelmdz=lonfirst+loninc*indgen(dimlonlmdz)
+endif
+
+if (lmdz eq 'zoom_128x89') then begin
+dimlonlmdz=128
+dimlatlmdz=89
+dimz=79
+output='/data/${USER}/CMIP6/VOLC/${lmdz}_${vv}/'
+latitudelmdz=[90., 86.85495, 83.8430634, 80.9792252, 78.2779922, 75.7526321, 73.4137878, 71.2680206, 69.3164978, $
+              67.5542145,   65.9702301, 64.5487747, 63.2709846, 62.1167107, 61.0660896, 60.1005592, 59.2034454, $ 
+              58.360157,    57.5581551, 56.7868576, 56.0374718, 55.302887,  54.5775719, 53.8574829, 53.1399384, $
+              52.4233971,   51.7071419, 50.9909439, 50.2747459, 49.5585518, 48.8423538, 48.1261597, 47.4099655, $
+              46.6937675,   45.9775734, 45.2613792, 44.5451813, 43.8289604, 43.1125679, 42.395504,  41.6765518, $ 
+              40.9534645,   40.2226982, 39.4792557, 38.7165909, 37.9265213, 37.099102,  36.2225075, 35.2828789, $ 
+              34.2642632,   33.1487236, 31.9167614, 30.5481758, 29.0234509, 27.32551,   25.4415455, 23.3643932, $ 
+              21.0930786,   18.6323872, 15.9917059, 13.1835403, 10.2221146, 7.12228632, 3.89880991, 0.565921485,$
+              -2.86286402, -6.37483072, -9.95808029, -13.6015415, -17.2949867, -21.029068, -24.7953758, -28.5864811, $ 
+             -32.3959961, -36.2186012, -40.0500298, -43.8870316, -47.7272758, -51.5692215, -55.4119453, -59.2549667, $
+             -63.0980797, -66.9412079, -70.7843399, -74.6274719, -78.4706039, -82.313736, -86.156868, -90.]
+endif
+
 NSW=6
 NLW=16
@@ -141 +170 @@
 ;
 ;---exact altitudes of LMDZ -- L79
-;--only exists for LR at the moment
+;--only exists for LR and zoom_128x89 at the moment
 if (dimz eq 79) then begin
-  filename='./zalt_zonmean_LR_l79_rev.nc'
+  filename='./zalt_zonmean_${lmdz}_l79_rev.nc'
   NETCDFREAD,filename,'GEOP',zz,dimzz
   NETCDFREAD,filename,'LAT',zzlat,dimzzlat
@@ -209 +238 @@
 ;
 for year=0, dimtime/month_in_year-1 do begin
-;for year=141, 142 do begin ;--Pinatubo
+;for year=141, 141 do begin ;--Pinatubo
 ;
 tau_550_lmdz(*,*,*)=1.e-15