Added new output variables to ocean.in and stations.in

[arango]

I updated several things:

• Added new input fields to output history NetCDF file using the Hout logical switch in ocean.in. Many thanks to Mathieu Dutour for bringing this to my attention:

  Hout(idWptp) == F                          ! wave surface period
  Hout(idWpbt) == F                          ! wave bottom period
  Hout(idWorb) == F                          ! wave bottom orbital velocity
  Hout(idWdis) == F                          ! wave dissipation
  

• Added new input fields to output stations NetCDF file using the Hout logical switch in stations.in. Many thanks to Mathieu Dutour for bringing this to my attention:

  Sout(idUbrs) == F                          ! bottom U-current stress
  Sout(idVbrs) == F                          ! bottom V-current stress
  Sout(idUbws) == F                          ! bottom U-wave stress
  Sout(idVbws) == F                          ! bottom V-wave stress
  Sout(idUbcs) == F                          ! bottom max wave-current U-stress
  Sout(idVbcs) == F                          ! bottom max wave-current V-stress
  
  Sout(idUbot) == F                          ! bed wave orbital U-velocity
  Sout(idVbot) == F                          ! bed wave orbital V-velocity
  Sout(idUbur) == F                          ! bottom U-velocity above bed
  Sout(idVbvr) == F                          ! bottom V-velocity above bed
  
  Sout(idW2xx) == F                          ! 2D radiation stress, Sxx component
  Sout(idW2xy) == F                          ! 2D radiation stress, Sxy component
  Sout(idW2yy) == F                          ! 2D radiation stress, Syy component
  Sout(idU2rs) == F                          ! 2D radiation U-stress
  Sout(idV2rs) == F                          ! 2D radiation V-stress
  Sout(idU2Sd) == F                          ! 2D U-Stokes velocity
  Sout(idV2Sd) == F                          ! 2D V-Stokes velocity
  
  Sout(idW3xx) == F                          ! 3D radiation stress, Sxx component
  Sout(idW3xy) == F                          ! 3D radiation stress, Sxy component
  Sout(idW3yy) == F                          ! 3D radiation stress, Syy component
  Sout(idW3zx) == F                          ! 3D radiation stress, Szx component
  Sout(idW3zy) == F                          ! 3D radiation stress, Szy component
  Sout(idU3rs) == F                          ! 3D U-radiation stress
  Sout(idV3rs) == F                          ! 3D V-radiation stress
  Sout(idU3Sd) == F                          ! 3D U-Stokes velocity
  Sout(idV3Sd) == F                          ! 3D V-Stokes velocity
  
  Sout(idWamp) == F                          ! wave height
  Sout(idWlen) == F                          ! wave length
  Sout(idWdir) == F                          ! wave direction
  Sout(idWptp) == F                          ! wave surface period
  Sout(idWpbt) == F                          ! wave bottom period
  Sout(idWorb) == F                          ! wave bottom orbital velocity
  Sout(idWdis) == F                          ! wave dissipation
  

• Update ROMS/Programs/types.F to include logic for REAL*16 type variables. Notice that there are several compilers that do not support REAL*16 data types. For example PGF90, G95, and GFORTRAN. However, other compilers like IFORT do have REAL*16 support.

• Corrected a compiler bug in the biology and sediment models include files *_mod.h. We need to use allocatable instead of pointer attributes during declaration. Many thanks to John Wilkin for reporting this problem. The PGI compiler didn't complained about. I am surprised how the PGI compiler relax the standard F90/F95 roules.