load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl" load "/scratch/y98/shared/ncl-functions/cordex_grid_plot_functions.ncl" begin x = 125 y = 160 domain = "d02" files_era = systemfunc("ls /scratch/y98/julia/WRF-CLIM/Era-Int/Processed/ERA_" + domain + "_clim_means.nc") files_bom = systemfunc("ls /group/y98/julia/BOM/Processed/BOM_clim_means_" + domain + ".nc") f_era = addfile(files_era,"r") f_bom = addfile(files_bom,"r") tmax_era = f_era->TMAX_mean_season tmin_era = f_era->TMIN_mean_season ppt_era = f_era->ppt_mean_season TXx_era = f_era->TXx_avg TXn_era = f_era->TXn_avg TNx_era = f_era->TNx_avg TNn_era = f_era->TNn_avg DTR_era = f_era->DTR_avg RX1d_era = f_era->RX1d_avg FD_era = f_era->FD SU_era = f_era->SU TR_era = f_era->TR R10mm_era = f_era->R10mm R20mm_era = f_era->R20mm tmax_bom = f_bom->TMAX_mean_season tmin_bom = f_bom->TMIN_mean_season ppt_bom = f_bom->ppt_mean_season TXx_bom = f_bom->TXx_avg TXn_bom = f_bom->TXn_avg TNx_bom = f_bom->TNx_avg TNn_bom = f_bom->TNn_avg DTR_bom = f_bom->DTR_avg RX1d_bom = f_bom->RX1d_avg FD_bom = f_bom->FD SU_bom = f_bom->SU TR_bom = f_bom->TR R10mm_bom = f_bom->R10mm R20mm_bom = f_bom->R20mm RX1d_annual_bom = dim_max_n(RX1d_bom,0) RX1d_annual_era = dim_max_n(RX1d_era,0) ; ---- calculate bias: model - obs --------- tmax_nomask = tmax_era - tmax_bom tmin_nomask = tmin_era - tmin_bom ppt_nomask = ppt_era - ppt_bom TXx_nomask = TXx_era - TXx_bom TXn_nomask = TXn_era - TXn_bom TNx_nomask = TNx_era - TNx_bom TNn_nomask = TNn_era - TNn_bom RX1d_nomask = (RX1d_era - RX1d_bom)/RX1d_bom*100 RX1d_annual_nomask = (RX1d_annual_era - RX1d_annual_bom)/RX1d_annual_bom*100 DTR_nomask = (DTR_era - DTR_bom)/DTR_bom*100 FD_nomask = FD_era SU_nomask = SU_era TR_nomask = TR_era R10mm_nomask = (R10mm_era - R10mm_bom)/R10mm_bom*100 R20mm_nomask = (R20mm_era - R20mm_bom)/R20mm_bom*100 delete(f_bom) ;get a landmask and some other stuff from a wrfout file ; load a wrfout file to get land-mask info out f = addfile("/scratch/y98/julia/WRF-CLIM/Era-Int/wrfout_" + domain + ".nc","r") landmask = wrf_user_getvar(f,"LANDMASK",0) ;dsizes_array = dimsizes(land_mask) ; also get some map info for later on truelat1 = f@TRUELAT1 truelat2 = f@TRUELAT2 cen_lon = f@CEN_LON cen_lat = f@CEN_LAT wrf_lat2d = f->XLAT(0,:,:) wrf_lon2d = f->XLONG(0,:,:) latr = wrf_lat2d lonr = wrf_lon2d ;mask water and remove relaxation zone dsizes_array = dimsizes(landmask) print(dsizes_array) landmask_nr = landmask(10:(dsizes_array(0)-11),10:(dsizes_array(1)-11)) lon = lonr(10:(dsizes_array(0)-11),10:(dsizes_array(1)-11)) lat = latr(10:(dsizes_array(0)-11),10:(dsizes_array(1)-11)) dsizes_landmask_nr = dimsizes(landmask_nr) print(dsizes_landmask_nr) tmax = mask(tmax_nomask,landmask_nr,1) tmin = mask(tmin_nomask,landmask_nr,1) ppt = mask(ppt_nomask,landmask_nr,1) FD = mask(FD_nomask,landmask_nr,1) SU = mask(SU_nomask,landmask_nr,1) dsizes_SU = dimsizes(SU) print(dsizes_SU) TR = mask(TR_nomask,landmask_nr,1) DTR = mask(DTR_nomask,landmask_nr,1) RX1d = mask(RX1d_nomask,landmask_nr,1) RX1d_annual = mask(RX1d_annual_nomask,landmask_nr,1) R10mm = mask(R10mm_nomask,landmask_nr,1) R20mm = mask(R20mm_nomask,landmask_nr,1) TXx = mask(TXx_nomask,landmask_nr,1) TXn = mask(TXn_nomask,landmask_nr,1) TNx = mask(TNx_nomask,landmask_nr,1) TNn = mask(TNn_nomask,landmask_nr,1) ; open a d03 file to get the lat long of the d03 boundary and apply this to the d02 data to get ij points for our contour corners f_d03 = addfile("/scratch/y98/julia/WRF-CLIM/Era-Int/wrfout_d03.nc","r") f_d02 = addfile("/scratch/y98/julia/WRF-CLIM/Era-Int/wrfout_d02.nc","r") wrf_lat2d_d03 = f_d03->XLAT(0,:,:) wrf_lon2d_d03 = f_d03->XLONG(0,:,:) ; get lat,lon of corners corner_lat_sw = wrf_lat2d_d03(9,9) print(corner_lat_sw) corner_lat_nw = wrf_lat2d_d03(9,99) corner_lat_ne = wrf_lat2d_d03(99,99) corner_lat_se = wrf_lat2d_d03(99,9) corner_lon_sw = wrf_lon2d_d03(9,9) print(corner_lon_sw) corner_lon_nw = wrf_lon2d_d03(9,99) corner_lon_ne = wrf_lon2d_d03(99,99) corner_lon_se = wrf_lon2d_d03(99,9) ; now work out the i,j points in the d02 domain that correspond to these lat long values opt = True loc_sw = wrf_user_ll_to_ij(f_d02,corner_lon_sw,corner_lat_sw,opt) loc_nw = wrf_user_ll_to_ij(f_d02,corner_lon_nw,corner_lat_nw,opt) loc_ne = wrf_user_ll_to_ij(f_d02,corner_lon_ne,corner_lat_ne,opt) loc_se = wrf_user_ll_to_ij(f_d02,corner_lon_se,corner_lat_se,opt) ; now we are converting to ncl subscripts and accounting fo the fact that the d02 data has already had the relaxation zone removed. loc_sw = loc_sw - 11 loc_nw = loc_nw - 11 loc_ne = loc_ne - 11 loc_se = loc_se - 11 print(loc_sw) print(loc_nw) print(loc_ne) print(loc_se) delete(f_d03) delete(wrf_lat2d_d03) delete(wrf_lon2d_d03) delete(corner_lat_sw) delete(corner_lat_nw) delete(corner_lat_ne) delete(corner_lat_se) delete(corner_lon_sw) delete(corner_lon_nw) delete(corner_lon_ne) delete(corner_lon_se) delete(loc_sw) delete(loc_nw) delete(loc_ne) delete(loc_se) ; now that we have our ij locations for overlaying the d02 grid onto the d03 grid, we are going to limit our variables to within these corners tmax_limit = tmax(:,45:94,71:120) tmin_limit = tmin(:,45:94,71:120) ppt_limit = ppt(:,45:94,71:120) FD_limit = FD(45:94,71:120) SU_limit = SU(45:94,71:120) TR_limit = TR(45:94,71:120) R10mm_limit = R10mm(45:94,71:120) R20mm_limit = R20mm(45:94,71:120) DTR_limit = DTR(:,45:94,71:120) RX1d_limit = RX1d(:,45:94,71:120) RX1d_annual_limit = RX1d_annual(45:94,71:120) TXx_limit = TXx(:,45:94,71:120) TXn_limit = TXn(:,45:94,71:120) TNx_limit = TNx(:,45:94,71:120) TNn_limit = TNn(:,45:94,71:120) lat_limit = lat(45:94,71:120) lon_limit = lon(45:94,71:120) dimsizes_FD_limit = dimsizes(FD_limit) print(dimsizes_FD_limit) ;*******Draw Indices Panel Plots************ ;fi = new((/1,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) ; ;fi(0,:,:) = FD_limit ;;fi(1,:,:) = SU_limit ;;fi(2,:,:) = TR_limit ;;fi(3,:,:) = R10mm_limit ;;fi(4,:,:) = R20mm_limit ; ;type = "pdf" ; choose x11 for pop-up window, or pdf, or eps ;file_save = "Plot/FD_count__" + domain ;tiMainString = (/"D02: Frost Days"/) ;tiMainString = (/"Frost Days","Summer Days","Tropical Nights","R10mm","R20mm"/) ;panel_x = 1 ;panel_y = 1 ;panel_title = "Bias: Annual Count of Frost Days" ;lbTitleString = "Bias (Days)" ;cnMinLevelValF = 0 ;cnMaxLevelValF = 30 ;cnLevelSpacingF = 1 ;plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) ; ;delete(fi) ;delete(file_save) ;delete(panel_title) ;delete(tiMainString) ;delete(plot) ;delete(cnMinLevelValF) ;delete(cnMaxLevelValF) ;delete(cnLevelSpacingF) ; ;fi = new((/1,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) ; ;fi(0,:,:) = SU_limit ;;fi(1,:,:) = SU_limit ;;fi(2,:,:) = TR_limit ;;fi(3,:,:) = R10mm_limit ;;fi(4,:,:) = R20mm_limit ; ;type = "pdf" ; choose x11 for pop-up window, or pdf, or eps ;file_save = "Plot/SU_count_" + domain ;tiMainString = (/"D0: Summer Days"/) ;;tiMainString = (/"Frost Days","Summer Days","Tropical Nights","R10mm","R20mm"/) ;panel_x = 1 ;panel_y = 1 ;panel_title = "Bias: Annual Count of Sumer Days" ;lbTitleString = "Bias (Days)" ;cnMinLevelValF = 0 ;cnMaxLevelValF = 50 ;cnLevelSpacingF = 10 ;plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) ; ;delete(fi) ;delete(file_save) ;delete(panel_title) ;delete(tiMainString) ;delete(plot) ;delete(cnMinLevelValF) ;delete(cnMaxLevelValF) ;delete(cnLevelSpacingF) ; ; ;fi = new((/1,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) ; ;fi(0,:,:) = TR_limit ;;fi(1,:,:) = SU_limit ;;fi(2,:,:) = TR_limit ;;fi(3,:,:) = R10mm_limit ;;fi(4,:,:) = R20mm_limit ; ;type = "pdf" ; choose x11 for pop-up window, or pdf, or eps ;file_save = "Plot/Bias_TR_count_Indices_" + domain ;tiMainString = (/"D02 Bias: Tropical Nights"/) ;;tiMainString = (/"Frost Days","Summer Days","Tropical Nights","R10mm","R20mm"/) ;panel_x = 1 ;panel_y = 1 ;panel_title = "Bias: Annual Count of Tropical Nights" ;lbTitleString = "Bias (Days)" ;cnMinLevelValF = -500 ;cnMaxLevelValF = 500 ;cnLevelSpacingF = 25 ;plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) ; ;delete(fi) ;delete(file_save) ;delete(panel_title) ;delete(tiMainString) ;delete(plot) ;delete(cnMinLevelValF) ;delete(cnMaxLevelValF) ;delete(cnLevelSpacingF) ; ; ;fi = new((/2,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) ; ;fi(0,:,:) = R10mm_limit ;fi(1,:,:) = R20mm_limit ;;fi(2,:,:) = TR_limit ;;fi(3,:,:) = R10mm_limit ;;fi(4,:,:) = R20mm_limit ; ;type = "pdf" ; choose x11 for pop-up window, or pdf, or eps ;file_save = "Plot/Bias_Rxmm_count_Indices_" + domain ;tiMainString = (/"D02 Bias: R10mm","D02 Bias: R20mm"/) ;;tiMainString = (/"Frost Days","Summer Days","Tropical Nights","R10mm","R20mm"/) ;panel_x = 1 ;panel_y = 2 ;panel_title = "Bias: Annual Count of Rainfall over Threshold" ;lbTitleString = "Bias (Days)" ;cnMinLevelValF = -100 ;cnMaxLevelValF = 100 ;cnLevelSpacingF = 10 ;plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) ; ;delete(fi) ;delete(file_save) ;delete(panel_title) ;delete(tiMainString) ;delete(plot) ;delete(cnMinLevelValF) ;delete(cnMaxLevelValF) ;delete(cnLevelSpacingF) ; ; ; fi = new((/4,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = tmax_limit(0,:,:) fi(1,:,:) = tmax_limit(1,:,:) fi(2,:,:) = tmax_limit(2,:,:) fi(3,:,:) = tmax_limit(3,:,:) ;fi(4,:,:) = tmin_limit(0,:,:) ;fi(5,:,:) = tmin_limit(1,:,:) ;fi(6,:,:) = tmin_limit(2,:,:) ;fi(7,:,:) = tmin_limit(3,:,:) ;fi(8,:,:) = ppt_limit(0,:,:) ;fi(9,:,:) = ppt_limit(1,:,:) ;fi(10,:,:) = ppt_limit(2,:,:) ;fi(11,:,:) = ppt_limit(3,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Seasonal_tmax_" + domain txString = "(a) Maximum Temperature W10k" tiMainString = (/"DJF", "MAM", "JJA", "SON"/) panel_x = 1 panel_y = 4 panel_title = "Seasonal TMAX: Bias D02" lbTitleString = "(~S~o~N~C)" cnMinLevelValF = -3 cnMaxLevelValF = 3 cnLevelSpacingF = 1 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) delete(cnLevelSpacingF) fi = new((/4,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = tmin_limit(0,:,:) fi(1,:,:) = tmin_limit(1,:,:) fi(2,:,:) = tmin_limit(2,:,:) fi(3,:,:) = tmin_limit(3,:,:) ;fi(4,:,:) = tmin_limit(0,:,:) ;fi(5,:,:) = tmin_limit(1,:,:) ;fi(6,:,:) = tmin_limit(2,:,:) ;fi(7,:,:) = tmin_limit(3,:,:) ;fi(8,:,:) = ppt_limit(0,:,:) ;fi(9,:,:) = ppt_limit(1,:,:) ;fi(10,:,:) = ppt_limit(2,:,:) ;fi(11,:,:) = ppt_limit(3,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Seasonal_tmin_" + domain txString = "(c) Minimum Temperature W10k" tiMainString = (/"DJF", "MAM", "JJA", "SON"/) panel_x = 1 panel_y = 4 panel_title = "Seasonal TMIN: Bias D02" lbTitleString = "(~S~o~N~C)" cnMinLevelValF = -3 cnMaxLevelValF = 3 cnLevelSpacingF = 1 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) fi = new((/4,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = ppt_limit(0,:,:) fi(1,:,:) = ppt_limit(1,:,:) fi(2,:,:) = ppt_limit(2,:,:) fi(3,:,:) = ppt_limit(3,:,:) ;fi(4,:,:) = tmin_limit(0,:,:) ;fi(5,:,:) = tmin_limit(1,:,:) ;fi(6,:,:) = tmin_limit(2,:,:) ;fi(7,:,:) = tmin_limit(3,:,:) ;fi(8,:,:) = ppt_limit(0,:,:) ;fi(9,:,:) = ppt_limit(1,:,:) ;fi(10,:,:) = ppt_limit(2,:,:) ;fi(11,:,:) = ppt_limit(3,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Seasonal_ppt_" + domain txString = "(a) W10k" tiMainString = (/"DJF", "MAM", "JJA", "SON"/) panel_x = 1 panel_y = 4 panel_title = "Seasonal PPT: Bias D02" lbTitleString = "(mm month~S~-1~N~)" cnMinLevelValF = -90 cnMaxLevelValF = 90 cnLevelSpacingF = 10 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) delete(tiMainString) ;********************************************** ;************ DTR Monthly Plots ************** ;********************************************** fi = new((/12,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = DTR_limit(0,:,:) fi(1,:,:) = DTR_limit(1,:,:) fi(2,:,:) = DTR_limit(2,:,:) fi(3,:,:) = DTR_limit(3,:,:) fi(4,:,:) = DTR_limit(4,:,:) fi(5,:,:) = DTR_limit(5,:,:) fi(6,:,:) = DTR_limit(6,:,:) fi(7,:,:) = DTR_limit(7,:,:) fi(8,:,:) = DTR_limit(8,:,:) fi(9,:,:) = DTR_limit(9,:,:) fi(10,:,:) = DTR_limit(10,:,:) fi(11,:,:) = DTR_limit(11,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Daily_Temperature_Range_" + domain tiMainString = (/"January D02", "February D02", "March D02", "April D02", "May D02", "June D02", "July D02", "August D02", "September D02", "October D02", "November D02","December D02"/) ;tiMainString = (/"CSIRO-1971", "BOM-1971"/) panel_x = 3 panel_y = 4 panel_title = "Bias D02: DTR" lbTitleString = "Bias(C)" cnMinLevelValF = -100 cnMaxLevelValF = 100 cnLevelSpacingF = 10 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) delete(cnLevelSpacingF) ;********************************************** ;************ RX1d Monthly Plots ************** ;********************************************** fi = new((/1,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = RX1d_annual_limit(:,:) ;fi(1,:,:) = RX1d_limit(1,:,:) ;fi(2,:,:) = RX1d_limit(2,:,:) ;fi(3,:,:) = RX1d_limit(3,:,:) ;fi(4,:,:) = RX1d_limit(4,:,:) ;fi(5,:,:) = RX1d_limit(5,:,:) ;fi(6,:,:) = RX1d_limit(6,:,:) ;fi(7,:,:) = RX1d_limit(7,:,:) ;fi(8,:,:) = RX1d_limit(8,:,:) ;fi(9,:,:) = RX1d_limit(9,:,:) ;fi(10,:,:) = RX1d_limit(10,:,:) ;fi(11,:,:) = RX1d_limit(11,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Max_1d_Precip_Rate_" + domain ;tiMainString = (/"January D02", "February D02", "March D02", "April D02", "May D02", "June D02", "July D02", "August D02", "September D02", "October D02", "November D02","December D02"/) tiMainString = (/"D02 RX1d Bias"/) ;tiMainString = (/"CSIRO-1971", "BOM-1971"/) panel_x = 1 panel_y = 1 panel_title = "RX1d" lbTitleString = "Bias(mm)" cnMinLevelValF = -100 cnMaxLevelValF = 100 cnLevelSpacingF = 10 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) delete(cnLevelSpacingF) ;********************************************** ;************ TXx Monthly Plots ************** ;********************************************** fi = new((/12,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = TXx_limit(0,:,:) fi(1,:,:) = TXx_limit(1,:,:) fi(2,:,:) = TXx_limit(2,:,:) fi(3,:,:) = TXx_limit(3,:,:) fi(4,:,:) = TXx_limit(4,:,:) fi(5,:,:) = TXx_limit(5,:,:) fi(6,:,:) = TXx_limit(6,:,:) fi(7,:,:) = TXx_limit(7,:,:) fi(8,:,:) = TXx_limit(8,:,:) fi(9,:,:) = TXx_limit(9,:,:) fi(10,:,:) = TXx_limit(10,:,:) fi(11,:,:) = TXx_limit(11,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Highest_Max_Temp_" + domain tiMainString = (/"January D02", "February D02", "March D02", "April D02", "May D02", "June D02", "July D02", "August D02", "September D02", "October D02", "November D02","December D02"/) ;tiMainString = (/"CSIRO-1971", "BOM-1971"/) panel_x = 3 panel_y = 4 panel_title = "TXx" lbTitleString = "Bias(C)" cnMinLevelValF = -100 cnMaxLevelValF = 100 cnLevelSpacingF = 10 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) delete(cnLevelSpacingF) ;********************************************** ;************ TXn Monthly Plots ************** ;********************************************** fi = new((/12,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = TXn_limit(0,:,:) fi(1,:,:) = TXn_limit(1,:,:) fi(2,:,:) = TXn_limit(2,:,:) fi(3,:,:) = TXn_limit(3,:,:) fi(4,:,:) = TXn_limit(4,:,:) fi(5,:,:) = TXn_limit(5,:,:) fi(6,:,:) = TXn_limit(6,:,:) fi(7,:,:) = TXn_limit(7,:,:) fi(8,:,:) = TXn_limit(8,:,:) fi(9,:,:) = TXn_limit(9,:,:) fi(10,:,:) = TXn_limit(10,:,:) fi(11,:,:) = TXn_limit(11,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Lowest_Max_Temp_" + domain tiMainString = (/"January D02", "February D02", "March D02", "April D02", "May D02", "June D02", "July D02", "August D02", "September D02", "October D02", "November D02","December D02"/) ;tiMainString = (/"CSIRO-1971", "BOM-1971"/) panel_x = 3 panel_y = 4 panel_title = "TXn" lbTitleString = "Bias(C)" cnMinLevelValF = -100 cnMaxLevelValF = 100 cnLevelSpacingF = 10 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) delete(cnLevelSpacingF) ;********************************************** ;************ TNx Monthly Plots ************** ;********************************************** fi = new((/12,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = TNx_limit(0,:,:) fi(1,:,:) = TNx_limit(1,:,:) fi(2,:,:) = TNx_limit(2,:,:) fi(3,:,:) = TNx_limit(3,:,:) fi(4,:,:) = TNx_limit(4,:,:) fi(5,:,:) = TNx_limit(5,:,:) fi(6,:,:) = TNx_limit(6,:,:) fi(7,:,:) = TNx_limit(7,:,:) fi(8,:,:) = TNx_limit(8,:,:) fi(9,:,:) = TNx_limit(9,:,:) fi(10,:,:) = TNx_limit(10,:,:) fi(11,:,:) = TNx_limit(11,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Highest_Min_Temp_" + domain tiMainString = (/"January D02", "February D02", "March D02", "April D02", "May D02", "June D02", "July D02", "August D02", "September D02", "October D02", "November D02","December D02"/) ;tiMainString = (/"CSIRO-1971", "BOM-1971"/) panel_x = 3 panel_y = 4 panel_title = "TNx" lbTitleString = "Bias(C)" cnMinLevelValF = -100 cnMaxLevelValF = 100 cnLevelSpacingF = 10 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) delete(cnLevelSpacingF) ;********************************************** ;************ TNn Monthly Plots ************** ;********************************************** fi = new((/12,dimsizes_FD_limit(0),dimsizes_FD_limit(1)/),float) fi(0,:,:) = TNn_limit(0,:,:) fi(1,:,:) = TNn_limit(1,:,:) fi(2,:,:) = TNn_limit(2,:,:) fi(3,:,:) = TNn_limit(3,:,:) fi(4,:,:) = TNn_limit(4,:,:) fi(5,:,:) = TNn_limit(5,:,:) fi(6,:,:) = TNn_limit(6,:,:) fi(7,:,:) = TNn_limit(7,:,:) fi(8,:,:) = TNn_limit(8,:,:) fi(9,:,:) = TNn_limit(9,:,:) fi(10,:,:) = TNn_limit(10,:,:) fi(11,:,:) = TNn_limit(11,:,:) type = "pdf" ; choose x11 for pop-up window, or pdf, or eps file_save = "Plot/Bias_Lowest_Min_Temp_" + domain tiMainString = (/"January D02", "February D02", "March D02", "April D02", "May D02", "June D02", "July D02", "August D02", "September D02", "October D02", "November D02","December D02"/) ;tiMainString = (/"CSIRO-1971", "BOM-1971"/) panel_x = 3 panel_y = 4 panel_title = "TNn" lbTitleString = "Bias(C)" cnMinLevelValF = -100 cnMaxLevelValF = 100 cnLevelSpacingF = 10 plot = plot_cordex_panel_scale(cen_lat,cen_lon,type,file_save,fi,lat_limit,lon_limit,txString,tiMainString,lbTitleString,panel_x,panel_y,cnMinLevelValF,cnMaxLevelValF,cnLevelSpacingF) delete(fi) delete(file_save) delete(panel_title) delete(plot) delete(cnMinLevelValF) delete(cnMaxLevelValF) delete(cnLevelSpacingF) end