subroutine ir00param4(ibasin,lat,rsst,vmax,shdc,delv,t200, + d200,tspd,rmiss,goes00,ierr) c c----------------------------------------------------------------------- c This is a subroutine of the SHIPS model. c c This program calculates the GOES parameters currently used in c SHIPS from various synoptic predictors. This is for use when c GOES satellite imagery is not available. c c c Inputs: ibasin - 1 for Atlantic, 2 for E. and C. Pacific, 3 for W. Pacific c 4 for N. Indian Ocean , 5 for S. Hemisphere c lat - Storm latitude (deg N) at t=0h c rsst - Reynolds SST (deg C) at t=0h c vmax - Maximum winds (kt) at t=0h c shdc - 850-200 mb corrected shear magnitude (kt) at t=0h c delv - Intensity change (kt) from t=-12h to 0h c t200 - 200 to 800km area average 200mb temperature (deg C) at t=0h c d200 - 200mb divergence for r=0-1000km (sec-1*10^7) at t=0h c tspd - Storm translational speed (kt) at t=0h c rmiss - default value for missing data (currently 9999.) c c Output: c goes00 - Array containing the 16 goes00 parameters (1-16) c ierr - error flag =0 for normal completion c =1 one or more predictors missing c c History: c Created April 2009 by A. Schumacher, CIRA/CSU c Updated June 2011 by A. Schumacher to add WP capabilities c Updated June 2013 by A. Schumacher and M. DeMaria to add IO/SH c capabilities c Updated Aug. 2015 by J. Knaff to add V500, R5 and Fr5 predictors to c the goes00 array for wind radii prediction (JHT), c added implicit none c Updated Jan. 2017 by K. Musgrave to fix subroutine name (3-->4) c Updated Feb. 2023 by K. Musgrave for NCO specifications c c----------------------------------------------------------------------- c implicit none c c Number of predictors integer nvar parameter (nvar=9) c c Input predictors real lat,rsst,vmax,shdc,delv,t200,d200,tspd real pvar(nvar) c c Output GOES parameters real goes00(19) c c general usage variables integer i,j,ibasin,ierr real sum01,sum02,sum03,sum04,sum05,sum06,sum07,sum08 real sum09,sum10,sum11,sum12,sum13,sum14,sum15,sum16 real sum17,sum18,sum19 real rmiss real sumv01,sumv02,sumv03,sumv04,sumv05,sumv06,sumv07 real sumv08,sumv09,sumv10,sumv11,sumv12,sumv13,sumv14 real sumv15,sumv16,sumv17,sumv18,sumv19 c c Atlantic arrays real mna01,mna02,mna03,mna04,mna05,mna06,mna07,mna08 real mna09,mna10,mna11,mna12,mna13,mna14,mna15,mna16 real mna17,mna18,mna19 c real sda01,sda02,sda03,sda04,sda05,sda06,sda07,sda08 real sda09,sda10,sda11,sda12,sda13,sda14,sda15,sda16 real sda17,sda18,sda19 c real mna(nvar) ! Predictor means real sda(nvar) ! Predictor standard deviations real cfa01(nvar) ! Normalized predictor coefficients, goes00(1) real cfa02(nvar) ! Normalized predictor coefficients, goes00(2) real cfa03(nvar) ! Normalized predictor coefficients, goes00(3) real cfa04(nvar) ! Normalized predictor coefficients, goes00(4) real cfa05(nvar) ! Normalized predictor coefficients, goes00(5) real cfa06(nvar) ! Normalized predictor coefficients, goes00(6) real cfa07(nvar) ! Normalized predictor coefficients, goes00(7) real cfa08(nvar) ! Normalized predictor coefficients, goes00(8) real cfa09(nvar) ! Normalized predictor coefficients, goes00(9) real cfa10(nvar) ! Normalized predictor coefficients, goes00(10) real cfa11(nvar) ! Normalized predictor coefficients, goes00(11) real cfa12(nvar) ! Normalized predictor coefficients, goes00(12) real cfa13(nvar) ! Normalized predictor coefficients, goes00(13) real cfa14(nvar) ! Normalized predictor coefficients, goes00(14) real cfa15(nvar) ! Normalized predictor coefficients, goes00(15) real cfa16(nvar) ! Normalized predictor coefficients, goes00(16) real cfa17(nvar) ! Normalized predictor coefficients, goes00(17) real cfa18(nvar) ! Normalized predictor coefficients, goes00(18) real cfa19(nvar) ! Normalized predictor coefficients, goes00(19) c c E Pacific arrays real mnp01,mnp02,mnp03,mnp04,mnp05,mnp06,mnp07,mnp08 real mnp09,mnp10,mnp11,mnp12,mnp13,mnp14,mnp15,mnp16 real mnp17,mnp18,mnp19 c real sdp01,sdp02,sdp03,sdp04,sdp05,sdp06,sdp07,sdp08 real sdp09,sdp10,sdp11,sdp12,sdp13,sdp14,sdp15,sdp16 real sdp17,sdp18,sdp19 c real mnp(nvar) ! Predictor means real sdp(nvar) ! Predictor standard deviations real cfp01(nvar) ! Normalized predictor coefficients, goes00(1) real cfp02(nvar) ! Normalized predictor coefficients, goes00(2) real cfp03(nvar) ! Normalized predictor coefficients, goes00(3) real cfp04(nvar) ! Normalized predictor coefficients, goes00(4) real cfp05(nvar) ! Normalized predictor coefficients, goes00(5) real cfp06(nvar) ! Normalized predictor coefficients, goes00(6) real cfp07(nvar) ! Normalized predictor coefficients, goes00(7) real cfp08(nvar) ! Normalized predictor coefficients, goes00(8) real cfp09(nvar) ! Normalized predictor coefficients, goes00(9) real cfp10(nvar) ! Normalized predictor coefficients, goes00(10) real cfp11(nvar) ! Normalized predictor coefficients, goes00(11) real cfp12(nvar) ! Normalized predictor coefficients, goes00(12) real cfp13(nvar) ! Normalized predictor coefficients, goes00(13) real cfp14(nvar) ! Normalized predictor coefficients, goes00(14) real cfp15(nvar) ! Normalized predictor coefficients, goes00(15) real cfp16(nvar) ! Normalized predictor coefficients, goes00(16) real cfp17(nvar) ! Normalized predictor coefficients, goes00(17) real cfp18(nvar) ! Normalized predictor coefficients, goes00(18) real cfp19(nvar) ! Normalized predictor coefficients, goes00(19) c c W Pacific arrays real mnw01,mnw02,mnw03,mnw04,mnw05,mnw06,mnw07,mnw08 real mnw09,mnw10,mnw11,mnw12,mnw13,mnw14,mnw15,mnw16 real mnw17,mnw18,mnw19 c real sdw01,sdw02,sdw03,sdw04,sdw05,sdw06,sdw07,sdw08 real sdw09,sdw10,sdw11,sdw12,sdw13,sdw14,sdw15,sdw16 real sdw17,sdw18,sdw19 c real mnw(nvar) ! Predictor means real sdw(nvar) ! Predictor standard deviations real cfw01(nvar) ! Normalized predictor coefficients, goes00(1) real cfw02(nvar) ! Normalized predictor coefficients, goes00(2) real cfw03(nvar) ! Normalized predictor coefficients, goes00(3) real cfw04(nvar) ! Normalized predictor coefficients, goes00(4) real cfw05(nvar) ! Normalized predictor coefficients, goes00(5) real cfw06(nvar) ! Normalized predictor coefficients, goes00(6) real cfw07(nvar) ! Normalized predictor coefficients, goes00(7) real cfw08(nvar) ! Normalized predictor coefficients, goes00(8) real cfw09(nvar) ! Normalized predictor coefficients, goes00(9) real cfw10(nvar) ! Normalized predictor coefficients, goes00(10) real cfw11(nvar) ! Normalized predictor coefficients, goes00(11) real cfw12(nvar) ! Normalized predictor coefficients, goes00(12) real cfw13(nvar) ! Normalized predictor coefficients, goes00(13) real cfw14(nvar) ! Normalized predictor coefficients, goes00(14) real cfw15(nvar) ! Normalized predictor coefficients, goes00(15) real cfw16(nvar) ! Normalized predictor coefficients, goes00(16) real cfw17(nvar) ! Normalized predictor coefficients, goes00(17) real cfw18(nvar) ! Normalized predictor coefficients, goes00(18) real cfw19(nvar) ! Normalized predictor coefficients, goes00(19) c c Indian Ocean arrays real mni01,mni02,mni03,mni04,mni05,mni06,mni07,mni08 real mni09,mni10,mni11,mni12,mni13,mni14,mni15,mni16 real mni17,mni18,mni19 c real sdi01,sdi02,sdi03,sdi04,sdi05,sdi06,sdi07,sdi08 real sdi09,sdi10,sdi11,sdi12,sdi13,sdi14,sdi15,sdi16 real sdi17,sdi18,sdi19 c real mni(nvar) ! Predictor means real sdi(nvar) ! Predictor standard deviations real cfi01(nvar) ! Normalized predictor coefficients, goes00(1) real cfi02(nvar) ! Normalized predictor coefficients, goes00(2) real cfi03(nvar) ! Normalized predictor coefficients, goes00(3) real cfi04(nvar) ! Normalized predictor coefficients, goes00(4) real cfi05(nvar) ! Normalized predictor coefficients, goes00(5) real cfi06(nvar) ! Normalized predictor coefficients, goes00(6) real cfi07(nvar) ! Normalized predictor coefficients, goes00(7) real cfi08(nvar) ! Normalized predictor coefficients, goes00(8) real cfi09(nvar) ! Normalized predictor coefficients, goes00(9) real cfi10(nvar) ! Normalized predictor coefficients, goes00(10) real cfi11(nvar) ! Normalized predictor coefficients, goes00(11) real cfi12(nvar) ! Normalized predictor coefficients, goes00(12) real cfi13(nvar) ! Normalized predictor coefficients, goes00(13) real cfi14(nvar) ! Normalized predictor coefficients, goes00(14) real cfi15(nvar) ! Normalized predictor coefficients, goes00(15) real cfi16(nvar) ! Normalized predictor coefficients, goes00(16) real cfi17(nvar) ! Normalized predictor coefficients, goes00(17) real cfi18(nvar) ! Normalized predictor coefficients, goes00(18) real cfi19(nvar) ! Normalized predictor coefficients, goes00(19) c c S. Hemispehre arrays real mns01,mns02,mns03,mns04,mns05,mns06,mns07,mns08 real mns09,mns10,mns11,mns12,mns13,mns14,mns15,mns16 real mns17,mns18,mns19 c real sds01,sds02,sds03,sds04,sds05,sds06,sds07,sds08 real sds09,sds10,sds11,sds12,sds13,sds14,sds15,sds16 real sds17,sds18,sds19 c real mns(nvar) ! Predictor means real sds(nvar) ! Predictor standard deviations real cfs01(nvar) ! Normalized predictor coefficients, goes00(1) real cfs02(nvar) ! Normalized predictor coefficients, goes00(2) real cfs03(nvar) ! Normalized predictor coefficients, goes00(3) real cfs04(nvar) ! Normalized predictor coefficients, goes00(4) real cfs05(nvar) ! Normalized predictor coefficients, goes00(5) real cfs06(nvar) ! Normalized predictor coefficients, goes00(6) real cfs07(nvar) ! Normalized predictor coefficients, goes00(7) real cfs08(nvar) ! Normalized predictor coefficients, goes00(8) real cfs09(nvar) ! Normalized predictor coefficients, goes00(9) real cfs10(nvar) ! Normalized predictor coefficients, goes00(10) real cfs11(nvar) ! Normalized predictor coefficients, goes00(11) real cfs12(nvar) ! Normalized predictor coefficients, goes00(12) real cfs13(nvar) ! Normalized predictor coefficients, goes00(13) real cfs14(nvar) ! Normalized predictor coefficients, goes00(14) real cfs15(nvar) ! Normalized predictor coefficients, goes00(15) real cfs16(nvar) ! Normalized predictor coefficients, goes00(16) real cfs17(nvar) ! Normalized predictor coefficients, goes00(17) real cfs18(nvar) ! Normalized predictor coefficients, goes00(18) real cfs19(nvar) ! Normalized predictor coefficients, goes00(19) c c----------------------------------------------------------------------- c ATLANTIC BASIN: c lsdiaga_1982_2012_rean_sat_nbc_ts.dat, 8555 cases w/ IR data c----------------------------------------------------------------------- c Predictor means and standard deviations c LAT RSST VMAX SHDC DELV T200 D200 TSPD VMX^2 data mna / 24.5, 27.5, 55.8, 16.7, 1.9,-53.1, 28.3, 11.1, 3772.4/ data sda / 8.2, 2.0, 25.7, 10.2, 9.6, 1.7, 35.1, 6.4, 3651.5/ c Model for IR00(01), explains 54.0% of the variance: data mna01 /-35.4/ data sda01 / 22.3/ data cfa01 / 0.27,-0.10,-1.30, 0.22,-0.07,-0.03,-0.08,-0.10, 0.83/ c Model for IR00(02), explains 23.7% of the variance: data mna02 / 15.9/ data sda02 / 6.8/ data cfa02 / 0.06, 0.17,-0.25, 0.14,-0.04,-0.13,-0.01,-0.01,-0.16/ c Model for IR00(03), explains 58.6% of the variance: data mna03 /-27.1/ data sda03 / 17.8/ data cfa03 / 0.23,-0.09,-0.94, 0.19,-0.05,-0.09,-0.18,-0.09, 0.44/ c Model for IR00(04), explains 29.3% of the variance: data mna04 / 21.0/ data sda04 / 7.0/ data cfa04 / 0.01, 0.27, 0.02, 0.22, 0.00,-0.17, 0.02, 0.03,-0.45/ c Model for IR00(05), explains 51.1% of the variance: data mna05 / 71.1/ data sda05 / 27.5/ data cfa05 /-0.21, 0.04, 1.30,-0.30, 0.04, 0.04, 0.10, 0.10,-0.85/ c Model for IR00(06), explains 54.1% of the variance: data mna06 / 64.2/ data sda06 / 28.9/ data cfa06 /-0.21, 0.08, 1.27,-0.26, 0.04, 0.02, 0.09, 0.10,-0.77/ c Model for IR00(07), explains 56.2% of the variance: data mna07 / 57.2/ data sda07 / 29.4/ data cfa07 /-0.21, 0.11, 1.18,-0.22, 0.05, 0.00, 0.08, 0.11,-0.65/ c Model for IR00(08), explains 57.6% of the variance: data mna08 / 49.5/ data sda08 / 29.3/ data cfa08 /-0.21, 0.15, 1.05,-0.18, 0.05, 0.00, 0.08, 0.11,-0.48/ c Model for IR00(09), explains 57.1% of the variance: data mna09 / 38.9/ data sda09 / 28.1/ data cfa09 /-0.22, 0.18, 0.79,-0.11, 0.06, 0.01, 0.08, 0.10,-0.22/ c Model for IR00(10), explains 51.0% of the variance: data mna10 / 24.2/ data sda10 / 23.9/ data cfa10 /-0.26, 0.16, 0.34,-0.02, 0.09, 0.06, 0.08, 0.05, 0.17/ c Model for IR00(11), explains 29.4% of the variance: data mna11 /-33.9/ data sda11 / 31.2/ data cfa11 / 0.28,-0.11,-1.61, 0.15,-0.05,-0.01,-0.05,-0.11, 1.48/ c Model for IR00(12), explains 32.2% of the variance: data mna12 /-38.0/ data sda12 / 30.2/ data cfa12 / 0.29,-0.11,-1.52, 0.17,-0.06, 0.00,-0.06,-0.11, 1.32/ c Model for IR00(13), explains 4.7% of the variance: data mna13 / 15.5/ data sda13 / 11.4/ data cfa13 /-0.07, 0.05, 0.59,-0.04, 0.01,-0.02, 0.03, 0.04,-0.72/ c Model for IR00(14), explains 47.9% of the variance: data mna14 /-46.3/ data sda14 / 25.2/ data cfa14 / 0.27,-0.10,-1.35, 0.22,-0.07, 0.01,-0.05,-0.09, 0.91/ c Model for IR00(15), explains 50.6% of the variance: data mna15 /-38.4/ data sda15 / 24.9/ data cfa15 / 0.28,-0.09,-1.31, 0.21,-0.07,-0.01,-0.05,-0.10, 0.86/ c Model for IR00(16), explains 10.7% of the variance: data mna16 / 57.5/ data sda16 / 36.2/ data cfa16 / 0.18,-0.08,-1.00, 0.10,-0.04, 0.06, 0.00,-0.06, 0.98/ c Model for IR00(17), explains 50.0% of the variance: data mna17 / 7.1/ data sda17 / 2.4/ data cfa17 / 0.56,-0.07, 0.07,-0.05, 0.03, 0.13, 0.29, 0.00, 0.15/ c Model for IR00(18), explains 50.0% of the variance: data mna18 / 11.9/ data sda18 / 3.8/ data cfa18 / 0.56,-0.07, 0.07,-0.05, 0.03, 0.13, 0.29, 0.00, 0.15/ c Model for IR00(19), explains 40.1% of the variance: data mna19 / 1.0/ data sda19 / 0.3/ data cfa19 / 0.61,-0.07,-0.75,-0.05, 0.04, 0.13, 0.31, 0.00, 0.66/ c----------------------------------------------------------------------- c EAST PACIFIC BASIN: c lsdiage_1982_2012_rean_sat_nbc_ts.dat, 10254 cases w/ IR data c----------------------------------------------------------------------- c Predictor means and standard deviations c LAT RSST VMAX SHDC DELV T200 D200 TSPD VMX^2 data mnp / 16.7, 27.3, 53.8, 13.1, 0.2,-52.3, 23.6, 8.9, 3608.5/ data sdp / 3.7, 1.9, 26.7, 8.0, 11.2, 1.0, 32.8, 4.1, 3707.2/ c Model for IR00(01), explains 65.1% of the variance: data mnp01 /-38.1/ data sdp01 / 24.2/ data cfp01 / 0.14,-0.37,-1.40, 0.12,-0.14,-0.02,-0.11,-0.08, 1.04/ c Model for IR00(02), explains 14.4% of the variance: data mnp02 / 13.6/ data sdp02 / 6.3/ data cfp02 /-0.01, 0.23, 0.03, 0.02,-0.16,-0.10,-0.05,-0.03,-0.29/ c Model for IR00(03), explains 64.3% of the variance: data mnp03 /-26.4/ data sdp03 / 20.0/ data cfp03 / 0.14,-0.35,-0.94, 0.07,-0.14,-0.07,-0.22,-0.10, 0.62/ c Model for IR00(04), explains 20.0% of the variance: data mnp04 / 17.7/ data sdp04 / 6.5/ data cfp04 /-0.03, 0.35, 0.37, 0.11,-0.11,-0.12,-0.01,-0.02,-0.58/ c Model for IR00(05), explains 62.1% of the variance: data mnp05 / 72.1/ data sdp05 / 31.2/ data cfp05 /-0.13, 0.36, 1.39,-0.16, 0.08, 0.03, 0.12, 0.08,-1.02/ c Model for IR00(06), explains 64.1% of the variance: data mnp06 / 65.1/ data sdp06 / 32.3/ data cfp06 /-0.14, 0.35, 1.26,-0.12, 0.12, 0.03, 0.13, 0.08,-0.86/ c Model for IR00(07), explains 64.5% of the variance: data mnp07 / 57.9/ data sdp07 / 32.3/ data cfp07 /-0.14, 0.34, 1.12,-0.09, 0.15, 0.02, 0.13, 0.09,-0.71/ c Model for IR00(08), explains 63.6% of the variance: data mnp08 / 50.1/ data sdp08 / 31.4/ data cfp08 /-0.13, 0.33, 0.96,-0.06, 0.17, 0.02, 0.14, 0.09,-0.55/ c Model for IR00(09), explains 59.3% of the variance: data mnp09 / 39.7/ data sdp09 / 29.0/ data cfp09 /-0.12, 0.32, 0.72,-0.01, 0.19, 0.02, 0.15, 0.09,-0.32/ c Model for IR00(10), explains 48.1% of the variance: data mnp10 / 25.6/ data sdp10 / 23.9/ data cfp10 /-0.10, 0.29, 0.41, 0.03, 0.21, 0.03, 0.15, 0.09,-0.09/ c Model for IR00(11), explains 44.5% of the variance: data mnp11 /-39.7/ data sdp11 / 32.2/ data cfp11 / 0.13,-0.35,-1.81, 0.11,-0.10, 0.00,-0.06,-0.06, 1.73/ c Model for IR00(12), explains 47.8% of the variance: data mnp12 /-43.4/ data sdp12 / 31.3/ data cfp12 / 0.13,-0.36,-1.72, 0.12,-0.10, 0.01,-0.07,-0.06, 1.57/ c Model for IR00(13), explains 7.3% of the variance: data mnp13 / 15.1/ data sdp13 / 10.7/ data cfp13 /-0.02, 0.15, 0.59,-0.06,-0.03, 0.00, 0.00, 0.01,-0.77/ c Model for IR00(14), explains 60.8% of the variance: data mnp14 /-50.1/ data sdp14 / 27.6/ data cfp14 / 0.14,-0.37,-1.49, 0.14,-0.11, 0.01,-0.07,-0.06, 1.12/ c Model for IR00(15), explains 63.1% of the variance: data mnp15 /-42.6/ data sdp15 / 26.7/ data cfp15 / 0.14,-0.36,-1.44, 0.13,-0.13, 0.00,-0.08,-0.07, 1.06/ c Model for IR00(16), explains 18.9% of the variance: data mnp16 / 52.6/ data sdp16 / 33.3/ data cfp16 / 0.08,-0.25,-1.19, 0.12,-0.04, 0.05, 0.01,-0.01, 1.15/ c Model for IR00(17), explains 40.6% of the variance: data mnp17 / 5.5/ data sdp17 / 2.1/ data cfp17 / 0.27, 0.20, 0.17, 0.02, 0.14, 0.11, 0.42, 0.11,-0.03/ c Model for IR00(18), explains 40.6% of the variance: data mnp18 / 9.3/ data sdp18 / 3.3/ data cfp18 / 0.27, 0.20, 0.17, 0.02, 0.14, 0.11, 0.42, 0.11,-0.03/ c Model for IR00(19), explains 37.6% of the variance: data mnp19 / 0.8/ data sdp19 / 0.3/ data cfp19 / 0.27, 0.21,-0.47, 0.02, 0.13, 0.11, 0.43, 0.12, 0.32/ c----------------------------------------------------------------------- c WEST PACIFIC BASIN: c lsdiagw_2000_2012_rean_sat_nbc_ts.dat, 2824 cases w/ IR data c----------------------------------------------------------------------- c Predictor means and standard deviations c LAT RSST VMAX SHDC DELV T200 D200 TSPD VMX^2 data mnw / 19.5, 28.9, 61.0, 13.8, 1.4,-50.5, 49.6, 10.0, 4690.7/ data sdw / 6.2, 1.1, 31.2, 8.0, 12.4, 1.6, 38.9, 5.0, 4697.1/ c Model for IR00(01), explains 45.6% of the variance: data mnw01 /-52.0/ data sdw01 / 18.8/ data cfw01 / 0.17,-0.24,-1.74, 0.13,-0.10, 0.09,-0.06,-0.08, 1.32/ c Model for IR00(02), explains 35.2% of the variance: data mnw02 / 14.0/ data sdw02 / 6.9/ data cfw02 / 0.08, 0.06,-1.29, 0.17,-0.06, 0.04,-0.01,-0.04, 0.81/ c Model for IR00(03), explains 53.2% of the variance: data mnw03 /-43.5/ data sdw03 / 16.1/ data cfw03 / 0.20,-0.22,-1.39, 0.15,-0.06,-0.09,-0.18,-0.07, 0.91/ c Model for IR00(04), explains 44.6% of the variance: data mnw04 / 19.2/ data sdw04 / 7.6/ data cfw04 / 0.09, 0.13,-0.93, 0.27,-0.01,-0.06,-0.05,-0.02, 0.40/ c Model for IR00(05), explains 41.7% of the variance: data mnw05 / 87.6/ data sdw05 / 19.2/ data cfw05 /-0.09, 0.15, 1.65,-0.29, 0.01, 0.02, 0.12, 0.08,-1.28/ c Model for IR00(06), explains 45.9% of the variance: data mnw06 / 82.6/ data sdw06 / 22.1/ data cfw06 /-0.10, 0.18, 1.69,-0.24, 0.03,-0.02, 0.11, 0.08,-1.25/ c Model for IR00(07), explains 49.2% of the variance: data mnw07 / 76.8/ data sdw07 / 24.4/ data cfw07 /-0.10, 0.20, 1.65,-0.19, 0.05,-0.04, 0.10, 0.08,-1.15/ c Model for IR00(08), explains 52.0% of the variance: data mnw08 / 70.4/ data sdw08 / 26.2/ data cfw08 /-0.10, 0.22, 1.58,-0.14, 0.07,-0.05, 0.09, 0.07,-1.04/ c Model for IR00(09), explains 53.8% of the variance: data mnw09 / 61.5/ data sdw09 / 27.6/ data cfw09 /-0.12, 0.25, 1.42,-0.08, 0.10,-0.06, 0.07, 0.06,-0.84/ c Model for IR00(10), explains 51.6% of the variance: data mnw10 / 47.7/ data sdw10 / 27.6/ data cfw10 /-0.17, 0.26, 1.06,-0.01, 0.13,-0.03, 0.05, 0.04,-0.49/ c Model for IR00(11), explains 26.5% of the variance: data mnw11 /-48.2/ data sdw11 / 29.6/ data cfw11 / 0.13,-0.21,-1.86, 0.05,-0.04, 0.18,-0.02,-0.08, 1.90/ c Model for IR00(12), explains 26.1% of the variance: data mnw12 /-53.0/ data sdw12 / 27.0/ data cfw12 / 0.13,-0.22,-1.80, 0.06,-0.06, 0.19,-0.01,-0.08, 1.77/ c Model for IR00(13), explains 7.4% of the variance: data mnw13 / 15.8/ data sdw13 / 11.8/ data cfw13 / 0.05, 0.09, 0.23, 0.00, 0.00,-0.09, 0.05, 0.01,-0.46/ c Model for IR00(14), explains 40.1% of the variance: data mnw14 /-62.5/ data sdw14 / 20.1/ data cfw14 / 0.12,-0.25,-1.57, 0.13,-0.08, 0.13,-0.05,-0.08, 1.17/ c Model for IR00(15), explains 41.5% of the variance: data mnw15 /-55.2/ data sdw15 / 21.0/ data cfw15 / 0.16,-0.23,-1.68, 0.12,-0.11, 0.13,-0.03,-0.08, 1.28/ c Model for IR00(16), explains 11.9% of the variance: data mnw16 / 54.4/ data sdw16 / 34.2/ data cfw16 / 0.03,-0.12,-0.87, 0.01,-0.02, 0.20, 0.03,-0.03, 1.01/ c Model for IR00(17), explains 54.1% of the variance: data mnw17 / 7.6/ data sdw17 / 2.0/ data cfw17 / 0.39, 0.10, 0.56,-0.04, 0.05, 0.29, 0.32, 0.01,-0.28/ c Model for IR00(18), explains 54.1% of the variance: data mnw18 / 12.7/ data sdw18 / 3.2/ data cfw18 / 0.39, 0.10, 0.56,-0.04, 0.05, 0.29, 0.32, 0.01,-0.28/ c Model for IR00(19), explains 41.1% of the variance: data mnw19 / 1.1/ data sdw19 / 0.3/ data cfw19 / 0.44, 0.12,-0.42,-0.04, 0.07, 0.31, 0.36, 0.01, 0.22/ c----------------------------------------------------------------------- c INDIAN OCEAN: c lsdiagi_1998_2012_rean_sat_nbc_ts.dat, 749 cases w/ IR data c----------------------------------------------------------------------- c Predictor means and standard deviations c LAT RSST VMAX SHDC DELV T200 D200 TSPD VMX^2 data mni / 14.7, 29.1, 48.5, 16.2, 2.5,-51.3, 50.9, 6.8, 3019.8/ data sdi / 4.7, 1.1, 25.9, 8.1, 11.2, 1.8, 43.0, 3.6, 3754.3/ c Model for IR00(01), explains 42.5% of the variance: data mni01 /-49.9/ data sdi01 / 21.8/ data cfi01 / 0.09,-0.22,-1.19, 0.16,-0.15, 0.00,-0.17,-0.12, 0.82/ c Model for IR00(02), explains 36.7% of the variance: data mni02 / 14.8/ data sdi02 / 7.5/ data cfi02 /-0.09,-0.06,-0.98, 0.27,-0.06, 0.04,-0.02,-0.14, 0.61/ c Model for IR00(03), explains 47.8% of the variance: data mni03 /-40.1/ data sdi03 / 17.9/ data cfi03 / 0.15,-0.19,-0.79, 0.15,-0.17,-0.13,-0.28,-0.12, 0.42/ c Model for IR00(04), explains 38.2% of the variance: data mni04 / 20.8/ data sdi04 / 7.7/ data cfi04 /-0.06, 0.05,-0.71, 0.35,-0.06,-0.03,-0.05,-0.12, 0.34/ c Model for IR00(05), explains 37.7% of the variance: data mni05 / 82.6/ data sdi05 / 23.3/ data cfi05 /-0.03, 0.11, 1.07,-0.27, 0.10, 0.09, 0.27, 0.06,-0.77/ c Model for IR00(06), explains 41.5% of the variance: data mni06 / 76.6/ data sdi06 / 25.5/ data cfi06 /-0.05, 0.14, 1.10,-0.26, 0.11, 0.07, 0.24, 0.08,-0.75/ c Model for IR00(07), explains 43.8% of the variance: data mni07 / 70.5/ data sdi07 / 26.9/ data cfi07 /-0.06, 0.16, 1.09,-0.23, 0.13, 0.05, 0.21, 0.09,-0.70/ c Model for IR00(08), explains 45.1% of the variance: data mni08 / 64.0/ data sdi08 / 27.9/ data cfi08 /-0.07, 0.18, 1.05,-0.19, 0.14, 0.03, 0.19, 0.10,-0.62/ c Model for IR00(09), explains 46.2% of the variance: data mni09 / 55.3/ data sdi09 / 28.5/ data cfi09 /-0.10, 0.20, 0.94,-0.13, 0.17, 0.01, 0.15, 0.11,-0.48/ c Model for IR00(10), explains 44.2% of the variance: data mni10 / 43.1/ data sdi10 / 27.9/ data cfi10 /-0.14, 0.22, 0.72,-0.05, 0.19, 0.00, 0.12, 0.13,-0.27/ c Model for IR00(11), explains 27.0% of the variance: data mni11 /-52.1/ data sdi11 / 29.5/ data cfi11 / 0.06,-0.25,-1.58, 0.12,-0.09, 0.08,-0.09,-0.12, 1.42/ c Model for IR00(12), explains 27.9% of the variance: data mni12 /-55.8/ data sdi12 / 28.4/ data cfi12 / 0.07,-0.25,-1.40, 0.12,-0.11, 0.08,-0.10,-0.12, 1.19/ c Model for IR00(13), explains 6.3% of the variance: data mni13 / 17.7/ data sdi13 / 11.8/ data cfi13 / 0.06, 0.09, 0.59, 0.04,-0.02,-0.16, 0.02, 0.08,-0.76/ c Model for IR00(14), explains 34.5% of the variance: data mni14 /-61.0/ data sdi14 / 23.5/ data cfi14 / 0.07,-0.24,-1.17, 0.14,-0.14, 0.04,-0.15,-0.10, 0.86/ c Model for IR00(15), explains 40.2% of the variance: data mni15 /-53.4/ data sdi15 / 23.9/ data cfi15 / 0.09,-0.23,-1.23, 0.16,-0.14, 0.02,-0.14,-0.11, 0.86/ c Model for IR00(16), explains 11.0% of the variance: data mni16 / 49.3/ data sdi16 / 33.6/ data cfi16 / 0.09,-0.17,-1.20, 0.02,-0.06, 0.11, 0.03,-0.07, 1.19/ c Model for IR00(17), explains 51.7% of the variance: data mni17 / 6.6/ data sdi17 / 2.0/ data cfi17 / 0.32, 0.11, 0.28,-0.02, 0.19, 0.30, 0.39, 0.02,-0.09/ c Model for IR00(18), explains 51.7% of the variance: data mni18 / 11.1/ data sdi18 / 3.2/ data cfi18 / 0.32, 0.11, 0.28,-0.02, 0.19, 0.30, 0.39, 0.02,-0.09/ c Model for IR00(19), explains 45.6% of the variance: data mni19 / 1.0/ data sdi19 / 0.3/ data cfi19 / 0.35, 0.12,-0.43,-0.02, 0.19, 0.32, 0.43, 0.02, 0.26/ c----------------------------------------------------------------------- c S HEMISPHERE: c lsdiags_1999_2012_rean_sat_nbc_ts.dat, 4752 cases w/ IR data c----------------------------------------------------------------------- c Predictor means and standard deviations c LAT RSST VMAX SHDC DELV T200 D200 TSPD VMX^2 data mns /-16.9, 28.2, 54.0, 14.4, 1.2,-51.5, 45.9, 8.4, 3661.5/ data sds / 4.7, 1.4, 27.4, 8.8, 11.9, 1.3, 39.4, 4.9, 3837.9/ c Model for IR00(01), explains 53.9% of the variance: data mns01 /-48.3/ data sds01 / 22.3/ data cfs01 /-0.12,-0.25,-1.49, 0.25,-0.11,-0.04,-0.11,-0.09, 1.13/ c Model for IR00(02), explains 22.7% of the variance: data mns02 / 14.1/ data sds02 / 6.6/ data cfs02 / 0.00,-0.01,-0.69, 0.06,-0.09,-0.01,-0.03,-0.03, 0.27/ c Model for IR00(03), explains 56.5% of the variance: data mns03 /-37.8/ data sds03 / 18.1/ data cfs03 /-0.13,-0.26,-1.16, 0.22,-0.10,-0.11,-0.22,-0.06, 0.81/ c Model for IR00(04), explains 23.8% of the variance: data mns04 / 19.6/ data sds04 / 7.0/ data cfs04 / 0.04, 0.06,-0.40, 0.13,-0.08,-0.04,-0.04,-0.01,-0.01/ c Model for IR00(05), explains 53.4% of the variance: data mns05 / 82.3/ data sds05 / 24.8/ data cfs05 / 0.04, 0.22, 1.45,-0.36, 0.06, 0.03, 0.16, 0.06,-1.12/ c Model for IR00(06), explains 55.8% of the variance: data mns06 / 76.6/ data sds06 / 26.9/ data cfs06 / 0.06, 0.23, 1.46,-0.31, 0.08, 0.02, 0.16, 0.07,-1.08/ c Model for IR00(07), explains 56.8% of the variance: data mns07 / 70.3/ data sds07 / 28.2/ data cfs07 / 0.08, 0.23, 1.43,-0.27, 0.10, 0.02, 0.15, 0.07,-1.01/ c Model for IR00(08), explains 56.4% of the variance: data mns08 / 63.3/ data sds08 / 28.7/ data cfs08 / 0.10, 0.23, 1.34,-0.22, 0.12, 0.02, 0.14, 0.07,-0.89/ c Model for IR00(09), explains 53.1% of the variance: data mns09 / 53.6/ data sds09 / 28.6/ data cfs09 / 0.13, 0.22, 1.17,-0.16, 0.14, 0.03, 0.12, 0.08,-0.71/ c Model for IR00(10), explains 44.0% of the variance: data mns10 / 39.7/ data sds10 / 26.9/ data cfs10 / 0.17, 0.22, 0.82,-0.08, 0.15, 0.07, 0.08, 0.08,-0.40/ c Model for IR00(11), explains 33.3% of the variance: data mns11 /-48.7/ data sds11 / 30.6/ data cfs11 /-0.14,-0.25,-1.74, 0.19,-0.06,-0.01,-0.08,-0.11, 1.65/ c Model for IR00(12), explains 36.9% of the variance: data mns12 /-53.0/ data sds12 / 29.1/ data cfs12 /-0.13,-0.25,-1.65, 0.22,-0.08, 0.00,-0.08,-0.11, 1.49/ c Model for IR00(13), explains 7.8% of the variance: data mns13 / 16.0/ data sds13 / 12.0/ data cfs13 / 0.04, 0.08, 0.62,-0.06,-0.02,-0.01, 0.03, 0.04,-0.83/ c Model for IR00(14), explains 48.7% of the variance: data mns14 /-59.7/ data sds14 / 24.5/ data cfs14 /-0.10,-0.24,-1.47, 0.27,-0.11,-0.01,-0.08,-0.09, 1.13/ c Model for IR00(15), explains 52.0% of the variance: data mns15 /-52.2/ data sds15 / 24.6/ data cfs15 /-0.11,-0.24,-1.48, 0.26,-0.12,-0.02,-0.08,-0.10, 1.10/ c Model for IR00(16), explains 13.1% of the variance: data mns16 / 50.5/ data sds16 / 32.7/ data cfs16 /-0.08,-0.14,-1.16, 0.14,-0.04, 0.07,-0.03,-0.08, 1.15/ c Model for IR00(17), explains 50.0% of the variance: data mns17 / 6.7/ data sds17 / 1.9/ data cfs17 /-0.39, 0.18, 0.49,-0.11, 0.09, 0.20, 0.38, 0.01,-0.27/ c Model for IR00(18), explains 50.0% of the variance: data mns18 / 11.3/ data sds18 / 3.1/ data cfs18 /-0.39, 0.18, 0.49,-0.11, 0.09, 0.20, 0.38, 0.01,-0.27/ c Model for IR00(19), explains 40.1% of the variance: data mns19 / 1.0/ data sds19 / 0.3/ data cfs19 /-0.42, 0.20,-0.30,-0.13, 0.09, 0.21, 0.41, 0.01, 0.11/ c c----------------------------------------------------------------------- c Scale factors for GOES predictors real sfac(19) data sfac /10.0,10.0,10.0,10.0, 1.0, 1.0, 1.0, 1.0, 1.0, 1.0, + 10.0,10.0, 1.0,10.0,10.0, 1.0, 10.0, 10.0, 100.0/ c c Max/Min caps - constrains predicted parameters to reasonable c values based on parameter properties and observed max/min real cmax(19) data cmax / 30.0, 45.0, 30.0, 45.0,100.0,100.0,100.0,100.0, + 100.0,100.0, 30.0, 30.0, 30.0, 30.0, 30.0,120.0, + 25.0,25.0, 3.0 / real cmin(19) data cmin /-90.0, 0.0,-90.0, 0.0, 0.0, 0.0, 0.0, 0.0, + 0.0, 0.0,-90.0,-90.0, 0.0,-90.0,-90.0, 20.0, + 2.0, 4.0, 0.3 / c c----------------------------------------------------------------------- c ierr = 0 c pvar(1) = lat pvar(2) = rsst pvar(3) = vmax pvar(4) = shdc pvar(5) = delv pvar(6) = t200 pvar(7) = d200 pvar(8) = tspd pvar(9) = vmax*vmax c do i=1,nvar if (pvar(i) .eq. rmiss) then c Error processing ierr = 1 do j=1,19 goes00(j) = rmiss enddo return endif enddo sumv01 = 0.0 sumv02 = 0.0 sumv03 = 0.0 sumv04 = 0.0 sumv05 = 0.0 sumv06 = 0.0 sumv07 = 0.0 sumv08 = 0.0 sumv09 = 0.0 sumv10 = 0.0 sumv11 = 0.0 sumv12 = 0.0 sumv13 = 0.0 sumv14 = 0.0 sumv15 = 0.0 sumv16 = 0.0 sumv17 = 0.0 sumv18 = 0.0 sumv19 = 0.0 if (ibasin .eq. 1) then do i=1,nvar sumv01 = sumv01 + cfa01(i)*((pvar(i) - mna(i))/sda(i)) sumv02 = sumv02 + cfa02(i)*((pvar(i) - mna(i))/sda(i)) sumv03 = sumv03 + cfa03(i)*((pvar(i) - mna(i))/sda(i)) sumv04 = sumv04 + cfa04(i)*((pvar(i) - mna(i))/sda(i)) sumv05 = sumv05 + cfa05(i)*((pvar(i) - mna(i))/sda(i)) sumv06 = sumv06 + cfa06(i)*((pvar(i) - mna(i))/sda(i)) sumv07 = sumv07 + cfa07(i)*((pvar(i) - mna(i))/sda(i)) sumv08 = sumv08 + cfa08(i)*((pvar(i) - mna(i))/sda(i)) sumv09 = sumv09 + cfa09(i)*((pvar(i) - mna(i))/sda(i)) sumv10 = sumv10 + cfa10(i)*((pvar(i) - mna(i))/sda(i)) sumv11 = sumv11 + cfa11(i)*((pvar(i) - mna(i))/sda(i)) sumv12 = sumv12 + cfa12(i)*((pvar(i) - mna(i))/sda(i)) sumv13 = sumv13 + cfa13(i)*((pvar(i) - mna(i))/sda(i)) sumv14 = sumv14 + cfa14(i)*((pvar(i) - mna(i))/sda(i)) sumv15 = sumv15 + cfa15(i)*((pvar(i) - mna(i))/sda(i)) sumv16 = sumv16 + cfa16(i)*((pvar(i) - mna(i))/sda(i)) sumv17 = sumv17 + cfa17(i)*((pvar(i) - mna(i))/sda(i)) sumv18 = sumv18 + cfa18(i)*((pvar(i) - mna(i))/sda(i)) sumv19 = sumv19 + cfa19(i)*((pvar(i) - mna(i))/sda(i)) enddo goes00(1) = mna01 + (sumv01*sda01) goes00(2) = mna02 + (sumv02*sda02) goes00(3) = mna03 + (sumv03*sda03) goes00(4) = mna04 + (sumv04*sda04) goes00(5) = mna05 + (sumv05*sda05) goes00(6) = mna06 + (sumv06*sda06) goes00(7) = mna07 + (sumv07*sda07) goes00(8) = mna08 + (sumv08*sda08) goes00(9) = mna09 + (sumv09*sda09) goes00(10) = mna10 + (sumv10*sda10) goes00(11) = mna11 + (sumv11*sda11) goes00(12) = mna12 + (sumv12*sda12) goes00(13) = mna13 + (sumv13*sda13) goes00(14) = mna14 + (sumv14*sda14) goes00(15) = mna15 + (sumv15*sda15) goes00(16) = mna16 + (sumv16*sda16) goes00(17) = mna17 + (sumv17*sda17) goes00(18) = mna18 + (sumv18*sda18) goes00(19) = mna19 + (sumv19*sda19) elseif (ibasin .eq. 2) then do i=1,nvar sumv01 = sumv01 + cfp01(i)*((pvar(i) - mnp(i))/sdp(i)) sumv02 = sumv02 + cfp02(i)*((pvar(i) - mnp(i))/sdp(i)) sumv03 = sumv03 + cfp03(i)*((pvar(i) - mnp(i))/sdp(i)) sumv04 = sumv04 + cfp04(i)*((pvar(i) - mnp(i))/sdp(i)) sumv05 = sumv05 + cfp05(i)*((pvar(i) - mnp(i))/sdp(i)) sumv06 = sumv06 + cfp06(i)*((pvar(i) - mnp(i))/sdp(i)) sumv07 = sumv07 + cfp07(i)*((pvar(i) - mnp(i))/sdp(i)) sumv08 = sumv08 + cfp08(i)*((pvar(i) - mnp(i))/sdp(i)) sumv09 = sumv09 + cfp09(i)*((pvar(i) - mnp(i))/sdp(i)) sumv10 = sumv10 + cfp10(i)*((pvar(i) - mnp(i))/sdp(i)) sumv11 = sumv11 + cfp11(i)*((pvar(i) - mnp(i))/sdp(i)) sumv12 = sumv12 + cfp12(i)*((pvar(i) - mnp(i))/sdp(i)) sumv13 = sumv13 + cfp13(i)*((pvar(i) - mnp(i))/sdp(i)) sumv14 = sumv14 + cfp14(i)*((pvar(i) - mnp(i))/sdp(i)) sumv15 = sumv15 + cfp15(i)*((pvar(i) - mnp(i))/sdp(i)) sumv16 = sumv16 + cfp16(i)*((pvar(i) - mnp(i))/sdp(i)) sumv17 = sumv17 + cfp17(i)*((pvar(i) - mnp(i))/sdp(i)) sumv18 = sumv18 + cfp18(i)*((pvar(i) - mnp(i))/sdp(i)) sumv19 = sumv19 + cfp19(i)*((pvar(i) - mnp(i))/sdp(i)) enddo goes00(1) = mnp01 + (sumv01*sdp01) goes00(2) = mnp02 + (sumv02*sdp02) goes00(3) = mnp03 + (sumv03*sdp03) goes00(4) = mnp04 + (sumv04*sdp04) goes00(5) = mnp05 + (sumv05*sdp05) goes00(6) = mnp06 + (sumv06*sdp06) goes00(7) = mnp07 + (sumv07*sdp07) goes00(8) = mnp08 + (sumv08*sdp08) goes00(9) = mnp09 + (sumv09*sdp09) goes00(10) = mnp10 + (sumv10*sdp10) goes00(11) = mnp11 + (sumv11*sdp11) goes00(12) = mnp12 + (sumv12*sdp12) goes00(13) = mnp13 + (sumv13*sdp13) goes00(14) = mnp14 + (sumv14*sdp14) goes00(15) = mnp15 + (sumv15*sdp15) goes00(16) = mnp16 + (sumv16*sdp16) goes00(17) = mnp17 + (sumv17*sdp17) goes00(18) = mnp18 + (sumv18*sdp18) goes00(19) = mnp19 + (sumv19*sdp19) elseif (ibasin .eq. 3) then do i=1,nvar sumv01 = sumv01 + cfw01(i)*((pvar(i) - mnw(i))/sdw(i)) sumv02 = sumv02 + cfw02(i)*((pvar(i) - mnw(i))/sdw(i)) sumv03 = sumv03 + cfw03(i)*((pvar(i) - mnw(i))/sdw(i)) sumv04 = sumv04 + cfw04(i)*((pvar(i) - mnw(i))/sdw(i)) sumv05 = sumv05 + cfw05(i)*((pvar(i) - mnw(i))/sdw(i)) sumv06 = sumv06 + cfw06(i)*((pvar(i) - mnw(i))/sdw(i)) sumv07 = sumv07 + cfw07(i)*((pvar(i) - mnw(i))/sdw(i)) sumv08 = sumv08 + cfw08(i)*((pvar(i) - mnw(i))/sdw(i)) sumv09 = sumv09 + cfw09(i)*((pvar(i) - mnw(i))/sdw(i)) sumv10 = sumv10 + cfw10(i)*((pvar(i) - mnw(i))/sdw(i)) sumv11 = sumv11 + cfw11(i)*((pvar(i) - mnw(i))/sdw(i)) sumv12 = sumv12 + cfw12(i)*((pvar(i) - mnw(i))/sdw(i)) sumv13 = sumv13 + cfw13(i)*((pvar(i) - mnw(i))/sdw(i)) sumv14 = sumv14 + cfw14(i)*((pvar(i) - mnw(i))/sdw(i)) sumv15 = sumv15 + cfw15(i)*((pvar(i) - mnw(i))/sdw(i)) sumv16 = sumv16 + cfw16(i)*((pvar(i) - mnw(i))/sdw(i)) sumv17 = sumv17 + cfw17(i)*((pvar(i) - mnw(i))/sdw(i)) sumv18 = sumv18 + cfw18(i)*((pvar(i) - mnw(i))/sdw(i)) sumv19 = sumv19 + cfw19(i)*((pvar(i) - mnw(i))/sdw(i)) enddo goes00(1) = mnw01 + (sumv01*sdw01) goes00(2) = mnw02 + (sumv02*sdw02) goes00(3) = mnw03 + (sumv03*sdw03) goes00(4) = mnw04 + (sumv04*sdw04) goes00(5) = mnw05 + (sumv05*sdw05) goes00(6) = mnw06 + (sumv06*sdw06) goes00(7) = mnw07 + (sumv07*sdw07) goes00(8) = mnw08 + (sumv08*sdw08) goes00(9) = mnw09 + (sumv09*sdw09) goes00(10) = mnw10 + (sumv10*sdw10) goes00(11) = mnw11 + (sumv11*sdw11) goes00(12) = mnw12 + (sumv12*sdw12) goes00(13) = mnw13 + (sumv13*sdw13) goes00(14) = mnw14 + (sumv14*sdw14) goes00(15) = mnw15 + (sumv15*sdw15) goes00(16) = mnw16 + (sumv16*sdw16) goes00(17) = mnw17 + (sumv17*sdw17) goes00(18) = mnw18 + (sumv18*sdw18) goes00(19) = mnw19 + (sumv19*sdw19) elseif (ibasin .eq. 4) then do i=1,nvar sumv01 = sumv01 + cfi01(i)*((pvar(i) - mni(i))/sdi(i)) sumv02 = sumv02 + cfi02(i)*((pvar(i) - mni(i))/sdi(i)) sumv03 = sumv03 + cfi03(i)*((pvar(i) - mni(i))/sdi(i)) sumv04 = sumv04 + cfi04(i)*((pvar(i) - mni(i))/sdi(i)) sumv05 = sumv05 + cfi05(i)*((pvar(i) - mni(i))/sdi(i)) sumv06 = sumv06 + cfi06(i)*((pvar(i) - mni(i))/sdi(i)) sumv07 = sumv07 + cfi07(i)*((pvar(i) - mni(i))/sdi(i)) sumv08 = sumv08 + cfi08(i)*((pvar(i) - mni(i))/sdi(i)) sumv09 = sumv09 + cfi09(i)*((pvar(i) - mni(i))/sdi(i)) sumv10 = sumv10 + cfi10(i)*((pvar(i) - mni(i))/sdi(i)) sumv11 = sumv11 + cfi11(i)*((pvar(i) - mni(i))/sdi(i)) sumv12 = sumv12 + cfi12(i)*((pvar(i) - mni(i))/sdi(i)) sumv13 = sumv13 + cfi13(i)*((pvar(i) - mni(i))/sdi(i)) sumv14 = sumv14 + cfi14(i)*((pvar(i) - mni(i))/sdi(i)) sumv15 = sumv15 + cfi15(i)*((pvar(i) - mni(i))/sdi(i)) sumv16 = sumv16 + cfi16(i)*((pvar(i) - mni(i))/sdi(i)) sumv17 = sumv17 + cfi17(i)*((pvar(i) - mni(i))/sdi(i)) sumv18 = sumv18 + cfi18(i)*((pvar(i) - mni(i))/sdi(i)) sumv19 = sumv19 + cfi19(i)*((pvar(i) - mni(i))/sdi(i)) enddo goes00(1) = mni01 + (sumv01*sdi01) goes00(2) = mni02 + (sumv02*sdi02) goes00(3) = mni03 + (sumv03*sdi03) goes00(4) = mni04 + (sumv04*sdi04) goes00(5) = mni05 + (sumv05*sdi05) goes00(6) = mni06 + (sumv06*sdi06) goes00(7) = mni07 + (sumv07*sdi07) goes00(8) = mni08 + (sumv08*sdi08) goes00(9) = mni09 + (sumv09*sdi09) goes00(10) = mni10 + (sumv10*sdi10) goes00(11) = mni11 + (sumv11*sdi11) goes00(12) = mni12 + (sumv12*sdi12) goes00(13) = mni13 + (sumv13*sdi13) goes00(14) = mni14 + (sumv14*sdi14) goes00(15) = mni15 + (sumv15*sdi15) goes00(16) = mni16 + (sumv16*sdi16) goes00(17) = mni17 + (sumv17*sdi17) goes00(18) = mni18 + (sumv18*sdi18) goes00(19) = mni19 + (sumv19*sdi19) elseif (ibasin .eq. 5) then do i=1,nvar sumv01 = sumv01 + cfs01(i)*((pvar(i) - mns(i))/sds(i)) sumv02 = sumv02 + cfs02(i)*((pvar(i) - mns(i))/sds(i)) sumv03 = sumv03 + cfs03(i)*((pvar(i) - mns(i))/sds(i)) sumv04 = sumv04 + cfs04(i)*((pvar(i) - mns(i))/sds(i)) sumv05 = sumv05 + cfs05(i)*((pvar(i) - mns(i))/sds(i)) sumv06 = sumv06 + cfs06(i)*((pvar(i) - mns(i))/sds(i)) sumv07 = sumv07 + cfs07(i)*((pvar(i) - mns(i))/sds(i)) sumv08 = sumv08 + cfs08(i)*((pvar(i) - mns(i))/sds(i)) sumv09 = sumv09 + cfs09(i)*((pvar(i) - mns(i))/sds(i)) sumv10 = sumv10 + cfs10(i)*((pvar(i) - mns(i))/sds(i)) sumv11 = sumv11 + cfs11(i)*((pvar(i) - mns(i))/sds(i)) sumv12 = sumv12 + cfs12(i)*((pvar(i) - mns(i))/sds(i)) sumv13 = sumv13 + cfs13(i)*((pvar(i) - mns(i))/sds(i)) sumv14 = sumv14 + cfs14(i)*((pvar(i) - mns(i))/sds(i)) sumv15 = sumv15 + cfs15(i)*((pvar(i) - mns(i))/sds(i)) sumv16 = sumv16 + cfs16(i)*((pvar(i) - mns(i))/sds(i)) sumv17 = sumv17 + cfs17(i)*((pvar(i) - mns(i))/sds(i)) sumv18 = sumv18 + cfs18(i)*((pvar(i) - mns(i))/sds(i)) sumv19 = sumv19 + cfs19(i)*((pvar(i) - mns(i))/sds(i)) enddo goes00(1) = mns01 + (sumv01*sds01) goes00(2) = mns02 + (sumv02*sds02) goes00(3) = mns03 + (sumv03*sds03) goes00(4) = mns04 + (sumv04*sds04) goes00(5) = mns05 + (sumv05*sds05) goes00(6) = mns06 + (sumv06*sds06) goes00(7) = mns07 + (sumv07*sds07) goes00(8) = mns08 + (sumv08*sds08) goes00(9) = mns09 + (sumv09*sds09) goes00(10) = mns10 + (sumv10*sds10) goes00(11) = mns11 + (sumv11*sds11) goes00(12) = mns12 + (sumv12*sds12) goes00(13) = mns13 + (sumv13*sds13) goes00(14) = mns14 + (sumv14*sds14) goes00(15) = mns15 + (sumv15*sds15) goes00(16) = mns16 + (sumv16*sds16) goes00(17) = mns17 + (sumv17*sds17) goes00(18) = mns18 + (sumv18*sds18) goes00(19) = mns19 + (sumv19*sds19) else c Error processing ierr = 1 do i=1,19 goes00(i) = rmiss enddo return endif c c c Apply max/min caps to parameters c do i=1,19 if (goes00(i) .lt. cmin(i)) goes00(i) = cmin(i) if (goes00(i) .gt. cmax(i)) goes00(i) = cmax(i) enddo c c Apply scale factors do i=1,19 goes00(i) = goes00(i)*sfac(i) enddo c return c end