--- MITgcm/model/src/packages_init_variables.F	2003/11/23 01:28:05	1.28
+++ MITgcm/model/src/packages_init_variables.F	2008/09/10 09:18:57	1.69
@@ -1,7 +1,8 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/packages_init_variables.F,v 1.28 2003/11/23 01:28:05 jmc Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/src/packages_init_variables.F,v 1.69 2008/09/10 09:18:57 mlosch Exp $
 C $Name:  $
 
 #include "PACKAGES_CONFIG.h"
+#include "AD_CONFIG.h"
 #include "CPP_OPTIONS.h"
 
 CBOP
@@ -11,40 +12,78 @@
 
 C     !DESCRIPTION: \bv
 C     *==========================================================*
-C     | SUBROUTINE PACKAGES_INIT_VARIABLES                        
-C     | o Does initialisation of package-related variable data    
+C     | SUBROUTINE PACKAGES_INIT_VARIABLES
+C     | o Does initialisation of package-related variable data
 C     *==========================================================*
 C     \ev
 
 C     !CALLING SEQUENCE:
 C     PACKAGES_INIT_VARIABLES
 C       |
+C       |-- DIAGNOSTICS_INIT_VARIA
+C       |
+C       |-- GAD_GAD_INIT_VARIA
+C       |
+C       |-- CD_CODE_INI_VARS
+C       |
 C       |-- GMREDI_INIT
 C       |
+C       |-- DWNSLP_INIT_VARIA
+C       |
 C       |-- KPP_INIT
 C       |-- KPP_OPEN_DIAGS
 C       |
-C       |-- SEAICE_INIT
-C       | 
-C       |-- OBCS_INIT_VARIABLES
+C       |-- PP81_INIT
+C       |
+C       |-- MY82_INIT_VARIA
+C       |
+C       |-- GGL90_INIT
 C       |
-C       |-- PTRACERS_INIT
+C       |-- SEAICE_INIT_VARIA
+C       |
+C       |-- SHELFICE_INIT_VARIA
+C       |
+C       |-- PTRACERS_INIT_VARIA
+C       |
+C       |-- GCHEM_INIT_VARI
 C       |
 C       |-- LAND_INI_VARS
 C       |
+C       |-- CTRL_INIT_VARIABLES
+C       |-- CTRL_MAP_INI_ECCO
 C       |-- CTRL_MAP_INI
 C       |
 C       |-- EXF_INIT
 C       |
-C       |-- INI_FORCING
+C       |-- EBM_INI_VARS
 C       |
-C       |-- ECCO_INIT
+C       |-- COST_INIT_VARIA
 C       |
-C       |-- COST_INIT
+C       |-- PROFILES_INIT_VARIA
 C       |
-C       |-- BULKF_INIT
+C       |-- OBCS_INIT_VARIABLES
+C       |
+C       |-- BULKF_INIT_VARIA
 C       |
 C       |-- THSICE_INI_VARS
+C       |
+C       |-- CPL_INI_VARS
+C       |
+C       |-- ATM2D_INIT_VARS
+C       |
+C       |-- FIZHI_INI_VARS
+C       |
+C       |-- MATRIX_INIT
+C       |
+C       |-- RBCS_INIT_VARIA
+C       |
+C       |-- REGRID_INIT_VARIA
+C       |
+C       |-- SALT_PLUME_INIT_VARIA
+C       |
+C       |-- CHEAPAML_INIT_VARIA
+C       |
+C       |-- MYPACKAGE_INIT_VARIA
 
 C     !USES:
       IMPLICIT NONE
@@ -52,6 +91,11 @@
 #include "SIZE.h"
 #include "EEPARAMS.h"
 #include "PARAMS.h"
+#ifdef ALLOW_AUTODIFF_TAMC
+# include "DYNVARS.h"
+# include "tamc.h"
+# include "tamc_keys.h"
+#endif
 
 C     !INPUT/OUTPUT PARAMETERS:
 C     === Routine arguments ===
@@ -59,23 +103,115 @@
       INTEGER myThid
 CEOP
 
+#ifdef ALLOW_DEBUG
+      IF (debugMode)
+     &     CALL DEBUG_ENTER('PACKAGES_INIT_VARIABLES',myThid)
+#endif /* ALLOW_DEBUG */
+
+#ifdef ALLOW_DIAGNOSTICS
+      IF ( useDiagnostics ) THEN
+        CALL DIAGNOSTICS_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_DIAGNOSTICS */
+
+#ifdef ALLOW_GENERIC_ADVDIFF
+      IF ( useGAD ) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('GAD_INIT_VARIA',myThid)
+# endif
+        CALL GAD_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_GENERIC_ADVDIFF */
+
+#ifdef ALLOW_CD_CODE
+C--   Initialize CD_CODE variables:
+C- note(jmc): before packaging CD_CODE, was done within ini_fields (=called before),
+C             therefore call CD-ini-vars before others pkg.
+      IF (useCDscheme) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('CD_CODE_INI_VARS',myThid)
+# endif
+        CALL CD_CODE_INI_VARS( myThid )
+      ENDIF
+#endif /* ALLOW_CD_CODE */
+
 #ifdef ALLOW_GMREDI
 C--   Initialize GM/Redi parameterization
-      IF (useGMRedi) CALL GMREDI_INIT( myThid )
-#endif
+      IF (useGMRedi) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('GMREDI_INIT',myThid)
+# endif
+        CALL GMREDI_INIT( myThid )
+      ENDIF
+#endif /* ALLOW_GMREDI */
+
+#ifdef ALLOW_DOWN_SLOPE
+      IF ( useDOWN_SLOPE ) THEN
+        CALL DWNSLP_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_DOWN_SLOPE */
 
 #ifdef ALLOW_KPP
 C--   Initialize KPP vertical mixing scheme.
       IF (useKPP) THEN
+# ifdef ALLOW_DEBUG
+       IF (debugMode)
+     &       CALL DEBUG_CALL('KPP_INIT + KPP_OPEN_DIAGS',myThid)
+# endif
        CALL KPP_INIT( myThid )
        CALL KPP_OPEN_DIAGS( myThid )
       ENDIF
-#endif
+#endif /* ALLOW_KPP */
+
+#ifdef ALLOW_PP81
+C--   Initialize PP81 vertical mixing scheme.
+      IF (usePP81) THEN
+# ifdef ALLOW_DEBUG
+       IF (debugMode) CALL DEBUG_CALL('PP81_INIT',myThid)
+# endif
+       CALL PP81_INIT( myThid )
+      ENDIF
+#endif /* ALLOW_PP81 */
+
+#ifdef ALLOW_MY82
+C--   Initialize MY82 vertical mixing scheme.
+      IF (useMY82) THEN
+       CALL MY82_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_MY82 */
+
+#ifdef ALLOW_GGL90
+C--   Initialize GGL90 vertical mixing scheme.
+      IF (useGGL90) THEN
+# ifdef ALLOW_DEBUG
+       IF (debugMode) CALL DEBUG_CALL('GGL90_INIT',myThid)
+# endif
+       CALL GGL90_INIT( myThid )
+      ENDIF
+#endif /* ALLOW_GGL90 */
 
 #ifdef ALLOW_SEAICE
 C--   Initialize SEAICE model.
-      IF (useSEAICE) CALL SEAICE_INIT( myThid )
-#endif
+cph# ifndef ALLOW_AUTODIFF_TAMC
+      IF (useSEAICE) THEN
+cph# endif
+# ifdef ALLOW_DEBUG
+       IF (debugMode) CALL DEBUG_CALL('SEAICE_INIT_VARIA',myThid)
+# endif
+       CALL SEAICE_INIT_VARIA( myThid )
+cph# ifndef ALLOW_AUTODIFF_TAMC
+      ENDIF
+cph# endif
+#endif /* ALLOW_SEAICE */
+
+#ifdef ALLOW_SHELFICE
+      IF (useShelfIce) THEN
+# ifdef ALLOW_DEBUG
+       IF (debugMode) CALL DEBUG_CALL('SHELFICE_INIT_VARIA',myThid)
+# endif
+       CALL SHELFICE_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_SHELFICE */
 
 cph-- moved further down
 cph#ifdef ALLOW_OBCS
@@ -89,53 +225,181 @@
 # ifndef ALLOW_AUTODIFF_TAMC
       IF ( usePTRACERS ) THEN
 # endif
-         CALL PTRACERS_INIT( mythid )
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('PTRACERS_INIT_VARIA',myThid)
+# endif
+        CALL PTRACERS_INIT_VARIA( myThid )
+# ifndef ALLOW_AUTODIFF_TAMC
+      ENDIF
+# endif
+#endif /* ALLOW_PTRACERS */
+
 #ifdef ALLOW_GCHEM
-         CALL GCHEM_INIT_VARI( mythid )
-#endif
+# ifndef ALLOW_AUTODIFF_TAMC
+      IF (useGCHEM) THEN
+# endif
+        CALL GCHEM_INIT_VARI( myThid )
 # ifndef ALLOW_AUTODIFF_TAMC
       ENDIF
 # endif
-#endif
+#endif /* ALLOW_GCHEM */
 
 #ifdef ALLOW_LAND
-      IF ( useLAND ) CALL LAND_INI_VARS( mythid )
-#endif
+      IF ( useLAND ) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('LAND_INI_VARS',myThid)
+# endif
+        CALL LAND_INI_VARS( myThid )
+      ENDIF
+#endif /* ALLOW_LAND */
+
+#ifdef ALLOW_SMOOTH
+      CALL SMOOTH_INIT_VARIA(myThid)
+#endif /* ALLOW_SMOOTH */
 
 #ifdef ALLOW_AUTODIFF
+CADJ STORE theta = tapelev_init, key = 1
 c--   Initialise auxiliary xx_ fields
-      CALL CTRL_INIT_VARIABLES ( mythid )
+      IF (debugMode) CALL DEBUG_CALL('CTRL_INIT_VARIABLES',myThid)
+      CALL CTRL_INIT_VARIABLES ( myThid )
 c--   Map the control variables onto the model state.
-      CALL CTRL_MAP_INI( mythid )
+# ifdef ALLOW_ECCO
+      IF (debugMode) CALL DEBUG_CALL('CTRL_MAP_INI_ECCO',myThid)
+      CALL CTRL_MAP_INI_ECCO( myThid )
+# else
+      IF (debugMode) CALL DEBUG_CALL('CTRL_MAP_INI',myThid)
+      CALL CTRL_MAP_INI( myThid )
+# endif
       _BARRIER
-#endif
+#endif /* ALLOW_AUTODIFF */
 
 #ifdef ALLOW_EXF
-      CALL EXF_INIT( mythid )
-#else
-      CALL INI_FORCING( mythid )
-#endif
+      IF (useEXF) THEN
+# ifdef ALLOW_DEBUG
+       IF (debugMode) CALL DEBUG_CALL('EXF_INIT',myThid)
+# endif
+       CALL EXF_INIT( myThid )
+      ENDIF
+#endif /* ALLOW_EXF */
+
+#ifdef ALLOW_EBM
+# ifdef ALLOW_AUTODIFF
+CADJ STORE theta = tapelev_init, key = 1
+# endif
+      IF (useEBM) THEN
+        CALL EBM_INI_VARS( myThid )
+      ENDIF
+#endif /* ALLOW_EBM */
 
 #ifdef ALLOW_COST
 c--   Initialise the cost function.
 ceh3 needs an IF ( useCOST ) THEN
-      CALL COST_INITVARIA( mythid )
+      CALL COST_INIT_VARIA( myThid )
       _BARRIER
-#endif
+#endif /* ALLOW_COST */
+
+#ifdef ALLOW_PROFILES
+c--   Initialise the cost function.
+      CALL PROFILES_INIT_VARIA( myThid )
+      _BARRIER
+#endif /* ALLOW_PROFILES */
 
 #ifdef ALLOW_OBCS
 C--   Open boundaries data
-      IF (useOBCS) CALL OBCS_INIT_VARIABLES( myThid )
-#endif
+      IF (useOBCS) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('OBCS_INIT_VARIABLES',myThid)
+# endif
+        CALL OBCS_INIT_VARIABLES( myThid )
+      ENDIF
+#endif /* ALLOW_OBCS */
 
 #ifdef ALLOW_BULK_FORCE
-      IF (useBulkForce) CALL BULKF_INIT( mythid)
-#endif
+      IF (useBulkForce) THEN
+        CALL BULKF_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_BULK_FORCE */
 
 #ifdef ALLOW_THSICE
-      IF (useThSIce) CALL THSICE_INI_VARS( mythid)
-#endif
+      IF (useThSIce) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('THSICE_INI_VARS',myThid)
+# endif
+        CALL THSICE_INI_VARS( myThid)
+      ENDIF
+#endif /* ALLOW_THSICE */
+
+#ifdef COMPONENT_MODULE
+      IF (useCoupler) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('CPL_INI_VARS',myThid)
+# endif
+        CALL CPL_INI_VARS( myThid )
+      ENDIF
+#endif /* COMPONENT_MODULE */
+
+#ifdef ALLOW_ATM2D
+      IF (useAtm2d) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('ATM2D_INIT_VARS',myThid)
+# endif
+        CALL ATM2D_INIT_VARS( myThid )
+      ENDIF
+#endif /* ALLOW_ATM2D */
+
+#ifdef ALLOW_FIZHI
+C Initialize FIZHI state variables
+      IF (useFIZHI) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('FIZHI_INIT_VARS',myThid)
+# endif
+        CALL FIZHI_INIT_VARS( myThid )
+      ENDIF
+#endif /* ALLOW_FIZHI */
+
+#ifdef ALLOW_MATRIX
+      IF ( useMATRIX ) THEN
+# ifdef ALLOW_DEBUG
+        IF (debugMode) CALL DEBUG_CALL('MATRIX_INIT',myThid)
+# endif
+        CALL MATRIX_INIT( myThid )
+      ENDIF
+#endif /* ALLOW_MATRIX */
+
+#ifdef ALLOW_RBCS
+      IF ( useRBCS ) THEN
+        CALL RBCS_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_RBCS */
+
+#ifdef ALLOW_REGRID
+      IF ( useREGRID ) THEN
+        CALL REGRID_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_REGRID */
+
+#ifdef ALLOW_SALT_PLUME
+      IF ( useSALT_PLUME ) THEN
+        CALL SALT_PLUME_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_SALT_PLUME */
+
+#ifdef ALLOW_CHEAPAML
+      IF (useCheapAML) THEN
+        CALL CHEAPAML_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_CHEAPAML */
+
+#ifdef ALLOW_MYPACKAGE
+      IF ( useMYPACKAGE ) THEN
+        CALL MYPACKAGE_INIT_VARIA( myThid )
+      ENDIF
+#endif /* ALLOW_MYPACKAGE */
 
+#ifdef ALLOW_DEBUG
+      IF (debugMode)
+     &     CALL DEBUG_LEAVE('PACKAGES_INIT_VARIABLES',myThid)
+#endif /* ALLOW_DEBUG */
 
       RETURN
       END