99 |
CEOP |
CEOP |
100 |
|
|
101 |
#ifdef ALLOW_DEBUG |
#ifdef ALLOW_DEBUG |
102 |
IF ( debugLevel .GE. debLevB ) |
IF ( debugLevel .GE. debLevB ) |
103 |
& CALL DEBUG_ENTER('DO_OCEANIC_PHYS',myThid) |
& CALL DEBUG_ENTER('DO_OCEANIC_PHYS',myThid) |
104 |
#endif |
#endif |
105 |
|
|
106 |
|
#ifdef ALLOW_THSICE |
107 |
|
IF ( useThSIce .AND. buoyancyRelation(1:7) .EQ. 'OCEANIC' ) THEN |
108 |
|
#ifdef ALLOW_DEBUG |
109 |
|
IF ( debugLevel .GE. debLevB ) |
110 |
|
& CALL DEBUG_CALL('THSICE_MAIN',myThid) |
111 |
|
#endif |
112 |
|
C-- Step forward Therm.Sea-Ice variables |
113 |
|
C and modify forcing terms including effects from ice |
114 |
|
CALL TIMER_START('THSICE_MAIN [DO_OCEANIC_PHYS]', myThid) |
115 |
|
CALL THSICE_MAIN( myTime, myIter, myThid ) |
116 |
|
CALL TIMER_STOP( 'THSICE_MAIN [DO_OCEANIC_PHYS]', myThid) |
117 |
|
ENDIF |
118 |
|
#endif /* ALLOW_THSICE */ |
119 |
|
|
120 |
|
C-- Freeze water at the surface |
121 |
|
#ifdef ALLOW_AUTODIFF_TAMC |
122 |
|
CADJ STORE theta = comlev1, key = ikey_dynamics |
123 |
|
#endif |
124 |
|
IF ( allowFreezing .AND. .NOT. useSEAICE |
125 |
|
& .AND. .NOT. useThSIce ) THEN |
126 |
|
CALL FREEZE_SURFACE( myTime, myIter, myThid ) |
127 |
|
ENDIF |
128 |
|
|
129 |
|
#ifdef COMPONENT_MODULE |
130 |
|
# ifndef ALLOW_AIM |
131 |
|
C-- Apply imported data (from coupled interface) to forcing fields |
132 |
|
C jmc: don't know precisely where to put this call (bf or af thSIce ?) |
133 |
|
IF ( useCoupler ) THEN |
134 |
|
CALL OCN_APPLY_IMPORT( .TRUE., myTime, myIter, myThid ) |
135 |
|
ENDIF |
136 |
|
# endif |
137 |
|
#endif /* COMPONENT_MODULE */ |
138 |
|
|
139 |
#ifdef ALLOW_AUTODIFF_TAMC |
#ifdef ALLOW_AUTODIFF_TAMC |
140 |
C-- dummy statement to end declaration part |
C-- dummy statement to end declaration part |
141 |
ikey = 1 |
ikey = 1 |