/[MITgcm]/MITgcm/model/src/external_forcing.F
ViewVC logotype

Diff of /MITgcm/model/src/external_forcing.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.10.2.2 by jmc, Fri Apr 6 16:02:52 2001 UTC revision 1.15 by adcroft, Mon Mar 25 16:17:31 2002 UTC
# Line 1  Line 1 
1  C     $Header$  C $Header$
2  C $Name$  C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6  CStartOfInterface  CBOP
7    C     !ROUTINE: EXTERNAL_FORCING_U
8    C     !INTERFACE:
9        SUBROUTINE EXTERNAL_FORCING_U(        SUBROUTINE EXTERNAL_FORCING_U(
10       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin, iMax, jMin, jMax,bi,bj,kLev,
11       I           myCurrentTime,myThid)       I           myCurrentTime,myThid)
12  C     /==========================================================\  C     !DESCRIPTION: \bv
13  C     | S/R EXTERNAL_FORCING_U                                   |  C     *==========================================================*
14  C     | o Contains problem specific forcing for zonal velocity.  |  C     | S/R EXTERNAL_FORCING_U                                    
15  C     |==========================================================|  C     | o Contains problem specific forcing for zonal velocity.  
16  C     | Adds terms to gU for forcing by external sources         |  C     *==========================================================*
17  C     | e.g. wind stress, bottom friction etc..................  |  C     | Adds terms to gU for forcing by external sources          
18  C     \==========================================================/  C     | e.g. wind stress, bottom friction etc..................  
19        IMPLICIT NONE  C     *==========================================================*
20    C     \ev
21    
22    C     !USES:
23          IMPLICIT NONE
24  C     == Global data ==  C     == Global data ==
25  #include "SIZE.h"  #include "SIZE.h"
26  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 23  C     == Global data == Line 28  C     == Global data ==
28  #include "GRID.h"  #include "GRID.h"
29  #include "DYNVARS.h"  #include "DYNVARS.h"
30  #include "FFIELDS.h"  #include "FFIELDS.h"
31    
32    C     !INPUT/OUTPUT PARAMETERS:
33  C     == Routine arguments ==  C     == Routine arguments ==
34  C     iMin - Working range of tile for applying forcing.  C     iMin - Working range of tile for applying forcing.
35  C     iMax  C     iMax
# Line 32  C     kLev Line 39  C     kLev
39        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
40        _RL myCurrentTime        _RL myCurrentTime
41        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
42    
43    C     !LOCAL VARIABLES:
44  C     == Local variables ==  C     == Local variables ==
45  C     Loop counters  C     Loop counters
46        INTEGER I, J        INTEGER I, J
47    CEOP
48    
49  C--   Forcing term  C--   Forcing term
50  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
# Line 50  C     Add windstress momentum impulse in Line 58  C     Add windstress momentum impulse in
58         ENDDO         ENDDO
59        ENDIF        ENDIF
60    
61    #ifdef ALLOW_OBCS
62    C     IF (useOBCS) THEN
63    C      CALL OBCS_SPONGE_U(
64    C    I           iMin, iMax, jMin, jMax,bi,bj,kLev,
65    C    I           myCurrentTime,myThid)
66    C     ENDIF
67    #endif
68    
69        RETURN        RETURN
70        END        END
71  CStartOfInterface  CBOP
72    C     !ROUTINE: EXTERNAL_FORCING_V
73    C     !INTERFACE:
74        SUBROUTINE EXTERNAL_FORCING_V(        SUBROUTINE EXTERNAL_FORCING_V(
75       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin, iMax, jMin, jMax,bi,bj,kLev,
76       I           myCurrentTime,myThid)       I           myCurrentTime,myThid)
77  C     /==========================================================\  C     !DESCRIPTION: \bv
78  C     | S/R EXTERNAL_FORCING_V                                   |  C     *==========================================================*
79  C     | o Contains problem specific forcing for merid velocity.  |  C     | S/R EXTERNAL_FORCING_V                                    
80  C     |==========================================================|  C     | o Contains problem specific forcing for merid velocity.  
81  C     | Adds terms to gV for forcing by external sources         |  C     *==========================================================*
82  C     | e.g. wind stress, bottom friction etc..................  |  C     | Adds terms to gV for forcing by external sources          
83  C     \==========================================================/  C     | e.g. wind stress, bottom friction etc..................  
84        IMPLICIT NONE  C     *==========================================================*
85    C     \ev
86    
87    C     !USES:
88          IMPLICIT NONE
89  C     == Global data ==  C     == Global data ==
90  #include "SIZE.h"  #include "SIZE.h"
91  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 73  C     == Global data == Line 94  C     == Global data ==
94  #include "DYNVARS.h"  #include "DYNVARS.h"
95  #include "FFIELDS.h"  #include "FFIELDS.h"
96    
97    C     !INPUT/OUTPUT PARAMETERS:
98  C     == Routine arguments ==  C     == Routine arguments ==
99  C     iMin - Working range of tile for applying forcing.  C     iMin - Working range of tile for applying forcing.
100  C     iMax  C     iMax
# Line 83  C     kLev Line 104  C     kLev
104        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
105        _RL myCurrentTime        _RL myCurrentTime
106        INTEGER myThid        INTEGER myThid
107  CEndOfInterface  
108    C     !LOCAL VARIABLES:
109  C     == Local variables ==  C     == Local variables ==
110  C     Loop counters  C     Loop counters
111        INTEGER I, J        INTEGER I, J
112    CEOP
113    
114  C--   Forcing term  C--   Forcing term
115  C     Add windstress momentum impulse into the top-layer  C     Add windstress momentum impulse into the top-layer
# Line 100  C     Add windstress momentum impulse in Line 123  C     Add windstress momentum impulse in
123         ENDDO         ENDDO
124        ENDIF        ENDIF
125    
126    #ifdef ALLOW_OBCS
127    C     IF (useOBCS) THEN
128    C      CALL OBCS_SPONGE_V(
129    C    I           iMin, iMax, jMin, jMax,bi,bj,kLev,
130    C    I           myCurrentTime,myThid)
131    C     ENDIF
132    #endif
133    
134        RETURN        RETURN
135        END        END
136  CStartOfInterface  CBOP
137    C     !ROUTINE: EXTERNAL_FORCING_T
138    C     !INTERFACE:
139        SUBROUTINE EXTERNAL_FORCING_T(        SUBROUTINE EXTERNAL_FORCING_T(
140       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin, iMax, jMin, jMax,bi,bj,kLev,
141       I           myCurrentTime,myThid)       I           myCurrentTime,myThid)
142  C     /==========================================================\  C     !DESCRIPTION: \bv
143  C     | S/R EXTERNAL_FORCING_T                                   |  C     *==========================================================*
144  C     | o Contains problem specific forcing for temperature.     |  C     | S/R EXTERNAL_FORCING_T                                    
145  C     |==========================================================|  C     | o Contains problem specific forcing for temperature.      
146  C     | Adds terms to gT for forcing by external sources         |  C     *==========================================================*
147  C     | e.g. heat flux, climatalogical relaxation..............  |  C     | Adds terms to gT for forcing by external sources          
148  C     \==========================================================/  C     | e.g. heat flux, climatalogical relaxation..............  
149        IMPLICIT NONE  C     *==========================================================*
150    C     \ev
151    
152    C     !USES:
153          IMPLICIT NONE
154  C     == Global data ==  C     == Global data ==
155  #include "SIZE.h"  #include "SIZE.h"
156  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 129  C     == Global data == Line 165  C     == Global data ==
165        _RL swfracb(two)        _RL swfracb(two)
166  #endif  #endif
167    
168    C     !INPUT/OUTPUT PARAMETERS:
169  C     == Routine arguments ==  C     == Routine arguments ==
170  C     iMin - Working range of tile for applying forcing.  C     iMin - Working range of tile for applying forcing.
171  C     iMax  C     iMax
# Line 141  C     kLev Line 177  C     kLev
177        INTEGER myThid        INTEGER myThid
178  CEndOfInterface  CEndOfInterface
179    
180    C     !LOCAL VARIABLES:
181  C     == Local variables ==  C     == Local variables ==
182  C     Loop counters  C     Loop counters
183        INTEGER I, J        INTEGER I, J
184    CEOP
185    
186  C--   Forcing term  C--   Forcing term
187  C     Add heat in top-layer  C     Add heat in top-layer
# Line 172  C Penetrating SW radiation Line 210  C Penetrating SW radiation
210         ENDDO         ENDDO
211        ENDDO        ENDDO
212  #endif  #endif
213    
214    #ifdef ALLOW_OBCS
215    C     IF (useOBCS) THEN
216    C      CALL OBCS_SPONGE_T(
217    C    I           iMin, iMax, jMin, jMax,bi,bj,kLev,
218    C    I           myCurrentTime,myThid)
219    C     ENDIF
220    #endif
221    
222        RETURN        RETURN
223        END        END
224  CStartOfInterface  CBOP
225    C     !ROUTINE: EXTERNAL_FORCING_S
226    C     !INTERFACE:
227        SUBROUTINE EXTERNAL_FORCING_S(        SUBROUTINE EXTERNAL_FORCING_S(
228       I           iMin, iMax, jMin, jMax,bi,bj,kLev,       I           iMin, iMax, jMin, jMax,bi,bj,kLev,
229       I           myCurrentTime,myThid)       I           myCurrentTime,myThid)
 C     /==========================================================\  
 C     | S/R EXTERNAL_FORCING_S                                   |  
 C     | o Contains problem specific forcing for merid velocity.  |  
 C     |==========================================================|  
 C     | Adds terms to gS for forcing by external sources         |  
 C     | e.g. fresh-water flux, climatalogical relaxation.......  |  
 C     \==========================================================/  
       IMPLICIT NONE  
230    
231    C     !DESCRIPTION: \bv
232    C     *==========================================================*
233    C     | S/R EXTERNAL_FORCING_S                                    
234    C     | o Contains problem specific forcing for merid velocity.  
235    C     *==========================================================*
236    C     | Adds terms to gS for forcing by external sources          
237    C     | e.g. fresh-water flux, climatalogical relaxation.......  
238    C     *==========================================================*
239    C     \ev
240    
241    C     !USES:
242          IMPLICIT NONE
243  C     == Global data ==  C     == Global data ==
244  #include "SIZE.h"  #include "SIZE.h"
245  #include "EEPARAMS.h"  #include "EEPARAMS.h"
# Line 195  C     == Global data == Line 248  C     == Global data ==
248  #include "DYNVARS.h"  #include "DYNVARS.h"
249  #include "FFIELDS.h"  #include "FFIELDS.h"
250    
251    C     !INPUT/OUTPUT PARAMETERS:
252  C     == Routine arguments ==  C     == Routine arguments ==
253  C     iMin - Working range of tile for applying forcing.  C     iMin - Working range of tile for applying forcing.
254  C     iMax  C     iMax
# Line 204  C     kLev Line 258  C     kLev
258        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj        INTEGER iMin, iMax, jMin, jMax, kLev, bi, bj
259        _RL myCurrentTime        _RL myCurrentTime
260        INTEGER myThid        INTEGER myThid
 CEndOfInterface  
261    
262    C     !LOCAL VARIABLES:
263  C     == Local variables ==  C     == Local variables ==
264  C     Loop counters  C     Loop counters
265        INTEGER I, J        INTEGER I, J
266    CEOP
267    
268  C--   Forcing term  C--   Forcing term
269  C     Add fresh-water in top-layer  C     Add fresh-water in top-layer
# Line 221  C     Add fresh-water in top-layer Line 276  C     Add fresh-water in top-layer
276         ENDDO         ENDDO
277        ENDIF        ENDIF
278    
279    #ifdef ALLOW_OBCS
280    C     IF (useOBCS) THEN
281    C      CALL OBCS_SPONGE_S(
282    C    I           iMin, iMax, jMin, jMax,bi,bj,kLev,
283    C    I           myCurrentTime,myThid)
284    C     ENDIF
285    #endif
286    
287        RETURN        RETURN
288        END        END
Legend:
Removed from v.1.10.2.2  
changed lines
  Added in v.1.15
  ViewVC Help
Powered by ViewVC 1.1.22