anchor.diff:         ALLOCATE (Bp      (CmEp1),                 STAT=estat ) ; errstat = MAX ( errstat, estat )
          ALLOCATE (Temperature (Cmlon, Cmlat, CmETA), STAT=estat ) ; errstat = MAX ( errstat, estat )
          ALLOCATE (Gas_profiles(Cmlon, Cmlat, CmETA), STAT=estat ) ; errstat = MAX ( errstat, estat )
  +       ALLOCATE (H2O_profiles(Cmlon, Cmlat, CmETA), STAT=estat ) ; errstat = MAX ( errstat, estat )
          ALLOCATE (Psurface(Cmlon,Cmlat),            STAT=estat ) ; errstat = MAX ( errstat, estat )
       CASE ('d')
          IF ( ALLOCATED ( latvals      ) )  DEALLOCATE ( latvals      )
  
positive.code:     ! Local variables
      ! ---------------
      INTEGER (KIND=i4) :: itimes, ixtrack, spix, epix, idx_lat, idx_lon, ilevel, n, n1
      REAL    (KIND=r8)                      :: rho1, rho2, grad, aircolumn
      REAL    (KIND=r8), DIMENSION (0:CmETA) :: lhgt, lpre, ltmp
  
      ! -----------------------
      ! Some physical constants
      ! -----------------------
      REAL (KIND=r8), PARAMETER ::                &
           rho_stand = 2.6867773e+19_r8,          & ! Loschmidt, in [cm^-3]
           pzero     = 1013.25_r8,                & ! P0 in [mb]
           tzero     = 273.15_r8,                 &
           rho_zero  = rho_stand * tzero / pzero, &
           du_to_cm2 = 2.6867773e+16_r8
  
      ! -------------------------------
      ! Air density conversion constant
      ! -------------------------------
      REAL (KIND=r8), PARAMETER :: km2cm  = 1.0E+05_r8
      
      ! ------------------------------
      ! Name of this module/subroutine
      ! ------------------------------
  
negative.code:     ! Local variables
      ! ---------------
      INTEGER   (KIND=i4)         :: nswath, locerrstat, swath_id, swath_file_id, swlen, he5stat, &
                                     ismonth, ndatafields, h2o_cli_idx
      INTEGER   (KIND=i4), DIMENSION(10) :: datafield_rank, datafield_type
      INTEGER   (KIND=C_LONG)     :: nswathcl, Cmlatcl, Cmloncl, CmETAcl, CmEp1cl
      CHARACTER (LEN=   maxchlen) :: swath_file, locswathname, gasdatafieldname, datafield_name
      CHARACTER (LEN=10*maxchlen) :: swath_name
  
      CHARACTER (LEN= 9), PARAMETER :: cli_lat_field         = 'Latitudes'
  