anchor.diff:           IF(LMFDUDV) THEN
               DO 224 JL=1,KLON
               IF(PMFD(JL,JK).LT.0.) THEN
  -               PUD(JL,JK)=0.5*(PUU(JL,JK)+PUEN(JL,JK-1))
  -               PVD(JL,JK)=0.5*(PVU(JL,JK)+PVEN(JL,JK-1))
  +				if (jk.eq.1) then !edg
  +	                PUD(JL,JK)=0.5*(PUU(JL,JK)+PUEN(JL,JK))
  +	                PVD(JL,JK)=0.5*(PVU(JL,JK)+PVEN(JL,JK))
  +				else
  +	               PUD(JL,JK)=0.5*(PUU(JL,JK)+PUEN(JL,JK-1))
  +	               PVD(JL,JK)=0.5*(PVU(JL,JK)+PVEN(JL,JK-1))
  +			   end if
               END IF
     224       CONTINUE
            END IF
  
positive.code:          PMFUQ(JL,JK)=0.
           PMFDQ(JL,JK)=0.
           PMFUL(JL,JK)=0.
           PDMFUP(JL,JK-1)=0.
           PDMFDP(JL,JK-1)=0.
           PLUDE(JL,JK-1)=0.
        END IF
    115 CONTINUE
    120 CONTINUE
  
negative.code:             PMFD(JL,JK)=0.
              PMFDS(JL,JK)=0.
              PMFDQ(JL,JK)=0.
  			if (jk.gt.1) then !edg with low model level top, things can get screwy and jk can endup = 1
  	            PDMFDP(JL,JK-1)=0.
  			end if
           END IF
        ELSE
           PMFU(JL,JK)=0.
  