trun.sh - pism - [fork] customized build of PISM, the parallel ice sheet model (tillflux branch)
 (HTM) git clone git://src.adamsgaard.dk/pism
 (DIR) Log
 (DIR) Files
 (DIR) Refs
 (DIR) LICENSE
       ---
       trun.sh (1653B)
       ---
            1 #!/bin/bash
            2 
            3 python create_inputs.py
            4 
            5 run() {
            6 suffix=$1
            7 warming=$2
            8 extra_vars=$3
            9 
           10 # bootstrap ice temperature (and so its enthalpy) using the mean-annual surface
           11 # temperature
           12 pismr -bootstrap -i input.nc \
           13       -regional \
           14       -regional.no_model_strip 0 \
           15       -Lz 200 -Mz 201 -grid.ice_vertical_spacing equal \
           16       -surface given \
           17       -y 0 \
           18       -o in.nc -verbose 1
           19 
           20 echo "running PISM (${suffix})..."
           21 
           22 # bootstrap again and bring in the enthalpy field computed by the run above
           23 pismr -bootstrap -i input.nc -regrid_file in.nc -regrid_vars enthalpy \
           24       -bootstrapping.defaults.geothermal_flux 0 \
           25       -regional \
           26       -regional.no_model_strip 0 \
           27       -Lz 200 -Mz 201 -grid.ice_vertical_spacing equal \
           28       -energy.ch_warming.enabled ${warming} \
           29       -energy.ch_warming.average_channel_spacing 20 \
           30       -energy.ch_warming.residual_water_fraction 0.005 \
           31       -energy.ch_warming.temperate_ice_thermal_conductivity_ratio 1.0 \
           32       -surface given,delta_T \
           33       -surface_delta_T_file input.nc -surface_delta_T_period 1 \
           34       -extra_file ex_${suffix}.nc \
           35       -extra_vars ${extra_vars} \
           36       -extra_times 10days \
           37       -y 10 \
           38       -calendar 360_day \
           39       -o o_${suffix}.nc -verbose 1
           40 
           41 rm in.nc
           42 }
           43 
           44 variables=temp,temp_pa,ice_surface_temp
           45 
           46 # run without cryo-hydrologic warming
           47 run no_warming False ${variables}
           48 
           49 # run with cryo-hydrologic warming
           50 run warming True ${variables},ch_heat_flux,ch_temp,ch_liqfrac
           51 
           52 ncpdq -a z,time,y,x -O ex_no_warming.nc ex_no_warming.nc
           53 ncpdq -a z,time,y,x -O ex_warming.nc ex_warming.nc
           54 
           55 # compare temperature fields
           56 ncdiff -O -v temp ex_warming.nc ex_no_warming.nc temp_difference.nc