! Station Extraction Parameters. ! !svn $Id: stations_test_head.in 804 2016-06-06 23:02:36Z arango $ !========================================================= Hernan G. Arango === ! Copyright (c) 2002-2016 The ROMS/TOMS Group ! ! Licensed under a MIT/X style license ! ! See License_ROMS.txt ! !============================================================================== ! ! ! Input parameters can be entered in ANY order, provided that the parameter ! ! KEYWORD (usually, upper case) is typed correctly followed by "=" or "==" ! ! symbols. Any comment lines are allowed and must begin with an exclamation ! ! mark (!) in column one. Comments may appear to the right of a parameter ! ! specification to improve documentation. Comments will be ignored during ! ! reading. Blank lines are also allowed and ignored. Continuation lines in ! ! a parameter specification are allowed and must be preceded by a backslash ! ! (\). In some instances, more than one value is required for a parameter. ! ! If fewer values are provided, the last value is assigned for the entire ! ! parameter array. The multiplication symbol (*), without blank spaces in ! ! between, is allowed for a parameter specification. For example, in a two ! ! grids nested application: ! ! ! ! AKT_BAK == 2*1.0d-6 2*5.0d-6 ! m2/s ! ! ! ! indicates that the first two entries of array AKT_BAK, in fortran column- ! ! major order, will have the same value of "1.0d-6" for grid 1, whereas the ! ! next two entries will have the same value of "5.0d-6" for grid 2. ! ! ! ! In multiple levels of nesting and/or multiple connected domains step-ups, ! ! "Ngrids" entries are expected for some of these parameters. In such case, ! ! the order of the entries for a parameter is extremely important. It must ! ! follow the same order (1:Ngrids) as in the state variable declaration. The ! ! USER may follow the above guidelines for specifying his/her values. These ! ! parameters are marked by "==" plural symbol after the KEYWORD. ! ! ! !============================================================================== ! ! Switch to control the writing of stations data within nested and/or multiple ! connected grids, [1:Ngrids]. Lstations == T ! Logical switches (TRUE/FALSE) to activate writing of fields in STATION ! output file, [Sout(:,ng), ng=1, Ngrids]. Sout(idUvel) == T ! u 3D U-velocity Sout(idVvel) == T ! v 3D V-velocity Sout(idu3dE) == F ! u_eastward 3D U-eastward at RHO-points Sout(idv3dN) == F ! v_northward 3D V-nortward at RHO-points Sout(idWvel) == T ! w 3D W-velocity Sout(idOvel) == F ! omega 3D omega vertical velocity Sout(idUbar) == F ! ubar 2D U-velocity Sout(idVbar) == F ! vbar 2D V-velocity Sout(idu2dE) == F ! ubar_eastward 2D U-eastward at RHO-points Sout(idv2dN) == F ! vbar_northward 2D V-northward at RHO-points Sout(idFsur) == T ! zeta free-surface Sout(idBath) == T ! bath time-dependent bathymetry Sout(idTvar) == T T ! temp, salt, ... all (NT) tracers Sout(idUsms) == F ! sustr surface U-stress Sout(idVsms) == F ! svstr surface V-stress Sout(idUbms) == F ! bustr bottom U-stress Sout(idVbms) == F ! bvstr bottom V-stress Sout(idUbrs) == F ! bustrc bottom U-current stress Sout(idVbrs) == F ! bvstrc bottom V-current stress Sout(idUbws) == F ! bustrw bottom U-wave stress Sout(idVbws) == F ! bvstrw bottom V-wave stress Sout(idUbcs) == F ! bustrcwmax bottom max wave-current U-stress Sout(idVbcs) == F ! bvstrcwmax bottom max wave-current V-stress Sout(idUbot) == F ! Ubot bed wave orbital U-velocity Sout(idVbot) == F ! Vbot bed wave orbital V-velocity Sout(idUbur) == F ! Ur bottom U-velocity above bed Sout(idVbvr) == F ! Vr bottom V-velocity above bed Sout(idW2xx) == F ! Sxx_bar 2D radiation stress, Sxx component Sout(idW2xy) == F ! Sxy_bar 2D radiation stress, Sxy component Sout(idW2yy) == F ! Syy_bar 2D radiation stress, Syy component Sout(idU2rs) == F ! Ubar_Rstress 2D radiation U-stress Sout(idV2rs) == F ! Vbar_Rstress 2D radiation V-stress Sout(idU2Sd) == F ! ubar_stokes 2D U-Stokes velocity Sout(idV2Sd) == F ! vbar_stokes 2D V-Stokes velocity Sout(idW3xx) == F ! Sxx 3D radiation stress, Sxx component Sout(idW3xy) == F ! Sxy 3D radiation stress, Sxy component Sout(idW3yy) == F ! Syy 3D radiation stress, Syy component Sout(idW3zx) == F ! Szx 3D radiation stress, Szx component Sout(idW3zy) == F ! Szy 3D radiation stress, Szy component Sout(idU3rs) == F ! u_Rstress 3D U-radiation stress Sout(idV3rs) == F ! v_Rstress 3D V-radiation stress Sout(idU3Sd) == F ! u_stokes 3D U-Stokes velocity Sout(idV3Sd) == F ! v_stokes 3D V-Stokes velocity Sout(idWamp) == F ! Hwave wave height Sout(idWlen) == F ! Lwave wave length Sout(idWdir) == F ! Dwave wave direction Sout(idWptp) == F ! Pwave_top wave surface period Sout(idWpbt) == F ! Pwave_bot wave bottom period Sout(idWorb) == F ! Ub_swan wave bottom orbital velocity Sout(idWdis) == F ! Wave_dissip wave dissipation Sout(idPair) == T ! Pair surface air pressure Sout(idUair) == T ! Uair surface U-wind component Sout(idVair) == T ! Vair surface V-wind component Sout(idTsur) == F F ! shflux, ssflux surface net heat and salt flux Sout(idLhea) == F ! latent latent heat flux Sout(idShea) == F ! sensible sensible heat flux Sout(idLrad) == F ! lwrad longwave radiation flux Sout(idSrad) == F ! swrad shortwave radiation flux Sout(idEmPf) == F ! EminusP E-P flux Sout(idevap) == F ! evaporation evaporation rate Sout(idrain) == F ! rain precipitation rate Sout(idDano) == F ! rho density anomaly Sout(idVvis) == F ! AKv vertical viscosity Sout(idTdif) == F ! AKt vertical T-diffusion Sout(idSdif) == F ! AKs vertical Salinity diffusion Sout(idHsbl) == F ! Hsbl depth of surface boundary layer Sout(idHbbl) == F ! Hbbl depth of bottom boundary layer Sout(idMtke) == F ! tke turbulent kinetic energy Sout(idMtls) == F ! gls turbulent length scale ! Logical switches (TRUE/FALSE) to activate writing of exposed sediment ! layer properties into STATIONS output file. Sout(isd50) == F ! grain_diameter mean grain diameter Sout(idens) == F ! grain_density mean grain density Sout(iwsed) == F ! settling_vel mean settling velocity Sout(itauc) == F ! erosion_stress critical erosion stress Sout(irlen) == F ! ripple_length ripple length Sout(irhgt) == F ! ripple_height ripple height Sout(ibwav) == F ! bed_wave_amp wave excursion amplitude Sout(izdef) == F ! Zo_def default bottom roughness Sout(izapp) == F ! Zo_app apparent bottom roughness Sout(izNik) == F ! Zo_Nik Nikuradse bottom roughness Sout(izbio) == F ! Zo_bio biological bottom roughness Sout(izbfm) == F ! Zo_bedform bed form bottom roughness Sout(izbld) == F ! Zo_bedload bed load bottom roughness Sout(izwbl) == F ! Zo_wbl wave bottom roughness Sout(iactv) == F ! active_layer_t... active layer thickness Sout(ishgt) == F ! saltation saltation height ! Number of stations to process in each nested grid. These values are ! essential because the station arrays are dynamically allocated using ! these values, [1:Ngrids]. NSTATION == 126 ! Station locations for all grids in any desired order. The horizontal ! location for a particular station may be specified in terms of fractional ! (I,J) grid pairs (FLAG=0) or (longitude,latitude) grid pairs (FLAG=1). ! Here, FLAG is a special switch and may be used for multiple purposes. ! The GRID column indicates nested grid number to process. This value must ! be one in non-nested applications. The COMMENT section is ignored during ! reading and may be used to help documentation. POS = GRID FLAG X-POS Y-POS COMMENT 1 0 779 742 "8410140,-66.982,44.903, (1) CO-OPS WL, station name: Eastport, Maine" 1 0 733 714 "8411060,-67.21,44.657, (2) CO-OPS WL, station name: Cutler Farris Wharf, Maine" 1 0 611 728 "8413320,-68.205,44.392, (3) CO-OPS WL, station name: Bar Harbor, Maine" 1 0 351 729 "8418150,-70.247,43.657, (4) CO-OPS WL, station name: Portland, Maine" 1 0 292 705 "8419317,-70.563,43.32, (5) CO-OPS WL, station name: Wells, Maine" 1 0 257 678 "8423898,-70.712,43.072, (6) CO-OPS WL, station name: Fort Pt., NH" 1 0 170 596 "8443970,-71.053,42.353, (7) CO-OPS WL, station name: Boston, Ma" 1 0 230 442 "8447435,-69.95,41.688, (9) CO-OPS WL, station name: Chatham, Ma" 1 0 141 461 "8447930,-70.672,41.523, (10) CO-OPS WL, station name: Woods Hole, Ma" 1 0 119 446 "8448725,-70.767,41.353, (11) CO-OPS WL, station name: Menemsha Harbor, Ma" 1 0 177 385 "8449130,-70.097,41.202, (12) CO-OPS WL, station name: Nantucket Island, Ma" 1 0 72 498 "8452660,-71.327,41.505, (13) CO-OPS WL, station name: Newport RI" 1 0 87 529 "8452944,-71.343,41.717, (14) CO-OPS WL, station name: Conimicut Light, RI" 1 0 779 742 "8410140,-66.982,44.903, (1) CO-OPS TS, station name: "(1):Eastport,Maine"" 1 0 733 714 "8411060,-67.21,44.657, (2) CO-OPS TS, station name: "(6):CutlerFarrisWharf,Maine"" 1 0 611 728 "8413320,-68.205,44.392, (3) CO-OPS TS, station name: "(9):BarHarbor,Maine"" 1 0 351 729 "8418150,-70.247,43.657, (4) CO-OPS TS, station name: "(27):Portland,Maine"" 1 0 292 705 "8419317,-70.563,43.32, (5) CO-OPS TS, station name: "(32):Wells,Maine"" 1 0 166 595 "8443970,-71.053,42.353, (6) CO-OPS TS, station name: "(52):Boston,Ma"" 1 0 142 464 "8447930,-70.672,41.523, (7) CO-OPS TS, station name: "(81):WoodsHole,Ma"" 1 0 177 385 "8449130,-70.097,41.202, (8) CO-OPS TS, station name: "(90):NantucketIsland,Ma"" 1 0 72 498 "8452660,-71.327,41.505, (9) CO-OPS TS, station name: "(93):NewportRI"" 1 0 230 595 "A01,-70.5652,42.5181, (1) NERACOOS TS, station name: A01" 1 0 295 679 "B01,-70.4277,43.1806, (2) NERACOOS TS, station name: B01" 1 0 444 694 "E01,-69.3558,43.7152, (3) NERACOOS TS, station name: E01" 1 0 506 723 "F01,-68.9982,44.055, (4) NERACOOS TS, station name: F01" 1 0 600 683 "I01,-68.1085,44.106, (5) NERACOOS TS, station name: I01" 1 0 608 573 "M01,-67.614,43.512, (6) NERACOOS TS, station name: M01" 1 0 704 315 "N01,-65.9069,42.3304, (7) NERACOOS TS, station name: N01" 1 0 53 414 "44097,-71.117,40.981, (1) NDBC met T, station name: Block Island, RI" 1 0 167 594 "bhbm3,-71.052,42.355, (2) NDBC met T, station name: Boston, Ma" 1 0 136 452 "44020,-70.672,41.443, (3) NDBC met T, station name: Nantucket Main Channel Lighted buoy" 1 0 141 461 "BZBM3-8447930,-70.672,41.523, (4) NDBC met T, station name: Woods Hole, Ma" 1 0 208 576 "44013,-70.651,42.346, (5) NDBC met T, station name: Boston, Ma" 1 0 229 595 "44029,-70.57,42.52, (6) NDBC met T, station name: Bouy A01, Ma" 1 0 296 679 "44030,-70.418,43.183, (7) NDBC met T, station name: Bouy A01, Ma-Western Maine Shelf" 1 0 292 612 "44098,-70.169,42.801, (8) NDBC met T, station name: Jeffrey's Ledge, NH" 1 0 350 712 "44007,-70.144,43.531, (9) NDBC met T, station name: 12 NM SE of Portland, Maine" 1 0 350 712 "WELM1-8419317,-70.144,43.531, (10) NDBC met T, station name: Wells Maine" 1 0 185 398 "NTKM3,-70.097,41.285, (11) NDBC met T, station name: Nantucket Island, Ma." 1 0 297 487 "44018,-69.63,42.126, (12) NDBC met T, station name: Cape Cod, Ma" 1 0 444 695 "44032,-69.358,43.715, (13) NDBC met T, station name: Buoy E01 Central Maine Shelf" 1 0 214 238 "44008,-69.247,40.502, (14) NDBC met T, station name: Nantucket" 1 0 429 610 "44005,-69.128,43.204, (15) NDBC met T, station name: Gulf of Maine" 1 0 506 724 "44033,-69,44.06, (16) NDBC met T, station name: Buoy F01-West Penobscot Bay" 1 0 611 728 "ATGM1-8413320,-68.205,44.392, (17) NDBC met T, station name: Bar Harbor, Maine" 1 0 601 684 "44034,-68.11,44.11, (18) NDBC met T, station name: Buoy Io1 Eastern Maine Shelf" 1 0 578 583 "44037,-67.883,43.484, (19) NDBC met T, station name: Buoy M01 Jordan Basin" 1 0 695 667 "44027,-67.307,44.287, (20) NDBC met T, station name: Jonesport Maine" 1 0 733 714 "CFWM1-8411060,-67.21,44.657, (21) NDBC met T, station name: Cutler Farris Wharf, Maine" 1 0 779 742 "PSBM1-8410140,-66.985,44.903, (22) NDBC met T, station name: Eastport, Maine" 1 0 543 176 "44011,-66.6,41.105, (23) NDBC met T, station name: Georges Bank" 1 0 700 313 "44024,-65.927,42.312, (24) NDBC met T, station name: Northeast Channel" 1 0 774 740 "8410140,-66.982,44.903, WL: 2012 - Jul 2013" 1 0 734 712 "8411060,-67.21,44.657, WL: 2012 - Jul 2013" 1 0 720 717 "8411250,-67.297,44.642, none" 1 0 657 728 "8412581,-67.875,44.54, none" 1 0 613 728 "8413320,-68.205,44.392, WL: 2012 - Jul 2013" 1 0 572 711 "8413825,-68.435,44.17, none" 1 0 555 719 "8414249,-68.62,44.192, none" 1 0 556 771 "8414692,-68.805,44.505, none X68" 1 0 556 771 "8414721,-68.813,44.472, none" 1 0 525 724 "8414856,-68.873,44.127, none" 1 0 525 733 "8414888,-68.887,44.157, none" 1 0 404 721 "8417177,-69.785,43.755, none" 1 0 366 740 "8417874,-70.14,43.752, none" 1 0 352 729 "8418150,-70.247,43.657, WL: 2012 - Jul 2013" 1 0 331 722 "8418445,-70.333,43.545, none" 1 0 322 715 "8418606,-70.382,43.462, none" 1 0 257 676 "8419688,-70.697,43.067, none" 1 0 256 677 "8423898,-70.712,43.072, WL: 2012 - Jul 2013" 1 0 228 650 "8440452,-70.82,42.817, none" 1 0 223 632 "8441241,-70.788,42.71, none" 1 0 236 617 "8441551,-70.615,42.658, none" 1 0 225 611 "8441841,-70.66,42.61, none" 1 0 201 613 "8442645,-70.877,42.523, none" 1 0 184 606 "8443187,-70.943,42.458, none" 1 0 182 587 "8444162,-70.892,42.328, none" 1 0 171 585 "8444525,-70.953,42.28, none" 1 0 190 558 "8445138,-70.727,42.202, none" 1 0 187 536 "8446009,-70.647,42.083, none" 1 0 234 505 "8446121,-70.182,42.048, none" 1 0 181 531 "8446166,-70.67,42.038, none" 1 0 177 524 "8446493,-70.662,41.96, none" 1 0 180 487 "8447180,-70.507,41.772, none" 1 0 214 465 "8447241,-70.155,41.752, none" 1 0 162 487 "8447295,-70.623,41.735, none" 1 0 153 491 "8447368,-70.715,41.712, none" 1 0 151 489 "8447416,-70.72,41.695, none" 1 0 229 441 "8447435,-69.95,41.688, WL: 2012 - Jul 2013" 1 0 217 446 "8447495,-70.057,41.668, none" 1 0 138 487 "8447531,-70.813,41.657, none" 1 0 189 455 "8447605,-70.3,41.632, none" 1 0 150 473 "8447685,-70.652,41.605, none" 1 0 124 485 "8447712,-70.9,41.593, none" 1 0 214 426 "8447819,-69.99,41.56, none" 1 0 116 479 "8447842,-70.928,41.538, none" 1 0 141 461 "8447930,-70.672,41.523, WL: 2012 - Jul 2013" 1 0 127 458 "8448151,-70.76,41.47, none" 1 0 147 451 "8448157,-70.6,41.458, none" 1 0 148 444 "8448208,-70.555,41.458, none" 1 0 108 465 "8448428,-70.922,41.45, none" 1 0 116 460 "8448251,-70.857,41.448, none" 1 0 127 449 "8448533,-70.712,41.41, none" 1 0 150 435 "8448558,-70.512,41.388, none" 1 0 119 445 "8448725,-70.767,41.353, none" 1 0 82 482 "8450768,-71.193,41.465, none" 1 0 89 513 "8451552,-71.255,41.637, none" 1 0 71 496 "8452660,-71.327,41.505, WL: 2012 - Jul 2013" 1 0 87 528 "8452944,-71.343,41.717, WL: 2012 - Jul 2013" 1 0 88 531 "8453033,-71.352,41.752, none" 1 0 86 535 "8453767,-71.388,41.762, none" 1 0 69 512 "8454099,-71.41,41.587, WL: 2012 - Jul 2013" 1 0 66 514 "8454538,-71.445,41.572, none" 1 0 51 492 "8454658,-71.455,41.422, none" 1 0 44 488 "8455083,-71.49,41.363, none" 1 0 13 499 "8458022,-71.762,41.328, none" 1 0 23 466 "8459338,-71.557,41.173, none" 1 0 17 461 "8459449,-71.575,41.148, none" 1 0 23 475 "8459479,-71.58,41.228, none" 1 0 15 468 "8459681,-71.61,41.163, none" 1 0 1 250 "IJ, 1 250" 1 0 1 1 "IJ, 1 1" 1 0 500 1 "IJ, 500 1" 1 0 1171 1 "IJ, 1171 1" 1 0 1171 250 "IJ, 1171 250" ! ! GLOSSARY: ! ========= ! ! This input file provides the positions of stations data to be extracted ! and written into STATIONS output file. This file is use to monitor model ! fields at frequent time intervals. It is also used to sample the model at ! observation locations in order to build the penalty function used in ! variational data assimilation. ! !------------------------------------------------------------------------------ ! Control Switch. !------------------------------------------------------------------------------ ! ! Lstations Switch to control the writing of station data within nested ! and/or multiple connected grids. Ngrid values are expected. ! By default this switch is set to TRUE in "mod_scalars" for ! all grids when the CPP option STATIONS is activated. The ! USER can control which grid to process by turning on/off ! this switch. ! !------------------------------------------------------------------------------ ! Logical switches (TRUE/FALSE) to activate writing of fields in STATIONS ! output file. Except for tracer type variables, Ngrids values are expected ! for all switches. !------------------------------------------------------------------------------ ! ! Sout(idUvel) Write out 3D U-velocity component. ! Sout(idVvel) Write out 3D V-velocity component. ! Sout(idu3dE) Write out 3D Eastward velocity component at RHO-points. ! Sout(idv3dN) Write out 3D Northward velocity component at RHO-points. ! Sout(idWvel) Write out 3D W-velocity component. ! Sout(idOvel) Write out 3D omega vertical velocity. ! Sout(idUbar) Write out 2D U-velocity component. ! Sout(idVbar) Write out 2D V-velocity component. ! Sout(idu2dE) Write out 2D Eastward velocity component at RHO-points. ! Sout(idv2dN) Write out 2D Northward velocity component at RHO-points. ! Sout(idFsur) Write out free-surface. ! Sout(idBath) Write out time-dependent bathymetry. ! ! Sout(idTvar) Write out tracer type variables, (NT,Ngrids) values are ! Expected. ! ! Sout(idUsms) Write out surface U-momentum stress. ! Sout(idVsms) Write out surface V-momentum stress. ! Sout(idUbms) Write out bottom U-momentum stress. ! Sout(idVbms) Write out bottom V-momentum stress. ! ! Sout(idUbrs) Write out current-induced, U-momentum stress. ! Sout(idVbrs) Write out current-induced, V-momentum stress. ! Sout(idUbws) Write out wind-induced, bottom U-wave stress. ! Sout(idVbws) Write out wind-induced, bottom V-wave stress. ! Sout(idUbcs) Write out bottom maximum wave and current U-stress. ! Sout(idVbcs) Write out bottom maximum wave and current V-stress. ! ! Sout(idUbot) Write out wind-induced, bed wave orbital U-velocity. ! Sout(idVbot) Write out wind-induced, bed wave orbital V-velocity. ! Sout(idUbur) Write out bottom U-velocity above bed. ! Sout(idVbvr) Write out bottom V-velocity above bed. ! ! Sout(idW2xx) Write out 2D radiation stress, Sxx component. ! Sout(idW2xy) Write out 2D radiation stress, Sxy component. ! Sout(idW2yy) Write out 2D radiation stress, Syy component. ! Sout(idU2rs) Write out 2D U-radiation stress. ! Sout(idV2rs) Write out 2D V-radiation stress. ! Sout(idU2Sd) Write out 2D U-Stokes velocity. ! Sout(idV2Sd) Write out 2D V-Stokes velocity. ! ! Sout(idW3xx) Write out 3D radiation stress, Sxx component. ! Sout(idW3xy) Write out 3D radiation stress, Sxy component. ! Sout(idW3yy) Write out 3D radiation stress, Syy component. ! Sout(idW3zx) Write out 3D radiation stress, Szx component. ! Sout(idW3zy) Write out 3D radiation stress, Szy component. ! Sout(idU3rs) Write out 3D U-radiation stress. ! Sout(idV3rs) Write out 3D V-radiation stress. ! Sout(idU3Sd) Write out 3D U-Stokes velocity. ! Sout(idV3Sd) Write out 3D V-Stokes velocity. ! ! Sout(idWamp) Write out wave height. ! Sout(idWlen) Write out wave length. ! Sout(idWdir) Write out wave direction. ! Sout(idWptp) Write out wave surface period. ! Sout(idWpbt) Write out wave bottom period. ! Sout(idWorb) Write out wave bottom orbital velocity. ! Sout(idWdis) Write out wave dissipation. ! ! Sout(idPair) Write out surface air pressure. ! Sout(idUair) Write out surface U-wind component. ! Sout(idVair) Write out surface V-wind component. ! ! Sout(idTsur) Write out surface net heat and salt flux. ! Sout(idLhea) Write out latent heat flux. ! Sout(idShea) Write out sensible heat flux. ! Sout(idLrad) Write out long-wave radiation flux. ! Sout(idSrad) Write out short-wave radiation flux. ! Sout(idEmPf) Write out E-P flux. ! Sout(idevap) Write out evaporation rate. ! Sout(idrain) Write out precipitation rate. ! ! Sout(idDano) Write out density anomaly. ! Sout(idVvis) Write out vertical viscosity coefficient. ! Sout(idTdif) Write out vertical diffusion coefficient of temperature. ! Sout(idSdif) Write out vertical diffusion coefficient of salinity. ! Sout(idHsbl) Write out depth of oceanic surface boundary layer. ! Sout(idHbbl) Write out depth of oceanic bottom boundary layer. ! Sout(idMtke) Write out turbulent kinetic energy. ! Sout(idMtls) Write out turbulent kinetic energy times length scale. ! ! Sout(isd50) Write out mean grain diameter. ! Sout(idens) Write out mean grain density. ! Sout(iwsed) Write out mean settling velocity. ! Sout(itauc) Write out critical erosion stress. ! Sout(irlen) Write out ripple length. ! Sout(irhgt) Write out ripple height. ! Sout(ibwav) Write out wave excursion amplitude. ! Sout(izdef) Write out default bottom roughness. ! Sout(izapp) Write out apparent bottom roughness. ! Sout(izNik) Write out Nikuradse bottom roughness. ! Sout(izbio) Write out biological bottom roughness. ! Sout(izbfm) Write out bed form bottom roughness. ! Sout(izbld) Write out bed load bottom roughness. ! Sout(izwbl) Write out wave bottom roughness. ! Sout(iactv) Write out active layer thickness. ! Sout(ishgt) Write out saltation height. ! !------------------------------------------------------------------------------ ! Station positions KEYWORDS. !------------------------------------------------------------------------------ ! ! NSTATION Number of stations to process in each nested grid. Ngrids ! values are expected. These values are essential because ! the station arrays are dynamically allocated using these ! values. ! ! POS Station locations for all grids in any desired order. The ! horizontal location for a particular station can be ! specified in terms of fractional (I,J) grid pairs (FLAG=0) ! or (longitude,latitude) grid pairs (FLAG=1). Here, FLAG is ! a special switch and may be used for multiple purposes. ! The GRID column indicates nested grid number to process. ! This value must be one in non-nested applications. The ! COMMENT section is ignored during reading and may be used ! to help documentation. This locations are used to ! interpolate model data before writting to output STATIONS ! file. The stations locations are read using the following ! statement: ! ! is(1:Ngrids)=0 ! DO WHILE (.true.) ! read (unit,*,ERR=10) igrid, flag, Xpos, Ypos ! ng=MAX(1,ABS(ng)) ! is(ng)=ic(ng)+1 ! SCALARS(ng)%Sflag(is(ng))=flag ! SCALARS(ng)%SposX(is(ng))=Xpos ! SCALARS(ng)%SposY(is(ng))=Ypos ! END DO ! ! Since there is no decoding during the reading of the stations ! locations, the POS keyword must be the last one to process. ! That is, it should be at the end of the input specification. ! !