pTitle = 'Pressure Pot.(p/rho)'

     I                 //pTitle//' Anomaly', myThid )

     I                 //pTitle//' Anomaly', myThid )

     I                           //pTitle//' Anomaly', myThid )

     I                           //pTitle//' Anomaly', myThid )

     I                           //pTitle//' Anomaly', myThid )

     I                           //pTitle//' Anomaly', myThid )

      diagTitle = DIAGS_MK_TITLE( 'Non-Hydrostatic '//pTitle, myThid )

      CHARACTER*(20) pTitle, sTitle

        pTitle = ' Geopotential       '