; Summarize daily data from WRF runs from Lake/No Lake runs and get differences load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl" begin getlatlonfname = "/media/rbehnke/MyBook1/NWS_GFS/20161201/wrfout_d02_2016-12-01_03:00:00_Lake" z = addfile(getlatlonfname,"r") lat2d = z->XLAT(0,:,:) lon2d = z->XLONG(0,:,:) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; Create netcdf file to write data to fname = "/home/rbehnke/Desktop/Test_Summary_d02.nc" setfileoption("nc","format","netcdf4") setfileoption("nc","headerReserveSpace",64000) setfileoption("nc","preFill",False) setfileoption("nc","defineMode",True) if (isfilepresent(fname)) then system("rm " + fname) end if ;Create the file f = addfile(fname,"c") ;Define dimensions no_days = 11 dim_names = (/"time","y","x"/) dim_sizes = (/no_days,120,120/) dim_unlimited = (/False,False,False/) filedimdef(f,dim_names,dim_sizes,dim_unlimited) ;Define file properties atts = True atts@description = "This file contains differences between WRF run with and without Flathead Lake for several variables." print("Defining file attributes") fileattdef(f,atts) ;Define time, lat, lon variables print("Defining file variables") filevardef(f,"time","string","time") filevardef(f,"lat","float",(/"y","x"/)) filevardef(f,"lon","float",(/"y","x"/)) timeatts = True timeatts@long_name = "Date" timeatts@units = "Date" filevarattdef(f,"time",timeatts) lonatts = True lonatts@long_name = "Longitude" lonatts@units = "degrees_east" filevarattdef(f,"lon",lonatts) latatts = True latatts@long_name = "Latitude" latatts@units = "degrees_north" filevarattdef(f,"lat",latatts) ;Define PHYSIOGRAPHIC Variables and their attributes (might need revision!) filevardef(f,"albedo","float",(/"y","x"/)); mean daily albedo - taken from no lake run filevardef(f,"terrain_height","float",(/"y","x"/)); terrain height (m) filevardef(f,"soil_type","integer",(/"y","x"/)); dominant soil category - taken from no lake run filevardef(f,"veg_type","integer",(/"y","x"/)); dominant vegetation type category - taken from no lake run filevardef(f,"lai","float",(/"y","x"/)); mean daily leaf area index - taken from no lake run filevardef(f,"lulc","float",(/"y","x"/)); land use category - taken from no lake run filevardef(f,"soil_moisture","float",(/"y","x"/)); mean daily soil moisture (top layer) - taken from no lake run filevardef(f,"sfc_skin_temp","float",(/"y","x"/)); mean daily surface skin temperature - taken from no lake run filevardef(f,"soil_temp","float",(/"y","x"/)); mean daily soil temperature (top layer) - taken from no lake run filevardef(f,"veg_frac","float",(/"y","x"/)); mean daily vegetation fraction - taken from no lake run albedo_atts = True albedo_atts@long_name = "Mean Daily Surface Albedo" albedo_atts@units = "fraction" filevarattdef(f,"albedo",albedo_atts) terrain_height_atts = True terrain_height_atts@long_name = "Terrain Height" terrain_height_atts@units = "m" filevarattdef(f,"terrain_height",terrain_height_atts) soil_type_atts = True soil_type_atts@long_name = "Dominant Soil Category" soil_type_atts@units = "" filevarattdef(f,"soil_type",soil_type_atts) veg_type_atts = True veg_type_atts@long_name = "Dominant Vegetation Category" veg_type_atts@units = "" filevarattdef(f,"veg_type",veg_type_atts) lai_atts = True lai_atts@long_name = "Mean Daily Leaf Area Index" lai_atts@units = "m2 m-2" filevarattdef(f,"lai",lai_atts) lulc_atts = True lulc_atts@long_name = "Land Use Category" lulc_atts@units = "" filevarattdef(f,"lulc",lulc_atts) soil_moisture_atts = True soil_moisture_atts@long_name = "Mean Daily Soil Moisture (Top Layer)" soil_moisture_atts@units = "m3 m-3" filevarattdef(f,"soil_moisture",soil_moisture_atts) sfc_skin_temp_atts = True sfc_skin_temp_atts@long_name = "Mean Daily Surface Skin Temperature" sfc_skin_temp_atts@units = "degrees Celsius" filevarattdef(f,"sfc_skin_temp",sfc_skin_temp_atts) soil_temp_atts = True soil_temp_atts@long_name = "Mean Daily Soil Temperature (Top Layer)" soil_temp_atts@units = "degrees Celsius" filevarattdef(f,"soil_temp",soil_temp_atts) veg_frac_atts = True veg_frac_atts@long_name = "Vegetation Fraction" veg_frac_atts@units = "" filevarattdef(f,"veg_frac",veg_frac_atts) ;Define precipitation variables filevardef(f,"precipsum_lk","float",(/"time","y","x"/)) filevardef(f,"precipaccum_lk","float",(/"time","y","x"/)) filevardef(f,"ncprecipsum_lk","float",(/"time","y","x"/)) filevardef(f,"conprecipsum_lk","float",(/"time","y","x"/)) filevardef(f,"tot_swe_lk","float",(/"time","y","x"/)) filevardef(f,"tot_snowdepth_lk","float",(/"time","y","x"/)) precipsum_lkatts = True precipsum_lkatts@long_name = "Daily Total Precipitation - Lake Run" precipsum_lkatts@units = "mm" filevarattdef(f,"precipsum_lk",precipsum_lkatts) precipaccum_lkatts = True precipaccum_lkatts@long_name = "Accumulated Daily Total Precipitation - Lake Run" precipaccum_lkatts@units = "mm" filevarattdef(f,"precipaccum_lk",precipaccum_lkatts) ncprecipsum_lkatts = True ncprecipsum_lkatts@long_name = "Daily Non-Convective Total Precipitation - Lake Run" ncprecipsum_lkatts@units = "mm" filevarattdef(f,"ncprecipsum_lk",ncprecipsum_lkatts) conprecipsum_lkatts = True conprecipsum_lkatts@long_name = "Daily Convective Total Precipitation - Lake Run" conprecipsum_lkatts@units = "mm" filevarattdef(f,"conprecipsum_lk",conprecipsum_lkatts) tot_swe_lkatts = True tot_swe_lkatts@long_name = "Total Change in Daily Snow Water Equivalent - Lake Run" tot_swe_lkatts@units = "mm" filevarattdef(f,"tot_swe_lk",tot_swe_lkatts) tot_snowdepth_lkatts = True tot_snowdepth_lkatts@long_name = "Total Change in Daily Snow Depth - Lake Run" tot_snowdepth_lkatts@units = "m" filevarattdef(f,"tot_snowdepth_lk",tot_snowdepth_lkatts) filevardef(f,"precipsum_nolk","float",(/"time","y","x"/)) filevardef(f,"precipaccum_nolk","float",(/"time","y","x"/)) filevardef(f,"ncprecipsum_nolk","float",(/"time","y","x"/)) filevardef(f,"conprecipsum_nolk","float",(/"time","y","x"/)) filevardef(f,"tot_swe_nolk","float",(/"time","y","x"/)) filevardef(f,"tot_snowdepth_nolk","float",(/"time","y","x"/)) precipsum_nolkatts = True precipsum_nolkatts@long_name = "Daily Total Precipitation - No Lake Run" precipsum_nolkatts@units = "mm" filevarattdef(f,"precipsum_nolk",precipsum_nolkatts) precipaccum_nolkatts = True precipaccum_nolkatts@long_name = "Accumulated Daily Total Precipitation - No Lake Run" precipaccum_nolkatts@units = "mm" filevarattdef(f,"precipaccum_nolk",precipaccum_nolkatts) ncprecipsum_nolkatts = True ncprecipsum_nolkatts@long_name = "Daily Non-Convective Total Precipitation - No Lake Run" ncprecipsum_nolkatts@units = "mm" filevarattdef(f,"ncprecipsum_nolk",ncprecipsum_nolkatts) conprecipsum_nolkatts = True conprecipsum_nolkatts@long_name = "Daily Convective Total Precipitation - No Lake Run" conprecipsum_nolkatts@units = "mm" filevarattdef(f,"conprecipsum_nolk",conprecipsum_nolkatts) tot_swe_nolkatts = True tot_swe_nolkatts@long_name = "Total Change in Daily Snow Water Equivalent - No Lake Run" tot_swe_nolkatts@units = "mm" filevarattdef(f,"tot_swe_nolk",tot_swe_nolkatts) tot_snowdepth_nolkatts = True tot_snowdepth_nolkatts@long_name = "Total Change in Daily Snow Depth - No Lake Run" tot_snowdepth_nolkatts@units = "m" filevarattdef(f,"tot_snowdepth_nolk",tot_snowdepth_nolkatts) filevardef(f,"precipsum_diff","float",(/"time","y","x"/)) filevardef(f,"precipaccum_diff","float",(/"time","y","x"/)) filevardef(f,"ncprecipsum_diff","float",(/"time","y","x"/)) filevardef(f,"conprecipsum_diff","float",(/"time","y","x"/)) filevardef(f,"swe_diff","float",(/"time","y","x"/)) filevardef(f,"snowdepth_diff","float",(/"time","y","x"/)) precipsum_diffatts = True precipsum_diffatts@long_name = "Daily Total Precipitation - No Lake Run minus Lake Run" precipsum_diffatts@units = "mm" filevarattdef(f,"precipsum_diff",precipsum_diffatts) precipaccum_diffatts = True precipaccum_diffatts@long_name = "Accumulated Daily Total Precipitation - No Lake Run minus Lake Run" precipaccum_diffatts@units = "mm" filevarattdef(f,"precipaccum_diff",precipaccum_diffatts) ncprecipsum_diffatts = True ncprecipsum_diffatts@long_name = "Daily Non-Convective Total Precipitation - No Lake Run minus Lake Run" ncprecipsum_diffatts@units = "mm" filevarattdef(f,"ncprecipsum_diff",ncprecipsum_diffatts) conprecipsum_diffatts = True conprecipsum_diffatts@long_name = "Daily Convective Total Precipitation - No Lake Run minus Lake Run" conprecipsum_diffatts@units = "mm" filevarattdef(f,"conprecipsum_diff",conprecipsum_diffatts) tot_swe_diffatts = True tot_swe_diffatts@long_name = "Total Change in Daily Snow Water Equivalent - No Lake Run minus Lake Run" tot_swe_diffatts@units = "mm" filevarattdef(f,"swe_diff",tot_swe_diffatts) tot_snowdepth_diffatts = True tot_snowdepth_diffatts@long_name = "Total Change in Daily Snow Depth - No Lake Run minus Lake Run" tot_snowdepth_diffatts@units = "m" filevarattdef(f,"snowdepth_diff",tot_snowdepth_diffatts) ;Define PBL variables filevardef(f,"min_pblh_lk","float",(/"time","y","x"/)) filevardef(f,"max_pblh_lk","float",(/"time","y","x"/)) filevardef(f,"avg_pblh_lk","float",(/"time","y","x"/)) min_pblh_lkatts = True min_pblh_lkatts@long_name = "Minimum Daily Planetary Boundary Layer Height - Lake Run" min_pblh_lkatts@units = "m" filevarattdef(f,"min_pblh_lk",min_pblh_lkatts) max_pblh_lkatts = True max_pblh_lkatts@long_name = "Maximum Daily Planetary Boundary Layer Height - Lake Run" max_pblh_lkatts@units = "m" filevarattdef(f,"max_pblh_lk",max_pblh_lkatts) avg_pblh_lkatts = True avg_pblh_lkatts@long_name = "Mean Daily Planetary Boundary Layer Height - Lake Run" avg_pblh_lkatts@units = "m" filevarattdef(f,"avg_pblh_lk",avg_pblh_lkatts) filevardef(f,"min_pblh_nolk","float",(/"time","y","x"/)) filevardef(f,"max_pblh_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_pblh_nolk","float",(/"time","y","x"/)) min_pblh_nolkatts = True min_pblh_nolkatts@long_name = "Minimum Daily Planetary Boundary Layer Height - No Lake Run" min_pblh_nolkatts@units = "m" filevarattdef(f,"min_pblh_nolk",min_pblh_nolkatts) max_pblh_nolkatts = True max_pblh_nolkatts@long_name = "Maximum Daily Planetary Boundary Layer Height - No Lake Run" max_pblh_nolkatts@units = "m" filevarattdef(f,"max_pblh_nolk",max_pblh_nolkatts) avg_pblh_nolkatts = True avg_pblh_nolkatts@long_name = "Mean Daily Planetary Boundary Layer Height - No Lake Run" avg_pblh_nolkatts@units = "m" filevarattdef(f,"avg_pblh_nolk",avg_pblh_nolkatts) filevardef(f,"min_pblh_diff","float",(/"time","y","x"/)) filevardef(f,"max_pblh_diff","float",(/"time","y","x"/)) filevardef(f,"avg_pblh_diff","float",(/"time","y","x"/)) min_pblh_diffatts = True min_pblh_diffatts@long_name = "Minimum Daily Planetary Boundary Layer Height - No Lake Run minus Lake Run" min_pblh_diffatts@units = "m" filevarattdef(f,"min_pblh_diff",min_pblh_diffatts) max_pblh_diffatts = True max_pblh_diffatts@long_name = "Maximum Daily Planetary Boundary Layer Height - No Lake Run minus Lake Run" max_pblh_diffatts@units = "m" filevarattdef(f,"max_pblh_diff",max_pblh_diffatts) avg_pblh_diffatts = True avg_pblh_diffatts@long_name = "Mean Daily Planetary Boundary Layer Height - No Lake Run minus Lake Run" avg_pblh_diffatts@units = "m" filevarattdef(f,"avg_pblh_diff",avg_pblh_diffatts) ;Define 2m Temperature variables filevardef(f,"min_t2_lk","float",(/"time","y","x"/)) filevardef(f,"max_t2_lk","float",(/"time","y","x"/)) filevardef(f,"avg_t2_lk","float",(/"time","y","x"/)) min_t2_lkatts = True min_t2_lkatts@long_name = "Minimum Daily 2m Temperature - Lake Run" min_t2_lkatts@units = "degrees Celsius" filevarattdef(f,"min_t2_lk",min_t2_lkatts) max_t2_lkatts = True max_t2_lkatts@long_name = "Maximum Daily 2m Temperature - Lake Run" max_t2_lkatts@units = "degrees Celsius" filevarattdef(f,"max_t2_lk",max_t2_lkatts) avg_t2_lkatts = True avg_t2_lkatts@long_name = "Mean Daily 2m Temperature - Lake Run" avg_t2_lkatts@units = "degrees Celsius" filevarattdef(f,"avg_t2_lk",avg_t2_lkatts) filevardef(f,"min_t2_nolk","float",(/"time","y","x"/)) filevardef(f,"max_t2_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_t2_nolk","float",(/"time","y","x"/)) min_t2_nolkatts = True min_t2_nolkatts@long_name = "Minimum Daily 2m Temperature - No Lake Run" min_t2_nolkatts@units = "degrees Celsius" filevarattdef(f,"min_t2_nolk",min_t2_nolkatts) max_t2_nolkatts = True max_t2_nolkatts@long_name = "Maximum Daily 2m Temperature - No Lake Run" max_t2_nolkatts@units = "degrees Celsius" filevarattdef(f,"max_t2_nolk",max_t2_nolkatts) avg_t2_nolkatts = True avg_t2_nolkatts@long_name = "Mean Daily 2m Temperature - No Lake Run" avg_t2_nolkatts@units = "degrees Celsius" filevarattdef(f,"avg_t2_nolk",avg_t2_nolkatts) filevardef(f,"min_t2_diff","float",(/"time","y","x"/)) filevardef(f,"max_t2_diff","float",(/"time","y","x"/)) filevardef(f,"avg_t2_diff","float",(/"time","y","x"/)) min_t2_diffatts = True min_t2_diffatts@long_name = "Minimum Daily 2m Temperature - No Lake Run minus Lake Run" min_t2_diffatts@units = "degrees Celsius" filevarattdef(f,"min_t2_diff",min_t2_diffatts) max_t2_diffatts = True max_t2_diffatts@long_name = "Maximum Daily 2m Temperature - No Lake Run minus Lake Run" max_t2_diffatts@units = "degrees Celsius" filevarattdef(f,"max_t2_diff",max_t2_diffatts) avg_t2_diffatts = True avg_t2_diffatts@long_name = "Mean Daily 2m Temperature - No Lake Run minus Lake Run" avg_t2_diffatts@units = "degrees Celsius" filevarattdef(f,"avg_t2_diff",avg_t2_diffatts) ;Define PSFC variables filevardef(f,"min_psfc_lk","float",(/"time","y","x"/)) filevardef(f,"max_psfc_lk","float",(/"time","y","x"/)) filevardef(f,"avg_psfc_lk","float",(/"time","y","x"/)) min_psfc_lkatts = True min_psfc_lkatts@long_name = "Minimum Daily Surface Pressure - Lake Run" min_psfc_lkatts@units = "Pascals" filevarattdef(f,"min_psfc_lk",min_psfc_lkatts) max_psfc_lkatts = True max_psfc_lkatts@long_name = "Maximum Daily Surface Pressure - Lake Run" max_psfc_lkatts@units = "Pascals" filevarattdef(f,"max_psfc_lk",max_psfc_lkatts) avg_psfc_lkatts = True avg_psfc_lkatts@long_name = "Mean Daily 2m Surface Pressure - Lake Run" avg_psfc_lkatts@units = "Pascals" filevarattdef(f,"avg_psfc_lk",avg_psfc_lkatts) filevardef(f,"min_psfc_nolk","float",(/"time","y","x"/)) filevardef(f,"max_psfc_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_psfc_nolk","float",(/"time","y","x"/)) min_psfc_nolkatts = True min_psfc_nolkatts@long_name = "Minimum Daily Surface Pressure - No Lake Run" min_psfc_nolkatts@units = "Pascals" filevarattdef(f,"min_psfc_nolk",min_psfc_nolkatts) max_psfc_nolkatts = True max_psfc_nolkatts@long_name = "Maximum Daily Surface Pressure - No Lake Run" max_psfc_nolkatts@units = "Pascals" filevarattdef(f,"max_psfc_nolk",max_psfc_nolkatts) avg_psfc_nolkatts = True avg_psfc_nolkatts@long_name = "Mean Daily 2m Surface Pressure - No Lake Run" avg_psfc_nolkatts@units = "Pascals" filevarattdef(f,"avg_psfc_nolk",avg_psfc_nolkatts) filevardef(f,"min_psfc_diff","float",(/"time","y","x"/)) filevardef(f,"max_psfc_diff","float",(/"time","y","x"/)) filevardef(f,"avg_psfc_diff","float",(/"time","y","x"/)) min_psfc_diffatts = True min_psfc_diffatts@long_name = "Minimum Daily Surface Pressure - No Lake Run minus Lake Run" min_psfc_diffatts@units = "Pascals" filevarattdef(f,"min_psfc_diff",min_psfc_diffatts) max_psfc_diffatts = True max_psfc_diffatts@long_name = "Maximum Daily Surface Pressure - No Lake Run minus Lake Run" max_psfc_diffatts@units = "Pascals" filevarattdef(f,"max_psfc_diff",max_psfc_diffatts) avg_psfc_diffatts = True avg_psfc_diffatts@long_name = "Mean Daily 2m Surface Pressure - No Lake Run minus Lake Run" avg_psfc_diffatts@units = "Pascals" filevarattdef(f,"avg_psfc_diff",avg_psfc_diffatts) ;Define Q2, RH, Td variables filevardef(f,"min_q2_lk","float",(/"time","y","x"/)) filevardef(f,"max_q2_lk","float",(/"time","y","x"/)) filevardef(f,"avg_q2_lk","float",(/"time","y","x"/)) min_q2_lkatts = True min_q2_lkatts@long_name = "Minimum 2m Daily Mixing Ratio - Lake Run" min_q2_lkatts@units = "grams per kilogram" filevarattdef(f,"min_q2_lk",min_q2_lkatts) max_q2_lkatts = True max_q2_lkatts@long_name = "Maximum 2m Daily Mixing Ratio - Lake Run" max_q2_lkatts@units = "grams per kilogram" filevarattdef(f,"max_q2_lk",max_q2_lkatts) avg_q2_lkatts = True avg_q2_lkatts@long_name = "Mean Daily 2m Mixing Ratio - Lake Run" avg_q2_lkatts@units = "grams per kilogram" filevarattdef(f,"avg_q2_lk",avg_q2_lkatts) filevardef(f,"min_rh2_lk","float",(/"time","y","x"/)) filevardef(f,"max_rh2_lk","float",(/"time","y","x"/)) filevardef(f,"avg_rh2_lk","float",(/"time","y","x"/)) min_rh2_lkatts = True min_rh2_lkatts@long_name = "Minimum 2m Daily Relative Humidity - Lake Run" min_rh2_lkatts@units = "percent" filevarattdef(f,"min_rh2_lk",min_rh2_lkatts) max_rh2_lkatts = True max_rh2_lkatts@long_name = "Maximum 2m Daily Relative Humidity - Lake Run" max_rh2_lkatts@units = "percent" filevarattdef(f,"max_rh2_lk",max_rh2_lkatts) avg_rh2_lkatts = True avg_rh2_lkatts@long_name = "Mean Daily 2m Relative Humidity - Lake Run" avg_rh2_lkatts@units = "percent" filevarattdef(f,"avg_rh2_lk",avg_rh2_lkatts) filevardef(f,"min_td2_lk","float",(/"time","y","x"/)) filevardef(f,"max_td2_lk","float",(/"time","y","x"/)) filevardef(f,"avg_td2_lk","float",(/"time","y","x"/)) min_td2_lkatts = True min_td2_lkatts@long_name = "Minimum 2m Daily Dew Point - Lake Run" min_td2_lkatts@units = "degrees Celsius" filevarattdef(f,"min_td2_lk",min_td2_lkatts) max_td2_lkatts = True max_td2_lkatts@long_name = "Maximum 2m Daily Dew Point - Lake Run" max_td2_lkatts@units = "degrees Celsius" filevarattdef(f,"max_td2_lk",max_td2_lkatts) avg_td2_lkatts = True avg_td2_lkatts@long_name = "Mean Daily 2m Dew Point - Lake Run" avg_td2_lkatts@units = "degrees Celsius" filevarattdef(f,"avg_td2_lk",avg_td2_lkatts) filevardef(f,"min_q2_nolk","float",(/"time","y","x"/)) filevardef(f,"max_q2_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_q2_nolk","float",(/"time","y","x"/)) min_q2_nolkatts = True min_q2_nolkatts@long_name = "Minimum 2m Daily Mixing Ratio - No Lake Run" min_q2_nolkatts@units = "grams per kilogram" filevarattdef(f,"min_q2_nolk",min_q2_nolkatts) max_q2_nolkatts = True max_q2_nolkatts@long_name = "Maximum 2m Daily Mixing Ratio - No Lake Run" max_q2_nolkatts@units = "grams per kilogram" filevarattdef(f,"max_q2_nolk",max_q2_nolkatts) avg_q2_nolkatts = True avg_q2_nolkatts@long_name = "Mean Daily 2m Mixing Ratio - No Lake Run" avg_q2_nolkatts@units = "grams per kilogram" filevarattdef(f,"avg_q2_nolk",avg_q2_nolkatts) filevardef(f,"min_rh2_nolk","float",(/"time","y","x"/)) filevardef(f,"max_rh2_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_rh2_nolk","float",(/"time","y","x"/)) min_rh2_nolkatts = True min_rh2_nolkatts@long_name = "Minimum 2m Daily Relative Humidity - No Lake Run" min_rh2_nolkatts@units = "percent" filevarattdef(f,"min_rh2_nolk",min_rh2_nolkatts) max_rh2_nolkatts = True max_rh2_nolkatts@long_name = "Maximum 2m Daily Relative Humidity - No Lake Run" max_rh2_nolkatts@units = "percent" filevarattdef(f,"max_rh2_nolk",max_rh2_nolkatts) avg_rh2_nolkatts = True avg_rh2_nolkatts@long_name = "Mean Daily 2m Relative Humidity - No Lake Run" avg_rh2_nolkatts@units = "percent" filevarattdef(f,"avg_rh2_nolk",avg_rh2_nolkatts) filevardef(f,"min_td2_nolk","float",(/"time","y","x"/)) filevardef(f,"max_td2_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_td2_nolk","float",(/"time","y","x"/)) min_td2_nolkatts = True min_td2_nolkatts@long_name = "Minimum 2m Daily Dew Point - No Lake Run" min_td2_nolkatts@units = "degrees Celsius" filevarattdef(f,"min_td2_nolk",min_td2_nolkatts) max_td2_nolkatts = True max_td2_nolkatts@long_name = "Maximum 2m Daily Dew Point - No Lake Run" max_td2_nolkatts@units = "degrees Celsius" filevarattdef(f,"max_td2_nolk",max_td2_nolkatts) avg_td2_nolkatts = True avg_td2_nolkatts@long_name = "Mean Daily 2m Dew Point - No Lake Run" avg_td2_nolkatts@units = "degrees Celsius" filevarattdef(f,"avg_td2_nolk",avg_td2_nolkatts) filevardef(f,"min_q2_diff","float",(/"time","y","x"/)) filevardef(f,"max_q2_diff","float",(/"time","y","x"/)) filevardef(f,"avg_q2_diff","float",(/"time","y","x"/)) min_q2_diffatts = True min_q2_diffatts@long_name = "Minimum 2m Daily Mixing Ratio - No Lake Run minus Lake Run" min_q2_diffatts@units = "grams per kilogram" filevarattdef(f,"min_q2_diff",min_q2_diffatts) max_q2_diffatts = True max_q2_diffatts@long_name = "Maximum 2m Daily Mixing Ratio - No Lake Run minus Lake Run" max_q2_diffatts@units = "grams per kilogram" filevarattdef(f,"max_q2_diff",max_q2_diffatts) avg_q2_diffatts = True avg_q2_diffatts@long_name = "Mean Daily 2m Mixing Ratio - No Lake Run minus Lake Run" avg_q2_diffatts@units = "grams per kilogram" filevarattdef(f,"avg_q2_diff",avg_q2_diffatts) filevardef(f,"min_rh2_diff","float",(/"time","y","x"/)) filevardef(f,"max_rh2_diff","float",(/"time","y","x"/)) filevardef(f,"avg_rh2_diff","float",(/"time","y","x"/)) min_rh2_diffatts = True min_rh2_diffatts@long_name = "Minimum 2m Daily Relative Humidity - No Lake Run minus Lake Run" min_rh2_diffatts@units = "percent" filevarattdef(f,"min_rh2_diff",min_rh2_diffatts) max_rh2_diffatts = True max_rh2_diffatts@long_name = "Maximum 2m Daily Relative Humidity - No Lake Run minus Lake Run" max_rh2_diffatts@units = "percent" filevarattdef(f,"max_rh2_diff",max_rh2_diffatts) avg_rh2_diffatts = True avg_rh2_diffatts@long_name = "Mean Daily 2m Relative Humidity - No Lake Run minus Lake Run" avg_rh2_diffatts@units = "percent" filevarattdef(f,"avg_rh2_diff",avg_rh2_diffatts) filevardef(f,"min_td2_diff","float",(/"time","y","x"/)) filevardef(f,"max_td2_diff","float",(/"time","y","x"/)) filevardef(f,"avg_td2_diff","float",(/"time","y","x"/)) min_td2_diffatts = True min_td2_diffatts@long_name = "Minimum 2m Daily Dew Point - No Lake Run minus Lake Run" min_td2_diffatts@units = "degrees Celsius" filevarattdef(f,"min_td2_diff",min_td2_diffatts) max_td2_diffatts = True max_td2_diffatts@long_name = "Maximum 2m Daily Dew Point - No Lake Run minus Lake Run" max_td2_diffatts@units = "degrees Celsius" filevarattdef(f,"max_td2_diff",max_td2_diffatts) avg_td2_diffatts = True avg_td2_diffatts@long_name = "Mean Daily 2m Dew Point - No Lake Run minus Lake Run" avg_td2_diffatts@units = "degrees Celsius" filevarattdef(f,"avg_td2_diff",avg_td2_diffatts) ;Define U10, V10, WS10 variables filevardef(f,"min_u10_lk","float",(/"time","y","x"/)) filevardef(f,"max_u10_lk","float",(/"time","y","x"/)) filevardef(f,"min_v10_lk","float",(/"time","y","x"/)) filevardef(f,"max_v10_lk","float",(/"time","y","x"/)) filevardef(f,"min_ws10_lk","float",(/"time","y","x"/)) filevardef(f,"max_ws10_lk","float",(/"time","y","x"/)) min_u10_lkatts = True min_u10_lkatts@long_name = "Minimum 10m Daily U Wind - Lake Run" min_u10_lkatts@units = "m/s" filevarattdef(f,"min_u10_lk",min_u10_lkatts) max_u10_lkatts = True max_u10_lkatts@long_name = "Maximum 10m Daily U Wind - Lake Run" max_u10_lkatts@units = "m/s" filevarattdef(f,"max_u10_lk",max_u10_lkatts) min_v10_lkatts = True min_v10_lkatts@long_name = "Minimum 10m Daily V Wind - Lake Run" min_v10_lkatts@units = "m/s" filevarattdef(f,"min_v10_lk",min_v10_lkatts) max_v10_lkatts = True max_v10_lkatts@long_name = "Maximum 10m Daily V Wind - Lake Run" max_v10_lkatts@units = "m/s" filevarattdef(f,"max_v10_lk",max_v10_lkatts) min_ws10_lkatts = True min_ws10_lkatts@long_name = "Minimum 10m Daily Wind Speed - Lake Run" min_ws10_lkatts@units = "m/s" filevarattdef(f,"min_ws10_lk",min_ws10_lkatts) max_ws10_lkatts = True max_ws10_lkatts@long_name = "Maximum 10m Daily Wind Speed - Lake Run" max_ws10_lkatts@units = "m/s" filevarattdef(f,"max_ws10_lk",max_ws10_lkatts) filevardef(f,"min_u10_nolk","float",(/"time","y","x"/)) filevardef(f,"max_u10_nolk","float",(/"time","y","x"/)) filevardef(f,"min_v10_nolk","float",(/"time","y","x"/)) filevardef(f,"max_v10_nolk","float",(/"time","y","x"/)) filevardef(f,"min_ws10_nolk","float",(/"time","y","x"/)) filevardef(f,"max_ws10_nolk","float",(/"time","y","x"/)) min_u10_nolkatts = True min_u10_nolkatts@long_name = "Minimum 10m Daily U Wind - No Lake Run" min_u10_nolkatts@units = "m/s" filevarattdef(f,"min_u10_nolk",min_u10_nolkatts) max_u10_nolkatts = True max_u10_nolkatts@long_name = "Maximum 10m Daily U Wind - No Lake Run" max_u10_nolkatts@units = "m/s" filevarattdef(f,"max_u10_nolk",max_u10_nolkatts) min_v10_nolkatts = True min_v10_nolkatts@long_name = "Minimum 10m Daily V Wind - No Lake Run" min_v10_nolkatts@units = "m/s" filevarattdef(f,"min_v10_nolk",min_v10_nolkatts) max_v10_nolkatts = True max_v10_nolkatts@long_name = "Maximum 10m Daily V Wind - No Lake Run" max_v10_nolkatts@units = "m/s" filevarattdef(f,"max_v10_nolk",max_v10_nolkatts) min_ws10_nolkatts = True min_ws10_nolkatts@long_name = "Minimum 10m Daily Wind Speed - No Lake Run" min_ws10_nolkatts@units = "m/s" filevarattdef(f,"min_ws10_nolk",min_ws10_nolkatts) max_ws10_nolkatts = True max_ws10_nolkatts@long_name = "Maximum 10m Daily Wind Speed - No Lake Run" max_ws10_nolkatts@units = "m/s" filevarattdef(f,"max_ws10_nolk",max_ws10_nolkatts) filevardef(f,"wspd10_lk_avg","float",(/"time","y","x"/)) filevardef(f,"wspd10_lk_std","float",(/"time","y","x"/)) filevardef(f,"wdir10_lk_avg","float",(/"time","y","x"/)) filevardef(f,"wdir10_lk_std","float",(/"time","y","x"/)) filevardef(f,"avgU10_lk","float",(/"time","y","x"/)) filevardef(f,"avgV10_lk","float",(/"time","y","x"/)) wspd10_lk_avgatts = True wspd10_lk_avgatts@long_name = "Average 10m Daily Wind Speed - Lake Run" wspd10_lk_avgatts@units = "m/s" filevarattdef(f,"wspd10_lk_avg",wspd10_lk_avgatts) wspd10_lk_stdatts = True wspd10_lk_stdatts@long_name = "Average 10m Daily Wind Speed Standard Deviation - Lake Run" wspd10_lk_stdatts@units = "m/s" filevarattdef(f,"wspd10_lk_std",wspd10_lk_stdatts) wdir10_lk_avgatts = True wdir10_lk_avgatts@long_name = "Average 10m Daily Wind Direction - Lake Run" wdir10_lk_avgatts@units = "degrees" filevarattdef(f,"wdir10_lk_avg",wdir10_lk_avgatts) wdir10_lk_stdatts = True wdir10_lk_stdatts@long_name = "Average 10m Daily Wind Direction Standard Deviation - Lake Run" wdir10_lk_stdatts@units = "degrees" filevarattdef(f,"wdir10_lk_std",wdir10_lk_stdatts) avgU10_lkatts = True avgU10_lkatts@long_name = "Average 10m Daily U Wind Speed - Lake Run" avgU10_lkatts@units = "m/s" filevarattdef(f,"avgU10_lk",avgU10_lkatts) avgV10_lkatts = True avgV10_lkatts@long_name = "Average 10m Daily V Wind Speed - Lake Run" avgV10_lkatts@units = "m/s" filevarattdef(f,"avgV10_lk",avgV10_lkatts) filevardef(f,"wspd10_nolk_avg","float",(/"time","y","x"/)) filevardef(f,"wspd10_nolk_std","float",(/"time","y","x"/)) filevardef(f,"wdir10_nolk_avg","float",(/"time","y","x"/)) filevardef(f,"wdir10_nolk_std","float",(/"time","y","x"/)) filevardef(f,"avgU10_nolk","float",(/"time","y","x"/)) filevardef(f,"avgV10_nolk","float",(/"time","y","x"/)) wspd10_nolk_avgatts = True wspd10_nolk_avgatts@long_name = "Average 10m Daily Wind Speed - No Lake Run" wspd10_nolk_avgatts@units = "m/s" filevarattdef(f,"wspd10_nolk_avg",wspd10_nolk_avgatts) wspd10_nolk_stdatts = True wspd10_nolk_stdatts@long_name = "Average 10m Daily Wind Speed Standard Deviation - No Lake Run" wspd10_nolk_stdatts@units = "m/s" filevarattdef(f,"wspd10_nolk_std",wspd10_nolk_stdatts) wdir10_nolk_avgatts = True wdir10_nolk_avgatts@long_name = "Average 10m Daily Wind Direction - No Lake Run" wdir10_nolk_avgatts@units = "degrees" filevarattdef(f,"wdir10_nolk_avg",wdir10_nolk_avgatts) wdir10_nolk_stdatts = True wdir10_nolk_stdatts@long_name = "Average 10m Daily Wind Direction Standard Deviation - No Lake Run" wdir10_nolk_stdatts@units = "degrees" filevarattdef(f,"wdir10_nolk_std",wdir10_nolk_stdatts) avgU10_nolkatts = True avgU10_nolkatts@long_name = "Average 10m Daily U Wind Speed - No Lake Run" avgU10_nolkatts@units = "m/s" filevarattdef(f,"avgU10_nolk",avgU10_nolkatts) avgV10_nolkatts = True avgV10_nolkatts@long_name = "Average 10m Daily V Wind Speed - No Lake Run" avgV10_nolkatts@units = "m/s" filevarattdef(f,"avgV10_nolk",avgV10_nolkatts) filevardef(f,"min_u10_diff","float",(/"time","y","x"/)) filevardef(f,"max_u10_diff","float",(/"time","y","x"/)) filevardef(f,"avg_u10_diff","float",(/"time","y","x"/)) filevardef(f,"min_v10_diff","float",(/"time","y","x"/)) filevardef(f,"max_v10_diff","float",(/"time","y","x"/)) filevardef(f,"avg_v10_diff","float",(/"time","y","x"/)) min_u10_diffatts = True min_u10_diffatts@long_name = "Minimum 10m Daily U Wind - No Lake Run minus Lake Run" min_u10_diffatts@units = "m/s" filevarattdef(f,"min_u10_diff",min_u10_diffatts) max_u10_diffatts = True max_u10_diffatts@long_name = "Maximum 10m Daily U Wind - No Lake Run minus Lake Run" max_u10_diffatts@units = "m/s" filevarattdef(f,"max_u10_diff",max_u10_diffatts) avg_u10_diffatts = True avg_u10_diffatts@long_name = "Average 10m Daily U Wind - No Lake Run minus Lake Run" avg_u10_diffatts@units = "m/s" filevarattdef(f,"avg_u10_diff",avg_u10_diffatts) min_v10_diffatts = True min_v10_diffatts@long_name = "Minimum 10m Daily V Wind - No Lake Run minus Lake Run" min_v10_diffatts@units = "m/s" filevarattdef(f,"min_v10_diff",min_v10_diffatts) max_v10_diffatts = True max_v10_diffatts@long_name = "Maximum 10m Daily V Wind - No Lake Run minus Lake Run" max_v10_diffatts@units = "m/s" filevarattdef(f,"max_v10_diff",max_v10_diffatts) avg_v10_diffatts = True avg_v10_diffatts@long_name = "Average 10m Daily V Wind - No Lake Run minus Lake Run" avg_v10_diffatts@units = "m/s" filevarattdef(f,"avg_v10_diff",avg_v10_diffatts) filevardef(f,"min_ws10_diff","float",(/"time","y","x"/)) filevardef(f,"max_ws10_diff","float",(/"time","y","x"/)) filevardef(f,"avg_ws10_diff","float",(/"time","y","x"/)) filevardef(f,"wspd10_std_diff","float",(/"time","y","x"/)) filevardef(f,"wdir10_std_diff","float",(/"time","y","x"/)) min_ws10_diffatts = True min_ws10_diffatts@long_name = "Minimum 10m Daily Wind Speed - No Lake Run minus Lake Run" min_ws10_diffatts@units = "m/s" filevarattdef(f,"min_ws10_diff",min_ws10_diffatts) max_ws10_diffatts = True max_ws10_diffatts@long_name = "Maximum 10m Daily Wind Speed - No Lake Run minus Lake Run" max_ws10_diffatts@units = "m/s" filevarattdef(f,"max_ws10_diff",max_ws10_diffatts) avg_ws10_diffatts = True avg_ws10_diffatts@long_name = "Average 10m Daily Wind Speed - No Lake Run minus Lake Run" avg_ws10_diffatts@units = "m/s" filevarattdef(f,"avg_ws10_diff",avg_ws10_diffatts) wspd10_std_diffatts = True wspd10_std_diffatts@long_name = "Wind Speed Standard Deviation- No Lake Run minus Lake Run" wspd10_std_diffatts@units = "m/s" filevarattdef(f,"wspd10_std_diff",wspd10_std_diffatts) wdir10_std_diffatts = True wdir10_std_diffatts@long_name = "Wind Direction Standard Deviation - No Lake Run minus Lake Run" wdir10_std_diffatts@units = "m/s" filevarattdef(f,"wdir10_std_diff",wdir10_std_diffatts) filevardef(f,"avg_U_700_lk","float",(/"time","y","x"/)) filevardef(f,"avg_V_700_lk","float",(/"time","y","x"/)) filevardef(f,"avg_wspd_700_lk","float",(/"time","y","x"/)) filevardef(f,"avg_wdir_700_lk","float",(/"time","y","x"/)) avg_U_700_lkatts = True avg_U_700_lkatts@long_name = "Average 700mb U Wind - Lake Run" avg_U_700_lkatts@units = "m/s" filevarattdef(f,"avg_U_700_lk",avg_U_700_lkatts) avg_V_700_lkatts = True avg_V_700_lkatts@long_name = "Average 700mb V Wind - Lake Run" avg_V_700_lkatts@units = "m/s" filevarattdef(f,"avg_V_700_lk",avg_V_700_lkatts) avg_wspd_700_lkatts = True avg_wspd_700_lkatts@long_name = "Average 700mb Wind Speed - Lake Run" avg_wspd_700_lkatts@units = "m/s" filevarattdef(f,"avg_wspd_700_lk",avg_wspd_700_lkatts) avg_wdir_700_lkatts = True avg_wdir_700_lkatts@long_name = "Average 700mb Wind Direction - Lake Run" avg_wdir_700_lkatts@units = "degrees" filevarattdef(f,"avg_wdir_700_lk",avg_wdir_700_lkatts) filevardef(f,"avg_U_800_lk","float",(/"time","y","x"/)) filevardef(f,"avg_V_800_lk","float",(/"time","y","x"/)) filevardef(f,"avg_wspd_800_lk","float",(/"time","y","x"/)) filevardef(f,"avg_wdir_800_lk","float",(/"time","y","x"/)) avg_U_800_lkatts = True avg_U_800_lkatts@long_name = "Average 800mb U Wind - Lake Run" avg_U_800_lkatts@units = "m/s" filevarattdef(f,"avg_U_800_lk",avg_U_800_lkatts) avg_V_800_lkatts = True avg_V_800_lkatts@long_name = "Average 800mb V Wind - Lake Run" avg_V_800_lkatts@units = "m/s" filevarattdef(f,"avg_V_800_lk",avg_V_800_lkatts) avg_wspd_800_lkatts = True avg_wspd_800_lkatts@long_name = "Average 800mb Wind Speed - Lake Run" avg_wspd_800_lkatts@units = "m/s" filevarattdef(f,"avg_wspd_800_lk",avg_wspd_800_lkatts) avg_wdir_800_lkatts = True avg_wdir_800_lkatts@long_name = "Average 800mb Wind Direction - Lake Run" avg_wdir_800_lkatts@units = "degrees" filevarattdef(f,"avg_wdir_800_lk",avg_wdir_800_lkatts) filevardef(f,"avg_U_700_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_V_700_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_wspd_700_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_wdir_700_nolk","float",(/"time","y","x"/)) avg_U_700_nolkatts = True avg_U_700_nolkatts@long_name = "Average 700mb U Wind - No Lake Run" avg_U_700_nolkatts@units = "m/s" filevarattdef(f,"avg_U_700_nolk",avg_U_700_nolkatts) avg_V_700_nolkatts = True avg_V_700_nolkatts@long_name = "Average 700mb V Wind - No Lake Run" avg_V_700_nolkatts@units = "m/s" filevarattdef(f,"avg_V_700_nolk",avg_V_700_nolkatts) avg_wspd_700_nolkatts = True avg_wspd_700_nolkatts@long_name = "Average 700mb Wind Speed - No Lake Run" avg_wspd_700_nolkatts@units = "m/s" filevarattdef(f,"avg_wspd_700_nolk",avg_wspd_700_nolkatts) avg_wdir_700_nolkatts = True avg_wdir_700_nolkatts@long_name = "Average 700mb Wind Direction - No Lake Run" avg_wdir_700_nolkatts@units = "degrees" filevarattdef(f,"avg_wdir_700_nolk",avg_wdir_700_nolkatts) filevardef(f,"avg_U_800_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_V_800_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_wspd_800_nolk","float",(/"time","y","x"/)) filevardef(f,"avg_wdir_800_nolk","float",(/"time","y","x"/)) avg_U_800_nolkatts = True avg_U_800_nolkatts@long_name = "Average 800mb U Wind - No Lake Run" avg_U_800_nolkatts@units = "m/s" filevarattdef(f,"avg_U_800_nolk",avg_U_800_nolkatts) avg_V_800_nolkatts = True avg_V_800_nolkatts@long_name = "Average 800mb V Wind - No Lake Run" avg_V_800_nolkatts@units = "m/s" filevarattdef(f,"avg_V_800_nolk",avg_V_800_nolkatts) avg_wspd_800_nolkatts = True avg_wspd_800_nolkatts@long_name = "Average 800mb Wind Speed - No Lake Run" avg_wspd_800_nolkatts@units = "m/s" filevarattdef(f,"avg_wspd_800_nolk",avg_wspd_800_nolkatts) avg_wdir_800_nolkatts = True avg_wdir_800_nolkatts@long_name = "Average 800mb Wind Direction - No Lake Run" avg_wdir_800_nolkatts@units = "degrees" filevarattdef(f,"avg_wdir_800_nolk",avg_wdir_800_nolkatts) filevardef(f,"avg_U_700_diff","float",(/"time","y","x"/)) filevardef(f,"avg_U_800_diff","float",(/"time","y","x"/)) filevardef(f,"avg_V_700_diff","float",(/"time","y","x"/)) filevardef(f,"avg_V_800_diff","float",(/"time","y","x"/)) filevardef(f,"avg_wspd_700_diff","float",(/"time","y","x"/)) filevardef(f,"avg_wspd_800_diff","float",(/"time","y","x"/)) avg_U_700_diffatts = True avg_U_700_diffatts@long_name = "Average 700mb U Wind - No Lake Run minus Lake Run" avg_U_700_diffatts@units = "m/s" filevarattdef(f,"avg_U_700_diff",avg_U_700_diffatts) avg_U_800_diffatts = True avg_U_800_diffatts@long_name = "Average 800mb U Wind - No Lake Run minus Lake Run" avg_U_800_diffatts@units = "m/s" filevarattdef(f,"avg_U_800_diff",avg_U_800_diffatts) avg_V_700_diffatts = True avg_V_700_diffatts@long_name = "Average 700mb V Wind - No Lake Run minus Lake Run" avg_V_700_diffatts@units = "m/s" filevarattdef(f,"avg_V_700_diff",avg_V_700_diffatts) avg_V_800_diffatts = True avg_V_800_diffatts@long_name = "Average 800mb V Wind - No Lake Run minus Lake Run" avg_V_800_diffatts@units = "m/s" filevarattdef(f,"avg_V_800_diff",avg_V_800_diffatts) avg_wspd_700_diffatts = True avg_wspd_700_diffatts@long_name = "Average 700mb Wind Speed - No Lake Run minus Lake Run" avg_wspd_700_diffatts@units = "m/s" filevarattdef(f,"avg_wspd_700_diff",avg_wspd_700_diffatts) avg_wspd_800_diffatts = True avg_wspd_800_diffatts@long_name = "Average 800mb Wind Speed - No Lake Run minus Lake Run" avg_wspd_800_diffatts@units = "m/s" filevarattdef(f,"avg_wspd_800_diff",avg_wspd_800_diffatts) ; Define Empty arrays for Accumulated Precipitation precipaccum_lk = new((/120, 120/),float,"No_FillValue") precipaccum_nolk = new((/120, 120/),float,"No_FillValue") precipaccum_lk(:,:) = 0 precipaccum_nolk(:,:) = 0 dates = (/"2016-12-01", "2016-12-02", "2016-12-03", "2016-12-04", "2016-12-05", "2016-12-06", "2016-12-07", "2016-12-08", "2016-12-09", "2016-12-10","2016-12-11", "2016-12-12", "2016-12-13", "2016-12-14", "2016-12-15", "2016-12-16", "2016-12-17", "2016-12-18", "2016-12-19", "2016-12-20","2016-12-21", "2016-12-22", "2016-12-23", "2016-12-24", "2016-12-25", "2016-12-26", "2016-12-27", "2016-12-28", "2016-12-29", "2016-12-30", "2016-12-31","2017-01-01", "2017-01-02", "2017-01-03", "2017-01-04", "2017-01-05", "2017-01-06", "2017-01-07", "2017-01-08", "2017-01-09", "2017-01-10", "2017-01-11", "2017-01-12", "2017-01-13", "2017-01-14", "2017-01-15", "2017-01-16", "2017-01-17", "2017-01-18", "2017-01-19", "2017-01-20", "2017-01-21", "2017-01-22", "2017-01-23", "2017-01-24", "2017-01-25", "2017-01-26", "2017-01-27", "2017-01-28", "2017-01-29", "2017-01-30", "2017-01-31", "2017-02-01", "2017-02-02", "2017-02-03", "2017-02-04", "2017-02-05", "2017-02-06", "2017-02-07", "2017-02-08", "2017-02-09", "2017-02-10", "2017-02-11", "2017-02-12", "2017-02-13", "2017-02-14", "2017-02-15", "2017-02-16", "2017-02-17", "2017-02-18", "2017-02-19", "2017-02-20", "2017-02-21", "2017-02-22", "2017-02-23", "2017-02-24", "2017-02-25", "2017-02-26", "2017-02-27", "2017-02-28", "2017-03-01", "2017-03-02", "2017-03-03", "2017-03-04", "2017-03-05", "2017-03-06", "2017-03-07", "2017-03-08", "2017-03-09", "2017-03-10", "2017-03-11", "2017-03-12", "2017-03-13", "2017-03-14", "2017-03-15", "2017-03-16", "2017-03-17", "2017-03-18", "2017-03-19", "2017-03-20", "2017-03-21", "2017-03-22", "2017-03-23", "2017-03-24", "2017-03-25", "2017-03-26", "2017-03-27", "2017-03-28", "2017-03-29", "2017-03-30", "2017-03-31", "2017-04-01", "2017-04-02", "2017-04-03", "2017-04-04", "2017-04-05", "2017-04-06", "2017-04-07", "2017-04-08", "2017-04-09", "2017-04-10", "2017-04-11", "2017-04-12", "2017-04-13", "2017-04-14", "2017-04-15", "2017-04-16", "2017-04-17", "2017-04-18", "2017-04-19", "2017-04-20", "2017-04-21", "2017-04-22", "2017-04-23", "2017-04-24", "2017-04-25", "2017-04-26", "2017-04-27", "2017-04-28", "2017-04-29", "2017-04-30", "2017-05-01", "2017-05-02", "2017-05-03", "2017-05-04", "2017-05-05", "2017-05-06", "2017-05-07", "2017-05-08", "2017-05-09", "2017-05-10", "2017-05-11", "2017-05-12", "2017-05-13", "2017-05-14", "2017-05-15", "2017-05-16", "2017-05-17", "2017-05-18", "2017-05-19", "2017-05-20", "2017-05-21", "2017-05-22", "2017-05-23", "2017-05-24", "2017-05-25", "2017-05-26", "2017-05-27", "2017-05-28", "2017-05-29", "2017-05-30", "2017-05-31"/) print("starting loop") do n=0,10 print(n) print("reading files") ; WITH LAKE lake_filename = "/media/rbehnke/MyBook1/NWS_GFS/" + str_sub_str(dates(n),"-","") + "/wrfout_d02_" + dates(n) + "_03:00:00_Lake" a = addfile(lake_filename,"r") latitude = wrf_user_getvar(a,"XLAT",0) longitude = wrf_user_getvar(a,"XLONG",0) nolake_filename = "/media/rbehnke/MyBook1/NWS_GFS/" + str_sub_str(dates(n),"-","") + "/wrfout_d02_" + dates(n) + "_03:00:00_NoLake" b = addfile(nolake_filename,"r") print("getting data") ; PHYSIOGRAPHIC Variables albedo_a = wrf_user_getvar(b,"ALBEDO",-1) albedo_b = albedo_a(4:27,:,:) albedo = dim_avg_n(albedo_b,0) albedo@lon2d = lon2d albedo@lat2d = lat2d terrain_height_a = wrf_user_getvar(b,"HGT",-1) terrain_height = terrain_height_a(4,:,:) terrain_height@lon2d = lon2d terrain_height@lat2d = lat2d soil_type_a = wrf_user_getvar(b,"ISLTYP",-1) soil_type = soil_type_a(4,:,:) soil_type@lon2d = lon2d soil_type@lat2d = lat2d veg_type_a = wrf_user_getvar(b,"IVGTYP",-1) veg_type = veg_type_a(4,:,:) veg_type@lon2d = lon2d veg_type@lat2d = lat2d lai_a = wrf_user_getvar(b,"LAI",-1) lai_b = lai_a(4:27,:,:) lai = dim_avg_n(lai_b,0) lai@lon2d = lon2d lai@lat2d = lat2d lulc_a = wrf_user_getvar(b,"LU_INDEX",-1) lulc = lulc_a(4,:,:) lulc@lon2d = lon2d lulc@lat2d = lat2d soil_moisture_a = wrf_user_getvar(b,"SMOIS",-1) soil_moisture_b = soil_moisture_a(4:27,0,:,:) soil_moisture = dim_avg_n(soil_moisture_b,0) soil_moisture@lon2d = lon2d soil_moisture@lat2d = lat2d sfc_skin_temp_a = wrf_user_getvar(b,"TSK",-1) sfc_skin_temp_b = sfc_skin_temp_a(4:27,:,:) sfc_skin_temp = dim_avg_n(sfc_skin_temp_b,0) sfc_skin_temp@lon2d = lon2d sfc_skin_temp@lat2d = lat2d soil_temp_a = wrf_user_getvar(b,"TSLB",-1) soil_temp_b = soil_temp_a(4:27,0,:,:) soil_temp = dim_avg_n(soil_temp_b,0) soil_temp@lon2d = lon2d soil_temp@lat2d = lat2d veg_frac_a = wrf_user_getvar(b,"VEGFRA",-1) veg_frac = veg_frac_a(4,:,:) veg_frac@lon2d = lon2d veg_frac@lat2d = lat2d ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; Precipitation LAKE precip_nc_lake = wrf_user_getvar(a,"RAINNC",-1); Convective Precipitation precip_nc_lk = precip_nc_lake(4:27,:,:) precip_con_lake = wrf_user_getvar(a,"RAINC",-1); Non-Convective Precipitation precip_con_lk = precip_con_lake(4:27,:,:) precip_tot_lk = precip_nc_lk + precip_con_lk; Total Precipitation (Convective and Non-Convective) swe_lake = wrf_user_getvar(a,"SNOW",-1); Snow Water Equivalent swe_lk = swe_lake(4:27,:,:) snow_depth_lake = wrf_user_getvar(a,"SNOWH",-1); Snow Depth snow_depth_lk = snow_depth_lake(4:27,:,:) ; Total Daily Precipitation LAKE precipsum_lk = (precip_nc_lake(27,:,:) - precip_nc_lake(3,:,:)) + (precip_con_lake(27,:,:) - precip_con_lake(3,:,:)) ;total precipitation - all kinds******* ncprecipsum_lk = precip_nc_lake(27,:,:) - precip_nc_lake(3,:,:) ;total non-convective precipitation - all kinds******* conprecipsum_lk = precip_con_lake(27,:,:) - precip_con_lake(3,:,:) ;total convective precipitation - all kinds******* tot_swe_lk = swe_lake(27,:,:) - swe_lake(3,:,:) ;Daily change in SWE******* tot_snowdepth_lk = snow_depth_lake(27,:,:) - snow_depth_lake(3,:,:) ;Daily change in Snow Depth******* precipaccum_lk = precipaccum_lk + precipsum_lk ;Calculate New Daily Accumulated Precipitation******* ; Precipitation NO LAKE precip_nc_nolake = wrf_user_getvar(b,"RAINNC",-1); Convective Precipitation precip_nc_nolk = precip_nc_nolake(4:27,:,:) precip_con_nolake = wrf_user_getvar(b,"RAINC",-1); Non-Convective Precipitation precip_con_nolk = precip_con_nolake(4:27,:,:) precip_tot_nolk = precip_nc_nolk + precip_con_nolk; Total Precipitation (Convective and Non-Convective) swe_nolake = wrf_user_getvar(b,"SNOW",-1); Snow Water Equivalent swe_nolk = swe_nolake(4:27,:,:) snow_depth_nolake = wrf_user_getvar(b,"SNOWH",-1); Snow Depth snow_depth_nolk = snow_depth_nolake(4:27,:,:) ; Total Daily Precipitation NO LAKE precipsum_nolk = (precip_nc_nolake(27,:,:) - precip_nc_nolake(3,:,:)) + (precip_con_nolake(27,:,:) - precip_con_nolake(3,:,:)) ;total precipitation - all kinds******* ncprecipsum_nolk = precip_nc_nolake(27,:,:) - precip_nc_nolake(3,:,:) ;total non-convective precipitation - all kinds******* conprecipsum_nolk = precip_con_nolake(27,:,:) - precip_con_nolake(3,:,:) ;total convective precipitation - all kinds******* tot_swe_nolk = swe_nolake(27,:,:) - swe_nolake(3,:,:) ;Daily change in SWE******* tot_snowdepth_nolk = snow_depth_nolake(27,:,:) - snow_depth_nolake(3,:,:) ;Daily change in Snow Depth******* precipaccum_nolk = precipaccum_nolk + precipsum_nolk ;Calculate New Daily Accumulated Precipitation******* precipsum_diff = precipsum_nolk - precipsum_lk ;diff in total precipitation******* ncprecipsum_diff = ncprecipsum_nolk - ncprecipsum_lk ;diff in non-convective precipitation******* conprecipsum_diff = conprecipsum_nolk - conprecipsum_lk ;diff in convective precipitation******* swe_diff = tot_swe_nolk - tot_swe_lk ;diff in swe******* snowdepth_diff = tot_snowdepth_nolk - tot_snowdepth_lk ;diff in snow depth******* precipaccum_diff = precipaccum_nolk - precipaccum_lk ;diff in accumulated total precipitation******* ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; Planetary Boundary Layer Height LAKE pblh_lake = wrf_user_getvar(a,"PBLH",-1) pblh_lk = pblh_lake(4:27,:,:) min_pblh_lk = dim_min_n(pblh_lk,0) max_pblh_lk = dim_max_n(pblh_lk,0) avg_pblh_lk = dim_avg_n(pblh_lk,0) ; Planetary Boundary Layer Height NO LAKE pblh_nolake = wrf_user_getvar(b,"PBLH",-1) pblh_nolk = pblh_nolake(4:27,:,:) min_pblh_nolk = dim_min_n(pblh_nolk,0) max_pblh_nolk = dim_max_n(pblh_nolk,0) avg_pblh_nolk = dim_avg_n(pblh_nolk,0) pblh_diff = pblh_nolk - pblh_lk ;diff in PBL Height (all times)******* min_pblh_diff = min_pblh_nolk - min_pblh_lk ;diff in minimum PBL Height******* max_pblh_diff = max_pblh_nolk - max_pblh_lk ;diff in maximum PBL Height******* avg_pblh_diff = avg_pblh_nolk - avg_pblh_lk ;diff in average PBL Height******* ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; 2m Temperature LAKE t2_lake = wrf_user_getvar(a,"T2",-1) - 273.15 t2_lk = t2_lake(4:27,:,:) min_t2_lk = dim_min_n(t2_lk,0) max_t2_lk = dim_max_n(t2_lk,0) avg_t2_lk = dim_avg_n(t2_lk,0) ; 2m Temperature NO LAKE t2_nolake = wrf_user_getvar(b,"T2",-1) - 273.15 t2_nolk = t2_nolake(4:27,:,:) min_t2_nolk = dim_min_n(t2_nolk,0) max_t2_nolk = dim_max_n(t2_nolk,0) avg_t2_nolk = dim_avg_n(t2_nolk,0) min_t2_diff = min_t2_nolk - min_t2_lk ;diff in minimum 2m Temperature******* max_t2_diff = max_t2_nolk - max_t2_lk ;diff in maximum 2m Temperature******* avg_t2_diff = avg_t2_nolk - avg_t2_lk ;diff in average 2m Temperature******* avg_t2_diff@lon2d = lon2d avg_t2_diff@lat2d = lat2d ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; Surface Pressure LAKE psfc_lake = wrf_user_getvar(a,"PSFC",-1) psfc_lk = psfc_lake(4:27,:,:) min_psfc_lk = dim_min_n(psfc_lk,0) max_psfc_lk = dim_max_n(psfc_lk,0) avg_psfc_lk = dim_avg_n(psfc_lk,0) ; Surface Pressure NO LAKE psfc_nolake = wrf_user_getvar(b,"PSFC",-1) psfc_nolk = psfc_nolake(4:27,:,:) min_psfc_nolk = dim_min_n(psfc_nolk,0) max_psfc_nolk = dim_max_n(psfc_nolk,0) avg_psfc_nolk = dim_avg_n(psfc_nolk,0) min_psfc_diff = min_psfc_nolk - min_psfc_lk ;diff in minimum Surface Pressure******* max_psfc_diff = max_psfc_nolk - max_psfc_lk ;diff in maximum Surface Pressure******* avg_psfc_diff = avg_psfc_nolk - avg_psfc_lk ;diff in average Surface Pressure******* ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; 2m Mixing Ratio, RH, and Td LAKE q2_lake = wrf_user_getvar(a,"Q2",-1)*1000 q2_lk = q2_lake(4:27,:,:) rh2_lk = wrf_rh(q2_lk/1000, psfc_lk, t2_lk+273.15) td2_lk = wrf_td(psfc_lk, q2_lk/1000) min_q2_lk = dim_min_n(q2_lk,0) max_q2_lk = dim_max_n(q2_lk,0) avg_q2_lk = dim_avg_n(q2_lk,0) min_rh2_lk = dim_min_n(rh2_lk,0) max_rh2_lk = dim_max_n(rh2_lk,0) avg_rh2_lk = dim_avg_n(rh2_lk,0) min_td2_lk = dim_min_n(td2_lk,0) max_td2_lk = dim_max_n(td2_lk,0) avg_td2_lk = dim_avg_n(td2_lk,0) ; 2m Mixing Ratio, RH, and Td NO LAKE q2_nolake = wrf_user_getvar(b,"Q2",-1)*1000 q2_nolk = q2_nolake(4:27,:,:) rh2_nolk = wrf_rh(q2_nolk/1000, psfc_nolk, t2_nolk+273.15) td2_nolk = wrf_td(psfc_nolk, q2_nolk/1000) min_q2_nolk = dim_min_n(q2_nolk,0) max_q2_nolk = dim_max_n(q2_nolk,0) avg_q2_nolk = dim_avg_n(q2_nolk,0) min_rh2_nolk = dim_min_n(rh2_nolk,0) max_rh2_nolk = dim_max_n(rh2_nolk,0) avg_rh2_nolk = dim_avg_n(rh2_nolk,0) min_td2_nolk = dim_min_n(td2_nolk,0) max_td2_nolk = dim_max_n(td2_nolk,0) avg_td2_nolk = dim_avg_n(td2_nolk,0) min_q2_diff = min_q2_nolk - min_q2_lk ;diff in minimum Mixing Ratio******* max_q2_diff = max_q2_nolk - max_q2_lk ;diff in maximum Mixing Ratio******* avg_q2_diff = avg_q2_nolk - avg_q2_lk ;diff in average Mixing Ratio******* min_rh2_diff = min_rh2_nolk - min_rh2_lk ;diff in minimum RH******* max_rh2_diff = max_rh2_nolk - max_rh2_lk ;diff in maximum RH******* avg_rh2_diff = avg_rh2_nolk - avg_rh2_lk ;diff in average RH******* min_td2_diff = min_td2_nolk - min_td2_lk ;diff in minimum Td******* max_td2_diff = max_td2_nolk - max_td2_lk ;diff in maximum Td******* avg_td2_diff = avg_td2_nolk - avg_td2_lk;diff in average Td******* ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ; SURFACE LAKE WIND u10_lake = wrf_user_getvar(a,"U10",-1) u10_lk = u10_lake(4:27,:,:) v10_lake = wrf_user_getvar(a,"V10",-1) v10_lk = v10_lake(4:27,:,:) wspd10_lk := wind_speed(u10_lk,v10_lk) wdir10_lk := wind_direction(u10_lk,v10_lk,0) min_u10_lk = dim_min_n(u10_lk,0) max_u10_lk = dim_max_n(u10_lk,0) min_v10_lk = dim_min_n(v10_lk,0) max_v10_lk = dim_max_n(v10_lk,0) min_ws10_lk = dim_min_n(wspd10_lk,0) max_ws10_lk = dim_max_n(wspd10_lk,0) wind10_lk := wind_stats(wspd10_lk, wdir10_lk, 0, False) ; variable of type 'list' wspd10_lk_avg = wind10_lk[0] ; average of 'wspd' wspd10_lk_std = wind10_lk[1] ; standard deviation of 'wspd' wdir10_lk_avg = wind10_lk[2] ; mean vector wind direction wdir10_lk_std = wind10_lk[3] ; standard deviation of the wind direction avgU10_lk = wind10_lk[4] ; average zonal component avgV10_lk = wind10_lk[5] ; average medidional wind component ; SURFACE NO LAKE WIND u10_nolake = wrf_user_getvar(b,"U10",-1) u10_nolk = u10_nolake(4:27,:,:) v10_nolake = wrf_user_getvar(b,"V10",-1) v10_nolk = v10_nolake(4:27,:,:) wspd10_nolk := wind_speed(u10_nolk,v10_nolk) wdir10_nolk := wind_direction(u10_nolk,v10_nolk,0) min_u10_nolk = dim_min_n(u10_nolk,0) max_u10_nolk = dim_max_n(u10_nolk,0) min_v10_nolk = dim_min_n(v10_nolk,0) max_v10_nolk = dim_max_n(v10_nolk,0) min_ws10_nolk = dim_min_n(wspd10_nolk,0) max_ws10_nolk = dim_max_n(wspd10_nolk,0) wind10_nolk := wind_stats(wspd10_nolk, wdir10_nolk, 0, False) ; variable of type 'list' wspd10_nolk_avg = wind10_nolk[0] ; average of 'wspd' wspd10_nolk_std = wind10_nolk[1] ; standard deviation of 'wspd' wdir10_nolk_avg = wind10_nolk[2] ; mean vector wind direction wdir10_nolk_std = wind10_nolk[3] ; standard deviation of the wind direction avgU10_nolk = wind10_nolk[4] ; average zonal component avgV10_nolk = wind10_nolk[5] ; average medidional wind component min_u10_diff = min_u10_nolk - min_u10_lk ;diff in minimum U10******* max_u10_diff = max_u10_nolk - max_u10_lk ;diff in maximum U10******* avg_u10_diff = avgU10_nolk - avgU10_lk ;diff in average U10******* avg_u10_diff@lon2d = lon2d avg_u10_diff@lat2d = lat2d min_v10_diff = min_v10_nolk - min_v10_lk ;diff in minimum V10******* max_v10_diff = max_v10_nolk - max_v10_lk ;diff in maximum V10******* avg_v10_diff = avgV10_nolk - avgV10_lk ;diff in average V10******* avg_v10_diff@lon2d = lon2d avg_v10_diff@lat2d = lat2d min_ws10_diff = min_ws10_nolk - min_ws10_lk ;diff in minimum WS10******* max_ws10_diff = max_ws10_nolk - max_ws10_lk ;diff in maximum WS10******* avg_ws10_diff = wspd10_nolk_avg - wspd10_lk_avg ;diff in average WS10******* wspd10_std_diff = wspd10_nolk_std - wspd10_lk_std; difference in wind speed standard deviation wdir10_std_diff = wdir10_nolk_std - wdir10_lk_std; difference in wind direction standard deviation ;700 and 800 mb LAKE WIND U_lake = wrf_user_getvar(a,"U",-1) U_lakeunstagger = wrf_user_unstagger(U_lake,U_lake@stagger) U_lk = U_lakeunstagger(4:27,:,:,:) V_lake = wrf_user_getvar(a,"V",-1) V_lakeunstagger = wrf_user_unstagger(V_lake,V_lake@stagger) V_lk = V_lakeunstagger(4:27,:,:,:) wspd_lk := wind_speed(U_lk,V_lk) wdir_lk := wind_direction(U_lk,V_lk,0) p_lk = wrf_user_getvar(a,"P",-1) pb_lk = wrf_user_getvar(a,"PB",-1) pf_lk = (p_lk(4:27,:,:,:) + pb_lk(4:27,:,:,:))*.01 ; Convert to hPa avg_U_700_lk = dim_avg_n(wrf_interp_3d_z(U_lk,pf_lk,700.),0) avg_V_700_lk = dim_avg_n(wrf_interp_3d_z(V_lk,pf_lk,700.),0) avg_wspd_700_lk = dim_avg_n(wrf_interp_3d_z(wspd_lk,pf_lk,700.),0) avg_wdir_700_lk = dim_avg_n(wrf_interp_3d_z(wdir_lk,pf_lk,700.),0) avg_U_800_lk = dim_avg_n(wrf_interp_3d_z(U_lk,pf_lk,800.),0) avg_V_800_lk = dim_avg_n(wrf_interp_3d_z(V_lk,pf_lk,800.),0) avg_wspd_800_lk = dim_avg_n(wrf_interp_3d_z(wspd_lk,pf_lk,800.),0) avg_wdir_800_lk = dim_avg_n(wrf_interp_3d_z(wdir_lk,pf_lk,800.),0) ;700 and 800 mb NO LAKE WIND U_nolake = wrf_user_getvar(b,"U",-1) U_nolakeunstagger = wrf_user_unstagger(U_nolake,U_nolake@stagger) U_nolk = U_nolakeunstagger(4:27,:,:,:) V_nolake = wrf_user_getvar(b,"V",-1) V_nolakeunstagger = wrf_user_unstagger(V_nolake,V_nolake@stagger) V_nolk = V_nolakeunstagger(4:27,:,:,:) wspd_nolk := wind_speed(U_nolk,V_nolk) wdir_nolk := wind_direction(U_nolk,V_nolk,0) p_nolk = wrf_user_getvar(b,"P",-1) pb_nolk = wrf_user_getvar(b,"PB",-1) pf_nolk = (p_nolk(4:27,:,:,:) + pb_nolk(4:27,:,:,:))*.01 ; Convert to hPa avg_U_700_nolk = dim_avg_n(wrf_interp_3d_z(U_nolk,pf_nolk,700.),0) avg_V_700_nolk = dim_avg_n(wrf_interp_3d_z(V_nolk,pf_nolk,700.),0) avg_wspd_700_nolk = dim_avg_n(wrf_interp_3d_z(wspd_nolk,pf_nolk,700.),0) avg_wdir_700_nolk = dim_avg_n(wrf_interp_3d_z(wdir_nolk,pf_nolk,700.),0) avg_U_800_nolk = dim_avg_n(wrf_interp_3d_z(U_nolk,pf_nolk,800.),0) avg_V_800_nolk = dim_avg_n(wrf_interp_3d_z(V_nolk,pf_nolk,800.),0) avg_wspd_800_nolk = dim_avg_n(wrf_interp_3d_z(wspd_nolk,pf_nolk,800.),0) avg_wdir_800_nolk = dim_avg_n(wrf_interp_3d_z(wdir_nolk,pf_nolk,800.),0) avg_U_700_diff = avg_U_700_nolk - avg_U_700_lk avg_U_800_diff = avg_U_800_nolk - avg_U_800_lk avg_V_700_diff = avg_V_700_nolk - avg_V_700_lk avg_V_800_diff = avg_V_800_nolk - avg_V_800_lk avg_wspd_700_diff = avg_wspd_700_nolk - avg_wspd_700_lk avg_wspd_800_diff = avg_wspd_800_nolk - avg_wspd_800_lk print("writing data") f->time = (/dates/) f->lat = (/latitude/) f->lon = (/longitude/) f->albedo = (/albedo/) f->terrain_height = (/terrain_height/) f->soil_type = (/soil_type/) f->veg_type = (/veg_type/) f->lai = (/lai/) f->lulc = (/lulc/) f->soil_moisture = (/soil_moisture/) f->sfc_skin_temp = (/sfc_skin_temp/) f->soil_temp = (/soil_temp/) f->veg_frac = (/veg_frac/) f->precipsum_lk(n,:,:) = (/precipsum_lk/) f->precipaccum_lk(n,:,:) = (/precipaccum_lk/) f->ncprecipsum_lk(n,:,:) = (/ncprecipsum_lk/) f->conprecipsum_lk(n,:,:) = (/conprecipsum_lk/) f->tot_swe_lk(n,:,:) = (/tot_swe_lk/) f->tot_snowdepth_lk(n,:,:) = (/tot_snowdepth_lk/) f->precipsum_nolk(n,:,:) = (/precipsum_nolk/) f->precipaccum_nolk(n,:,:) = (/precipaccum_nolk/) f->ncprecipsum_nolk(n,:,:) = (/ncprecipsum_nolk/) f->conprecipsum_nolk(n,:,:) = (/conprecipsum_nolk/) f->tot_swe_nolk(n,:,:) = (/tot_swe_nolk/) f->tot_snowdepth_nolk(n,:,:) = (/tot_snowdepth_nolk/) f->precipsum_diff(n,:,:) = (/precipsum_diff/) f->precipaccum_diff(n,:,:) = (/precipaccum_diff/) f->ncprecipsum_diff(n,:,:) = (/ncprecipsum_diff/) f->conprecipsum_diff(n,:,:) = (/conprecipsum_diff/) f->swe_diff(n,:,:) = (/swe_diff/) f->snowdepth_diff(n,:,:) = (/snowdepth_diff/) f->min_pblh_lk(n,:,:) = (/min_pblh_lk/) f->max_pblh_lk(n,:,:) = (/max_pblh_lk/) f->avg_pblh_lk(n,:,:) = (/avg_pblh_lk/) f->min_pblh_nolk(n,:,:) = (/min_pblh_nolk/) f->max_pblh_nolk(n,:,:) = (/max_pblh_nolk/) f->avg_pblh_nolk(n,:,:) = (/avg_pblh_nolk/) f->min_pblh_diff(n,:,:) = (/min_pblh_diff/) f->max_pblh_diff(n,:,:) = (/max_pblh_diff/) f->avg_pblh_diff(n,:,:) = (/avg_pblh_diff/) f->min_t2_lk(n,:,:) = (/min_t2_lk/) f->max_t2_lk(n,:,:) = (/max_t2_lk/) f->avg_t2_lk(n,:,:) = (/avg_t2_lk/) f->min_t2_nolk(n,:,:) = (/min_t2_nolk/) f->max_t2_nolk(n,:,:) = (/max_t2_nolk/) f->avg_t2_nolk(n,:,:) = (/avg_t2_nolk/) f->min_t2_diff(n,:,:) = (/min_t2_diff/) f->max_t2_diff(n,:,:) = (/max_t2_diff/) f->avg_t2_diff(n,:,:) = (/avg_t2_diff/) f->min_psfc_lk(n,:,:) = (/min_psfc_lk/) f->max_psfc_lk(n,:,:) = (/max_psfc_lk/) f->avg_psfc_lk(n,:,:) = (/avg_psfc_lk/) f->min_psfc_nolk(n,:,:) = (/min_psfc_nolk/) f->max_psfc_nolk(n,:,:) = (/max_psfc_nolk/) f->avg_psfc_nolk(n,:,:) = (/avg_psfc_nolk/) f->min_psfc_diff(n,:,:) = (/min_psfc_diff/) f->max_psfc_diff(n,:,:) = (/max_psfc_diff/) f->avg_psfc_diff(n,:,:) = (/avg_psfc_diff/) f->min_q2_lk(n,:,:) = (/min_q2_lk/) f->max_q2_lk(n,:,:) = (/max_q2_lk/) f->avg_q2_lk(n,:,:) = (/avg_q2_lk/) f->min_q2_nolk(n,:,:) = (/min_q2_nolk/) f->max_q2_nolk(n,:,:) = (/max_q2_nolk/) f->avg_q2_nolk(n,:,:) = (/avg_q2_nolk/) f->min_q2_diff(n,:,:) = (/min_q2_diff/) f->max_q2_diff(n,:,:) = (/max_q2_diff/) f->avg_q2_diff(n,:,:) = (/avg_q2_diff/) f->min_td2_lk(n,:,:) = (/min_td2_lk/) f->max_td2_lk(n,:,:) = (/max_td2_lk/) f->avg_td2_lk(n,:,:) = (/avg_td2_lk/) f->min_td2_nolk(n,:,:) = (/min_td2_nolk/) f->max_td2_nolk(n,:,:) = (/max_td2_nolk/) f->avg_td2_nolk(n,:,:) = (/avg_td2_nolk/) f->min_td2_diff(n,:,:) = (/min_td2_diff/) f->max_td2_diff(n,:,:) = (/max_td2_diff/) f->avg_td2_diff(n,:,:) = (/avg_td2_diff/) f->min_rh2_lk(n,:,:) = (/min_rh2_lk/) f->max_rh2_lk(n,:,:) = (/max_rh2_lk/) f->avg_rh2_lk(n,:,:) = (/avg_rh2_lk/) f->min_rh2_nolk(n,:,:) = (/min_rh2_nolk/) f->max_rh2_nolk(n,:,:) = (/max_rh2_nolk/) f->avg_rh2_nolk(n,:,:) = (/avg_td2_nolk/) f->min_rh2_diff(n,:,:) = (/min_rh2_diff/) f->max_rh2_diff(n,:,:) = (/max_rh2_diff/) f->avg_rh2_diff(n,:,:) = (/avg_rh2_diff/) f->min_u10_lk(n,:,:) = (/min_u10_lk/) f->max_u10_lk(n,:,:) = (/max_u10_lk/) f->min_v10_lk(n,:,:) = (/min_v10_lk/) f->max_v10_lk(n,:,:) = (/max_v10_lk/) f->min_ws10_lk(n,:,:) = (/min_ws10_lk/) f->max_ws10_lk(n,:,:) = (/max_ws10_lk/) f->wspd10_lk_avg(n,:,:) = (/wspd10_lk_avg/) f->wspd10_lk_std(n,:,:) = (/wspd10_lk_std/) f->wdir10_lk_avg(n,:,:) = (/wdir10_lk_avg/) f->wdir10_lk_std(n,:,:) = (/wdir10_lk_std/) f->avgU10_lk(n,:,:) = (/avgU10_lk/) f->avgV10_lk(n,:,:) = (/avgV10_lk/) f->min_u10_nolk(n,:,:) = (/min_u10_nolk/) f->max_u10_nolk(n,:,:) = (/max_u10_nolk/) f->min_v10_nolk(n,:,:) = (/min_v10_nolk/) f->max_v10_nolk(n,:,:) = (/max_v10_nolk/) f->min_ws10_nolk(n,:,:) = (/min_ws10_nolk/) f->max_ws10_nolk(n,:,:) = (/max_ws10_nolk/) f->wspd10_nolk_avg(n,:,:) = (/wspd10_nolk_avg/) f->wspd10_nolk_std(n,:,:) = (/wspd10_nolk_std/) f->wdir10_nolk_avg(n,:,:) = (/wdir10_nolk_avg/) f->wdir10_nolk_std(n,:,:) = (/wdir10_nolk_std/) f->avgU10_nolk(n,:,:) = (/avgU10_nolk/) f->avgV10_nolk(n,:,:) = (/avgV10_nolk/) f->min_u10_diff(n,:,:) = (/min_u10_diff/) f->max_u10_diff(n,:,:) = (/max_u10_diff/) f->avg_u10_diff(n,:,:) = (/avg_u10_diff/) f->min_v10_diff(n,:,:) = (/min_v10_diff/) f->max_v10_diff(n,:,:) = (/max_v10_diff/) f->avg_v10_diff(n,:,:) = (/avg_v10_diff/) f->min_ws10_diff(n,:,:) = (/min_ws10_diff/) f->max_ws10_diff(n,:,:) = (/max_ws10_diff/) f->avg_ws10_diff(n,:,:) = (/avg_ws10_diff/) f->wspd10_std_diff(n,:,:) = (/wspd10_std_diff/) f->wdir10_std_diff(n,:,:) = (/wdir10_std_diff/) f->avg_U_700_lk(n,:,:) = (/avg_U_700_lk/) f->avg_V_700_lk(n,:,:) = (/avg_V_700_lk/) f->avg_wspd_700_lk(n,:,:) = (/avg_wspd_700_lk/) f->avg_wdir_700_lk(n,:,:) = (/avg_wdir_700_lk/) f->avg_U_800_lk(n,:,:) = (/avg_U_800_lk/) f->avg_V_800_lk(n,:,:) = (/avg_V_800_lk/) f->avg_wspd_800_lk(n,:,:) = (/avg_wspd_800_lk/) f->avg_wdir_800_lk(n,:,:) = (/avg_wdir_800_lk/) f->avg_U_700_nolk(n,:,:) = (/avg_U_700_nolk/) f->avg_V_700_nolk(n,:,:) = (/avg_V_700_nolk/) f->avg_wspd_700_nolk(n,:,:) = (/avg_wspd_700_nolk/) f->avg_wdir_700_nolk(n,:,:) = (/avg_wdir_700_nolk/) f->avg_U_800_nolk(n,:,:) = (/avg_U_800_nolk/) f->avg_V_800_nolk(n,:,:) = (/avg_V_800_nolk/) f->avg_wspd_800_nolk(n,:,:) = (/avg_wspd_800_nolk/) f->avg_wdir_800_nolk(n,:,:) = (/avg_wdir_800_nolk/) f->avg_U_700_diff(n,:,:) = (/avg_U_700_diff/) f->avg_U_800_diff(n,:,:) = (/avg_U_800_diff/) f->avg_V_700_diff(n,:,:) = (/avg_V_700_diff/) f->avg_V_800_diff(n,:,:) = (/avg_V_800_diff/) f->avg_wspd_700_diff(n,:,:) = (/avg_wspd_700_diff/) f->avg_wspd_800_diff(n,:,:) = (/avg_wspd_800_diff/) end do end