calmet.inp

Changes made to calmet.inp

• GEODAT – Set to data/geo_masaya.dat.
• LCFILES – Set to T
• NUSTA – Set to 0 (we don’t have any upper-air stations)
• NOWSTA – Set to 0 (we don’t have any overwater stations)
• UPDAT – Comment out all three lines (we don’t have upper air data)
• SEADAT – Comment out (we don’t have overwater stations data)
• M3DDAT – Set to data/out_caleta.m3d
• IBYR, IBMO, IBDY, IBHR – Set to required simulation start date/- time (e.g. 2017, 12, 04, 0)
• IEYR, IEMO, IEDY, IEHR – Set to required simulation end date/- time (e.g. 2017, 12, 06, 0)
• ABTZ – Set to UTC+0000 (we want times in UTC)
• LCALGRD – Set to F (don’t compute special data fields)
• IUTMZN – Set to 16
• DATUM – Set to WGS-84
• NX – Set to 91 (1km res) or 901 (100m res)
• NY – Set to 55 (1km res) or 541 (100m res)
• DGRIDKM – Set to 1.0 (1km res) or 0.1 (100m res)
• XORIGKM – Set to 525.09417
• YORIGKM – Set to 1294.23926
• NZ – Set to 24
• ZFACE – Set to "0, 20, 60, 140, 260, 540, 1000, 2000, 3000, 4000, 5000, 6000, 7000, 8000, 9000, 10000, 11000, 12000, 13000, 14000, 15000, 16000, 17000, 18000, 19000". These are the heights of the output grid
• NOOBS – Set to 2 (use 3D data for surface, overwater and upper air data)
• NSSTA – Set to 0 (no surface stations)
• NPSTA – Set to -1 (flag to use 3D precip data)
• MCLOUD – Set to 4 (use prognostic 3D RH fields to get gridded cloud cover at all levels)
• BIAS – Set to 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 (i.e. zero bias at each Z level)
• IPROG – Set to 14 (i.e. use 3D winds as initial guess field)
• ISTEPPGS – Set to 10800 (timestep of the prognostic met model) This was set to 3 in Mark’s input file, but this caused an error for me, plus 10800 seems correct.
• RMAX1-RMAX3 – Set all to 5.0 (radii of influence in km for stations)
• RPROG – Set to 1.0 (weighting parameter of the 3D data).
• NSMTH – Set to 2, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4, 4 ,4 (number of passes in the smoothing procedure for each Z level)
• NINTR2 – (not needed?) Set to 5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5,5 (max stations used per Z level for interpolation)
• FEXTR2 – (not needed?) Set to 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0 (not used but need NZ values)
• ISURFT – Set to -2 (use domain-average prognostic lapse rate)
• IUPT – Set to -2 (use domain-average prognostic lapse rate)
• ZUPT – Set to 100. (depth through which the domain-scale lapse rate is computed)
• IMIXH – Set to 2 (Batchvarova and Gryning mixing height method used over all surfaces)
• DZZI – Set to 300. (depth, m, of layer above current convective mixing height through which lapse rate is computed)
• ZIMIN – Set to 250. (minimum overland mixing height)
• ITPROG – Set to 1 (use prognostic RH)
• ITPROG – Set to 2 (no surface or upper-air temperature observations; use 3D data).