Title: SST indices regression to SLP
Description: |
    SST indices included NAO, AMO, subpolar, subtropical gyre and NINO3.4 index.
COMMENT: |
    SST indices included AMO, subpolar, subtropical gyre and NINO3.4 index.
    NAO is defined as the difference between the monthly SLP anomalies averaged over the domains of (50°W-10°E, 25°-55°N) and (40°W-20°E, 55-85°N).
    AMO, defined as 10-year low-pass filtered area-averaged SST anomalies over the North Atlantic basin (0°N-65°N, 80°W-0°E). 
    Mean SST of subpolar grye region,  40°W-10°W, 50°N-60°N
    Mean SST of subtropical grye region, 70°W-20°W, 20°N-40°N
    Mean DJF SST of NINO3.4 region, 170°W-120°W, 5°S-5°N

#Thumbnail: 


## plots setting
FIGFILES: $(cat figfilelist.txt) # for mk_html_pages.sh
#### NAO
Plot_NAO:
    plotScript: mpiexm/plot_mpiexm_NAO.ncl
    TITLE: "NAO index NAME "
    FIGFILENAME: CASE_slp_nmodNMOD_nao_DJF
    MONTHS: 12,13,14
    COMPONENT:  echam6
    VARIABLE:  slp
    PLOTRES: |
        res@tiMainString = title
        res@gsnLeftString = "model nyear: NYR"
        res@trYMinF = -3.
        res@trYMaxF =  3.
        res@vpHeightF = .3
        res@tiYAxisString = "std. SLP dif."
        res@gsnYRefLine = 0.
NAOts_forRegMap:
    TSNETCDFFILE: CASE_slp_nmodNMOD_nao_DJF
    TSVARNAME: tsdif
    TSNA: NAO
    MONTHS: 12,13,14
    RANGE: |
        res@cnLevelSelectionMode = "ExplicitLevels"
        res@cnLevels = fspan(-.6,.6,13)
        

#### AMO
Plot_AMO:
    plotScript: mpiexm/plot_mpiexm_ts.ncl
    TITLE: "AMO index NAME "
    FIGFILENAME: CASE_sst_amo_nmodNMOD_ANN
    COMPONENT:  mpiom
    VARIABLE: sst
    LATBE: 0.,65.
    LONBE: -80.,0.
    MONTHS: 1,2,3,4,5,6,7,8,9,10,11,12
    TSPOSTPROC: |
        ;; low pass
        wgt = filwgts_lanczos(21,0,1./10.,-999,1.)
        ts := wgt_runave(ts,wgt,0)
        x = ispan(min(ybe),max(ybe),1)
    PLOTRES: |
        res@tiMainString = title
        res@gsnLeftString = "model nyear: NYR"
        res@vpHeightF = .3
        res@trYMinF = -0.4
        res@trYMaxF =  0.4
        res@tmXBFormat = "f3.0"
        ts_mean = avg(ts)
        ts := ts - ts_mean ;; remove mean
        res@gsnRightString = "avg="+sprintf("%3.2f",ts_mean)+" K"
AMOts_forRegMap:
    TSNETCDFFILE: CASE_sst_amo_nmodNMOD_ANN
    TSVARNAME: ts
    TSNA: AMO
    MONTHS: 1,2,3,4,5,6,7,8,9,10,11,12
    RANGE: |
        res@cnLevelSelectionMode = "ExplicitLevels"
        res@cnLevels = fspan(-600.,600.,13)
    LATBE: 0.,65.
    LONBE: -80.,0.
    POSTPLOT: |  ## add box
        latbe = (/LATBE/)
        lonbe = (/LONBE/)
        boxx = (/min(lonbe),max(lonbe),max(lonbe),min(lonbe),min(lonbe)/)
        boxy = (/min(latbe),min(latbe),max(latbe),max(latbe),min(latbe)/)
        lres = True
            lres@gsLineThicknessF = 3.0
        b = gsn_add_polyline(wks,plot,boxx,boxy,lres)

#### subtropical mean SST
Plot_sst_subT:
    plotScript: mpiexm/plot_mpiexm_ts.ncl
    TITLE: "SubT SST NAME "
    FIGFILENAME: CASE_sst_nmodNMOD_subT_ANN
    COMPONENT:  mpiom
    VARIABLE: sst
    LATBE: 20.,40.
    LONBE: -70.,-20.
    MONTHS: 1,2,3,4,5,6,7,8,9,10,11,12
    TSPOSTPROC: |
        ;; low pass
        ;;;wgt = filwgts_lanczos(21,0,1./10.,-999,1.)
        ;;;ts := wgt_runave(ts,wgt,0)
        x = ispan(min(ybe),max(ybe),1)
    PLOTRES: |
        res@tiMainString = title
        res@gsnLeftString = "model nyear: NYR"
        res@vpHeightF = .3
        res@trYMinF = -0.4
        res@trYMaxF =  0.4
        res@tmXBFormat = "f3.0"
        ts_mean = avg(ts)
        ts := ts - ts_mean ;; remove mean
        res@gsnRightString = "avg="+sprintf("%3.2f",ts_mean)+" K"
sst_subTts_forRegMap:
    TSNETCDFFILE: CASE_sst_nmodNMOD_subT_ANN
    TSVARNAME: ts
    TSNA: subT 
    MONTHS: 1,2,3,4,5,6,7,8,9,10,11,12
    RANGE: |
        res@cnLevelSelectionMode = "ExplicitLevels"
        res@cnLevels = fspan(-510.,510.,13)
    LATBE: 20.,40.
    LONBE: -70.,-20.
    POSTPLOT: |  ## add box
        latbe = (/LATBE/)
        lonbe = (/LONBE/)
        boxx = (/min(lonbe),max(lonbe),max(lonbe),min(lonbe),min(lonbe)/)
        boxy = (/min(latbe),min(latbe),max(latbe),max(latbe),min(latbe)/)
        lres = True
            lres@gsLineThicknessF = 3.0
        b = gsn_add_polyline(wks,plot,boxx,boxy,lres)

#### subpolar mean SST
Plot_sst_subP:
    plotScript: mpiexm/plot_mpiexm_ts.ncl
    TITLE: "SubP SST NAME "
    FIGFILENAME: CASE_sst_nmodNMOD_subP_ANN
    COMPONENT:  mpiom
    VARIABLE: sst
    LATBE: 50.,60.
    LONBE: -40.,-10.
    MONTHS: 1,2,3,4,5,6,7,8,9,10,11,12
    TSPOSTPROC: |
        ;; low pass
        ;;;wgt = filwgts_lanczos(21,0,1./10.,-999,1.)
        ;;;ts := wgt_runave(ts,wgt,0)
        x = ispan(min(ybe),max(ybe),1)
    PLOTRES: |
        res@tiMainString = title
        res@gsnLeftString = "model nyear: NYR"
        res@vpHeightF = .3
        res@trYMinF = -1.0
        res@trYMaxF =  1.0
        res@tmXBFormat = "f3.0"
        ts_mean = avg(ts)
        ts := ts - ts_mean ;; remove mean
        res@gsnRightString = "avg="+sprintf("%3.2f",ts_mean)+" K"
sst_subPts_forRegMap:
    TSNETCDFFILE: CASE_sst_nmodNMOD_subP_ANN
    TSVARNAME: ts
    TSNA: subP 
    MONTHS: 1,2,3,4,5,6,7,8,9,10,11,12
    RANGE: |
        res@cnLevelSelectionMode = "ExplicitLevels"
        res@cnLevels = fspan(-210.,210.,13)
    LATBE: 50.,60.
    LONBE: -40.,-10.
    POSTPLOT: |  ## add box
        latbe = (/LATBE/)
        lonbe = (/LONBE/)
        boxx = (/min(lonbe),max(lonbe),max(lonbe),min(lonbe),min(lonbe)/)
        boxy = (/min(latbe),min(latbe),max(latbe),max(latbe),min(latbe)/)
        lres = True
            lres@gsLineThicknessF = 3.0
        b = gsn_add_polyline(wks,plot,boxx,boxy,lres)

#### NINO3.4 mean SST
Plot_sst_nino34:
    plotScript: mpiexm/plot_mpiexm_ts.ncl
    TITLE: "NINO3.4 SST NAME "
    FIGFILENAME: CASE_sst_nmodNMOD_nino34_DJF
    COMPONENT:  mpiom
    VARIABLE: sst
    LATBE: -5.,5.
    LONBE: -170.,-120.
    MONTHS: 12,13,14
    TSPOSTPROC: |
        ;; low pass
        ;;;wgt = filwgts_lanczos(21,0,1./10.,-999,1.)
        ;;;ts := wgt_runave(ts,wgt,0)
        x = ispan(min(ybe),max(ybe),1)
    PLOTRES: |
        res@tiMainString = title
        res@gsnLeftString = "model nyear: NYR"
        res@vpHeightF = .3
        res@trYMinF = -3.1
        res@trYMaxF =  -res@trYMinF
        res@tmXBFormat = "f3.0"
        ts_mean = avg(ts)
        ts := ts - ts_mean ;; remove mean
        res@gsnRightString = "avg="+sprintf("%3.2f",ts_mean)+" K"
sst_nino34ts_forRegMap:
    TSNETCDFFILE: CASE_sst_nmodNMOD_nino34_DJF
    TSVARNAME: ts
    TSNA: NINO3.4
    MONTHS: 12,13,14
    RANGE: |
        res@cnLevelSelectionMode = "ExplicitLevels"
        res@cnLevels = fspan(-210.,210.,13)
    LATBE: -5.,5.
    LONBE: -170.,-120.
    POSTPLOT: |  ## add box
        latbe = (/LATBE/)
        lonbe = (/LONBE/)
        boxx = (/min(lonbe),max(lonbe),max(lonbe),min(lonbe),min(lonbe)/)
        boxy = (/min(latbe),min(latbe),max(latbe),max(latbe),min(latbe)/)
        lres = True
            lres@gsLineThicknessF = 3.0
        b = gsn_add_polyline(wks,plot,boxx,boxy,lres)

#### Regression to SLP
Plot_SLPRegMap:
    plotScript: mpiexm/plot_mpiexm_regmap.ncl
    TITLE: "TSNA index reg. NAME (model year: NYR)"
    FIGFILENAME: CASE_TSNA_nmodNMOD_slp_regmap
    COMPONENT:  echam6
    VARIABLE:  slp
    PLOTRES: |
        res@mpLimitMode = "LatLon"
        res@mpMinLatF = -90.
        res@mpMaxLatF = 90.
        res@mpMinLonF = -330.
        res@mpMaxLonF =   30.
        res@pmTickMarkDisplayMode = "Always"
        res@mpCenterLonF =   -150.
        res@cnLinesOn = False
        res@cnFillOn = True
        res@lbOrientation = "Vertical"
        res@pmLabelBarWidthF = 0.1
        res@pmLabelBarOrthogonalPosF   = 0.01
        RANGE

## case common setting    
CASEDIR: /projects/NS9207K/pgchiu/CASE
    ## CASE define later
YEARBE: 1101,1250
NYR: 150

## case setting
Case_CTL:
    CASE: f08_1AGCM
    NAME: CTL
    NMOD: 0
Case_CTL2:
    CASE: f08_5AGCMs_NASPG_sample
    NAME: CTL2
    MODELN: 1
    NMOD: 0
Case_AVG:
    CASE: f08_5AGCMs_NASPG_average
    NAME: IEave
    MODELN: 1
    NMOD: 1
Case_AVG5:
    CASE: f08_5AGCMs_NASPG_average
    NAME: IEave5
    MODELN: (/1,2,3,4,5/)
    NMOD: 5
Case_SAM:
    CASE: f09_5AGCMs_NASPG_sample
    NAME: IEsam
    MODELN: 1
    NMOD: 1
Case_SAM5:
    CASE: f09_5AGCMs_NASPG_sample
    NAME: IEsam5
    MODELN: (/1,2,3,4,5/)
    NMOD: 5


Scripts:
-  # Remove figfilelist.txt which is for mk_html_page.sh at bottom
    plotScript: mk_runcode.sh
    RUNTHESECODES: |
        rm -f figfilelist.txt
## NAO index
- varPack: Plot_NAO, Case_CTL
- varPack: Plot_NAO, Case_SAM
- varPack: Plot_NAO, Case_AVG 
- varPack: Plot_NAO, Case_CTL2
- varPack: Plot_NAO, Case_SAM5
- varPack: Plot_NAO, Case_AVG5
- varPack: Plot_SLPRegMap, NAOts_forRegMap, Case_CTL
- varPack: Plot_SLPRegMap, NAOts_forRegMap, Case_SAM
- varPack: Plot_SLPRegMap, NAOts_forRegMap, Case_AVG
- varPack: Plot_SLPRegMap, NAOts_forRegMap, Case_CTL2
- varPack: Plot_SLPRegMap, NAOts_forRegMap, Case_SAM5
- varPack: Plot_SLPRegMap, NAOts_forRegMap, Case_AVG5

## AMO index
- varPack: Plot_AMO, Case_CTL
- varPack: Plot_AMO, Case_SAM
- varPack: Plot_AMO, Case_AVG 
- varPack: Plot_AMO, Case_CTL2
- varPack: Plot_AMO, Case_SAM5
- varPack: Plot_AMO, Case_AVG5
- varPack: Plot_SLPRegMap, AMOts_forRegMap, Case_CTL
- varPack: Plot_SLPRegMap, AMOts_forRegMap, Case_SAM
- varPack: Plot_SLPRegMap, AMOts_forRegMap, Case_AVG
- varPack: Plot_SLPRegMap, AMOts_forRegMap, Case_CTL2
- varPack: Plot_SLPRegMap, AMOts_forRegMap, Case_SAM5
- varPack: Plot_SLPRegMap, AMOts_forRegMap, Case_AVG5

## sst_subT index
- varPack: Plot_sst_subT, Case_CTL
- varPack: Plot_sst_subT, Case_SAM
- varPack: Plot_sst_subT, Case_AVG 
- varPack: Plot_sst_subT, Case_CTL2
- varPack: Plot_sst_subT, Case_SAM5
- varPack: Plot_sst_subT, Case_AVG5
- varPack: Plot_SLPRegMap, sst_subTts_forRegMap, Case_CTL
- varPack: Plot_SLPRegMap, sst_subTts_forRegMap, Case_SAM
- varPack: Plot_SLPRegMap, sst_subTts_forRegMap, Case_AVG
- varPack: Plot_SLPRegMap, sst_subTts_forRegMap, Case_CTL2
- varPack: Plot_SLPRegMap, sst_subTts_forRegMap, Case_SAM5
- varPack: Plot_SLPRegMap, sst_subTts_forRegMap, Case_AVG5

## sst_subP index
- varPack: Plot_sst_subP, Case_CTL
- varPack: Plot_sst_subP, Case_SAM
- varPack: Plot_sst_subP, Case_AVG 
- varPack: Plot_sst_subP, Case_CTL2
- varPack: Plot_sst_subP, Case_SAM5
- varPack: Plot_sst_subP, Case_AVG5
- varPack: Plot_SLPRegMap, sst_subPts_forRegMap, Case_CTL
- varPack: Plot_SLPRegMap, sst_subPts_forRegMap, Case_SAM
- varPack: Plot_SLPRegMap, sst_subPts_forRegMap, Case_AVG
- varPack: Plot_SLPRegMap, sst_subPts_forRegMap, Case_CTL2
- varPack: Plot_SLPRegMap, sst_subPts_forRegMap, Case_SAM5
- varPack: Plot_SLPRegMap, sst_subPts_forRegMap, Case_AVG5

## sst_nino34 index
- varPack: Plot_sst_nino34, Case_CTL
- varPack: Plot_sst_nino34, Case_SAM
- varPack: Plot_sst_nino34, Case_AVG 
- varPack: Plot_sst_nino34, Case_CTL2
- varPack: Plot_sst_nino34, Case_SAM5
- varPack: Plot_sst_nino34, Case_AVG5
- varPack: Plot_SLPRegMap, sst_nino34ts_forRegMap, Case_CTL
- varPack: Plot_SLPRegMap, sst_nino34ts_forRegMap, Case_SAM
- varPack: Plot_SLPRegMap, sst_nino34ts_forRegMap, Case_AVG
- varPack: Plot_SLPRegMap, sst_nino34ts_forRegMap, Case_CTL2
- varPack: Plot_SLPRegMap, sst_nino34ts_forRegMap, Case_SAM5
- varPack: Plot_SLPRegMap, sst_nino34ts_forRegMap, Case_AVG5


-  
    plotScript: mk_html_page.sh
    COLUMN: >
        'CTL' 'IEsam(1st AGCM)' 'IEavg(1st AGCM)' 'CTL2' tr 
        'NAO index DJF' tr 
        4 
        'Ensemble ->' 2
        'NAO index DJF SLP regression map' tr 
        4
        'Ensemble ->' 2
        'AMO index' tr 
        4 
        'Ensemble ->' 2
        'AMO index SLP regression map' tr 
        4
        'Ensemble ->' 2
        'Mean subtropical gyre annual SST' tr 
        4 
        'Ensemble ->' 2
        'Mean subtropical gyre annual SST to SLP regression map' tr 
        4 
        'Ensemble ->' 2
        'Mean subpolar gyre annual SST' tr 
        4 
        'Ensemble ->' 2
        'Mean subpolar gyre annual SST to SLP regression map' tr 
        4
        'Ensemble ->' 2
        'Mean NINO3.4 DJF SST' tr 
        4 
        'Ensemble ->' 2
        'Mean NINO3.4 DJF SST to SLP regression map' tr 
        4
        'Ensemble ->' 2