--- MITgcm/pkg/dic/dic_surfforcing.F	2003/10/06 20:11:10	1.3
+++ MITgcm/pkg/dic/dic_surfforcing.F	2004/07/13 18:03:31	1.6
@@ -1,20 +1,24 @@
-#include "CPP_OPTIONS.h"
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/dic/dic_surfforcing.F,v 1.6 2004/07/13 18:03:31 jmc Exp $
+C $Name:  $
+
+#include "DIC_OPTIONS.h"
 #include "PTRACERS_OPTIONS.h"
 #include "GCHEM_OPTIONS.h"
 
-CStartOfInterFace
+CBOP
+C !ROUTINE: DIC_SURFFORCING
+
+C !INTERFACE: ==========================================================
       SUBROUTINE DIC_SURFFORCING( PTR_CO2 , GDC, 
      I           bi,bj,imin,imax,jmin,jmax,
      I           myIter,myTime,myThid)
 
-C     /==========================================================\
-C     | SUBROUTINE DIC_SURFFORCING                               |
-C     | o Calculate the carbon air-sea flux terms                |
-C     | o following external_forcing_dic.F from Mick             |
-C     |==========================================================|
-      IMPLICIT NONE
+C !DESCRIPTION:
+C  Calculate the carbon air-sea flux terms              
+C  following external_forcing_dic.F (OCMIP run) from Mick            
 
-C     == GLobal variables ==
+C !USES: ===============================================================
+      IMPLICIT NONE
 #include "SIZE.h"
 #include "DYNVARS.h"
 #include "EEPARAMS.h"
@@ -23,18 +27,27 @@
 #include "FFIELDS.h"
 #include "DIC_ABIOTIC.h"
 #ifdef DIC_BIOTIC
+#include "PTRACERS_SIZE.h"
 #include "PTRACERS.h"
 #endif
 
-C     == Routine arguments ==
+C !INPUT PARAMETERS: ===================================================
+C  myThid               :: thread number
+C  myIter               :: current timestep
+C  myTime               :: current time
+c  PTR_CO2              :: DIC tracer field
       INTEGER myIter, myThid
       _RL myTime
       _RL  PTR_CO2(1-OLx:sNx+OLx,1-OLy:sNy+OLy,Nr)
-      _RL  GDC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
       INTEGER iMin,iMax,jMin,jMax, bi, bj
 
+C !OUTPUT PARAMETERS: ===================================================
+c GDC                   :: tendency term due to air-sea exchange
+      _RL  GDC(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+
 #ifdef ALLOW_PTRACERS
-C     == Local variables ==
+
+C !LOCAL VARIABLES: ====================================================
        INTEGER I,J, kLev, it
 C Number of iterations for pCO2 solvers...
 C Solubility relation coefficients
@@ -46,6 +59,7 @@
       _RL surfphos(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
       _RL surfsi(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
       _RL VirtualFlux(1-OLx:sNx+OLx,1-OLy:sNy+OLy)
+CEOP
 
 cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc