--- MITgcm/pkg/thsice/thsice_ini_vars.F	2006/05/25 18:03:24	1.6
+++ MITgcm/pkg/thsice/thsice_ini_vars.F	2006/06/25 22:35:17	1.7
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/thsice/thsice_ini_vars.F,v 1.6 2006/05/25 18:03:24 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/thsice/thsice_ini_vars.F,v 1.7 2006/06/25 22:35:17 jmc Exp $
 C $Name:  $
 
 #include "THSICE_OPTIONS.h"
@@ -23,7 +23,7 @@
 #include "EEPARAMS.h"
 #include "PARAMS.h"
 #include "GRID.h"
-c #include "DYNVARS.h"
+#include "FFIELDS.h"
 #include "THSICE_PARAMS.h"
 #include "THSICE_VARS.h"
 
@@ -135,6 +135,23 @@
         _EXCH_XY_R8(dFdT, myThid)
 #endif
 
+#ifdef ATMOSPHERIC_LOADING
+C--     Initialise Sea-Ice Loading for SeaIce-Dynamics :
+        IF ( useSEAICE ) THEN
+          DO bj = myByLo(myThid), myByHi(myThid)
+           DO bi = myBxLo(myThid), myBxHi(myThid)
+            DO j=1-OLy,sNy+OLy
+             DO i=1-OLx,sNx+OLx
+               sIceLoad(i,j,bi,bj) = ( snowHeight(i,j,bi,bj)*rhos
+     &                               + iceHeight(i,j,bi,bj)*rhoi
+     &                               )*iceMask(i,j,bi,bj)
+             ENDDO
+            ENDDO
+           ENDDO
+          ENDDO
+        ENDIF
+#endif /* ATMOSPHERIC_LOADING */
+
 #endif /* ALLOW_THSICE */
 
       RETURN