Changeset 4605 for TOOLS/CMIP6_FORCING/AER_STRAT/volc.sh

Timestamp:

07/04/19 18:56:39 (5 years ago)

Author:

oboucher

Message:

netcdf files are now on /bdd/input4MIPs/ZALT
volc.sh is modified accordingly
zalt files are now compulsory

File:

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

TOOLS/CMIP6_FORCING/AER_STRAT/volc.sh

@@ -1 +1 @@
 dirpwd=`pwd`
 
+#original CMIP6 files
+vv='v3'
+#corrected CMIP6 files + extended in time
 vv='v4'
 
 #--choose the resolution
 #--only LR and zoom works fully for now !
-#lmdz='VVLR'
-#lmdz='VLR_L79'
-#lmdz='VLR'
-lmdz='LR'
-#lmdz='MR'
+#lmdz='VVLR_l19'
+lmdz='VLR_l39'
+#lmdz='VLR_l79'
+#lmdz='LR_l79'
+#lmdz='MR_l79'
 #--added for Frederique Cheruy
-#lmdz='zoom_128x89'
+#lmdz='zoom_128x89_l79'
 
 dirout='/data/'${USER}'/CMIP6/VOLC/'${lmdz}'_'${vv}'/'
@@ -26 +29 @@
 lmdz='${lmdz}'
 ;
-if (lmdz eq 'VVLR') then begin
+if (lmdz eq 'VVLR_l19') then begin
 dimlonlmdz=48
 dimlatlmdz=37
 dimz=19
-output='/data/${USER}/CMIP6/VOLC/VVLR_${vv}/'
+output='/data/${USER}/CMIP6/VOLC/${lmdz}_${vv}/'
 latfirst=90.
 latinc=-180./float(dimlatlmdz-1)
@@ -36 +39 @@
 endif
 
-if (lmdz eq 'VLR') then begin
+if (lmdz eq 'VLR_l39') then begin
 dimlonlmdz=96
 dimlatlmdz=96
 dimz=39
-output='/data/${USER}/CMIP6/VOLC/VLR_${vv}/'
+output='/data/${USER}/CMIP6/VOLC/${lmdz}_${vv}/'
 latfirst=90.
 latinc=-180./float(dimlatlmdz-1)
@@ -46 +49 @@
 endif
 
-if (lmdz eq 'VLR_L79') then begin
+if (lmdz eq 'VLR_l79') then begin
 dimlonlmdz=96
 dimlatlmdz=96
 dimz=79
-output='/data/${USER}/CMIP6/VOLC/VLR_L79_${vv}/'
+output='/data/${USER}/CMIP6/VOLC/${lmdz}_${vv}/'
 latfirst=90.
 latinc=-180./float(dimlatlmdz-1)
@@ -56 +59 @@
 endif
 
-if (lmdz eq 'LR') then begin
+if (lmdz eq 'LR_l79') then begin
 dimlonlmdz=144
 dimlatlmdz=143
 dimz=79
-output='/data/${USER}/CMIP6/VOLC/LR_${vv}/'
+output='/data/${USER}/CMIP6/VOLC/${lmdz}_${vv}/'
 latfirst=90.
 latinc=-180./float(dimlatlmdz-1)
@@ -66 +69 @@
 endif
 
-if (lmdz eq 'MR') then begin
+if (lmdz eq 'MR_l79') then begin
 dimlonlmdz=256
 dimlatlmdz=257
 dimz=79
-output='/data/${USER}/CMIP6/VOLC/MR_${vv}/'
+output='/data/${USER}/CMIP6/VOLC/${lmdz}_${vv}/'
 latfirst=90.
 latinc=-180./float(dimlatlmdz-1)
@@ -76 +79 @@
 endif
 
-if (lmdz eq 'zoom_128x89') then begin
+if (lmdz eq 'zoom_128x89_l79') then begin
 dimlonlmdz=128
 dimlatlmdz=89
@@ -152 +155 @@
 zz=fltarr(dimz)
 ;
-;---approximate altitudes of LMDZ - L19 and L39
-;---NOT RECOMMENDED
-;
-if (dimz eq 19) then begin
-zz=[ 0.07125416, 0.2402302, 0.4865242, 0.8627044, 1.426484, 2.244528,    $
-    3.394085, 4.952366, 6.950302, 9.26943, 11.59338, 13.71214, 15.83651, $
-    18.27338, 21.16397, 24.63156, 28.9039, 34.57184, 44.51565] ;
-endif
-;
-if (dimz eq 39) then begin
-zz=[0.0338988, 0.1106602, 0.2139233, 0.3609663, 0.5685416, 0.8534464,     $
-    1.233232, 1.726868, 2.355076, 3.139813, 4.101937, 5.25541, 6.596076,  $
-    8.085103, 9.636482, 11.13441, 12.50023, 13.75765, 15.00424, 16.32207, $
-    17.74403, 19.27899, 20.93219, 22.70941, 24.61736, 26.66389, 28.85831, $
-    31.21178, 33.73778, 36.45278, 39.37709, 42.53607, 45.96209, 49.69786, $
-    53.80379, 58.37727, 63.61497, 70.07092, 80.52708]
-endif
-;
 ;---exact altitudes of LMDZ -- L79
-;--only exists for LR, MR and zoom_128x89 at the moment
-if (dimz eq 79) then begin
-  filename='./zalt_zonmean_${lmdz}_l79_rev.nc'
-  NETCDFREAD,filename,'GEOP',zz,dimzz
-  NETCDFREAD,filename,'LAT',zzlat,dimzzlat
-  NETCDFREAD,filename,'TIME_COUNTER',zztime,dimzztime
-  dimzzilat=dimzzlat[0]
-  dimzzitime=dimzztime[0]
-  ;--becareful zz comes with four dimensions
-  ;--lon lat k time 
-  print, 'GEOP size=', size(zz)
-  ;--becareful zzlat comes with South Pole first
-  print, 'LAT from zalt field=', zzlat
-  if (dimzzilat ne dimlatlmdz) then begin
-    print , 'PB dimension latitude'
-  endif
-endif
-;
-if (dimz eq 39) then begin
-  filename='./zalt_zonmean_${lmdz}_l39_rev.nc'
-  NETCDFREAD,filename,'GEOP',zz,dimzz
-  NETCDFREAD,filename,'LAT',zzlat,dimzzlat
-  NETCDFREAD,filename,'TIME_COUNTER',zztime,dimzztime
-  dimzzilat=dimzzlat[0]
-  dimzzitime=dimzztime[0]
-  ;--becareful zz comes with four dimensions
-  ;--lon lat k time
-  print, 'GEOP size=', size(zz)
-  ;--becareful zzlat comes with South Pole first
-  print, 'LAT from zalt field=', zzlat
-  if (dimzzilat ne dimlatlmdz) then begin
-    print , 'PB dimension latitude'
-  endif
+;--only exists for VLR, LR, MR and zoom_128x89 at the moment
+;
+filename='/bdd/input4MIPs/ZALT/zalt_zonmean_${lmdz}_rev.nc'
+NETCDFREAD,filename,'GEOP',zz,dimzz
+NETCDFREAD,filename,'LAT',zzlat,dimzzlat
+NETCDFREAD,filename,'TIME_COUNTER',zztime,dimzztime
+dimzzilat=dimzzlat[0]
+dimzzitime=dimzztime[0]
+;--becareful zz comes with four dimensions
+;--lon lat k time 
+print, 'GEOP size=', size(zz)
+;--becareful zzlat comes with South Pole first
+print, 'LAT from zalt field=', zzlat
+if (dimzzilat ne dimlatlmdz) then begin
+  print , 'PB dimension latitude'
 endif
 ;
@@ -578 +545 @@
 cd $dirpwd
 
-module load idl
-
 idl << eof
 .r netcdf