anchor.diff:    CALL roundoff_2darr_r8 ( n_roff_dig, nx, nt, colloc(1:nx,0:nt-1) )
     locerrstat = HE5_SWWRFLD ( pge_swath_id, TRIM(ADJUSTL(dcol_field)), &
          he5_start_2d, he5_stride_2d, he5_edge_2d, colloc(1:nx,0:nt-1) )
  +  errstat = MAX ( errstat, locerrstat )
  +
  +  ! ---------------------------
  +  ! (6) AMF profile shape error
  +  ! ---------------------------
  +  amfloc = REAL ( amf_pro_err, KIND=r4 )
  +  CALL roundoff_2darr_r4 ( n_roff_dig, nx, nt, amfloc(1:nx,0:nt-1) )
  +  locerrstat = HE5_SWWRFLD ( pge_swath_id, TRIM(ADJUSTL(amferrpr_field)), &
  +       he5_start_2d, he5_stride_2d, he5_edge_2d, amfloc(1:nx,0:nt-1) )
     errstat = MAX ( errstat, locerrstat )
   
     RETURN
  
positive.code:   ! ---------------
    INTEGER (KIND=i4),                          INTENT (IN) :: pge_idx, nx, nt
    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1), INTENT (IN) :: saocol, saodco
    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1), INTENT (IN) :: amfmol, amfgeo
    REAL    (KIND=r8), DIMENSION (1:nx,0:nt-1), INTENT (IN) :: amfcfr, amfctp
    INTEGER (KIND=i2), DIMENSION (1:nx,0:nt-1), INTENT (IN) :: amfdiag
  
  
negative.code: 
  SUBROUTINE he5_amf_write ( &
           pge_idx, nx, nt, saocol, saodco, amfmol, amfgeo, amfdiag, &
           amfcfr, amfctp, amf_pro_err, errstat )
  
    USE OMSAO_precision_module, ONLY: i2, i4, r8
    USE OMSAO_he5_module
  