This package of codes is to fitting statistical modelling for ERAint Hs ~ ERAint SLP using 6hly ERAint data, and further apply the model fitting on other simulation runs, to predict wave height from SLP data. There are 2 sets of codes, one set is working on ERAint Hs ~ SLP, to get model fitting parameters, all those code locate at ERAint/, some of the data dir in the code need to be adjust if users want to run them. The data steps for the model fitting using ERAint Hs and SLP as: Prepare Pt, trGt, PtOnHs, trGtOnHs slp_get_seasonal.f -- split 1x1 SLP data into 4 seasons slp_get_avg.f -- get base period (1981-2000) mean for SLP fields slp_get_std.f -- get base period standard deviation for SLP fields slp_pt_standardized.f -- standardized SLP 1x1 by deduct base period mean and divided by base period std, get global Pt fields ready pt_surroundings.f -- Prepare Pt's for EOF/PC calculation for each area/season combination slp_get_gt.f -- get global Gt fields from SLP gt_boxcox.f -- get trGt fields from Gt fields gt_get_avg.f -- get base period mean for trGt fields gt_get_std.f -- get base period standard deviation fro trGt fields gt_mask.f -- get regional trGt fields gt_surroundings.f -- get trGt on Hs grid locations get_PtPC.r -- calculate EOF/PCs for Pt fields get_trGtPC.r -- calculate EOF/PCs for trGt fields Prepare Hs for model fitting: getHt.f -- get seasonal Hs from original ERAint Hs data file get_Ht_lambda.r -- get coxbox coefficient for ERAint Hs get_Ht_maxFromStd.r -- get max Hs for each grid point (for all available original ERAint Hs) for each season get_Ht_mean.r -- get Hs base period mean for each season get_Ht_std.r -- get Hs base period standard deviation Model fitting: get_coef_PM_30yr.r -- Model fitting, get all model fitting parameters using Hs/PtPCs/trGtPCs, mean and std fields of Hs/PtPCs/trGtPCs from above steps Part II, apply above model on SLP fields to get projected Hs fields, all codes and data steps locate at Proj/, Proj/readme.txt explains how the data steps be arranged.