*IDENT SETFILwmc
*DECLARE SETFIL1A
*D GPB3F403.115,GPB3F403.121
C WMC. Include prints of area filtered: very useful info!
C      If all is well, results should be 25-33% ish.
      WRITE(6,*) 'SETFIL : Updating filtering area to ',                  GPB3F403.115
     &           '2 - > ',NORTHERN_FILTERED_P_ROW,' and ',                GPB3F403.116
     &           SOUTHERN_FILTERED_P_ROW,' -> ',glsize(2)-1               GPB3F403.117
      WRITE(6,*) 'SETFIL: Filtering ',                                    GPB3F403.118
     &  REAL((glsize(2)-2)-                                               GPB3F403.119
     &  (SOUTHERN_FILTERED_P_ROW-NORTHERN_FILTERED_P_ROW-1))/             GPB3F403.120
     &  REAL(glsize(2)-2),' frac of grid (not counting polar rows)'       GPB3F403.121