--- MITgcm_contrib/jscott/igsm/src/surface.F	2006/08/22 20:25:52	1.2
+++ MITgcm_contrib/jscott/igsm/src/surface.F	2007/04/23 21:20:18	1.3
@@ -44,7 +44,7 @@
 #include "BD2G04.COM"
 
 #if ( defined CLM )
-#include "CLM.COM"
+#include "CLM.h"
 #endif
 
       COMMON/SPEC2/KM,KINC,COEK,C3LAND(IO0,JM0),C3OICE(IO0,JM0)         5808.1  
@@ -63,7 +63,8 @@
       DIMENSION WMGMINO(JM0)
       LOGICAL POLE,PRNT,HPRNT
       common/conprn/HPRNT          
-      common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)
+!     common/TSUR/TSURFC(JM0,0:13),TSURFT(JM0),TSURFD(JM0),DTEMSR(JM0)
+#include "TSRF.COM"
       common/SURRAD/TRSURF(JM0,4),SRSURF(JM0,4)
 c     REAL*8 B,TGV,TKV,TSV0,TSV1,TSV                                    5818.   
       COMMON/CWMG/WMGEA(JM0),NWMGEA(JM0),CHAVER(JM0),DTAV(JM0),DQAV(JM0)
@@ -323,12 +324,13 @@
       IM1=IM                                                            5969.   
 #if ( defined CLM )
          if(NS.eq.1)then
-           tsl4clm(j)=0.0
-           qs4clm(j)=0.0
-           ps4clm(j)=0.0
-           ws4clm(j)=0.0
-           us4clm(j)=0.0
-           vs4clm(j)=0.0
+           i=1
+           tsl4clm(i,j)=0.0
+           qs4clm(i,j)=0.0
+           ps4clm(i,j)=0.0
+           ws4clm(i,j)=0.0
+           us4clm(i,j)=0.0
+           vs4clm(i,j)=0.0
          endif
 #endif
       DO 6000 I=1,IMAX                                                  5970.   
@@ -595,6 +597,13 @@
         print *,'TG=',TG,' QG=',QG
         print *,'RVX=',RVX,' TG1=',TG1
       endif
+!      if(TAU.gt.3623.0.and.TAU.lt.3744.0) then
+!        if(ITYPE.eq.4)then
+!         if(J.eq.29)then
+!           write (329) ,TAU,TOFDAY,TG1
+!         endif
+!         endif
+!       endif
 
 #if ( defined CPL_OCEANCO2  || defined OCEAN_3D)
       if(ITYPE.eq.1)then
@@ -677,12 +686,12 @@
 #if ( defined CLM )
          if(NS.eq.1)then
           if(ITYPE.EQ.4.or.ITYPE.EQ.3)then
-           tsl4clm(j)=tsl4clm(j)+TS*PTYPE/PLAND
-           qs4clm(j)=qs4clm(j)+QS*PTYPE/PLAND
-           ps4clm(j)=ps4clm(j)+PS*PTYPE/PLAND
-           ws4clm(j)=ws4clm(j)+WS*PTYPE/PLAND
-           us4clm(j)=us4clm(j)+US*PTYPE/PLAND
-           vs4clm(j)=vs4clm(j)+VS*PTYPE/PLAND
+           tsl4clm(i,j)=tsl4clm(i,j)+TS*PTYPE/PLAND
+           qs4clm(i,j)=qs4clm(i,j)+QS*PTYPE/PLAND
+           ps4clm(i,j)=ps4clm(i,j)+PS*PTYPE/PLAND
+           ws4clm(i,j)=ws4clm(i,j)+WS*PTYPE/PLAND
+           us4clm(i,j)=us4clm(i,j)+US*PTYPE/PLAND
+           vs4clm(i,j)=vs4clm(i,j)+VS*PTYPE/PLAND
          endif
         endif
 #endif