## Postprocessing of files of coupled CESM version CESM version 1.2 runs for use in LIVVkit
(note: a test version for E3SM data 'main_e3sm.sh' is also here, but its not production yet)

This directory holds subdirectories with scripts that will postprocess model data for use in creating the plots and data used by LIVVkit. 

### What is needed to run postproc scripts:
	*Create a working directory (with write permissions) to put all target files
	*Refer to check list below to ensure all the files are accessable 
	*Run on OLCFs rhea system, to ensure enough memory for postprocessing CISM-Albany data

### Files needed to run coupled model postprocessing (right now a copy exists in the $WORLDWORK/cli115//4ue space)

1. b.e10.BG20TRCN.f09_g16.002
	* /lnd/hist
	* /atm/hist
	* /glc/hist
2. surfdata_0.9x1.25_simyr1850_c110725.nc
3. /obs_data
	* /CLOUDSAT_*	for all months, DJF, JJA, and ANN
	* /ISCCP_*    for all months, DJF, JJA, and ANN	
4. /racmo23_GRN_monthly (all this data is already separated by variable )
5. Greenland_5km_v1.1_SacksRev_c110629.nc 

### Main script for analysis 

#### Usage:
	*Copy over main.sh to separate file (e.g. main_deal.sh) 
	*Set parameters at the top of main_deal.sh to match your environment
	*Below that, opt out, if desired, to postprocessing, remapping, and plots/tables

```sh
 bash main_deal.sh 
```

#### Postprocessing
* loads computing environment for rhea, with hooks to add other machines
* makes necessary mask files if they do not already exist
* if 'postprocess_obs' is selected
	* calls atm_obs.sh
		* remaps ISCCP and CLOUDSAT climo data to model grid for selected variables for plots
		* makes ice sheet area averaged climos from original climo data
	* calls lnd_obs.sh 
		* makes climos from original racmo climo data
		* makes yearly time series from original racmo monthly files 
		* makes ice sheet area average (but not ice sheet %) of yearly time series of racmo data
* calls atm.sh
	* makes climos and ice sheet area averages
	* calls monthly_atm_ts.sh 
		* make monthly and yearly averaged time series 
		* make ice sheet area average of time series 
* calls lnd.sh
	* makes climos and ice area average of climos
 	* makes daily avergages over time record 
	* submits rhea_tseries_lnd.pbs  
		* makes time series of daily avg data
		* makes ice sheet average of daily time series	
	* makes annualized time series of daily avg data
	* makes ice sheet average of annualized daily time series	
	* calls tseries_lnd_h0.sh
		* makes time series from monthly averaged data
		* makes yearly averages of monthly averaged data
		* makes ice sheet average of monthly averaged data 
* calls glc.sh 
	* select flag to compute velnorm (turn off if velocities do not exist)
	* makes annual climos
	* makes ice sheet average of climos
	* makes time series of yearly averaged data
	* area average time series (it is not masked right now)


#### Make plots and tables
* contour plots of 2m T MODEL, RACMO and difference for DJF and JJA
* plots of ANN low, med, and high clouds for MODEL, ISCCP, and CLOUDSAT
* plots of the annual cycle of low, med, and high clouds for MODEL, ISCCP, and CLOUDSAT
* plot of time series of 2m T MODEL and RACMO
* print ice sheet are and time averages of suite of variables for model and RACMO
* time series of downscaled CISM SMB