diff --git a/util/convert_inflow.com b/util/convert_inflow.com index 2d910786363e9acfd6aaf05230211662ca7d3403..32bffded9e3f663af1325ccb79d6ee922de933d1 100755 --- a/util/convert_inflow.com +++ b/util/convert_inflow.com @@ -133,7 +133,8 @@ fi # # *** Is mask file for the mouths on the source grid provided. if [ "$dn_src_mouth" = "" ]; then - echo " No mouth mask on source grid provided - defintion in script expected " + echo " No mouth mask on source grid provided " + echo " Mask is obtained from HD parameter variable FDIR, or it must be defined in the script" else if [ -s $dn_ocean ] ; then echo " Mouth mask on source grid: " $dn_src_mouth @@ -514,7 +515,8 @@ esac # *** Allocate HD mouths ${F90} ${DSRC}/mo_grid.f90 ${DSRC}/mo_time.f90 ${DSRC}/mo_flow_inout.f90 ${DSRC}/mo_interpol.f90 ${DSRC}/mo_convert.f90 ${DSRC}/convert_discharge.f90 -o convert.exe $NC_INCLUDE $NC_LIB # -case $ISRC in +if [ "$dn_src_mouth" = "" ]; then + case $ISRC in [123] ) cdo setctomiss,0. -eqc,0. -selvar,FDIR $DNHDPARA hdmouth_mask.nc cdo setvar,FMOUTH hdmouth_mask.nc rivmouth_source.nc rm hdmouth_mask.nc ;; @@ -533,7 +535,18 @@ case $ISRC in rm t[12345].nc ;; [67] ) cdo setvar,FMOUTH $DNHDPARA rivmouth_source.nc ;; * ) echo 'Treatment of mouth mask needs to be specified -> STOP' ; exit ;; -esac + esac +else + case $ISRC in + [1234] ) cdo -selvar,FDIR $DNHDPARA ts.nc + cdo setvar,FMOUTH -setgrid,ts.nc $dn_src_mouth rivmouth_source.nc + rm ts.nc ;; + 5 ) cdo -selvar,lon $DNHDPARA ts.nc + cdo setvar,FMOUTH -setgrid,ts.nc $dn_src_mouth rivmouth_source.nc + rm ts.nc ;; + * ) cdo setvar,FMOUTH $dn_src_mouth rivmouth_source.nc ;; + esac +fi # # *** Run the Program to create data that will be used for the remapping of mouth points # *** --> Creates files hd_to_ocean_mouth.nc, hdmouth_on_oceangrid.nc @@ -549,7 +562,18 @@ fi # # Add script settings to coupling file ncatted -O -h -a Settings,global,o,c,"IMODE = $IMODE, ISRC=$ISRC, IOCEAN=$IOCEAN, IBGC=$IBGC" $DNOUT - +if [ "$dn_src_mouth" != "" ]; then + ncatted -O -h -a Source_mask,global,o,c,"River mouth mask on source grid: $dn_src_mouth" $DNOUT +fi +YEAR=$YBEG +while [ $YEAR -le $YEND ] ; do + DNOUTFLOW=${DATA_OUT}/${CFLOW}_${CDIS}_on_${OM}${OVS}_${YEAR}.nc + ncatted -O -h -a Settings,global,o,c,"IMODE = $IMODE, ISRC=$ISRC, IOCEAN=$IOCEAN, IBGC=$IBGC" $DNOUTFLOW + if [ "$dn_src_mouth" != "" ]; then + ncatted -O -h -a Source_mask,global,o,c,"River mouth mask on source grid: $dn_src_mouth" $DNOUTFLOW + fi + YEAR=`expr $YEAR + 1` +done # if [ $ISRC -eq 6 ] ; then mv $DNOUT ${DIN}/${OM}/${EXP}_$DNOUT