Commit 86169053 authored by Neven S. Fuckar's avatar Neven S. Fuckar
Browse files

more general approach to siasiesiv_Arcticreg1

parent 7089a230
......@@ -63,8 +63,8 @@ ln -sf ${CON_FILES}/mesh_mask_nemo.${NEMOVERSION}.nc mesh_zgr.nc
ln -sf ${CON_FILES}/mesh_mask_nemo.${NEMOVERSION}.nc mask.nc
ln -sf ${CON_FILES}/new_maskglo.${NEMOVERSION}.nc new_maskglo.nc
if [ -e ${CON_FILES}/mask.ArcticSeas.${NEMOVERSION}.nc ] ; then
cp ${CON_FILES}/mask.ArcticSeas.${NEMOVERSION}.nc mask_regions1.nc
if [ -e ${CON_FILES}/mask.regions.${NEMOVERSION}.nc ] ; then
cp ${CON_FILES}/mask.regions.${NEMOVERSION}.nc mask_regions.nc
fi
###############################################################################
......@@ -770,9 +770,9 @@ fi
# Created in April 2012 Author : vguemas@ic3.cat #
# Modified in June 2014 Author : neven.fuckar@ic3.cat #
# #
# Computation of the properties in 34 Arctic regions according to #
# mask.ArcticSeas.ORCA1.nc (mask_regions1.nc) is based on modification of #
# mask.regions.ORCA1.noverticalinfo.Matt.nc from Matthieu Chevallier. #
# Computation of the properties in various selected regions according to #
# mask.regions.${NEMOVERSION}.nc (mask_regions.nc) is based on modification #
# of mask.regions.ORCA1.noverticalinfo.Matt.nc from Matthieu Chevallier. #
# #
###############################################################################
......@@ -796,7 +796,7 @@ typeset var jt
if [ ! -z "$3" ] ; then
mv mask.nc mask_tmp.nc
mv mask_regions1.nc mask.nc
mv mask_regions.nc mask.nc
ncrename -h -v $3,tmask mask.nc
fi
......@@ -828,7 +828,7 @@ setminmax ice_S_${2} sia sie siv
if [ ! -z "$3" ] ; then
ncrename -h -v tmask,$3 mask.nc
mv mask.nc mask_regions1.nc
mv mask.nc mask_regions.nc
mv mask_tmp.nc mask.nc
fi
......
......@@ -6,6 +6,7 @@ listpost=( )
# Beware that the max_moc diagnostics can not be computed if you don't process complete years (because that's a diagnostic computed from annual means)
raw_vars_ocean=( 'default' ) # If listpost contains "ext_raw_oce" option, this is the list ocean variables you want to extract. If nothing is specified, all variables present in input file will be treated. If raw_vars_ocean='default', sosstsst, sosaline, somixhgt and somxl010 will be extracted.
raw_vars_ice=( 'default' ) # If listpost contains "ext_raw_ice" option, this is the list of seaice variables you want to extract. If nothing is specified, all variables will be treated. If raw_vars_ice='default', isnowthi, iicethic, ileadfra, iicetemp, and ice_pres will be extracted.
raw_regions_ice=( 'default' ) # If listpost contains "siasiesiv_Arcticreg1" option, this is the list of seaice regional variables you want to extract. If nothing is specified, all regions will be treated. If raw_regions_ice='default', YET TO BE DECIDED (NOW ONLY THREE TEST SEAS) will be preocessed.
expid=i00k # expid or nemovar_s4 / nemovar_combine / glorys2v1
mod='ecearth' # nemo / ecearth
typeoutput='MMO' # diags / MMO
......@@ -18,7 +19,7 @@ intsdate=1 # interval between start dates
chunklen=4 # length of the chunks (in months)
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ltime0=1 # first leadtime to post-process
ltimef=24 # last leadtime to postprocess
ltimef=16 # last leadtime to postprocess
# Fill up either ltime0/ltimef or year0/yearf
year0= # first year to post-process in the fist start date
yearf= # last year to post-process in the fist start date
......
......@@ -347,7 +347,17 @@ for ((yeari=$syeari;yeari<=$syearf;yeari=$(($yeari+intsdate)))) ; do
'MMO' ) pref='icemod' ;;
'diags') pref='ice' ;;
esac
lstseas=$( cdo showvar mask_regions1.nc )
if [[ $raw_regions_ice == '' ]] ; then
lstseas=$( cdo showvar mask_regions.nc )
else
if [[ $raw_regions_ice == 'default' ]] ; then
# defaults options needs to be updated properly, these three seas are only for a test case
lstseas="Nordic_Seas East_Siberian_Sea Barents_Sea"
else
lstseas=$( echo ${raw_regions_ice[@]} )
fi
fi
for sea in $lstseas ; do
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment