diff --git a/diag/pbldia_new_sfx.f90 b/diag/pbldia_new_sfx.f90
index 308fd8fcfe412a310596450e52158754a387b175..3773d0674761774a1e662d0c95ba75b66c0f6ca5 100755
--- a/diag/pbldia_new_sfx.f90
+++ b/diag/pbldia_new_sfx.f90
@@ -38,8 +38,8 @@ module pbldia_new_sfx
       !C*    UFTEMP     - UNIVERSAL FUNCTION FOR SCALARS
       
       
-      subroutine pbldia_new_log(AR2,ARDIN,ARDOUT,lat)
-            real,intent(in):: AR2(11),ARDIN(6),lat
+      subroutine pbldia_new_log(AR2,ARDIN,ARDOUT)
+            real,intent(in):: AR2(11),ARDIN(6)
             real,intent(out)::ARDOUT(3)
             real hwind,htemp,ustar,tstar,qstar
             real dteta,dq,wind,HIN,zeta
@@ -66,10 +66,10 @@ module pbldia_new_sfx
             return
         end subroutine pbldia_new_log
       
-      subroutine pbldia_new_sheba(AR2,ARDIN,ARDOUT,lat)
+      subroutine pbldia_new_sheba(AR2,ARDIN,ARDOUT)
       use sfx_sheba, only: get_psi_sheba => get_psi
       
-            real,intent(in):: AR2(11),ARDIN(6),lat
+            real,intent(in):: AR2(11),ARDIN(6)
             real,intent(out)::ARDOUT(3)
             real,parameter::zetalim = 1. !maximum value of z/L for 
       !      stable SL for wind
@@ -122,12 +122,12 @@ module pbldia_new_sfx
             return
         end subroutine pbldia_new_sheba
       
-      subroutine pbldia_new_sheba_coare(AR2,ARDIN,ARDOUT,lat)
+      subroutine pbldia_new_sheba_coare(AR2,ARDIN,ARDOUT)
       use sfx_sheba_coare, only: &
                       get_psi_coare => get_psi_a, &
                       get_psi_stable_sheba => get_psi_stable
       
-            real,intent(in):: AR2(11),ARDIN(6),lat
+            real,intent(in):: AR2(11),ARDIN(6)
             real,intent(out)::ARDOUT(3)
             real,parameter::zetalim = 1. !maximum value of z/L for stable SL
             real psi_m,psi_h,psi_m_hs,psi_h_hs
@@ -185,11 +185,11 @@ module pbldia_new_sfx
         end subroutine pbldia_new_sheba_coare
       
       
-      subroutine pbldia_new_sheba_noit(AR2,ARDIN,ARDOUT,lat)
+      subroutine pbldia_new_sheba_noit(AR2,ARDIN,ARDOUT)
       use sfx_sheba_noniterative, only: &
                       get_psi_stable_sheba => get_psi_stable
       
-            real,intent(in):: AR2(11),ARDIN(6),lat
+            real,intent(in):: AR2(11),ARDIN(6)
             real,intent(out)::ARDOUT(3)
             real,parameter::zetalim = 1. !maximum value of z/L for stable SL
             real psi_m,psi_h,psi_m_hs,psi_h_hs
@@ -252,11 +252,11 @@ module pbldia_new_sfx
             return
         end subroutine pbldia_new_sheba_noit
       
-      subroutine pbldia_new_most(AR2,ARDIN,ARDOUT,lat)
+      subroutine pbldia_new_most(AR2,ARDIN,ARDOUT)
       use sfx_most, only: &
                       get_psi_most => get_psi
       
-            real,intent(in):: AR2(11),ARDIN(6),lat
+            real,intent(in):: AR2(11),ARDIN(6)
             real,intent(out)::ARDOUT(3)
             real,parameter::zetalim = 1. !maximum value of z/L for stable SL
             real psi_m,psi_h,psi_m_hs,psi_h_hs
@@ -305,8 +305,8 @@ module pbldia_new_sfx
             return
         end subroutine pbldia_new_most
       
-      subroutine pbldia_new_esm(AR2,ARDIN,ARDOUT,lat)
-            real,intent(in):: AR2(11),ARDIN(6),lat
+      subroutine pbldia_new_esm(AR2,ARDIN,ARDOUT)
+            real,intent(in):: AR2(11),ARDIN(6)
             real,intent(out)::ARDOUT(3)
             real,parameter::zetalim = 1. !maximum value of z/L for stable SL
             real psi_m,psi_h,psi_m_hs,psi_h_hs