{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"# instructions\n",
"* Wenchang Yang (wenchang@prnceton.edu)\n",
"* Department of Geosciences, Princeton University"
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {
"ExecuteTime": {
"end_time": "2020-02-25T15:54:12.459756Z",
"start_time": "2020-02-25T15:54:12.448734Z"
},
"code_folding": []
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"\n",
"[imported]: os.path, sys, os\n",
"[imported]: xarray as xr, numpy as np, matplotlib.pyplot as plt, pandas as pd\n",
"[imported]: from matplotlib.pyplot import figure, close\n",
"[executed]: plt.ion()\n",
"[config]: plt.rcParams['figure.dpi'] = 128\n",
"[iPython config]: InlineBackend.figure_format ='retina'\n"
]
}
],
"source": [
"# init\n",
"%matplotlib inline\n",
"import sys\n",
"wython = '/tigress/wenchang/wython'\n",
"if wython not in sys.path:\n",
" sys.path.append(wython)\n",
"%run -im wystart\n",
"# https://www.dataquest.io/blog/jupyter-notebook-tips-tricks-shortcuts/\n",
"# from IPython.core.interactiveshell import InteractiveShell\n",
"# InteractiveShell.ast_node_interactivity = 'last' # options: 'all', 'none'\n",
"# https://stackoverflow.com/questions/41125690/matplotlib-notebook-showing-a-blank-histogram\n",
"# in case to switch back to the notebook backend, run the 2 lines below, \n",
"# reload(plt)\n",
"# %matplotlib notebook\n",
"# import xaddon, xfilter, xtc"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## data info"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
" ### 00010101.ice_month.nc"
]
},
{
"cell_type": "code",
"execution_count": 6,
"metadata": {
"ExecuteTime": {
"end_time": "2020-02-25T15:58:23.029008Z",
"start_time": "2020-02-25T15:58:22.883519Z"
},
"scrolled": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"/tigress/wenchang/MODEL_OUT/CTL1860_noleap_tigercpu_intelmpi_18_576PE/POSTP/00010101.ice_month.nc\n",
"netcdf \\00010101.ice_month {\r\n",
"dimensions:\r\n",
"\ttime = UNLIMITED ; // (12 currently)\r\n",
"\tyt = 200 ;\r\n",
"\txt = 360 ;\r\n",
"\tct = 5 ;\r\n",
"\tyv = 200 ;\r\n",
"\txv = 360 ;\r\n",
"\tnv = 2 ;\r\n",
"\txb = 361 ;\r\n",
"\tyb = 201 ;\r\n",
"variables:\r\n",
"\tfloat ALB(time, yt, xt) ;\r\n",
"\t\tALB:long_name = \"surface albedo\" ;\r\n",
"\t\tALB:units = \"0-1\" ;\r\n",
"\t\tALB:missing_value = -1.e+34f ;\r\n",
"\t\tALB:_FillValue = -1.e+34f ;\r\n",
"\t\tALB:cell_methods = \"time: mean\" ;\r\n",
"\t\tALB:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat BHEAT(time, yt, xt) ;\r\n",
"\t\tBHEAT:long_name = \"ocean to ice heat flux\" ;\r\n",
"\t\tBHEAT:units = \"W/m^2\" ;\r\n",
"\t\tBHEAT:missing_value = -1.e+34f ;\r\n",
"\t\tBHEAT:_FillValue = -1.e+34f ;\r\n",
"\t\tBHEAT:cell_methods = \"time: mean\" ;\r\n",
"\t\tBHEAT:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat BMELT(time, yt, xt) ;\r\n",
"\t\tBMELT:long_name = \"bottom surface melting energy flux\" ;\r\n",
"\t\tBMELT:units = \"W/m^2\" ;\r\n",
"\t\tBMELT:missing_value = -1.e+34f ;\r\n",
"\t\tBMELT:_FillValue = -1.e+34f ;\r\n",
"\t\tBMELT:cell_methods = \"time: mean\" ;\r\n",
"\t\tBMELT:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat BSNK(time, yt, xt) ;\r\n",
"\t\tBSNK:long_name = \"frozen water local bottom sink\" ;\r\n",
"\t\tBSNK:units = \"kg/(m^2*yr)\" ;\r\n",
"\t\tBSNK:missing_value = -1.e+34f ;\r\n",
"\t\tBSNK:_FillValue = -1.e+34f ;\r\n",
"\t\tBSNK:cell_methods = \"time: mean\" ;\r\n",
"\t\tBSNK:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat CALVING(time, yt, xt) ;\r\n",
"\t\tCALVING:long_name = \"frozen runoff\" ;\r\n",
"\t\tCALVING:units = \"kg/(m^2*s)\" ;\r\n",
"\t\tCALVING:missing_value = -1.e+34f ;\r\n",
"\t\tCALVING:_FillValue = -1.e+34f ;\r\n",
"\t\tCALVING:cell_methods = \"time: mean\" ;\r\n",
"\t\tCALVING:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat CELL_AREA(yt, xt) ;\r\n",
"\t\tCELL_AREA:long_name = \"cell area\" ;\r\n",
"\t\tCELL_AREA:units = \"sphere\" ;\r\n",
"\t\tCELL_AREA:missing_value = 1.e+20f ;\r\n",
"\t\tCELL_AREA:_FillValue = 1.e+20f ;\r\n",
"\t\tCELL_AREA:cell_methods = \"time: point\" ;\r\n",
"\tfloat CN(time, ct, yt, xt) ;\r\n",
"\t\tCN:long_name = \"ice concentration\" ;\r\n",
"\t\tCN:units = \"0-1\" ;\r\n",
"\t\tCN:missing_value = -1.e+34f ;\r\n",
"\t\tCN:_FillValue = -1.e+34f ;\r\n",
"\t\tCN:cell_methods = \"time: mean\" ;\r\n",
"\t\tCN:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat COSROT(yt, xt) ;\r\n",
"\t\tCOSROT:long_name = \"COSROT,SINROT points east\" ;\r\n",
"\t\tCOSROT:units = \"none\" ;\r\n",
"\t\tCOSROT:missing_value = 1.e+20f ;\r\n",
"\t\tCOSROT:_FillValue = 1.e+20f ;\r\n",
"\t\tCOSROT:cell_methods = \"time: point\" ;\r\n",
"\tfloat EXT(time, yt, xt) ;\r\n",
"\t\tEXT:long_name = \"ice modeled\" ;\r\n",
"\t\tEXT:units = \"0 or 1\" ;\r\n",
"\t\tEXT:missing_value = -1.e+34f ;\r\n",
"\t\tEXT:_FillValue = -1.e+34f ;\r\n",
"\t\tEXT:cell_methods = \"time: mean\" ;\r\n",
"\t\tEXT:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat EXT_max(time, yt, xt) ;\r\n",
"\t\tEXT_max:long_name = \"ice modeled\" ;\r\n",
"\t\tEXT_max:units = \"0 or 1\" ;\r\n",
"\t\tEXT_max:missing_value = -1.e+34f ;\r\n",
"\t\tEXT_max:_FillValue = -1.e+34f ;\r\n",
"\t\tEXT_max:cell_methods = \"time: max\" ;\r\n",
"\t\tEXT_max:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat EXT_min(time, yt, xt) ;\r\n",
"\t\tEXT_min:long_name = \"ice modeled\" ;\r\n",
"\t\tEXT_min:units = \"0 or 1\" ;\r\n",
"\t\tEXT_min:missing_value = -1.e+34f ;\r\n",
"\t\tEXT_min:_FillValue = -1.e+34f ;\r\n",
"\t\tEXT_min:cell_methods = \"time: min\" ;\r\n",
"\t\tEXT_min:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat FA_X(time, yv, xv) ;\r\n",
"\t\tFA_X:long_name = \"air stress on ice - x component\" ;\r\n",
"\t\tFA_X:units = \"Pa\" ;\r\n",
"\t\tFA_X:missing_value = -1.e+34f ;\r\n",
"\t\tFA_X:_FillValue = -1.e+34f ;\r\n",
"\t\tFA_X:cell_methods = \"time: mean\" ;\r\n",
"\t\tFA_X:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat FA_Y(time, yv, xv) ;\r\n",
"\t\tFA_Y:long_name = \"air stress on ice - y component\" ;\r\n",
"\t\tFA_Y:units = \"Pa\" ;\r\n",
"\t\tFA_Y:missing_value = -1.e+34f ;\r\n",
"\t\tFA_Y:_FillValue = -1.e+34f ;\r\n",
"\t\tFA_Y:cell_methods = \"time: mean\" ;\r\n",
"\t\tFA_Y:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat FI_X(time, yv, xv) ;\r\n",
"\t\tFI_X:long_name = \"ice internal stress - x component\" ;\r\n",
"\t\tFI_X:units = \"Pa\" ;\r\n",
"\t\tFI_X:missing_value = -1.e+34f ;\r\n",
"\t\tFI_X:_FillValue = -1.e+34f ;\r\n",
"\t\tFI_X:cell_methods = \"time: mean\" ;\r\n",
"\t\tFI_X:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat FI_Y(time, yv, xv) ;\r\n",
"\t\tFI_Y:long_name = \"ice internal stress - y component\" ;\r\n",
"\t\tFI_Y:units = \"Pa\" ;\r\n",
"\t\tFI_Y:missing_value = -1.e+34f ;\r\n",
"\t\tFI_Y:_FillValue = -1.e+34f ;\r\n",
"\t\tFI_Y:cell_methods = \"time: mean\" ;\r\n",
"\t\tFI_Y:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat FRAZIL(time, yt, xt) ;\r\n",
"\t\tFRAZIL:long_name = \"energy flux of frazil formation\" ;\r\n",
"\t\tFRAZIL:units = \"W/m^2\" ;\r\n",
"\t\tFRAZIL:missing_value = -1.e+34f ;\r\n",
"\t\tFRAZIL:_FillValue = -1.e+34f ;\r\n",
"\t\tFRAZIL:cell_methods = \"time: mean\" ;\r\n",
"\t\tFRAZIL:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat GEOLAT(yt, xt) ;\r\n",
"\t\tGEOLAT:long_name = \"latitude\" ;\r\n",
"\t\tGEOLAT:units = \"degrees\" ;\r\n",
"\t\tGEOLAT:missing_value = 1.e+20f ;\r\n",
"\t\tGEOLAT:_FillValue = 1.e+20f ;\r\n",
"\t\tGEOLAT:cell_methods = \"time: point\" ;\r\n",
"\tfloat GEOLON(yt, xt) ;\r\n",
"\t\tGEOLON:long_name = \"longitude\" ;\r\n",
"\t\tGEOLON:units = \"degrees\" ;\r\n",
"\t\tGEOLON:missing_value = 1.e+20f ;\r\n",
"\t\tGEOLON:_FillValue = 1.e+20f ;\r\n",
"\t\tGEOLON:cell_methods = \"time: point\" ;\r\n",
"\tfloat HI(time, yt, xt) ;\r\n",
"\t\tHI:long_name = \"ice thickness\" ;\r\n",
"\t\tHI:units = \"m-ice\" ;\r\n",
"\t\tHI:missing_value = -1.e+34f ;\r\n",
"\t\tHI:_FillValue = -1.e+34f ;\r\n",
"\t\tHI:cell_methods = \"time: mean\" ;\r\n",
"\t\tHI:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat HS(time, yt, xt) ;\r\n",
"\t\tHS:long_name = \"snow thickness\" ;\r\n",
"\t\tHS:units = \"m-snow\" ;\r\n",
"\t\tHS:missing_value = -1.e+34f ;\r\n",
"\t\tHS:_FillValue = -1.e+34f ;\r\n",
"\t\tHS:cell_methods = \"time: mean\" ;\r\n",
"\t\tHS:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat IX_TRANS(time, yt, xv) ;\r\n",
"\t\tIX_TRANS:long_name = \"x-direction ice transport\" ;\r\n",
"\t\tIX_TRANS:units = \"kg/s\" ;\r\n",
"\t\tIX_TRANS:missing_value = -1.e+34f ;\r\n",
"\t\tIX_TRANS:_FillValue = -1.e+34f ;\r\n",
"\t\tIX_TRANS:cell_methods = \"time: mean\" ;\r\n",
"\t\tIX_TRANS:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat IY_TRANS(time, yv, xt) ;\r\n",
"\t\tIY_TRANS:long_name = \"y-direction ice transport\" ;\r\n",
"\t\tIY_TRANS:units = \"kg/s\" ;\r\n",
"\t\tIY_TRANS:missing_value = -1.e+34f ;\r\n",
"\t\tIY_TRANS:_FillValue = -1.e+34f ;\r\n",
"\t\tIY_TRANS:cell_methods = \"time: mean\" ;\r\n",
"\t\tIY_TRANS:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat LH(time, yt, xt) ;\r\n",
"\t\tLH:long_name = \"latent heat flux\" ;\r\n",
"\t\tLH:units = \"W/m^2\" ;\r\n",
"\t\tLH:missing_value = -1.e+34f ;\r\n",
"\t\tLH:_FillValue = -1.e+34f ;\r\n",
"\t\tLH:cell_methods = \"time: mean\" ;\r\n",
"\t\tLH:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat LSNK(time, yt, xt) ;\r\n",
"\t\tLSNK:long_name = \"frozen water local sink\" ;\r\n",
"\t\tLSNK:units = \"kg/(m^2*yr)\" ;\r\n",
"\t\tLSNK:missing_value = -1.e+34f ;\r\n",
"\t\tLSNK:_FillValue = -1.e+34f ;\r\n",
"\t\tLSNK:cell_methods = \"time: mean\" ;\r\n",
"\t\tLSNK:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat LSRC(time, yt, xt) ;\r\n",
"\t\tLSRC:long_name = \"frozen water local source\" ;\r\n",
"\t\tLSRC:units = \"kg/(m^2*yr)\" ;\r\n",
"\t\tLSRC:missing_value = -1.e+34f ;\r\n",
"\t\tLSRC:_FillValue = -1.e+34f ;\r\n",
"\t\tLSRC:cell_methods = \"time: mean\" ;\r\n",
"\t\tLSRC:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat LW(time, yt, xt) ;\r\n",
"\t\tLW:long_name = \"long wave heat flux over ice\" ;\r\n",
"\t\tLW:units = \"W/m^2\" ;\r\n",
"\t\tLW:missing_value = -1.e+34f ;\r\n",
"\t\tLW:_FillValue = -1.e+34f ;\r\n",
"\t\tLW:cell_methods = \"time: mean\" ;\r\n",
"\t\tLW:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat MI(time, yt, xt) ;\r\n",
"\t\tMI:long_name = \"ice mass\" ;\r\n",
"\t\tMI:units = \"kg/m^2\" ;\r\n",
"\t\tMI:missing_value = -1.e+34f ;\r\n",
"\t\tMI:_FillValue = -1.e+34f ;\r\n",
"\t\tMI:cell_methods = \"time: mean\" ;\r\n",
"\t\tMI:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat RAIN(time, yt, xt) ;\r\n",
"\t\tRAIN:long_name = \"rate of rain fall\" ;\r\n",
"\t\tRAIN:units = \"kg/(m^2*s)\" ;\r\n",
"\t\tRAIN:missing_value = -1.e+34f ;\r\n",
"\t\tRAIN:_FillValue = -1.e+34f ;\r\n",
"\t\tRAIN:cell_methods = \"time: mean\" ;\r\n",
"\t\tRAIN:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat RUNOFF(time, yt, xt) ;\r\n",
"\t\tRUNOFF:long_name = \"liquid runoff\" ;\r\n",
"\t\tRUNOFF:units = \"kg/(m^2*s)\" ;\r\n",
"\t\tRUNOFF:missing_value = -1.e+34f ;\r\n",
"\t\tRUNOFF:_FillValue = -1.e+34f ;\r\n",
"\t\tRUNOFF:cell_methods = \"time: mean\" ;\r\n",
"\t\tRUNOFF:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat SALTF(time, yt, xt) ;\r\n",
"\t\tSALTF:long_name = \"ice to ocean salt flux\" ;\r\n",
"\t\tSALTF:units = \"kg/(m^2*s)\" ;\r\n",
"\t\tSALTF:missing_value = -1.e+34f ;\r\n",
"\t\tSALTF:_FillValue = -1.e+34f ;\r\n",
"\t\tSALTF:cell_methods = \"time: mean\" ;\r\n",
"\t\tSALTF:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat SH(time, yt, xt) ;\r\n",
"\t\tSH:long_name = \"sensible heat flux\" ;\r\n",
"\t\tSH:units = \"W/m^2\" ;\r\n",
"\t\tSH:missing_value = -1.e+34f ;\r\n",
"\t\tSH:_FillValue = -1.e+34f ;\r\n",
"\t\tSH:cell_methods = \"time: mean\" ;\r\n",
"\t\tSH:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat SINROT(yt, xt) ;\r\n",
"\t\tSINROT:long_name = \"-SINROT,COSROT points north\" ;\r\n",
"\t\tSINROT:units = \"none\" ;\r\n",
"\t\tSINROT:missing_value = 1.e+20f ;\r\n",
"\t\tSINROT:_FillValue = 1.e+20f ;\r\n",
"\t\tSINROT:cell_methods = \"time: point\" ;\r\n",
"\tfloat SN2IC(time, yt, xt) ;\r\n",
"\t\tSN2IC:long_name = \"rate of snow to ice conversion\" ;\r\n",
"\t\tSN2IC:units = \"kg/(m^2*s)\" ;\r\n",
"\t\tSN2IC:missing_value = -1.e+34f ;\r\n",
"\t\tSN2IC:_FillValue = -1.e+34f ;\r\n",
"\t\tSN2IC:cell_methods = \"time: mean\" ;\r\n",
"\t\tSN2IC:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat SNOWFL(time, yt, xt) ;\r\n",
"\t\tSNOWFL:long_name = \"rate of snow fall\" ;\r\n",
"\t\tSNOWFL:units = \"kg/(m^2*s)\" ;\r\n",
"\t\tSNOWFL:missing_value = -1.e+34f ;\r\n",
"\t\tSNOWFL:_FillValue = -1.e+34f ;\r\n",
"\t\tSNOWFL:cell_methods = \"time: mean\" ;\r\n",
"\t\tSNOWFL:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat SSH(time, yt, xt) ;\r\n",
"\t\tSSH:long_name = \"sea surface height\" ;\r\n",
"\t\tSSH:units = \"m\" ;\r\n",
"\t\tSSH:missing_value = -1.e+34f ;\r\n",
"\t\tSSH:_FillValue = -1.e+34f ;\r\n",
"\t\tSSH:cell_methods = \"time: mean\" ;\r\n",
"\t\tSSH:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat SSS(time, yt, xt) ;\r\n",
"\t\tSSS:long_name = \"sea surface salinity\" ;\r\n",
"\t\tSSS:units = \"psu\" ;\r\n",
"\t\tSSS:missing_value = -1.e+34f ;\r\n",
"\t\tSSS:_FillValue = -1.e+34f ;\r\n",
"\t\tSSS:cell_methods = \"time: mean\" ;\r\n",
"\t\tSSS:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat SST(time, yt, xt) ;\r\n",
"\t\tSST:long_name = \"sea surface temperature\" ;\r\n",
"\t\tSST:units = \"deg-C\" ;\r\n",
"\t\tSST:missing_value = -1.e+34f ;\r\n",
"\t\tSST:_FillValue = -1.e+34f ;\r\n",
"\t\tSST:cell_methods = \"time: mean\" ;\r\n",
"\t\tSST:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat SW(time, yt, xt) ;\r\n",
"\t\tSW:long_name = \"short wave heat flux\" ;\r\n",
"\t\tSW:units = \"W/m^2\" ;\r\n",
"\t\tSW:missing_value = -1.e+34f ;\r\n",
"\t\tSW:_FillValue = -1.e+34f ;\r\n",
"\t\tSW:cell_methods = \"time: mean\" ;\r\n",
"\t\tSW:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat T1(time, yt, xt) ;\r\n",
"\t\tT1:long_name = \"upper ice layer temperature\" ;\r\n",
"\t\tT1:units = \"C\" ;\r\n",
"\t\tT1:missing_value = -1.e+34f ;\r\n",
"\t\tT1:_FillValue = -1.e+34f ;\r\n",
"\t\tT1:cell_methods = \"time: mean\" ;\r\n",
"\t\tT1:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat T2(time, yt, xt) ;\r\n",
"\t\tT2:long_name = \"lower ice layer temperature\" ;\r\n",
"\t\tT2:units = \"C\" ;\r\n",
"\t\tT2:missing_value = -1.e+34f ;\r\n",
"\t\tT2:_FillValue = -1.e+34f ;\r\n",
"\t\tT2:cell_methods = \"time: mean\" ;\r\n",
"\t\tT2:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat TMELT(time, yt, xt) ;\r\n",
"\t\tTMELT:long_name = \"upper surface melting energy flux\" ;\r\n",
"\t\tTMELT:units = \"W/m^2\" ;\r\n",
"\t\tTMELT:missing_value = -1.e+34f ;\r\n",
"\t\tTMELT:_FillValue = -1.e+34f ;\r\n",
"\t\tTMELT:cell_methods = \"time: mean\" ;\r\n",
"\t\tTMELT:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat TS(time, yt, xt) ;\r\n",
"\t\tTS:long_name = \"surface temperature\" ;\r\n",
"\t\tTS:units = \"C\" ;\r\n",
"\t\tTS:missing_value = -1.e+34f ;\r\n",
"\t\tTS:_FillValue = -1.e+34f ;\r\n",
"\t\tTS:cell_methods = \"time: mean\" ;\r\n",
"\t\tTS:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat UI(time, yv, xv) ;\r\n",
"\t\tUI:long_name = \"ice velocity - x component\" ;\r\n",
"\t\tUI:units = \"m/s\" ;\r\n",
"\t\tUI:missing_value = -1.e+34f ;\r\n",
"\t\tUI:_FillValue = -1.e+34f ;\r\n",
"\t\tUI:cell_methods = \"time: mean\" ;\r\n",
"\t\tUI:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat UO(time, yv, xv) ;\r\n",
"\t\tUO:long_name = \"surface current - x component\" ;\r\n",
"\t\tUO:units = \"m/s\" ;\r\n",
"\t\tUO:missing_value = -1.e+34f ;\r\n",
"\t\tUO:_FillValue = -1.e+34f ;\r\n",
"\t\tUO:cell_methods = \"time: mean\" ;\r\n",
"\t\tUO:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat VI(time, yv, xv) ;\r\n",
"\t\tVI:long_name = \"ice velocity - y component\" ;\r\n",
"\t\tVI:units = \"m/s\" ;\r\n",
"\t\tVI:missing_value = -1.e+34f ;\r\n",
"\t\tVI:_FillValue = -1.e+34f ;\r\n",
"\t\tVI:cell_methods = \"time: mean\" ;\r\n",
"\t\tVI:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat VO(time, yv, xv) ;\r\n",
"\t\tVO:long_name = \"surface current - y component\" ;\r\n",
"\t\tVO:units = \"m/s\" ;\r\n",
"\t\tVO:missing_value = -1.e+34f ;\r\n",
"\t\tVO:_FillValue = -1.e+34f ;\r\n",
"\t\tVO:cell_methods = \"time: mean\" ;\r\n",
"\t\tVO:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tfloat XPRT(time, yt, xt) ;\r\n",
"\t\tXPRT:long_name = \"frozen water transport convergence\" ;\r\n",
"\t\tXPRT:units = \"kg/(m^2*yr)\" ;\r\n",
"\t\tXPRT:missing_value = -1.e+34f ;\r\n",
"\t\tXPRT:_FillValue = -1.e+34f ;\r\n",
"\t\tXPRT:cell_methods = \"time: mean\" ;\r\n",
"\t\tXPRT:time_avg_info = \"average_T1,average_T2,average_DT\" ;\r\n",
"\tdouble average_DT(time) ;\r\n",
"\t\taverage_DT:long_name = \"Length of average period\" ;\r\n",
"\t\taverage_DT:units = \"days\" ;\r\n",
"\t\taverage_DT:missing_value = 1.e+20 ;\r\n",
"\t\taverage_DT:_FillValue = 1.e+20 ;\r\n",
"\tdouble average_T1(time) ;\r\n",
"\t\taverage_T1:long_name = \"Start time for average period\" ;\r\n",
"\t\taverage_T1:units = \"days since 0001-01-01 00:00:00\" ;\r\n",
"\t\taverage_T1:missing_value = 1.e+20 ;\r\n",
"\t\taverage_T1:_FillValue = 1.e+20 ;\r\n",
"\tdouble average_T2(time) ;\r\n",
"\t\taverage_T2:long_name = \"End time for average period\" ;\r\n",
"\t\taverage_T2:units = \"days since 0001-01-01 00:00:00\" ;\r\n",
"\t\taverage_T2:missing_value = 1.e+20 ;\r\n",
"\t\taverage_T2:_FillValue = 1.e+20 ;\r\n",
"\tdouble ct(ct) ;\r\n",
"\t\tct:long_name = \"thickness\" ;\r\n",
"\t\tct:units = \"meters\" ;\r\n",
"\t\tct:cartesian_axis = \"Z\" ;\r\n",
"\tdouble nv(nv) ;\r\n",
"\t\tnv:long_name = \"vertex number\" ;\r\n",
"\t\tnv:units = \"none\" ;\r\n",
"\t\tnv:cartesian_axis = \"N\" ;\r\n",
"\tdouble time(time) ;\r\n",
"\t\ttime:long_name = \"time\" ;\r\n",
"\t\ttime:units = \"days since 0001-01-01 00:00:00\" ;\r\n",
"\t\ttime:cartesian_axis = \"T\" ;\r\n",
"\t\ttime:calendar_type = \"NOLEAP\" ;\r\n",
"\t\ttime:calendar = \"NOLEAP\" ;\r\n",
"\t\ttime:bounds = \"time_bounds\" ;\r\n",
"\tdouble time_bounds(time, nv) ;\r\n",
"\t\ttime_bounds:long_name = \"time axis boundaries\" ;\r\n",
"\t\ttime_bounds:units = \"days\" ;\r\n",
"\t\ttime_bounds:missing_value = 1.e+20 ;\r\n",
"\t\ttime_bounds:_FillValue = 1.e+20 ;\r\n",
"\tdouble xb(xb) ;\r\n",
"\t\txb:long_name = \"longitude\" ;\r\n",
"\t\txb:units = \"degrees_E\" ;\r\n",
"\t\txb:cartesian_axis = \"X\" ;\r\n",
"\tdouble xt(xt) ;\r\n",
"\t\txt:long_name = \"longitude\" ;\r\n",
"\t\txt:units = \"degrees_E\" ;\r\n",
"\t\txt:cartesian_axis = \"X\" ;\r\n",
"\t\txt:edges = \"xb\" ;\r\n",
"\tdouble xv(xv) ;\r\n",
"\t\txv:long_name = \"longitude\" ;\r\n",
"\t\txv:units = \"degrees_E\" ;\r\n",
"\t\txv:cartesian_axis = \"X\" ;\r\n",
"\tdouble yb(yb) ;\r\n",
"\t\tyb:long_name = \"latitude\" ;\r\n",
"\t\tyb:units = \"degrees_N\" ;\r\n",
"\t\tyb:cartesian_axis = \"Y\" ;\r\n",
"\tdouble yt(yt) ;\r\n",
"\t\tyt:long_name = \"latitude\" ;\r\n",
"\t\tyt:units = \"degrees_N\" ;\r\n",
"\t\tyt:cartesian_axis = \"Y\" ;\r\n",
"\t\tyt:edges = \"yb\" ;\r\n",
"\tdouble yv(yv) ;\r\n",
"\t\tyv:long_name = \"latitude\" ;\r\n",
"\t\tyv:units = \"degrees_N\" ;\r\n",
"\t\tyv:cartesian_axis = \"Y\" ;\r\n",
"\r\n",
"// global attributes:\r\n",
"\t\t:filename = \"ice_month.nc\" ;\r\n",
"\t\t:title = \"CM2.5_A_Control-1990_FLOR_B01\" ;\r\n",
"\t\t:grid_type = \"regular\" ;\r\n",
"\t\t:grid_tile = \"N/A\" ;\r\n",
"\t\t:history = \"Sun Jan 12 20:12:14 2020: ncks -7 -L 1 POSTP/00010101.ice_month.nc POSTP_deflate1/00010101.ice_month.nc\" ;\r\n",
"\t\t:NCO = \"netCDF Operators version 4.8.1 (Homepage = http://nco.sf.net, Code = http://github.com/nco/nco)\" ;\r\n",
"}\r\n"
]
}
],
"source": [
"# \n",
"expname = 'CTL1860_noleap_tigercpu_intelmpi_18_576PE'\n",
"ifile = f'/tigress/wenchang/MODEL_OUT/{expname}/POSTP/00010101.ice_month.nc'\n",
"print(ifile)\n",
"!ncdump -h $ifile"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"### 00010101.atmos_month.nc"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {
"ExecuteTime": {
"end_time": "2020-02-25T18:11:53.702776Z",
"start_time": "2020-02-25T18:11:52.948129Z"
},
"scrolled": false
},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"/tigress/wenchang/MODEL_OUT/CTL1860_noleap_tigercpu_intelmpi_18_576PE/POSTP/00010101.atmos_month.nc\n",
"netcdf \\00010101.atmos_month {\n",
"dimensions:\n",
"\ttime = UNLIMITED ; // (12 currently)\n",
"\tgrid_yt = 360 ;\n",
"\tgrid_xt = 576 ;\n",
"\tpfull = 32 ;\n",
"\tbnds = 2 ;\n",
"\tphalf = 33 ;\n",
"variables:\n",
"\tfloat IWP(time, grid_yt, grid_xt) ;\n",
"\t\tIWP:long_name = \"Ice water path\" ;\n",
"\t\tIWP:units = \"kg/m2\" ;\n",
"\t\tIWP:missing_value = 1.e+20f ;\n",
"\t\tIWP:_FillValue = 1.e+20f ;\n",
"\t\tIWP:cell_methods = \"time: mean\" ;\n",
"\t\tIWP:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tIWP:interp_method = \"conserve_order1\" ;\n",
"\tfloat LWP(time, grid_yt, grid_xt) ;\n",
"\t\tLWP:long_name = \"Liquid water path\" ;\n",
"\t\tLWP:units = \"kg/m2\" ;\n",
"\t\tLWP:missing_value = 1.e+20f ;\n",
"\t\tLWP:_FillValue = 1.e+20f ;\n",
"\t\tLWP:cell_methods = \"time: mean\" ;\n",
"\t\tLWP:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tLWP:interp_method = \"conserve_order1\" ;\n",
"\tfloat WVP(time, grid_yt, grid_xt) ;\n",
"\t\tWVP:long_name = \"Column integrated water vapor\" ;\n",
"\t\tWVP:units = \"kg/m2\" ;\n",
"\t\tWVP:missing_value = 1.e+20f ;\n",
"\t\tWVP:_FillValue = 1.e+20f ;\n",
"\t\tWVP:cell_methods = \"time: mean\" ;\n",
"\t\tWVP:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tWVP:interp_method = \"conserve_order1\" ;\n",
"\tfloat alb_sfc(time, grid_yt, grid_xt) ;\n",
"\t\talb_sfc:long_name = \"surface albedo\" ;\n",
"\t\talb_sfc:units = \"percent\" ;\n",
"\t\talb_sfc:missing_value = -999.f ;\n",
"\t\talb_sfc:_FillValue = -999.f ;\n",
"\t\talb_sfc:cell_methods = \"time: mean\" ;\n",
"\t\talb_sfc:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\talb_sfc:interp_method = \"conserve_order1\" ;\n",
"\tfloat area(grid_yt, grid_xt) ;\n",
"\t\tarea:long_name = \"cell area\" ;\n",
"\t\tarea:units = \"m**2\" ;\n",
"\t\tarea:missing_value = 1.e+20f ;\n",
"\t\tarea:_FillValue = 1.e+20f ;\n",
"\t\tarea:cell_methods = \"time: point\" ;\n",
"\t\tarea:interp_method = \"conserve_order1\" ;\n",
"\tdouble average_DT(time) ;\n",
"\t\taverage_DT:long_name = \"Length of average period\" ;\n",
"\t\taverage_DT:units = \"days\" ;\n",
"\t\taverage_DT:missing_value = 1.e+20 ;\n",
"\t\taverage_DT:_FillValue = 1.e+20 ;\n",
"\tdouble average_T1(time) ;\n",
"\t\taverage_T1:long_name = \"Start time for average period\" ;\n",
"\t\taverage_T1:units = \"days since 0001-01-01 00:00:00\" ;\n",
"\t\taverage_T1:missing_value = 1.e+20 ;\n",
"\t\taverage_T1:_FillValue = 1.e+20 ;\n",
"\tdouble average_T2(time) ;\n",
"\t\taverage_T2:long_name = \"End time for average period\" ;\n",
"\t\taverage_T2:units = \"days since 0001-01-01 00:00:00\" ;\n",
"\t\taverage_T2:missing_value = 1.e+20 ;\n",
"\t\taverage_T2:_FillValue = 1.e+20 ;\n",
"\tfloat cat15(time, grid_yt, grid_xt) ;\n",
"\t\tcat15:long_name = \"de-pression < 1000\" ;\n",
"\t\tcat15:units = \"mb\" ;\n",
"\t\tcat15:missing_value = -1.e+10f ;\n",
"\t\tcat15:_FillValue = -1.e+10f ;\n",
"\t\tcat15:cell_methods = \"time: mean\" ;\n",
"\t\tcat15:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tcat15:interp_method = \"conserve_order1\" ;\n",
"\tfloat cat15_max(time, grid_yt, grid_xt) ;\n",
"\t\tcat15_max:long_name = \"de-pression < 1000\" ;\n",
"\t\tcat15_max:units = \"mb\" ;\n",
"\t\tcat15_max:missing_value = -1.e+10f ;\n",
"\t\tcat15_max:_FillValue = -1.e+10f ;\n",
"\t\tcat15_max:cell_methods = \"time: max\" ;\n",
"\t\tcat15_max:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tcat15_max:interp_method = \"conserve_order1\" ;\n",
"\tfloat cld_amt(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tcld_amt:long_name = \"cloud fraction\" ;\n",
"\t\tcld_amt:units = \"none\" ;\n",
"\t\tcld_amt:missing_value = -1.e+10f ;\n",
"\t\tcld_amt:_FillValue = -1.e+10f ;\n",
"\t\tcld_amt:cell_methods = \"time: mean\" ;\n",
"\t\tcld_amt:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tcld_amt:interp_method = \"conserve_order1\" ;\n",
"\tfloat drag_heat(time, grid_yt, grid_xt) ;\n",
"\t\tdrag_heat:long_name = \"drag coeff for heat\" ;\n",
"\t\tdrag_heat:units = \"none\" ;\n",
"\t\tdrag_heat:missing_value = 1.e+20f ;\n",
"\t\tdrag_heat:_FillValue = 1.e+20f ;\n",
"\t\tdrag_heat:cell_methods = \"time: mean\" ;\n",
"\t\tdrag_heat:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tdrag_heat:interp_method = \"conserve_order1\" ;\n",
"\tfloat drag_moist(time, grid_yt, grid_xt) ;\n",
"\t\tdrag_moist:long_name = \"drag coeff for moisture\" ;\n",
"\t\tdrag_moist:units = \"none\" ;\n",
"\t\tdrag_moist:missing_value = 1.e+20f ;\n",
"\t\tdrag_moist:_FillValue = 1.e+20f ;\n",
"\t\tdrag_moist:cell_methods = \"time: mean\" ;\n",
"\t\tdrag_moist:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tdrag_moist:interp_method = \"conserve_order1\" ;\n",
"\tfloat drag_mom(time, grid_yt, grid_xt) ;\n",
"\t\tdrag_mom:long_name = \"drag coeff for momentum\" ;\n",
"\t\tdrag_mom:units = \"none\" ;\n",
"\t\tdrag_mom:missing_value = 1.e+20f ;\n",
"\t\tdrag_mom:_FillValue = 1.e+20f ;\n",
"\t\tdrag_mom:cell_methods = \"time: mean\" ;\n",
"\t\tdrag_mom:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tdrag_mom:interp_method = \"conserve_order1\" ;\n",
"\tfloat evap(time, grid_yt, grid_xt) ;\n",
"\t\tevap:long_name = \"evaporation rate\" ;\n",
"\t\tevap:units = \"kg/m2/s\" ;\n",
"\t\tevap:missing_value = 1.e+20f ;\n",
"\t\tevap:_FillValue = 1.e+20f ;\n",
"\t\tevap:cell_methods = \"time: mean\" ;\n",
"\t\tevap:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tevap:interp_method = \"conserve_order1\" ;\n",
"\tfloat f15(time, grid_yt, grid_xt) ;\n",
"\t\tf15:long_name = \"Cat15 frequency\" ;\n",
"\t\tf15:units = \"none\" ;\n",
"\t\tf15:missing_value = -1.e+10f ;\n",
"\t\tf15:_FillValue = -1.e+10f ;\n",
"\t\tf15:cell_methods = \"time: mean\" ;\n",
"\t\tf15:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tf15:interp_method = \"conserve_order1\" ;\n",
"\tdouble grid_xt(grid_xt) ;\n",
"\t\tgrid_xt:long_name = \"T-cell longitude\" ;\n",
"\t\tgrid_xt:units = \"degrees_E\" ;\n",
"\t\tgrid_xt:cartesian_axis = \"X\" ;\n",
"\t\tgrid_xt:bounds = \"grid_xt_bnds\" ;\n",
"\tdouble grid_xt_bnds(grid_xt, bnds) ;\n",
"\t\tgrid_xt_bnds:long_name = \"T-cell longitude\" ;\n",
"\t\tgrid_xt_bnds:units = \"degrees_E\" ;\n",
"\t\tgrid_xt_bnds:cartesian_axis = \"X\" ;\n",
"\tdouble grid_yt(grid_yt) ;\n",
"\t\tgrid_yt:long_name = \"T-cell latitude\" ;\n",
"\t\tgrid_yt:units = \"degrees_N\" ;\n",
"\t\tgrid_yt:cartesian_axis = \"Y\" ;\n",
"\t\tgrid_yt:bounds = \"grid_yt_bnds\" ;\n",
"\tdouble grid_yt_bnds(grid_yt, bnds) ;\n",
"\t\tgrid_yt_bnds:long_name = \"T-cell latitude\" ;\n",
"\t\tgrid_yt_bnds:units = \"degrees_N\" ;\n",
"\t\tgrid_yt_bnds:cartesian_axis = \"Y\" ;\n",
"\tfloat high_cld_amt(time, grid_yt, grid_xt) ;\n",
"\t\thigh_cld_amt:long_name = \"high cloud amount\" ;\n",
"\t\thigh_cld_amt:units = \"percent\" ;\n",
"\t\thigh_cld_amt:missing_value = 1.e+20f ;\n",
"\t\thigh_cld_amt:_FillValue = 1.e+20f ;\n",
"\t\thigh_cld_amt:cell_methods = \"time: mean\" ;\n",
"\t\thigh_cld_amt:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\thigh_cld_amt:interp_method = \"conserve_order1\" ;\n",
"\tfloat ice_mask(time, grid_yt, grid_xt) ;\n",
"\t\tice_mask:long_name = \"fractional amount of sea ice\" ;\n",
"\t\tice_mask:units = \"none\" ;\n",
"\t\tice_mask:valid_range = -0.01f, 1.01f ;\n",
"\t\tice_mask:missing_value = 1.e+20f ;\n",
"\t\tice_mask:_FillValue = 1.e+20f ;\n",
"\t\tice_mask:cell_methods = \"time: mean\" ;\n",
"\t\tice_mask:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tice_mask:interp_method = \"conserve_order1\" ;\n",
"\tfloat ice_wat(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tice_wat:long_name = \"cloud ice water specific humidity\" ;\n",
"\t\tice_wat:units = \"kg/kg\" ;\n",
"\t\tice_wat:missing_value = -1.e+10f ;\n",
"\t\tice_wat:_FillValue = -1.e+10f ;\n",
"\t\tice_wat:cell_methods = \"time: mean\" ;\n",
"\t\tice_wat:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tice_wat:interp_method = \"conserve_order1\" ;\n",
"\tfloat land_mask(grid_yt, grid_xt) ;\n",
"\t\tland_mask:long_name = \"fractional amount of land\" ;\n",
"\t\tland_mask:units = \"none\" ;\n",
"\t\tland_mask:valid_range = -0.01f, 1.01f ;\n",
"\t\tland_mask:missing_value = 1.e+20f ;\n",
"\t\tland_mask:_FillValue = 1.e+20f ;\n",
"\t\tland_mask:cell_methods = \"time: point\" ;\n",
"\t\tland_mask:interp_method = \"conserve_order1\" ;\n",
"\tfloat liq_wat(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tliq_wat:long_name = \"cloud liquid specific humidity\" ;\n",
"\t\tliq_wat:units = \"kg/kg\" ;\n",
"\t\tliq_wat:missing_value = -1.e+10f ;\n",
"\t\tliq_wat:_FillValue = -1.e+10f ;\n",
"\t\tliq_wat:cell_methods = \"time: mean\" ;\n",
"\t\tliq_wat:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tliq_wat:interp_method = \"conserve_order1\" ;\n",
"\tfloat low_cld_amt(time, grid_yt, grid_xt) ;\n",
"\t\tlow_cld_amt:long_name = \"low cloud amount\" ;\n",
"\t\tlow_cld_amt:units = \"percent\" ;\n",
"\t\tlow_cld_amt:missing_value = 1.e+20f ;\n",
"\t\tlow_cld_amt:_FillValue = 1.e+20f ;\n",
"\t\tlow_cld_amt:cell_methods = \"time: mean\" ;\n",
"\t\tlow_cld_amt:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tlow_cld_amt:interp_method = \"conserve_order1\" ;\n",
"\tfloat lwdn_sfc(time, grid_yt, grid_xt) ;\n",
"\t\tlwdn_sfc:long_name = \"LW flux down at surface\" ;\n",
"\t\tlwdn_sfc:units = \"watts/m2\" ;\n",
"\t\tlwdn_sfc:missing_value = -999.f ;\n",
"\t\tlwdn_sfc:_FillValue = -999.f ;\n",
"\t\tlwdn_sfc:cell_methods = \"time: mean\" ;\n",
"\t\tlwdn_sfc:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tlwdn_sfc:interp_method = \"conserve_order1\" ;\n",
"\tfloat lwdn_sfc_clr(time, grid_yt, grid_xt) ;\n",
"\t\tlwdn_sfc_clr:long_name = \"clear skyLW flux down at surface\" ;\n",
"\t\tlwdn_sfc_clr:units = \"watts/m2\" ;\n",
"\t\tlwdn_sfc_clr:missing_value = -999.f ;\n",
"\t\tlwdn_sfc_clr:_FillValue = -999.f ;\n",
"\t\tlwdn_sfc_clr:cell_methods = \"time: mean\" ;\n",
"\t\tlwdn_sfc_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tlwdn_sfc_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat lwflx(time, grid_yt, grid_xt) ;\n",
"\t\tlwflx:long_name = \"net (down-up) longwave flux\" ;\n",
"\t\tlwflx:units = \"w/m2\" ;\n",
"\t\tlwflx:missing_value = 1.e+20f ;\n",
"\t\tlwflx:_FillValue = 1.e+20f ;\n",
"\t\tlwflx:cell_methods = \"time: mean\" ;\n",
"\t\tlwflx:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tlwflx:interp_method = \"conserve_order1\" ;\n",
"\tfloat lwup_sfc(time, grid_yt, grid_xt) ;\n",
"\t\tlwup_sfc:long_name = \"LW flux up at surface\" ;\n",
"\t\tlwup_sfc:units = \"watts/m2\" ;\n",
"\t\tlwup_sfc:missing_value = -999.f ;\n",
"\t\tlwup_sfc:_FillValue = -999.f ;\n",
"\t\tlwup_sfc:cell_methods = \"time: mean\" ;\n",
"\t\tlwup_sfc:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tlwup_sfc:interp_method = \"conserve_order1\" ;\n",
"\tfloat lwup_sfc_clr(time, grid_yt, grid_xt) ;\n",
"\t\tlwup_sfc_clr:long_name = \"clear skyLW flux up at surface\" ;\n",
"\t\tlwup_sfc_clr:units = \"watts/m2\" ;\n",
"\t\tlwup_sfc_clr:missing_value = -999.f ;\n",
"\t\tlwup_sfc_clr:_FillValue = -999.f ;\n",
"\t\tlwup_sfc_clr:cell_methods = \"time: mean\" ;\n",
"\t\tlwup_sfc_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tlwup_sfc_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat mc(time, phalf, grid_yt, grid_xt) ;\n",
"\t\tmc:long_name = \"Cumulus Mass Flux from RAS\" ;\n",
"\t\tmc:units = \"kg/m2/s\" ;\n",
"\t\tmc:missing_value = -999.f ;\n",
"\t\tmc:_FillValue = -999.f ;\n",
"\t\tmc:cell_methods = \"time: mean\" ;\n",
"\t\tmc:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tmc:interp_method = \"conserve_order1\" ;\n",
"\tfloat mid_cld_amt(time, grid_yt, grid_xt) ;\n",
"\t\tmid_cld_amt:long_name = \"mid cloud amount\" ;\n",
"\t\tmid_cld_amt:units = \"percent\" ;\n",
"\t\tmid_cld_amt:missing_value = 1.e+20f ;\n",
"\t\tmid_cld_amt:_FillValue = 1.e+20f ;\n",
"\t\tmid_cld_amt:cell_methods = \"time: mean\" ;\n",
"\t\tmid_cld_amt:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tmid_cld_amt:interp_method = \"conserve_order1\" ;\n",
"\tfloat netlw_200hPa(time, grid_yt, grid_xt) ;\n",
"\t\tnetlw_200hPa:long_name = \"net LW flux at 200 hPa\" ;\n",
"\t\tnetlw_200hPa:units = \"watts/m2\" ;\n",
"\t\tnetlw_200hPa:missing_value = -999.f ;\n",
"\t\tnetlw_200hPa:_FillValue = -999.f ;\n",
"\t\tnetlw_200hPa:cell_methods = \"time: mean\" ;\n",
"\t\tnetlw_200hPa:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tnetlw_200hPa:interp_method = \"conserve_order1\" ;\n",
"\tfloat netlw_200hPa_clr(time, grid_yt, grid_xt) ;\n",
"\t\tnetlw_200hPa_clr:long_name = \"clear skynet LW flux at 200 hPa\" ;\n",
"\t\tnetlw_200hPa_clr:units = \"watts/m2\" ;\n",
"\t\tnetlw_200hPa_clr:missing_value = -999.f ;\n",
"\t\tnetlw_200hPa_clr:_FillValue = -999.f ;\n",
"\t\tnetlw_200hPa_clr:cell_methods = \"time: mean\" ;\n",
"\t\tnetlw_200hPa_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tnetlw_200hPa_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat netlw_lin_trop(time, grid_yt, grid_xt) ;\n",
"\t\tnetlw_lin_trop:long_name = \"net LW flux at linear tropopause\" ;\n",
"\t\tnetlw_lin_trop:units = \"watts/m2\" ;\n",
"\t\tnetlw_lin_trop:missing_value = -999.f ;\n",
"\t\tnetlw_lin_trop:_FillValue = -999.f ;\n",
"\t\tnetlw_lin_trop:cell_methods = \"time: mean\" ;\n",
"\t\tnetlw_lin_trop:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tnetlw_lin_trop:interp_method = \"conserve_order1\" ;\n",
"\tfloat netlw_lin_trop_clr(time, grid_yt, grid_xt) ;\n",
"\t\tnetlw_lin_trop_clr:long_name = \"clear skynet LW flux at linear tropopause\" ;\n",
"\t\tnetlw_lin_trop_clr:units = \"watts/m2\" ;\n",
"\t\tnetlw_lin_trop_clr:missing_value = -999.f ;\n",
"\t\tnetlw_lin_trop_clr:_FillValue = -999.f ;\n",
"\t\tnetlw_lin_trop_clr:cell_methods = \"time: mean\" ;\n",
"\t\tnetlw_lin_trop_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tnetlw_lin_trop_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat netrad_toa(time, grid_yt, grid_xt) ;\n",
"\t\tnetrad_toa:long_name = \"net radiation (lw + sw) at toa\" ;\n",
"\t\tnetrad_toa:units = \"watts/m2\" ;\n",
"\t\tnetrad_toa:missing_value = -999.f ;\n",
"\t\tnetrad_toa:_FillValue = -999.f ;\n",
"\t\tnetrad_toa:cell_methods = \"time: mean\" ;\n",
"\t\tnetrad_toa:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tnetrad_toa:interp_method = \"conserve_order1\" ;\n",
"\tfloat netrad_toa_clr(time, grid_yt, grid_xt) ;\n",
"\t\tnetrad_toa_clr:long_name = \"clear skynet radiation (lw + sw) at toa\" ;\n",
"\t\tnetrad_toa_clr:units = \"watts/m2\" ;\n",
"\t\tnetrad_toa_clr:missing_value = -999.f ;\n",
"\t\tnetrad_toa_clr:_FillValue = -999.f ;\n",
"\t\tnetrad_toa_clr:cell_methods = \"time: mean\" ;\n",
"\t\tnetrad_toa_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tnetrad_toa_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat olr(time, grid_yt, grid_xt) ;\n",
"\t\tolr:long_name = \"outgoing longwave radiation\" ;\n",
"\t\tolr:units = \"watts/m2\" ;\n",
"\t\tolr:missing_value = -999.f ;\n",
"\t\tolr:_FillValue = -999.f ;\n",
"\t\tolr:cell_methods = \"time: mean\" ;\n",
"\t\tolr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tolr:interp_method = \"conserve_order1\" ;\n",
"\tfloat olr_clr(time, grid_yt, grid_xt) ;\n",
"\t\tolr_clr:long_name = \"clear skyoutgoing longwave radiation\" ;\n",
"\t\tolr_clr:units = \"watts/m2\" ;\n",
"\t\tolr_clr:missing_value = -999.f ;\n",
"\t\tolr_clr:_FillValue = -999.f ;\n",
"\t\tolr_clr:cell_methods = \"time: mean\" ;\n",
"\t\tolr_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tolr_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat omega(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tomega:long_name = \"omega\" ;\n",
"\t\tomega:units = \"Pa/s\" ;\n",
"\t\tomega:missing_value = -1.e+10f ;\n",
"\t\tomega:_FillValue = -1.e+10f ;\n",
"\t\tomega:cell_methods = \"time: mean\" ;\n",
"\t\tomega:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tomega:interp_method = \"conserve_order1\" ;\n",
"\tdouble pfull(pfull) ;\n",
"\t\tpfull:long_name = \"ref full pressure level\" ;\n",
"\t\tpfull:units = \"mb\" ;\n",
"\t\tpfull:cartesian_axis = \"Z\" ;\n",
"\t\tpfull:positive = \"down\" ;\n",
"\tdouble phalf(phalf) ;\n",
"\t\tphalf:long_name = \"ref half pressure level\" ;\n",
"\t\tphalf:units = \"mb\" ;\n",
"\t\tphalf:cartesian_axis = \"Z\" ;\n",
"\t\tphalf:positive = \"down\" ;\n",
"\tfloat prec_conv(time, grid_yt, grid_xt) ;\n",
"\t\tprec_conv:long_name = \"Precipitation rate from convection\" ;\n",
"\t\tprec_conv:units = \"kg(h2o)/m2/s\" ;\n",
"\t\tprec_conv:missing_value = 1.e+20f ;\n",
"\t\tprec_conv:_FillValue = 1.e+20f ;\n",
"\t\tprec_conv:cell_methods = \"time: mean\" ;\n",
"\t\tprec_conv:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tprec_conv:interp_method = \"conserve_order1\" ;\n",
"\tfloat prec_ls(time, grid_yt, grid_xt) ;\n",
"\t\tprec_ls:long_name = \"Precipitation rate from strat cloud\" ;\n",
"\t\tprec_ls:units = \"kg/m2/s\" ;\n",
"\t\tprec_ls:missing_value = 1.e+20f ;\n",
"\t\tprec_ls:_FillValue = 1.e+20f ;\n",
"\t\tprec_ls:cell_methods = \"time: mean\" ;\n",
"\t\tprec_ls:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tprec_ls:interp_method = \"conserve_order1\" ;\n",
"\tfloat precip(time, grid_yt, grid_xt) ;\n",
"\t\tprecip:long_name = \"Total precipitation rate\" ;\n",
"\t\tprecip:units = \"kg/m2/s\" ;\n",
"\t\tprecip:missing_value = 1.e+20f ;\n",
"\t\tprecip:_FillValue = 1.e+20f ;\n",
"\t\tprecip:cell_methods = \"time: mean\" ;\n",
"\t\tprecip:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tprecip:interp_method = \"conserve_order1\" ;\n",
"\tfloat ps(time, grid_yt, grid_xt) ;\n",
"\t\tps:long_name = \"surface pressure\" ;\n",
"\t\tps:units = \"Pa\" ;\n",
"\t\tps:missing_value = -1.e+10f ;\n",
"\t\tps:_FillValue = -1.e+10f ;\n",
"\t\tps:cell_methods = \"time: mean\" ;\n",
"\t\tps:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tps:interp_method = \"conserve_order1\" ;\n",
"\tfloat qdt_conv(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tqdt_conv:long_name = \"Spec humidity tendency from convection\" ;\n",
"\t\tqdt_conv:units = \"kg/kg/s\" ;\n",
"\t\tqdt_conv:missing_value = -999.f ;\n",
"\t\tqdt_conv:_FillValue = -999.f ;\n",
"\t\tqdt_conv:cell_methods = \"time: mean\" ;\n",
"\t\tqdt_conv:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tqdt_conv:interp_method = \"conserve_order1\" ;\n",
"\tfloat qdt_ls(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tqdt_ls:long_name = \"Spec humidity tendency from strat cloud\" ;\n",
"\t\tqdt_ls:units = \"kg/kg/s\" ;\n",
"\t\tqdt_ls:missing_value = -999.f ;\n",
"\t\tqdt_ls:_FillValue = -999.f ;\n",
"\t\tqdt_ls:cell_methods = \"time: mean\" ;\n",
"\t\tqdt_ls:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tqdt_ls:interp_method = \"conserve_order1\" ;\n",
"\tfloat qo3(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tqo3:long_name = \"ozone mixing ratio\" ;\n",
"\t\tqo3:units = \"kg/kg\" ;\n",
"\t\tqo3:missing_value = -999.f ;\n",
"\t\tqo3:_FillValue = -999.f ;\n",
"\t\tqo3:cell_methods = \"time: mean\" ;\n",
"\t\tqo3:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tqo3:interp_method = \"conserve_order1\" ;\n",
"\tfloat qo3_col(time, grid_yt, grid_xt) ;\n",
"\t\tqo3_col:long_name = \"ozone column\" ;\n",
"\t\tqo3_col:units = \"DU\" ;\n",
"\t\tqo3_col:missing_value = -999.f ;\n",
"\t\tqo3_col:_FillValue = -999.f ;\n",
"\t\tqo3_col:cell_methods = \"time: mean\" ;\n",
"\t\tqo3_col:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tqo3_col:interp_method = \"conserve_order1\" ;\n",
"\tfloat rh(time, pfull, grid_yt, grid_xt) ;\n",
"\t\trh:long_name = \"relative humidity\" ;\n",
"\t\trh:units = \"percent\" ;\n",
"\t\trh:missing_value = -999.f ;\n",
"\t\trh:_FillValue = -999.f ;\n",
"\t\trh:cell_methods = \"time: mean\" ;\n",
"\t\trh:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\trh:interp_method = \"conserve_order1\" ;\n",
"\tfloat rh_ref(time, grid_yt, grid_xt) ;\n",
"\t\trh_ref:long_name = \"relative humidity at 2 m\" ;\n",
"\t\trh_ref:units = \"percent\" ;\n",
"\t\trh_ref:missing_value = 1.e+20f ;\n",
"\t\trh_ref:_FillValue = 1.e+20f ;\n",
"\t\trh_ref:cell_methods = \"time: mean\" ;\n",
"\t\trh_ref:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\trh_ref:interp_method = \"conserve_order1\" ;\n",
"\tfloat shflx(time, grid_yt, grid_xt) ;\n",
"\t\tshflx:long_name = \"sensible heat flux\" ;\n",
"\t\tshflx:units = \"w/m2\" ;\n",
"\t\tshflx:missing_value = 1.e+20f ;\n",
"\t\tshflx:_FillValue = 1.e+20f ;\n",
"\t\tshflx:cell_methods = \"time: mean\" ;\n",
"\t\tshflx:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tshflx:interp_method = \"conserve_order1\" ;\n",
"\tfloat slp(time, grid_yt, grid_xt) ;\n",
"\t\tslp:long_name = \"sea-level pressure\" ;\n",
"\t\tslp:units = \"mb\" ;\n",
"\t\tslp:valid_range = 800.f, 1200.f ;\n",
"\t\tslp:missing_value = -1.e+10f ;\n",
"\t\tslp:_FillValue = -1.e+10f ;\n",
"\t\tslp:cell_methods = \"time: mean\" ;\n",
"\t\tslp:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tslp:interp_method = \"conserve_order1\" ;\n",
"\tfloat snow_conv(time, grid_yt, grid_xt) ;\n",
"\t\tsnow_conv:long_name = \"Frozen precip rate from convection\" ;\n",
"\t\tsnow_conv:units = \"kg(h2o)/m2/s\" ;\n",
"\t\tsnow_conv:missing_value = 1.e+20f ;\n",
"\t\tsnow_conv:_FillValue = 1.e+20f ;\n",
"\t\tsnow_conv:cell_methods = \"time: mean\" ;\n",
"\t\tsnow_conv:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tsnow_conv:interp_method = \"conserve_order1\" ;\n",
"\tfloat snow_ls(time, grid_yt, grid_xt) ;\n",
"\t\tsnow_ls:long_name = \"Frozen precip rate from strat cloud\" ;\n",
"\t\tsnow_ls:units = \"kg/m2/s\" ;\n",
"\t\tsnow_ls:missing_value = 1.e+20f ;\n",
"\t\tsnow_ls:_FillValue = 1.e+20f ;\n",
"\t\tsnow_ls:cell_methods = \"time: mean\" ;\n",
"\t\tsnow_ls:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tsnow_ls:interp_method = \"conserve_order1\" ;\n",
"\tfloat snow_tot(time, grid_yt, grid_xt) ;\n",
"\t\tsnow_tot:long_name = \"Frozen precip rate from all sources\" ;\n",
"\t\tsnow_tot:units = \"kg(h2o)/m2/s\" ;\n",
"\t\tsnow_tot:missing_value = 1.e+20f ;\n",
"\t\tsnow_tot:_FillValue = 1.e+20f ;\n",
"\t\tsnow_tot:cell_methods = \"time: mean\" ;\n",
"\t\tsnow_tot:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tsnow_tot:interp_method = \"conserve_order1\" ;\n",
"\tfloat sphum(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tsphum:long_name = \"specific humidity\" ;\n",
"\t\tsphum:units = \"kg/kg\" ;\n",
"\t\tsphum:missing_value = -1.e+10f ;\n",
"\t\tsphum:_FillValue = -1.e+10f ;\n",
"\t\tsphum:cell_methods = \"time: mean\" ;\n",
"\t\tsphum:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tsphum:interp_method = \"conserve_order1\" ;\n",
"\tfloat swdn_200hPa(time, grid_yt, grid_xt) ;\n",
"\t\tswdn_200hPa:long_name = \"SW flux down at 200 hPa\" ;\n",
"\t\tswdn_200hPa:units = \"watts/m2\" ;\n",
"\t\tswdn_200hPa:missing_value = -999.f ;\n",
"\t\tswdn_200hPa:_FillValue = -999.f ;\n",
"\t\tswdn_200hPa:cell_methods = \"time: mean\" ;\n",
"\t\tswdn_200hPa:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswdn_200hPa:interp_method = \"conserve_order1\" ;\n",
"\tfloat swdn_200hPa_clr(time, grid_yt, grid_xt) ;\n",
"\t\tswdn_200hPa_clr:long_name = \"clear skySW flux down at 200 hPa\" ;\n",
"\t\tswdn_200hPa_clr:units = \"watts/m2\" ;\n",
"\t\tswdn_200hPa_clr:missing_value = -999.f ;\n",
"\t\tswdn_200hPa_clr:_FillValue = -999.f ;\n",
"\t\tswdn_200hPa_clr:cell_methods = \"time: mean\" ;\n",
"\t\tswdn_200hPa_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswdn_200hPa_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat swdn_lin_trop(time, grid_yt, grid_xt) ;\n",
"\t\tswdn_lin_trop:long_name = \"SW flux down at linear tropopause\" ;\n",
"\t\tswdn_lin_trop:units = \"watts/m2\" ;\n",
"\t\tswdn_lin_trop:missing_value = -999.f ;\n",
"\t\tswdn_lin_trop:_FillValue = -999.f ;\n",
"\t\tswdn_lin_trop:cell_methods = \"time: mean\" ;\n",
"\t\tswdn_lin_trop:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswdn_lin_trop:interp_method = \"conserve_order1\" ;\n",
"\tfloat swdn_lin_trop_clr(time, grid_yt, grid_xt) ;\n",
"\t\tswdn_lin_trop_clr:long_name = \"clear skySW flux down at linear tropopause\" ;\n",
"\t\tswdn_lin_trop_clr:units = \"watts/m2\" ;\n",
"\t\tswdn_lin_trop_clr:missing_value = -999.f ;\n",
"\t\tswdn_lin_trop_clr:_FillValue = -999.f ;\n",
"\t\tswdn_lin_trop_clr:cell_methods = \"time: mean\" ;\n",
"\t\tswdn_lin_trop_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswdn_lin_trop_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat swdn_sfc(time, grid_yt, grid_xt) ;\n",
"\t\tswdn_sfc:long_name = \"SW flux down at surface\" ;\n",
"\t\tswdn_sfc:units = \"watts/m2\" ;\n",
"\t\tswdn_sfc:missing_value = -999.f ;\n",
"\t\tswdn_sfc:_FillValue = -999.f ;\n",
"\t\tswdn_sfc:cell_methods = \"time: mean\" ;\n",
"\t\tswdn_sfc:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswdn_sfc:interp_method = \"conserve_order1\" ;\n",
"\tfloat swdn_sfc_clr(time, grid_yt, grid_xt) ;\n",
"\t\tswdn_sfc_clr:long_name = \"clear skySW flux down at surface\" ;\n",
"\t\tswdn_sfc_clr:units = \"watts/m2\" ;\n",
"\t\tswdn_sfc_clr:missing_value = -999.f ;\n",
"\t\tswdn_sfc_clr:_FillValue = -999.f ;\n",
"\t\tswdn_sfc_clr:cell_methods = \"time: mean\" ;\n",
"\t\tswdn_sfc_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswdn_sfc_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat swdn_toa(time, grid_yt, grid_xt) ;\n",
"\t\tswdn_toa:long_name = \"SW flux down at TOA\" ;\n",
"\t\tswdn_toa:units = \"watts/m2\" ;\n",
"\t\tswdn_toa:missing_value = -999.f ;\n",
"\t\tswdn_toa:_FillValue = -999.f ;\n",
"\t\tswdn_toa:cell_methods = \"time: mean\" ;\n",
"\t\tswdn_toa:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswdn_toa:interp_method = \"conserve_order1\" ;\n",
"\tfloat swdn_toa_clr(time, grid_yt, grid_xt) ;\n",
"\t\tswdn_toa_clr:long_name = \"clear skySW flux down at TOA\" ;\n",
"\t\tswdn_toa_clr:units = \"watts/m2\" ;\n",
"\t\tswdn_toa_clr:missing_value = -999.f ;\n",
"\t\tswdn_toa_clr:_FillValue = -999.f ;\n",
"\t\tswdn_toa_clr:cell_methods = \"time: mean\" ;\n",
"\t\tswdn_toa_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswdn_toa_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat swup_200hPa(time, grid_yt, grid_xt) ;\n",
"\t\tswup_200hPa:long_name = \"SW flux up at 200 hPa\" ;\n",
"\t\tswup_200hPa:units = \"watts/m2\" ;\n",
"\t\tswup_200hPa:missing_value = -999.f ;\n",
"\t\tswup_200hPa:_FillValue = -999.f ;\n",
"\t\tswup_200hPa:cell_methods = \"time: mean\" ;\n",
"\t\tswup_200hPa:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswup_200hPa:interp_method = \"conserve_order1\" ;\n",
"\tfloat swup_200hPa_clr(time, grid_yt, grid_xt) ;\n",
"\t\tswup_200hPa_clr:long_name = \"clear skySW flux up at 200 hPa\" ;\n",
"\t\tswup_200hPa_clr:units = \"watts/m2\" ;\n",
"\t\tswup_200hPa_clr:missing_value = -999.f ;\n",
"\t\tswup_200hPa_clr:_FillValue = -999.f ;\n",
"\t\tswup_200hPa_clr:cell_methods = \"time: mean\" ;\n",
"\t\tswup_200hPa_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswup_200hPa_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat swup_lin_trop(time, grid_yt, grid_xt) ;\n",
"\t\tswup_lin_trop:long_name = \"SW flux up at linear tropopause\" ;\n",
"\t\tswup_lin_trop:units = \"watts/m2\" ;\n",
"\t\tswup_lin_trop:missing_value = -999.f ;\n",
"\t\tswup_lin_trop:_FillValue = -999.f ;\n",
"\t\tswup_lin_trop:cell_methods = \"time: mean\" ;\n",
"\t\tswup_lin_trop:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswup_lin_trop:interp_method = \"conserve_order1\" ;\n",
"\tfloat swup_lin_trop_clr(time, grid_yt, grid_xt) ;\n",
"\t\tswup_lin_trop_clr:long_name = \"clear skySW flux up at linear tropopause\" ;\n",
"\t\tswup_lin_trop_clr:units = \"watts/m2\" ;\n",
"\t\tswup_lin_trop_clr:missing_value = -999.f ;\n",
"\t\tswup_lin_trop_clr:_FillValue = -999.f ;\n",
"\t\tswup_lin_trop_clr:cell_methods = \"time: mean\" ;\n",
"\t\tswup_lin_trop_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswup_lin_trop_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat swup_sfc(time, grid_yt, grid_xt) ;\n",
"\t\tswup_sfc:long_name = \"SW flux up at surface\" ;\n",
"\t\tswup_sfc:units = \"watts/m2\" ;\n",
"\t\tswup_sfc:missing_value = -999.f ;\n",
"\t\tswup_sfc:_FillValue = -999.f ;\n",
"\t\tswup_sfc:cell_methods = \"time: mean\" ;\n",
"\t\tswup_sfc:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswup_sfc:interp_method = \"conserve_order1\" ;\n",
"\tfloat swup_sfc_clr(time, grid_yt, grid_xt) ;\n",
"\t\tswup_sfc_clr:long_name = \"clear skySW flux up at surface\" ;\n",
"\t\tswup_sfc_clr:units = \"watts/m2\" ;\n",
"\t\tswup_sfc_clr:missing_value = -999.f ;\n",
"\t\tswup_sfc_clr:_FillValue = -999.f ;\n",
"\t\tswup_sfc_clr:cell_methods = \"time: mean\" ;\n",
"\t\tswup_sfc_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswup_sfc_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat swup_toa(time, grid_yt, grid_xt) ;\n",
"\t\tswup_toa:long_name = \"SW flux up at TOA\" ;\n",
"\t\tswup_toa:units = \"watts/m2\" ;\n",
"\t\tswup_toa:missing_value = -999.f ;\n",
"\t\tswup_toa:_FillValue = -999.f ;\n",
"\t\tswup_toa:cell_methods = \"time: mean\" ;\n",
"\t\tswup_toa:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswup_toa:interp_method = \"conserve_order1\" ;\n",
"\tfloat swup_toa_clr(time, grid_yt, grid_xt) ;\n",
"\t\tswup_toa_clr:long_name = \"clear skySW flux up at TOA\" ;\n",
"\t\tswup_toa_clr:units = \"watts/m2\" ;\n",
"\t\tswup_toa_clr:missing_value = -999.f ;\n",
"\t\tswup_toa_clr:_FillValue = -999.f ;\n",
"\t\tswup_toa_clr:cell_methods = \"time: mean\" ;\n",
"\t\tswup_toa_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tswup_toa_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat t_ref(time, grid_yt, grid_xt) ;\n",
"\t\tt_ref:long_name = \"temperature at 2 m\" ;\n",
"\t\tt_ref:units = \"deg_k\" ;\n",
"\t\tt_ref:valid_range = 100.f, 400.f ;\n",
"\t\tt_ref:missing_value = 1.e+20f ;\n",
"\t\tt_ref:_FillValue = 1.e+20f ;\n",
"\t\tt_ref:cell_methods = \"time: mean\" ;\n",
"\t\tt_ref:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tt_ref:interp_method = \"conserve_order1\" ;\n",
"\tfloat t_ref_max(time, grid_yt, grid_xt) ;\n",
"\t\tt_ref_max:long_name = \"temperature at 2 m\" ;\n",
"\t\tt_ref_max:units = \"deg_k\" ;\n",
"\t\tt_ref_max:valid_range = 100.f, 400.f ;\n",
"\t\tt_ref_max:missing_value = 1.e+20f ;\n",
"\t\tt_ref_max:_FillValue = 1.e+20f ;\n",
"\t\tt_ref_max:cell_methods = \"time: max\" ;\n",
"\t\tt_ref_max:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tt_ref_max:interp_method = \"conserve_order1\" ;\n",
"\tfloat t_ref_min(time, grid_yt, grid_xt) ;\n",
"\t\tt_ref_min:long_name = \"temperature at 2 m\" ;\n",
"\t\tt_ref_min:units = \"deg_k\" ;\n",
"\t\tt_ref_min:valid_range = 100.f, 400.f ;\n",
"\t\tt_ref_min:missing_value = 1.e+20f ;\n",
"\t\tt_ref_min:_FillValue = 1.e+20f ;\n",
"\t\tt_ref_min:cell_methods = \"time: min\" ;\n",
"\t\tt_ref_min:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tt_ref_min:interp_method = \"conserve_order1\" ;\n",
"\tfloat t_surf(time, grid_yt, grid_xt) ;\n",
"\t\tt_surf:long_name = \"surface temperature\" ;\n",
"\t\tt_surf:units = \"deg_k\" ;\n",
"\t\tt_surf:valid_range = 100.f, 400.f ;\n",
"\t\tt_surf:missing_value = 1.e+20f ;\n",
"\t\tt_surf:_FillValue = 1.e+20f ;\n",
"\t\tt_surf:cell_methods = \"time: mean\" ;\n",
"\t\tt_surf:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tt_surf:interp_method = \"conserve_order1\" ;\n",
"\tfloat tau_x(time, grid_yt, grid_xt) ;\n",
"\t\ttau_x:long_name = \"zonal wind stress\" ;\n",
"\t\ttau_x:units = \"pa\" ;\n",
"\t\ttau_x:missing_value = 1.e+20f ;\n",
"\t\ttau_x:_FillValue = 1.e+20f ;\n",
"\t\ttau_x:cell_methods = \"time: mean\" ;\n",
"\t\ttau_x:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttau_x:interp_method = \"conserve_order1\" ;\n",
"\tfloat tau_y(time, grid_yt, grid_xt) ;\n",
"\t\ttau_y:long_name = \"meridional wind stress\" ;\n",
"\t\ttau_y:units = \"pa\" ;\n",
"\t\ttau_y:missing_value = 1.e+20f ;\n",
"\t\ttau_y:_FillValue = 1.e+20f ;\n",
"\t\ttau_y:cell_methods = \"time: mean\" ;\n",
"\t\ttau_y:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttau_y:interp_method = \"conserve_order1\" ;\n",
"\tfloat taubx(time, grid_yt, grid_xt) ;\n",
"\t\ttaubx:long_name = \"x base flux for grav wave drag\" ;\n",
"\t\ttaubx:units = \"N/m2\" ;\n",
"\t\ttaubx:missing_value = -999.f ;\n",
"\t\ttaubx:_FillValue = -999.f ;\n",
"\t\ttaubx:cell_methods = \"time: mean\" ;\n",
"\t\ttaubx:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttaubx:interp_method = \"conserve_order1\" ;\n",
"\tfloat tauby(time, grid_yt, grid_xt) ;\n",
"\t\ttauby:long_name = \"y base flux for grav wave drag\" ;\n",
"\t\ttauby:units = \"N/m2\" ;\n",
"\t\ttauby:missing_value = -999.f ;\n",
"\t\ttauby:_FillValue = -999.f ;\n",
"\t\ttauby:cell_methods = \"time: mean\" ;\n",
"\t\ttauby:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttauby:interp_method = \"conserve_order1\" ;\n",
"\tfloat tdt_conv(time, pfull, grid_yt, grid_xt) ;\n",
"\t\ttdt_conv:long_name = \"Temperature tendency from convection\" ;\n",
"\t\ttdt_conv:units = \"deg_K/s\" ;\n",
"\t\ttdt_conv:missing_value = -999.f ;\n",
"\t\ttdt_conv:_FillValue = -999.f ;\n",
"\t\ttdt_conv:cell_methods = \"time: mean\" ;\n",
"\t\ttdt_conv:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttdt_conv:interp_method = \"conserve_order1\" ;\n",
"\tfloat tdt_ls(time, pfull, grid_yt, grid_xt) ;\n",
"\t\ttdt_ls:long_name = \"Temperature tendency from strat cloud\" ;\n",
"\t\ttdt_ls:units = \"deg_K/s\" ;\n",
"\t\ttdt_ls:missing_value = -999.f ;\n",
"\t\ttdt_ls:_FillValue = -999.f ;\n",
"\t\ttdt_ls:cell_methods = \"time: mean\" ;\n",
"\t\ttdt_ls:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttdt_ls:interp_method = \"conserve_order1\" ;\n",
"\tfloat tdt_lw(time, pfull, grid_yt, grid_xt) ;\n",
"\t\ttdt_lw:long_name = \"temperature tendency for LW radiation\" ;\n",
"\t\ttdt_lw:units = \"deg_K/sec\" ;\n",
"\t\ttdt_lw:missing_value = -999.f ;\n",
"\t\ttdt_lw:_FillValue = -999.f ;\n",
"\t\ttdt_lw:cell_methods = \"time: mean\" ;\n",
"\t\ttdt_lw:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttdt_lw:interp_method = \"conserve_order1\" ;\n",
"\tfloat tdt_lw_clr(time, pfull, grid_yt, grid_xt) ;\n",
"\t\ttdt_lw_clr:long_name = \"clear skytemperature tendency for LW radiation\" ;\n",
"\t\ttdt_lw_clr:units = \"deg_K/sec\" ;\n",
"\t\ttdt_lw_clr:missing_value = -999.f ;\n",
"\t\ttdt_lw_clr:_FillValue = -999.f ;\n",
"\t\ttdt_lw_clr:cell_methods = \"time: mean\" ;\n",
"\t\ttdt_lw_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttdt_lw_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat tdt_sw(time, pfull, grid_yt, grid_xt) ;\n",
"\t\ttdt_sw:long_name = \"temperature tendency for SW radiation\" ;\n",
"\t\ttdt_sw:units = \"deg_K/sec\" ;\n",
"\t\ttdt_sw:missing_value = -999.f ;\n",
"\t\ttdt_sw:_FillValue = -999.f ;\n",
"\t\ttdt_sw:cell_methods = \"time: mean\" ;\n",
"\t\ttdt_sw:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttdt_sw:interp_method = \"conserve_order1\" ;\n",
"\tfloat tdt_sw_clr(time, pfull, grid_yt, grid_xt) ;\n",
"\t\ttdt_sw_clr:long_name = \"clear skytemperature tendency for SW radiation\" ;\n",
"\t\ttdt_sw_clr:units = \"deg_K/sec\" ;\n",
"\t\ttdt_sw_clr:missing_value = -999.f ;\n",
"\t\ttdt_sw_clr:_FillValue = -999.f ;\n",
"\t\ttdt_sw_clr:cell_methods = \"time: mean\" ;\n",
"\t\ttdt_sw_clr:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttdt_sw_clr:interp_method = \"conserve_order1\" ;\n",
"\tfloat temp(time, pfull, grid_yt, grid_xt) ;\n",
"\t\ttemp:long_name = \"temperature\" ;\n",
"\t\ttemp:units = \"K\" ;\n",
"\t\ttemp:valid_range = 100.f, 350.f ;\n",
"\t\ttemp:missing_value = -1.e+10f ;\n",
"\t\ttemp:_FillValue = -1.e+10f ;\n",
"\t\ttemp:cell_methods = \"time: mean\" ;\n",
"\t\ttemp:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttemp:interp_method = \"conserve_order1\" ;\n",
"\tdouble time(time) ;\n",
"\t\ttime:long_name = \"time\" ;\n",
"\t\ttime:units = \"days since 0001-01-01 00:00:00\" ;\n",
"\t\ttime:cartesian_axis = \"T\" ;\n",
"\t\ttime:calendar_type = \"NOLEAP\" ;\n",
"\t\ttime:calendar = \"NOLEAP\" ;\n",
"\t\ttime:bounds = \"time_bounds\" ;\n",
"\tdouble time_bounds(time, bnds) ;\n",
"\t\ttime_bounds:long_name = \"time axis boundaries\" ;\n",
"\t\ttime_bounds:units = \"days\" ;\n",
"\t\ttime_bounds:missing_value = 1.e+20 ;\n",
"\t\ttime_bounds:_FillValue = 1.e+20 ;\n",
"\tfloat tot_cld_amt(time, grid_yt, grid_xt) ;\n",
"\t\ttot_cld_amt:long_name = \"total cloud amount\" ;\n",
"\t\ttot_cld_amt:units = \"percent\" ;\n",
"\t\ttot_cld_amt:missing_value = 1.e+20f ;\n",
"\t\ttot_cld_amt:_FillValue = 1.e+20f ;\n",
"\t\ttot_cld_amt:cell_methods = \"time: mean\" ;\n",
"\t\ttot_cld_amt:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\ttot_cld_amt:interp_method = \"conserve_order1\" ;\n",
"\tfloat u_ref(time, grid_yt, grid_xt) ;\n",
"\t\tu_ref:long_name = \"zonal wind component at 10 m\" ;\n",
"\t\tu_ref:units = \"m/s\" ;\n",
"\t\tu_ref:valid_range = -400.f, 400.f ;\n",
"\t\tu_ref:missing_value = 1.e+20f ;\n",
"\t\tu_ref:_FillValue = 1.e+20f ;\n",
"\t\tu_ref:cell_methods = \"time: mean\" ;\n",
"\t\tu_ref:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tu_ref:interp_method = \"conserve_order1\" ;\n",
"\tfloat ucomp(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tucomp:long_name = \"zonal wind\" ;\n",
"\t\tucomp:units = \"m/sec\" ;\n",
"\t\tucomp:valid_range = -330.f, 330.f ;\n",
"\t\tucomp:missing_value = -1.e+10f ;\n",
"\t\tucomp:_FillValue = -1.e+10f ;\n",
"\t\tucomp:cell_methods = \"time: mean\" ;\n",
"\t\tucomp:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tucomp:interp_method = \"conserve_order1\" ;\n",
"\tfloat v_ref(time, grid_yt, grid_xt) ;\n",
"\t\tv_ref:long_name = \"meridional wind component at 10 m\" ;\n",
"\t\tv_ref:units = \"m/s\" ;\n",
"\t\tv_ref:valid_range = -400.f, 400.f ;\n",
"\t\tv_ref:missing_value = 1.e+20f ;\n",
"\t\tv_ref:_FillValue = 1.e+20f ;\n",
"\t\tv_ref:cell_methods = \"time: mean\" ;\n",
"\t\tv_ref:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tv_ref:interp_method = \"conserve_order1\" ;\n",
"\tfloat vcomp(time, pfull, grid_yt, grid_xt) ;\n",
"\t\tvcomp:long_name = \"meridional wind\" ;\n",
"\t\tvcomp:units = \"m/sec\" ;\n",
"\t\tvcomp:valid_range = -330.f, 330.f ;\n",
"\t\tvcomp:missing_value = -1.e+10f ;\n",
"\t\tvcomp:_FillValue = -1.e+10f ;\n",
"\t\tvcomp:cell_methods = \"time: mean\" ;\n",
"\t\tvcomp:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tvcomp:interp_method = \"conserve_order1\" ;\n",
"\tfloat vort850(time, grid_yt, grid_xt) ;\n",
"\t\tvort850:long_name = \"850-mb vorticity\" ;\n",
"\t\tvort850:units = \"1/s\" ;\n",
"\t\tvort850:missing_value = -1.e+10f ;\n",
"\t\tvort850:_FillValue = -1.e+10f ;\n",
"\t\tvort850:cell_methods = \"time: mean\" ;\n",
"\t\tvort850:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\tvort850:interp_method = \"conserve_order1\" ;\n",
"\tfloat wind(time, grid_yt, grid_xt) ;\n",
"\t\twind:long_name = \"wind speed for flux calculations\" ;\n",
"\t\twind:units = \"m/s\" ;\n",
"\t\twind:valid_range = 0.f, 400.f ;\n",
"\t\twind:missing_value = 1.e+20f ;\n",
"\t\twind:_FillValue = 1.e+20f ;\n",
"\t\twind:cell_methods = \"time: mean\" ;\n",
"\t\twind:time_avg_info = \"average_T1,average_T2,average_DT\" ;\n",
"\t\twind:interp_method = \"conserve_order1\" ;\n",
"\tfloat zsurf(grid_yt, grid_xt) ;\n",
"\t\tzsurf:long_name = \"surface height\" ;\n",
"\t\tzsurf:units = \"m\" ;\n",
"\t\tzsurf:missing_value = 1.e+20f ;\n",
"\t\tzsurf:_FillValue = 1.e+20f ;\n",
"\t\tzsurf:cell_methods = \"time: point\" ;\n",
"\t\tzsurf:interp_method = \"conserve_order1\" ;\n",
"\n",
"// global attributes:\n",
"\t\t:filename = \"atmos_month.tile1.nc\" ;\n",
"\t\t:title = \"CM2.5_A_Control-1990_FLOR_B01\" ;\n",
"\t\t:grid_type = \"regular\" ;\n",
"\t\t:grid_tile = \"N/A\" ;\n",
"\t\t:history = \"Sun Jan 12 18:40:03 2020: ncks -7 -L 1 POSTP/00010101.atmos_month.nc POSTP_deflate1/00010101.atmos_month.nc\\n\",\n",
"\t\t\t\"/tigress/wenchang/CM2.5_tiger2_wy/CM2.5/bin/fregrid_parallel.tigercpu_intelmpi_18 --input_mosaic C180_mosaic.nc --input_file 00010101.atmos_month --interp_method conserve_order1 --remap_file fregrid_remap_file --nlon 576 --nlat 360 --scalar_field (**please see the field list in this file**)\" ;\n",
"\t\t:code_version = \"$Name: tikal_201409 $\" ;\n",
"\t\t:NCO = \"netCDF Operators version 4.8.1 (Homepage = http://nco.sf.net, Code = http://github.com/nco/nco)\" ;\n",
"}\n"
]
}
],
"source": [
"# \n",
"expname = 'CTL1860_noleap_tigercpu_intelmpi_18_576PE'\n",
"ifile = f'/tigress/wenchang/MODEL_OUT/{expname}/POSTP/00010101.atmos_month.nc'\n",
"print(ifile)\n",
"!ncdump -h $ifile"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## WY"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {
"ExecuteTime": {
"end_time": "2019-12-19T19:42:25.970633Z",
"start_time": "2019-12-19T19:42:25.967123Z"
}
},
"outputs": [
{
"data": {
"text/html": [
"\n"
],
"text/plain": [
""
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"%%html\n",
""
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.7.3"
},
"toc": {
"base_numbering": 1,
"nav_menu": {},
"number_sections": true,
"sideBar": true,
"skip_h1_title": false,
"title_cell": "Table of Contents",
"title_sidebar": "Contents",
"toc_cell": false,
"toc_position": {},
"toc_section_display": true,
"toc_window_display": true
}
},
"nbformat": 4,
"nbformat_minor": 2
}