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/csm/contributed.ncl" load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl" load "taylor_diagram_lh.ncl" ; load "taylor_diagram.ncl" ; load "taylor_metrics_table.ncl" load "./taylor_stats.ncl" ;load "./taylor_diagram_cam.ncl" begin ; ============================================================== ; ENSEMBLES temperature and precipitation without bias correction ; ============================================================== ; yrStrt = 1961 ; modification from 1970 ; yrLast = 1990 ; modification from 1999 ; ============================================================== ; Open the file: Read only the user specified period first observations then model ; ============================================================== f = addfile("E-OBS_1961-1990_tas_sm.nc", "r") ; TO = f->tas(:,:,:) ;printVarSummary(TO) ; (time, lat,lon) ; ============================================================== ; compute desired regional seasonal climatology (4,nlat,mlon) ;============================================================model2 deg temp f1= addfile("CNRM_control_tas_sm.nc", "r") ; T1 = f1->tas(:,:,:) ;printVarSummary(T1) ; (time, lat,lon) ;============================================================model1 deg cam4temp f2= addfile("DMI_control_tas_sm.nc", "r") ; T2 = f2->tas(:,:,:) ;printVarSummary(T2) ;============================================================model1deg cam5 temp f3= addfile("ICTP_control_tas_sm.nc", "r") ; T3 = f3->tas(:,:,:) ;printVarSummary(T3) ; (time, lat,lon) ;============================================================model1deg cam5 temp f4= addfile("KNMI_control_tas_sm.nc", "r") ; T4 = f4->tas(:,:,:) ;printVarSummary(T4) ;============================================================model1deg cam5 temp f5= addfile("SMHI_control_tas_sm.nc", "r") ; T5 = f5->tas(:,:,:) ;printVarSummary(T5) ;==================================================================precip ;================================================================== f6 = addfile ("E-OBS_1961-1990_pr_sm.nc", "r") pm = f6->pr(:,:,:) ;printVarSummary(pm) ; (time, lat,lon) ;=================================================================== f7= addfile("CNRM_control_pr_sm.nc", "r") ; P1 = f7->pr(:,:,:) ;printVarSummary(P1) ; (time, lat,lon) ;============================================================ f8= addfile("DMI_control_pr_sm.nc", "r") ; P2 = f8->pr(:,:,:) ;printVarSummary(P2) ;============================================================ f9= addfile("ICTP_control_pr_sm.nc", "r") ; P3 = f9->pr(:,:,:) ;printVarSummary(P3) ; (time, lat,lon) ;============================================================ f10= addfile("KNMI_control_pr_sm.nc", "r") ; P4 = f10->pr(:,:,:) ;printVarSummary(P4) ;============================================================ f11= addfile("SMHI_control_pr_sm.nc", "r") ; P5 = f11->pr(:,:,:) ;printVarSummary(P5) ;================================================================== ;Taylor diagram calculations ;=================================================================== case = (/ "ALADIN", "HIRHAM" , "RegCM", "RACMO2", "RCA"/) nCase = dimsizes(case ) ; # of Cases [Cases] ;var = (/ "AL", "BI", "EA", "FR", "IP", "MD", "ME", "SC" /) ; nVar = dimsizes(var) ; # of Variables ; more info to be added [all are bogus] var = (/ "Temp(~S~o~N~C)", "Precip(mm)"/) nVar = dimsizes(var) SEASONS = (/ "DJF", "JJA" /) nSeason = dimsizes( SEASONS ) ;================================================================== ; create weights: sqrt(cos(lat)) [or sqrt(gw) ] ; ================================================================= ; rad = 4.*atan(1.)/180. ; lat = f->latitude ; if (typeof(lat).eq."double") then ; clat = sqrt( cos(rad*tofloat(lat)) ) ; else ; clat = sqrt( cos(rad*lat) ) ; end if ; copy_VarCoords(lat, clat) ; contributed ; printVarSummary(clat) ;==============================taylorstats============= ;DJF AL ;====================================================== tA = TO(0,75:91,182:222) t1A = T1(0,75:91,182:222) t2A = T2(0,75:91,182:222) t3A = T3(0,75:91,182:222) t4A = T4(0,75:91,182:222) t5A = T5(0,75:91,182:222) PA = pm(0,75:91,182:222) P1A = P1(0,75:91,182:222) P2A = P2(0,75:91,182:222) P3A = P3(0,75:91,182:222) P4A = P4(0,75:91,182:222) P5A = P5(0,75:91,182:222) w = 1.0 stat_taylor1A = taylor_stats(t1A, tA, w, 0) stat_taylor2A = taylor_stats(t2A, tA, w, 0) stat_taylor3A = taylor_stats(t3A, tA, w, 0) stat_taylor4A = taylor_stats(t4A, tA, w, 0) stat_taylor5A = taylor_stats(t5A, tA, w, 0) stat_taylor1Ap = taylor_stats(P1A, PA, w, 0) stat_taylor2Ap = taylor_stats(P2A, PA, w, 0) stat_taylor3Ap = taylor_stats(P3A, PA, w, 0) stat_taylor4Ap = taylor_stats(P4A, PA, w, 0) stat_taylor5Ap = taylor_stats(P5A, PA, w, 0) ;==============================taylorstats============= ;DJF BI ;====================================================== tB = TO(0,99:135,122:170) t1B = T1(0,99:135,122:170) t2B = T2(0,99:135,122:170) t3B = T3(0,99:135,122:170) t4B = T4(0,99:135,122:170) t5B = T5(0,99:135,122:170) PB = pm(0,99:135,122:170) P1B = P1(0,99:135,122:170) P2B = P2(0,99:135,122:170) P3B = P3(0,99:135,122:170) P4B = P4(0,99:135,122:170) P5B = P5(0,99:135,122:170) stat_taylor1B = taylor_stats(t1B, tB, w, 0) stat_taylor2B = taylor_stats(t2B, tB, w, 0) stat_taylor3B = taylor_stats(t3B, tB, w, 0) stat_taylor4B = taylor_stats(t4B, tB, w, 0) stat_taylor5B = taylor_stats(t5B, tB, w, 0) stat_taylor1Bp = taylor_stats(P1B, PB, w, 0) stat_taylor2Bp = taylor_stats(P2B, PB, w, 0) stat_taylor3Bp = taylor_stats(P3B, PB, w, 0) stat_taylor4Bp = taylor_stats(P4B, PB, w, 0) stat_taylor5Bp = taylor_stats(P5B, PB, w, 0) ;==============================taylorstats============= ;DJF EA ;====================================================== tE = TO(0,76:119,227:282) t1E = T1(0,76:119,227:282) t2E = T2(0,76:119,227:282) t3E = T3(0,76:119,227:282) t4E = T4(0,76:119,227:282) t5E = T5(0,76:119,227:282) PE = pm(0,76:119,227:282) P1E = P1(0,76:119,227:282) P2E = P2(0,76:119,227:282) P3E = P3(0,76:119,227:282) P4E = P4(0,76:119,227:282) P5E = P5(0,76:119,227:282) stat_taylor1E = taylor_stats(t1E, tE, w, 0) stat_taylor2E = taylor_stats(t2E, tE, w, 0) stat_taylor3E = taylor_stats(t3E, tE, w, 0) stat_taylor4E = taylor_stats(t4E, tE, w, 0) stat_taylor5E = taylor_stats(t5E, tE, w, 0) stat_taylor1Ep = taylor_stats(P1E, PE, w, 0) stat_taylor2Ep = taylor_stats(P2E, PE, w, 0) stat_taylor3Ep = taylor_stats(P3E, PE, w, 0) stat_taylor4Ep = taylor_stats(P4E, PE, w, 0) stat_taylor5Ep = taylor_stats(P5E, PE, w, 0) ;==============================taylorstats============= ;DJF FR ;====================================================== tF = TO(0,75:99,142:182) t1F = T1(0,75:99,142:182) t2F = T2(0,75:99,142:182) t3F = T3(0,75:99,142:182) t4F = T4(0,75:99,142:182) t5F = T5(0,75:99,142:182) PF = pm(0,75:99,142:182) P1F = P1(0,75:99,142:182) P2F = P2(0,75:99,142:182) P3F = P3(0,75:99,142:182) P4F = P4(0,75:99,142:182) P5F = P5(0,75:99,142:182) stat_taylor1F = taylor_stats(t1F, tF, w, 0) stat_taylor2F = taylor_stats(t2F, tF, w, 0) stat_taylor3F = taylor_stats(t3F, tF, w, 0) stat_taylor4F = taylor_stats(t4F, tF, w, 0) stat_taylor5F = taylor_stats(t5F, tF, w, 0) stat_taylor1Fp = taylor_stats(P1F, PF, w, 0) stat_taylor2Fp = taylor_stats(P2F, PF, w, 0) stat_taylor3Fp = taylor_stats(P3F, PF, w, 0) stat_taylor4Fp = taylor_stats(P4F, PF, w, 0) stat_taylor5Fp = taylor_stats(P5F, PF, w, 0) ;==============================taylorstats============= ;DJF IP ;====================================================== tI = TO(0,44:75,123:147) t1I = T1(0,44:75,123:147) t2I = T2(0,44:75,123:147) t3I = T3(0,44:75,123:147) t4I = T4(0,44:75,123:147) t5I = T5(0,44:75,123:147) PI = pm(0,44:75,123:147) P1I = P1(0,44:75,123:147) P2I = P2(0,44:75,123:147) P3I = P3(0,44:75,123:147) P4I = P4(0,44:75,123:147) P5I = P5(0,44:75,123:147) stat_taylor1I = taylor_stats(t1I, tI, w, 0) stat_taylor2I = taylor_stats(t2I, tI, w, 0) stat_taylor3I = taylor_stats(t3I, tI, w, 0) stat_taylor4I = taylor_stats(t4I, tI, w, 0) stat_taylor5I = taylor_stats(t5I, tI, w, 0) stat_taylor1Ip = taylor_stats(P1I, PI, w, 0) stat_taylor2Ip = taylor_stats(P2I, PI, w, 0) stat_taylor3Ip = taylor_stats(P3I, PI, w, 0) stat_taylor4Ip = taylor_stats(P4I, PI, w, 0) stat_taylor5Ip = taylor_stats(P5I, PI, w, 0) ;==============================taylorstats============= ;DJF MD ;====================================================== tMD = TO(0,44:75,175:262) t1MD = T1(0,44:75,175:262) t2MD = T2(0,44:75,175:262) t3MD = T3(0,44:75,175:262) t4MD = T4(0,44:75,175:262) t5MD = T5(0,44:75,175:262) PMD = pm(0,44:75,175:262) P1MD = P1(0,44:75,175:262) P2MD = P2(0,44:75,175:262) P3MD = P3(0,44:75,175:262) P4MD = P4(0,44:75,175:262) P5MD = P5(0,44:75,175:262) stat_taylor1MD = taylor_stats(t1MD, tMD, w, 0) stat_taylor2MD = taylor_stats(t2MD, tMD, w, 0) stat_taylor3MD = taylor_stats(t3MD, tMD, w, 0) stat_taylor4MD = taylor_stats(t4MD, tMD, w, 0) stat_taylor5MD = taylor_stats(t5MD, tMD, w, 0) stat_taylor1MDp = taylor_stats(P1MD, PMD, w, 0) stat_taylor2MDp = taylor_stats(P2MD, PMD, w, 0) stat_taylor3MDp = taylor_stats(P3MD, PMD, w, 0) stat_taylor4MDp = taylor_stats(P4MD, PMD, w, 0) stat_taylor5MDp = taylor_stats(P5MD, PMD, w, 0) ;==============================taylorstats============= ;DJF ME ;====================================================== tME = TO(0,92:119,171:276) t1ME = T1(0,92:119,171:276) t2ME = T2(0,92:119,171:276) t3ME = T3(0,92:119,171:276) t4ME = T4(0,92:119,171:276) t5ME = T5(0,92:119,171:276) PME = pm(0,92:119,171:276) P1ME = P1(0,92:119,171:276) P2ME = P2(0,92:119,171:276) P3ME = P3(0,92:119,171:276) P4ME = P4(0,92:119,171:276) P5ME = P5(0,92:119,171:276) stat_taylor1ME = taylor_stats(t1ME, tME, w, 0) stat_taylor2ME = taylor_stats(t2ME, tME, w, 0) stat_taylor3ME = taylor_stats(t3ME, tME, w, 0) stat_taylor4ME = taylor_stats(t4ME, tME, w, 0) stat_taylor5ME = taylor_stats(t5ME, tME, w, 0) stat_taylor1MEp = taylor_stats(P1ME, PME, w, 0) stat_taylor2MEp = taylor_stats(P2ME, PME, w, 0) stat_taylor3MEp = taylor_stats(P3ME, PME, w, 0) stat_taylor4MEp = taylor_stats(P4ME, PME, w, 0) stat_taylor5MEp = taylor_stats(P5ME, PME, w, 0) ;==============================taylorstats============= ;DJF SC ;====================================================== tS = TO(0,120:179,183:282) t1S = T1(0,120:179,183:282) t2S = T2(0,120:179,183:282) t3S = T3(0,120:179,183:282) t4S = T4(0,120:179,183:282) t5S = T5(0,120:179,183:282) PS = pm(0,120:179,183:282) P1S = P1(0,120:179,183:282) P2S = P2(0,120:179,183:282) P3S = P3(0,120:179,183:282) P4S = P4(0,120:179,183:282) P5S = P5(0,120:179,183:282) stat_taylor1S = taylor_stats(t1S, tS, w, 0) stat_taylor2S = taylor_stats(t2S, tS, w, 0) stat_taylor3S = taylor_stats(t3S, tS, w, 0) stat_taylor4S = taylor_stats(t4S, tS, w, 0) stat_taylor5S = taylor_stats(t5S, tS, w, 0) stat_taylor1Sp = taylor_stats(P1S, PS, w, 0) stat_taylor2Sp = taylor_stats(P2S, PS, w, 0) stat_taylor3Sp = taylor_stats(P3S, PS, w, 0) stat_taylor4Sp = taylor_stats(P4S, PS, w, 0) stat_taylor5Sp = taylor_stats(P5S, PS, w, 0) ;==============================taylorstats============= ;JJA AL ;====================================================== ttA = TO(2,75:91,182:222) t11A = T1(2,75:91,182:222) t12A = T2(2,75:91,182:222) t13A = T3(2,75:91,182:222) t14A = T4(2,75:91,182:222) t15A = T5(2,75:91,182:222) PPA = pm(2,75:91,182:222) P11A = P1(2,75:91,182:222) P12A = P2(2,75:91,182:222) P13A = P3(2,75:91,182:222) P14A = P4(2,75:91,182:222) P15A = P5(2,75:91,182:222) stat_taylor11A = taylor_stats(t11A, ttA, w, 0) stat_taylor12A = taylor_stats(t12A, ttA, w, 0) stat_taylor13A = taylor_stats(t13A, ttA, w, 0) stat_taylor14A = taylor_stats(t14A, ttA, w, 0) stat_taylor15A = taylor_stats(t15A, ttA, w, 0) stat_taylor11Ap = taylor_stats(P11A, PPA, w, 0) stat_taylor12Ap = taylor_stats(P12A, PPA, w, 0) stat_taylor13Ap = taylor_stats(P13A, PPA, w, 0) stat_taylor14Ap = taylor_stats(P14A, PPA, w, 0) stat_taylor15Ap = taylor_stats(P15A, PPA, w, 0) ;==============================taylorstats============= ;JJA BI ;====================================================== ttB = TO(2,99:135,122:170) ;JJA t11B = T1(2,99:135,122:170) t12B = T2(2,99:135,122:170) t13B = T3(2,99:135,122:170) t14B = T4(2,99:135,122:170) t15B = T5(2,99:135,122:170) PPB = pm(2,99:135,122:170) ;JJA P11B = P1(2,99:135,122:170) P12B = P2(2,99:135,122:170) P13B = P3(2,99:135,122:170) P14B = P4(2,99:135,122:170) P15B = P5(2,99:135,122:170) stat_taylor11B = taylor_stats(t11B, ttB, w, 0) stat_taylor12B = taylor_stats(t12B, ttB, w, 0) stat_taylor13B = taylor_stats(t13B, ttB, w, 0) stat_taylor14B = taylor_stats(t14B, ttB, w, 0) stat_taylor15B = taylor_stats(t15B, ttB, w, 0) stat_taylor11Bp = taylor_stats(P11B, PPB, w, 0) stat_taylor12Bp = taylor_stats(P12B, PPB, w, 0) stat_taylor13Bp = taylor_stats(P13B, PPB, w, 0) stat_taylor14Bp = taylor_stats(P14B, PPB, w, 0) stat_taylor15Bp = taylor_stats(P15B, PPB, w, 0) ;==============================taylorstats============= ;JJA EA ;====================================================== ttE = TO(2,76:119,227:282) ;JJA t11E = T1(2,76:119,227:282) t12E = T2(2,76:119,227:282) t13E = T3(2,76:119,227:282) t14E = T4(2,76:119,227:282) t15E = T5(2,76:119,227:282) PPE = pm(2,76:119,227:282) ;JJA P11E = P1(2,76:119,227:282) P12E = P2(2,76:119,227:282) P13E = P3(2,76:119,227:282) P14E = P4(2,76:119,227:282) P15E = P5(2,76:119,227:282) stat_taylor11E = taylor_stats(t11E, ttE, w, 0) stat_taylor12E = taylor_stats(t12E, ttE, w, 0) stat_taylor13E = taylor_stats(t13E, ttE, w, 0) stat_taylor14E = taylor_stats(t14E, ttE, w, 0) stat_taylor15E = taylor_stats(t15E, ttE, w, 0) stat_taylor11Ep = taylor_stats(P11E, PPE, w, 0) stat_taylor12Ep = taylor_stats(P12E, PPE, w, 0) stat_taylor13Ep = taylor_stats(P13E, PPE, w, 0) stat_taylor14Ep = taylor_stats(P14E, PPE, w, 0) stat_taylor15Ep = taylor_stats(P15E, PPE, w, 0) ;==============================taylorstats============= ;JJA FR ;====================================================== ttF = TO(2,75:99,142:182) t11F = T1(2,75:99,142:182) t12F = T2(2,75:99,142:182) t13F = T3(2,75:99,142:182) t14F = T4(2,75:99,142:182) t15F = T5(2,75:99,142:182) PPF = pm(2,75:99,142:182) P11F = P1(2,75:99,142:182) P12F = P2(2,75:99,142:182) P13F = P3(2,75:99,142:182) P14F = P4(2,75:99,142:182) P15F = P5(2,75:99,142:182) stat_taylor11Fp = taylor_stats(P11F, PPF, w, 0) stat_taylor12Fp = taylor_stats(P12F, PPF, w, 0) stat_taylor13Fp = taylor_stats(P13F, PPF, w, 0) stat_taylor14Fp = taylor_stats(P14F, PPF, w, 0) stat_taylor15Fp = taylor_stats(P15F, PPF, w, 0) ;==============================taylorstats============= ;JJA IP ;====================================================== ttI = TO(2,44:75,123:147) t11I = T1(2,44:75,123:147) t12I = T2(2,44:75,123:147) t13I = T3(2,44:75,123:147) t14I = T4(2,44:75,123:147) t15I = T5(2,44:75,123:147) PPI = pm(2,44:75,123:147) P11I = P1(2,44:75,123:147) P12I = P2(2,44:75,123:147) P13I = P3(2,44:75,123:147) P14I = P4(2,44:75,123:147) P15I = P5(2,44:75,123:147) stat_taylor11I = taylor_stats(t11I, ttI, w, 0) stat_taylor12I = taylor_stats(t12I, ttI, w, 0) stat_taylor13I = taylor_stats(t13I, ttI, w, 0) stat_taylor14I = taylor_stats(t14I, ttI, w, 0) stat_taylor15I = taylor_stats(t15I, ttI, w, 0) stat_taylor11Ip = taylor_stats(P11I, PPI, w, 0) stat_taylor12Ip = taylor_stats(P12I, PPI, w, 0) stat_taylor13Ip = taylor_stats(P13I, PPI, w, 0) stat_taylor14Ip = taylor_stats(P14I, PPI, w, 0) stat_taylor15Ip = taylor_stats(P15I, PPI, w, 0) ;==============================taylorstats============= ;JJA MD ;====================================================== ttMD = TO(2,44:75,175:262) t11MD = T1(2,44:75,175:262) t12MD = T2(2,44:75,175:262) t13MD = T3(2,44:75,175:262) t14MD = T4(2,44:75,175:262) t15MD = T5(2,44:75,175:262) PPMD = pm(2,44:75,175:262) P11MD = P1(2,44:75,175:262) P12MD = P2(2,44:75,175:262) P13MD = P3(2,44:75,175:262) P14MD = P4(2,44:75,175:262) P15MD = P5(2,44:75,175:262) stat_taylor11MD = taylor_stats(t11MD, ttMD, w, 0) stat_taylor12MD = taylor_stats(t12MD, ttMD, w, 0) stat_taylor13MD = taylor_stats(t13MD, ttMD, w, 0) stat_taylor14MD = taylor_stats(t14MD, ttMD, w, 0) stat_taylor15MD = taylor_stats(t15MD, ttMD, w, 0) stat_taylor11MDp = taylor_stats(P11MD, PPMD, w, 0) stat_taylor12MDp = taylor_stats(P12MD, PPMD, w, 0) stat_taylor13MDp = taylor_stats(P13MD, PPMD, w, 0) stat_taylor14MDp = taylor_stats(P14MD, PPMD, w, 0) stat_taylor15MDp = taylor_stats(P15MD, PPMD, w, 0) ;==============================taylorstats============= ;JJA ME ;====================================================== ttME = TO(2,92:119,171:276) t11ME = T1(2,92:119,171:276) t12ME = T2(2,92:119,171:276) t13ME = T3(2,92:119,171:276) t14ME = T4(2,92:119,171:276) t15ME = T5(2,92:119,171:276) PPME = pm(2,92:119,171:276) P11ME = P1(2,92:119,171:276) P12ME = P2(2,92:119,171:276) P13ME = P3(2,92:119,171:276) P14ME = P4(2,92:119,171:276) P15ME = P5(2,92:119,171:276) stat_taylor11ME = taylor_stats(t11ME, ttME, w, 0) stat_taylor12ME = taylor_stats(t12ME, ttME, w, 0) stat_taylor13ME = taylor_stats(t13ME, ttME, w, 0) stat_taylor14ME = taylor_stats(t14ME, ttME, w, 0) stat_taylor15ME = taylor_stats(t15ME, ttME, w, 0) stat_taylor11MEp = taylor_stats(P11ME, PPME, w, 0) stat_taylor12MEp = taylor_stats(P12ME, PPME, w, 0) stat_taylor13MEp = taylor_stats(P13ME, PPME, w, 0) stat_taylor14MEp = taylor_stats(P14ME, PPME, w, 0) stat_taylor15MEp = taylor_stats(P15ME, PPME, w, 0) ;==============================taylorstats============= ;JJA SC ;====================================================== ttS = TO(2,120:179,183:282) t11S = T1(2,120:179,183:282) t12S = T2(2,120:179,183:282) t13S = T3(2,120:179,183:282) t14S = T4(2,120:179,183:282) t15S = T5(2,120:179,183:282) PPS = pm(2,120:179,183:282) P11S = P1(2,120:179,183:282) P12S = P2(2,120:179,183:282) P13S = P3(2,120:179,183:282) P14S = P4(2,120:179,183:282) P15S = P5(2,120:179,183:282) stat_taylor11S = taylor_stats(t11S, ttS, w, 0) stat_taylor12S = taylor_stats(t12S, ttS, w, 0) stat_taylor13S = taylor_stats(t13S, ttS, w, 0) stat_taylor14S = taylor_stats(t14S, ttS, w, 0) stat_taylor15S = taylor_stats(t15S, ttS, w, 0) stat_taylor11Sp = taylor_stats(P11S, PPS, w, 0) stat_taylor12Sp = taylor_stats(P12S, PPS, w, 0) stat_taylor13Sp = taylor_stats(P13S, PPS, w, 0) stat_taylor14Sp = taylor_stats(P14S, PPS, w, 0) stat_taylor15Sp = taylor_stats(P15S, PPS, w, 0) ;=====DJF================================================================ CA_cc_AL = (/stat_taylor1A(0), stat_taylor1Ap(0)/) CA_cc1_AL = (/stat_taylor2A(0), stat_taylor2Ap(0)/) CA_cc2_AL = (/stat_taylor3A(0), stat_taylor3Ap(0)/) CA_cc3_AL = (/stat_taylor4A(0), stat_taylor4Ap(0)/) CA_cc4_AL = (/stat_taylor5A(0), stat_taylor5Ap(0)/) CA_cc_BI = (/stat_taylor1B(0), stat_taylor1Bp(0)/) CA_cc1_BI = (/stat_taylor2B(0), stat_taylor2Bp(0)/) CA_cc2_BI = (/stat_taylor3B(0), stat_taylor3Bp(0)/) CA_cc3_BI = (/stat_taylor4B(0), stat_taylor4Bp(0)/) CA_cc4_BI = (/stat_taylor5B(0), stat_taylor5Bp(0)/) CA_cc_EA = (/stat_taylor1E(0), stat_taylor1Ep(0)/) CA_cc1_EA = (/stat_taylor2E(0), stat_taylor2Ep(0)/) CA_cc2_EA = (/stat_taylor3E(0), stat_taylor3Ep(0)/) CA_cc3_EA = (/stat_taylor4E(0), stat_taylor4Ep(0)/) CA_cc4_EA = (/stat_taylor5E(0), stat_taylor5Ep(0)/) CA_cc_FR = (/stat_taylor1F(0), stat_taylor1Fp(0)/) CA_cc1_FR = (/stat_taylor2F(0), stat_taylor2Fp(0)/) CA_cc2_FR = (/stat_taylor3F(0), stat_taylor3Fp(0)/) CA_cc3_FR = (/stat_taylor4F(0), stat_taylor4Fp(0)/) CA_cc4_FR = (/stat_taylor5F(0), stat_taylor5Fp(0)/) CA_cc_IP = (/stat_taylor1I(0), stat_taylor1Ip(0)/) CA_cc1_IP = (/stat_taylor2I(0), stat_taylor2Ip(0)/) CA_cc2_IP = (/stat_taylor3I(0), stat_taylor3Ip(0)/) CA_cc3_IP = (/stat_taylor4I(0), stat_taylor4Ip(0)/) CA_cc4_IP = (/stat_taylor5I(0), stat_taylor5Ip(0)/) CA_cc_MD = (/stat_taylor1MD(0), stat_taylor1MDp(0)/) CA_cc1_MD = (/stat_taylor2MD(0), stat_taylor2MDp(0)/) CA_cc2_MD = (/stat_taylor3MD(0), stat_taylor3MDp(0)/) CA_cc3_MD = (/stat_taylor4MD(0), stat_taylor4MDp(0)/) CA_cc4_MD = (/stat_taylor5MD(0), stat_taylor5MDp(0)/) CA_cc_ME = (/stat_taylor1ME(0), stat_taylor1MEp(0)/) CA_cc1_ME = (/stat_taylor2ME(0), stat_taylor2MEp(0)/) CA_cc2_ME = (/stat_taylor3ME(0), stat_taylor3MEp(0)/) CA_cc3_ME = (/stat_taylor4ME(0), stat_taylor4MEp(0)/) CA_cc4_ME = (/stat_taylor5ME(0), stat_taylor5MEp(0)/) CA_cc_SC = (/stat_taylor1S(0), stat_taylor1Sp(0)/) CA_cc1_SC = (/stat_taylor2S(0), stat_taylor2Sp(0)/) CA_cc2_SC = (/stat_taylor3S(0), stat_taylor3Sp(0)/) CA_cc3_SC = (/stat_taylor4S(0), stat_taylor4Sp(0)/) CA_cc4_SC = (/stat_taylor5S(0), stat_taylor5Sp(0)/) CA_rat_AL = (/stat_taylor1A(1), stat_taylor1Ap(1)/) CA_rat1_AL = (/stat_taylor2A(1), stat_taylor2Ap(1)/) CA_rat2_AL = (/stat_taylor3A(1), stat_taylor3Ap(1)/) CA_rat3_AL = (/stat_taylor4A(1), stat_taylor4Ap(1)/) CA_rat4_AL = (/stat_taylor5A(1), stat_taylor5Ap(1)/) CA_rat_BI = (/stat_taylor1B(1), stat_taylor1Bp(1)/) CA_rat1_BI = (/stat_taylor2B(1), stat_taylor2Bp(1)/) CA_rat2_BI = (/stat_taylor3B(1), stat_taylor3Bp(1)/) CA_rat3_BI = (/stat_taylor4B(1), stat_taylor4Bp(1)/) CA_rat4_BI = (/stat_taylor5B(1), stat_taylor5Bp(1)/) CA_rat_EA = (/stat_taylor1E(1), stat_taylor1Ep(1)/) CA_rat1_EA = (/stat_taylor2E(1), stat_taylor2Ep(1)/) CA_rat2_EA = (/stat_taylor3E(1), stat_taylor3Ep(1)/) CA_rat3_EA = (/stat_taylor4E(1), stat_taylor4Ep(1)/) CA_rat4_EA = (/stat_taylor5E(1), stat_taylor5Ep(1)/) CA_rat_FR = (/stat_taylor1F(1), stat_taylor1Fp(1)/) CA_rat1_FR = (/stat_taylor2F(1), stat_taylor2Fp(1)/) CA_rat2_FR = (/stat_taylor3F(1), stat_taylor3Fp(1)/) CA_rat3_FR = (/stat_taylor4F(1), stat_taylor4Fp(1)/) CA_rat4_FR = (/stat_taylor5F(1), stat_taylor5Fp(1)/) CA_rat_IP = (/stat_taylor1I(1), stat_taylor1Ip(1)/) CA_rat1_IP = (/stat_taylor2I(1), stat_taylor2Ip(1)/) CA_rat2_IP = (/stat_taylor3I(1), stat_taylor3Ip(1)/) CA_rat3_IP = (/stat_taylor4I(1), stat_taylor4Ip(1)/) CA_rat4_IP = (/stat_taylor5I(1), stat_taylor5Ip(1)/) CA_rat_MD = (/stat_taylor1MD(1), stat_taylor1MDp(1)/) CA_rat1_MD = (/stat_taylor2MD(1), stat_taylor2MDp(1)/) CA_rat2_MD = (/stat_taylor3MD(1), stat_taylor3MDp(1)/) CA_rat3_MD = (/stat_taylor4MD(1), stat_taylor4MDp(1)/) CA_rat4_MD = (/stat_taylor5MD(1), stat_taylor5MDp(1)/) CA_rat_ME = (/stat_taylor1ME(1), stat_taylor1MEp(1)/) CA_rat1_ME = (/stat_taylor2ME(1), stat_taylor2MEp(1)/) CA_rat2_ME = (/stat_taylor3ME(1), stat_taylor3MEp(1)/) CA_rat3_ME = (/stat_taylor4ME(1), stat_taylor4MEp(1)/) CA_rat4_ME = (/stat_taylor5ME(1), stat_taylor5MEp(1)/) CA_rat_SC = (/stat_taylor1S(1), stat_taylor1Sp(1)/) CA_rat1_SC = (/stat_taylor2S(1), stat_taylor2Sp(1)/) CA_rat2_SC = (/stat_taylor3S(1), stat_taylor3Sp(1)/) CA_rat3_SC = (/stat_taylor4S(1), stat_taylor4Sp(1)/) CA_rat4_SC = (/stat_taylor5S(1), stat_taylor5Sp(1)/) ;********************************** ; Put the ratios and pattern correlations into ; arrays for plotting ;********************************** ; arrays to be passed to taylor_diagram. It will calculate the x xnd y coordinates. ratio_AL = new ((/nCase, nVar/), typeof(CA_cc_AL) ) ratio_BI = new ((/nCase, nVar/), typeof(CA_cc_AL) ) ratio_EA = new ((/nCase, nVar/), typeof(CA_cc_AL) ) ratio_FR = new ((/nCase, nVar/), typeof(CA_cc_AL) ) ratio_IP = new ((/nCase, nVar/), typeof(CA_cc_AL) ) ratio_MD = new ((/nCase, nVar/), typeof(CA_cc_AL) ) ratio_ME = new ((/nCase, nVar/), typeof(CA_cc_AL) ) ratio_SC = new ((/nCase, nVar/), typeof(CA_cc_AL) ) cc_AL = new ((/nCase, nVar/), typeof(CA_cc_AL) ) cc_BI = new ((/nCase, nVar/), typeof(CA_cc_AL) ) cc_EA = new ((/nCase, nVar/), typeof(CA_cc_AL) ) cc_FR = new ((/nCase, nVar/), typeof(CA_cc_AL) ) cc_IP = new ((/nCase, nVar/), typeof(CA_cc_AL) ) cc_MD = new ((/nCase, nVar/), typeof(CA_cc_AL) ) cc_ME = new ((/nCase, nVar/), typeof(CA_cc_AL) ) cc_SC = new ((/nCase, nVar/), typeof(CA_cc_AL) ) ratio_AL(0,:) = CA_rat_AL ratio_AL(1,:) = CA_rat1_AL ratio_AL(2,:) = CA_rat2_AL ratio_AL(3,:) = CA_rat3_AL ratio_AL(4,:) = CA_rat4_AL ratio_BI(0,:) = CA_rat_BI ratio_BI(1,:) = CA_rat1_BI ratio_BI(2,:) = CA_rat2_BI ratio_BI(3,:) = CA_rat3_BI ratio_BI(4,:) = CA_rat4_BI ratio_EA(0,:) = CA_rat_EA ratio_EA(1,:) = CA_rat1_EA ratio_EA(2,:) = CA_rat2_EA ratio_EA(3,:) = CA_rat3_EA ratio_EA(4,:) = CA_rat4_EA ratio_FR(0,:) = CA_rat_FR ratio_FR(1,:) = CA_rat1_FR ratio_FR(2,:) = CA_rat2_FR ratio_FR(3,:) = CA_rat3_FR ratio_FR(4,:) = CA_rat4_FR ratio_IP(0,:) = CA_rat_IP ratio_IP(1,:) = CA_rat1_IP ratio_IP(2,:) = CA_rat2_IP ratio_IP(3,:) = CA_rat3_IP ratio_IP(4,:) = CA_rat4_IP ratio_MD(0,:) = CA_rat_MD ratio_MD(1,:) = CA_rat1_MD ratio_MD(2,:) = CA_rat2_MD ratio_MD(3,:) = CA_rat3_MD ratio_MD(4,:) = CA_rat4_MD ratio_ME(0,:) = CA_rat_ME ratio_ME(1,:) = CA_rat1_ME ratio_ME(2,:) = CA_rat2_ME ratio_ME(3,:) = CA_rat3_ME ratio_ME(4,:) = CA_rat4_ME ratio_SC(0,:) = CA_rat_SC ratio_SC(1,:) = CA_rat1_SC ratio_SC(2,:) = CA_rat2_SC ratio_SC(3,:) = CA_rat3_SC ratio_SC(4,:) = CA_rat4_SC cc_AL(0,:) = CA_cc_AL cc_AL(1,:) = CA_cc1_AL cc_AL(2,:) = CA_cc2_AL cc_AL(3,:) = CA_cc3_AL cc_AL(4,:) = CA_cc4_AL cc_BI(0,:) = CA_cc_BI cc_BI(1,:) = CA_cc1_BI cc_BI(2,:) = CA_cc2_BI cc_BI(3,:) = CA_cc3_BI cc_BI(4,:) = CA_cc4_BI cc_EA(0,:) = CA_cc_EA cc_EA(1,:) = CA_cc1_EA cc_EA(2,:) = CA_cc2_EA cc_EA(3,:) = CA_cc3_EA cc_EA(4,:) = CA_cc4_EA cc_FR(0,:) = CA_cc_FR cc_FR(1,:) = CA_cc1_FR cc_FR(2,:) = CA_cc2_FR cc_FR(3,:) = CA_cc3_FR cc_FR(4,:) = CA_cc4_FR cc_IP(0,:) = CA_cc_IP cc_IP(1,:) = CA_cc1_IP cc_IP(2,:) = CA_cc2_IP cc_IP(3,:) = CA_cc3_IP cc_IP(4,:) = CA_cc4_IP cc_MD(0,:) = CA_cc_MD cc_MD(1,:) = CA_cc1_MD cc_MD(2,:) = CA_cc2_MD cc_MD(3,:) = CA_cc3_MD cc_MD(4,:) = CA_cc4_MD cc_ME(0,:) = CA_cc_ME cc_ME(1,:) = CA_cc1_ME cc_ME(2,:) = CA_cc2_ME cc_ME(3,:) = CA_cc3_ME cc_ME(4,:) = CA_cc4_ME cc_SC(0,:) = CA_cc_SC cc_SC(1,:) = CA_cc1_SC cc_SC(2,:) = CA_cc2_SC cc_SC(3,:) = CA_cc3_SC cc_SC(4,:) = CA_cc4_SC ;print(bias) ;print(cc) ;print(ratio) ;================================ ; PLOTS ;================================ res = True ; default taylor diagram wks = gsn_open_wks("pdf","taylor_pr_tas_ysm_DJF") ; send graphics to PNG file plot = new(8,graphic) ; varSource = var +"_"+ source ; add extra info [*not* required] varSource = var ; add extra info [*not* required] ty_opt = True ; taylor diagram with options ty_opt@Markers = (/16, 16, 16, 16, 16/) ; make all solid fill ty_opt@Colors = (/"red", "blue", "green", "orange", "magenta" /) ty_opt@varLabels = varSource ty_opt@caseLabels = case ty_opt@varLabelsYloc = 0.7 ; Move location of variable labels [default 0.45] ty_opt@caseLabelsFontHeightF = 0.03 ; make slight larger [default=0.12 ] ty_opt@varLabelsFontHeightF = 0.02 ; make slight smaller [default=0.013] ty_opt@tiMainString = "DJF" ; title ty_opt@stnRad = (/ 0.5,1.5/) ; additional standard radii ty_opt@ccRays = (/ 0.6, 0.9 /) ; correllation rays ty_opt@centerDiffRMS = True ; RMS 'circles' ty_opt@vpWidthF = 0.10 ty_opt@vpHeightF = 0.02*nCase ; ty_opt@tiMainString = "DJF" plot(0) = taylor_diagram(wks, ratio_AL, cc_AL, ty_opt) plot(1) = taylor_diagram(wks, ratio_BI, cc_BI, ty_opt) plot(2) = taylor_diagram(wks, ratio_EA, cc_EA, ty_opt) plot(3) = taylor_diagram(wks, ratio_FR, cc_FR, ty_opt) plot(4) = taylor_diagram(wks, ratio_IP, cc_IP, ty_opt) plot(5) = taylor_diagram(wks, ratio_MD, cc_MD, ty_opt) plot(6) = taylor_diagram(wks, ratio_ME, cc_ME, ty_opt) plot(7) = taylor_diagram(wks, ratio_SC, cc_SC, ty_opt) ;*************************************************** ; Panel plot ;*************************************************** resP = True ; resP@txString = "Comarison of bias corrections" resP@gsnPanelYWhiteSpacePercent = 3 ; add extra white space resP@gsnPanelFigureStringsPerimOn = False resP@amJust = "TopLeft" ; resP@gsnPanelFigureStrings = (/"AL","BI","EA","FR","IP",\ ; "MD","ME","SC"/) resP@gsnPanelFigureStrings = (/"AL","BI","EA","FR"/) resP@gsnPanelFigureStringsFontHeightF = 0.02 gsn_panel(wks,plot(0:3),(/2,2/),resP) end