Update documentation for MSFPS.f, NGMFLD.f, OTLFT.f, OTLIFT.f, and ZE…

…NSUN.f (#681)

* generate html docs

* fix warnings on GFSPOSTSIG.F

* fix formatting

* add data types for 3 vars

* update ZENSUN.f

* fix file warnings for OTLIFT.f, OTLFT.f, NGMFLD.f, and MSFPS.f

* revert accidental changes to compile_upp.sh

* edit typos

---------

Co-authored-by: gspetro <gillian.s.petro@gmail.com>

sorc/ncep_post.fd/MSFPS.f

@@ -13,6 +13,16 @@
 !> 2006-11-01 | Rozumalski | Swiped from WRF si package
 !>
 !> @author Rozumalski @date 2006-11-01
+
+! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+!> msfps() computes the map scale factor for a polar stereographic grid at a give latitude.
+!>
+!> This subroutine computes the map scale factor for a polar stereographic grid at a give latitude.
+!>
+!> @param[in] LAT Latitude at which map factor is valid.
+!> @param[in] TRUELAT1 TRUELAT 1.
+!> @param[out] MSF Map scale factor.
+
       SUBROUTINE MSFPS(LAT,TRUELAT1,MSF)
 
 

sorc/ncep_post.fd/NGMFLD.f

@@ -4,7 +4,7 @@
 !> This routine computes a handful of NGM layer mean 
 !> fields.  This is done to provide a fully complete 
 !> ETA NGM look-alike output file.
-!> ### The sigma (layer) fields computed bu this routine are tabulated below.
+!> ### The sigma (layer) fields computed by this routine are tabulated below.
 !> Sigma (layer) | Field(s) | 
 !> --------------|----------|
 !> 0.47191 - 1.00000 |  RH    |
@@ -43,6 +43,19 @@
 !> 2021-09-30 | JESSE MENG   | 2D DECOMPOSITION
 !>
 !> @author Russ Treadon W/NP2 @date 1992-12-22
+
+!- - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+!> ngmfld() computes layer mean NGM fields
+!>
+!> This routine computes a handful of NGM layer mean 
+!> fields. 
+!>
+!> @param[out] RH4710 Sigma layer 0.47-1.00 mean relative humidity.
+!> @param[out] RH4796 Sigma layer 0.47-0.96 mean relative humidity.
+!> @param[out] RH1847 Sigma layer 0.18-0.47 mean relative humidity.
+!> @param[out] RH8498 Sigma layer 0.84-0.98 mean relative humidity.
+!> @param[out] QM8510 Sigma layer 0.85-1.00 mean moisture convergence.
+
       SUBROUTINE NGMFLD(RH4710,RH4796,RH1847,RH8498,QM8510)
 
 !     

sorc/ncep_post.fd/OTLFT.f

@@ -1,9 +1,9 @@
 !> @file
 !> @brief otlft() computes lifted index.
 !>
-!> This routine computes lifts a parcel specified by the
+!> This routine lifts a parcel specified by the
 !> passed pressure, temperature, and specific humidity to
-!> 500mb and then computes a lifted index.  This lifted
+!> 500mb and then computes a lifted index. This
 !> lifted index is the difference between the lifted
 !> parcel's temperature at 500mb and the ambient 500mb
 !> temperature.
@@ -25,6 +25,22 @@
 !> 2021-09-30 | JESSE MENG    | 2D DECOMPOSITION
 !>
 !> @author Russ Treadon W/NP2 @date 1993-03-10
+
+! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+!> otlft() computes lifted index.
+!>
+!> This routine lifts a parcel specified by the
+!> passed pressure, temperature, and specific humidity to
+!> 500mb and then computes a lifted index. This 
+!> lifted index is the difference between the lifted
+!> parcel's temperature at 500mb and the ambient 500mb
+!> temperature.
+!>
+!> @param[in] PBND Parcel pressure.
+!> @param[in] TBND Parcel temperature.
+!> @param[in] QBND Parcel specific humidity.
+!> @param[out] SLINDX Lifted index.
+
       SUBROUTINE OTLFT(PBND,TBND,QBND,SLINDX)
 
 !     

sorc/ncep_post.fd/OTLIFT.f

@@ -2,7 +2,7 @@
 !> @brief otlift() computes SFC to 500mb lifted index.
 !>
 !> This routine computes a surface to 500mb lifted index.
-!> The lifted parcel is from the first atmpspheric ETA
+!> The lifted parcel is from the first atmospheric ETA
 !> layer (ie, the ETA layer closest to the model ground).
 !> The lifted index is the difference between this parcel's
 !> temperature at 500mb and the ambient 500mb temperature.
@@ -23,6 +23,18 @@
 !> 2021-09-30 | JESSE MENG    | 2D DECOMPOSITION
 !>
 !> @author Russ Treadon W/NP2 @date 1993-03-10
+
+! - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - -
+!> otlift() computes SFC to 500mb lifted index.
+!>
+!> This routine computes a surface to 500mb lifted index.
+!> The lifted parcel is from the first atmospheric ETA
+!> layer (ie, the ETA layer closest to the model ground).
+!> The lifted index is the difference between this parcel's
+!> temperature at 500mb and the ambient 500mb temperature.
+!>
+!> @param[out] SLINDX lifted index.
+
       SUBROUTINE OTLIFT(SLINDX)
 
 !     

sorc/ncep_post.fd/ZENSUN.f

@@ -45,12 +45,13 @@
 !>  This data is characterized by 74 points.
 !>
 !>
-!> @param[in] day Julian day (positive scalar or vector), (spring equinox =  80), (summer solstice= 171), (fall equinox   = 266), (winter solstice= 356).
-!> @param[in] time Universal Time in hours (scalar or vector).
-!> @param[in] lat Geographic latitude of point on earth's surface (degrees).
-!> @param[in] lon Geographic longitude of point on earth's surface (degrees).
-!> @param[out] sun_zenith  - solar zenith angle.
-!> @param[out] sun_azimuth - solar azimuth angle.
+!> @param[in] day integer Julian day (positive scalar or vector), (spring equinox =  80), (summer solstice= 171), (fall equinox   = 266), (winter solstice= 356).
+!> @param[in] time real Universal Time in hours (scalar or vector).
+!> @param[in] lat real Geographic latitude of point on earth's surface (degrees).
+!> @param[in] lon real Geographic longitude of point on earth's surface (degrees).
+!> @param[in] pi real The mathematical constant pi.
+!> @param[out] sun_zenith real Solar zenith angle.
+!> @param[out] sun_azimuth real Solar azimuth angle.
 !>
 !> ### Program history log:
 !> Date | Programmer | Comments