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

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

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

revision 1.7 by jmc, Wed Jan 21 14:35:16 2015 UTC revision 1.8 by jmc, Sun Feb 14 00:08:49 2016 UTC
# Line 109  C--   Forcing term Line 109  C--   Forcing term
109        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_U(        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_U(
110       U                       gU_arr,       U                       gU_arr,
111       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
112       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
113  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
114    
115  #ifdef ALLOW_ATM_PHYS  #ifdef ALLOW_ATM_PHYS
116        IF ( useAtm_Phys ) CALL ATM_PHYS_TENDENCY_APPLY_U(        IF ( useAtm_Phys ) CALL ATM_PHYS_TENDENCY_APPLY_U(
117       U                       gU_arr,       U                       gU_arr,
118       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
119       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
120  #endif /* ALLOW_ATM_PHYS */  #endif /* ALLOW_ATM_PHYS */
121    
122  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
123        IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_U(        IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_U(
124       U                       gU_arr,       U                       gU_arr,
125       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
126       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
127  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
128    
129  C     Ocean: Add momentum surface forcing (e.g., wind-stress) in surface level  C     Ocean: Add momentum surface forcing (e.g., wind-stress) in surface level
# Line 153  C-jmc: Without CD-scheme, this is OK ; b Line 153  C-jmc: Without CD-scheme, this is OK ; b
153           CALL TAUEDDY_TENDENCY_APPLY_U(           CALL TAUEDDY_TENDENCY_APPLY_U(
154       U                 gU_arr,       U                 gU_arr,
155       I                 iMin,iMax,jMin,jMax, k, bi,bj,       I                 iMin,iMax,jMin,jMax, k, bi,bj,
156       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
157  #endif  #endif
158    
159  #ifdef ALLOW_RBCS  #ifdef ALLOW_RBCS
# Line 161  C-jmc: Without CD-scheme, this is OK ; b Line 161  C-jmc: Without CD-scheme, this is OK ; b
161          CALL RBCS_ADD_TENDENCY(          CALL RBCS_ADD_TENDENCY(
162       U                 gU_arr,       U                 gU_arr,
163       I                 k, bi, bj, -1,       I                 k, bi, bj, -1,
164       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
165    
166        ENDIF        ENDIF
167  #endif /* ALLOW_RBCS */  #endif /* ALLOW_RBCS */
# Line 171  C-jmc: Without CD-scheme, this is OK ; b Line 171  C-jmc: Without CD-scheme, this is OK ; b
171          CALL OBCS_SPONGE_U(          CALL OBCS_SPONGE_U(
172       U                   gU_arr,       U                   gU_arr,
173       I                   iMin,iMax,jMin,jMax, k, bi,bj,       I                   iMin,iMax,jMin,jMax, k, bi,bj,
174       I                   myTime, 0, myThid )       I                   myTime, myIter, myThid )
175        ENDIF        ENDIF
176  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
177    
# Line 180  C-jmc: Without CD-scheme, this is OK ; b Line 180  C-jmc: Without CD-scheme, this is OK ; b
180          CALL MYPACKAGE_TENDENCY_APPLY_U(          CALL MYPACKAGE_TENDENCY_APPLY_U(
181       U                 gU_arr,       U                 gU_arr,
182       I                 iMin,iMax,jMin,jMax, k, bi,bj,       I                 iMin,iMax,jMin,jMax, k, bi,bj,
183       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
184        ENDIF        ENDIF
185  #endif /* ALLOW_MYPACKAGE */  #endif /* ALLOW_MYPACKAGE */
186    
# Line 287  C--   Forcing term Line 287  C--   Forcing term
287        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_V(        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_V(
288       U                       gV_arr,       U                       gV_arr,
289       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
290       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
291  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
292    
293  #ifdef ALLOW_ATM_PHYS  #ifdef ALLOW_ATM_PHYS
294        IF ( useAtm_Phys ) CALL ATM_PHYS_TENDENCY_APPLY_V(        IF ( useAtm_Phys ) CALL ATM_PHYS_TENDENCY_APPLY_V(
295       U                       gV_arr,       U                       gV_arr,
296       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
297       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
298  #endif /* ALLOW_ATM_PHYS */  #endif /* ALLOW_ATM_PHYS */
299    
300  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
301        IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_V(        IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_V(
302       U                       gV_arr,       U                       gV_arr,
303       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
304       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
305  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
306    
307  C     Ocean: Add momentum surface forcing (e.g., wind-stress) in surface level  C     Ocean: Add momentum surface forcing (e.g., wind-stress) in surface level
# Line 331  C-jmc: Without CD-scheme, this is OK ; b Line 331  C-jmc: Without CD-scheme, this is OK ; b
331           CALL TAUEDDY_TENDENCY_APPLY_V(           CALL TAUEDDY_TENDENCY_APPLY_V(
332       U                 gV_arr,       U                 gV_arr,
333       I                 iMin,iMax,jMin,jMax, k, bi,bj,       I                 iMin,iMax,jMin,jMax, k, bi,bj,
334       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
335  #endif  #endif
336    
337  #ifdef ALLOW_RBCS  #ifdef ALLOW_RBCS
# Line 339  C-jmc: Without CD-scheme, this is OK ; b Line 339  C-jmc: Without CD-scheme, this is OK ; b
339          CALL RBCS_ADD_TENDENCY(          CALL RBCS_ADD_TENDENCY(
340       U                 gV_arr,       U                 gV_arr,
341       I                 k, bi, bj, -2,       I                 k, bi, bj, -2,
342       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
343        ENDIF        ENDIF
344  #endif /* ALLOW_RBCS */  #endif /* ALLOW_RBCS */
345    
# Line 348  C-jmc: Without CD-scheme, this is OK ; b Line 348  C-jmc: Without CD-scheme, this is OK ; b
348          CALL OBCS_SPONGE_V(          CALL OBCS_SPONGE_V(
349       U                   gV_arr,       U                   gV_arr,
350       I                   iMin,iMax,jMin,jMax, k, bi,bj,       I                   iMin,iMax,jMin,jMax, k, bi,bj,
351       I                   myTime, 0, myThid )       I                   myTime, myIter, myThid )
352        ENDIF        ENDIF
353  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
354    
# Line 357  C-jmc: Without CD-scheme, this is OK ; b Line 357  C-jmc: Without CD-scheme, this is OK ; b
357          CALL MYPACKAGE_TENDENCY_APPLY_V(          CALL MYPACKAGE_TENDENCY_APPLY_V(
358       U                 gV_arr,       U                 gV_arr,
359       I                 iMin,iMax,jMin,jMax, k, bi,bj,       I                 iMin,iMax,jMin,jMax, k, bi,bj,
360       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
361        ENDIF        ENDIF
362  #endif /* ALLOW_MYPACKAGE */  #endif /* ALLOW_MYPACKAGE */
363    
# Line 469  C--   Forcing term Line 469  C--   Forcing term
469        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_T(        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_T(
470       U                       gT_arr,       U                       gT_arr,
471       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
472       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
473  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
474    
475  #ifdef ALLOW_ATM_PHYS  #ifdef ALLOW_ATM_PHYS
476        IF ( useAtm_Phys ) CALL ATM_PHYS_TENDENCY_APPLY_T(        IF ( useAtm_Phys ) CALL ATM_PHYS_TENDENCY_APPLY_T(
477       U                       gT_arr,       U                       gT_arr,
478       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
479       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
480  #endif /* ALLOW_ATM_PHYS */  #endif /* ALLOW_ATM_PHYS */
481    
482  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
483        IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_T(        IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_T(
484       U                       gT_arr,       U                       gT_arr,
485       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
486       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
487  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
488    
489  #ifdef ALLOW_ADDFLUID  #ifdef ALLOW_ADDFLUID
# Line 677  c     ENDIF Line 677  c     ENDIF
677       &     CALL FRAZIL_TENDENCY_APPLY_T(       &     CALL FRAZIL_TENDENCY_APPLY_T(
678       U                 gT_arr,       U                 gT_arr,
679       I                 iMin,iMax,jMin,jMax, k, bi,bj,       I                 iMin,iMax,jMin,jMax, k, bi,bj,
680       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
681  #endif /* ALLOW_FRAZIL */  #endif /* ALLOW_FRAZIL */
682    
683  #ifdef ALLOW_SHELFICE  #ifdef ALLOW_SHELFICE
# Line 685  c     ENDIF Line 685  c     ENDIF
685       &     CALL SHELFICE_FORCING_T(       &     CALL SHELFICE_FORCING_T(
686       U                   gT_arr,       U                   gT_arr,
687       I                   iMin,iMax,jMin,jMax, k, bi,bj,       I                   iMin,iMax,jMin,jMax, k, bi,bj,
688       I                   myTime, 0, myThid )       I                   myTime, myIter, myThid )
689  #endif /* ALLOW_SHELFICE */  #endif /* ALLOW_SHELFICE */
690    
691  #ifdef ALLOW_ICEFRONT  #ifdef ALLOW_ICEFRONT
692        IF ( useICEFRONT )        IF ( useICEFRONT )
693       &     CALL ICEFRONT_TENDENCY_APPLY_T(       &     CALL ICEFRONT_TENDENCY_APPLY_T(
694       U                   gT_arr,       U                   gT_arr,
695       I                   k, bi, bj, myTime, 0, myThid )       I                   k, bi, bj, myTime, myIter, myThid )
696  #endif /* ALLOW_ICEFRONT */  #endif /* ALLOW_ICEFRONT */
697    
698  #ifdef ALLOW_SALT_PLUME  #ifdef ALLOW_SALT_PLUME
# Line 700  c     ENDIF Line 700  c     ENDIF
700       &     CALL SALT_PLUME_TENDENCY_APPLY_T(       &     CALL SALT_PLUME_TENDENCY_APPLY_T(
701       U                     gT_arr,       U                     gT_arr,
702       I                     iMin,iMax,jMin,jMax, k, bi,bj,       I                     iMin,iMax,jMin,jMax, k, bi,bj,
703       I                     myTime, 0, myThid )       I                     myTime, myIter, myThid )
704  #endif /* ALLOW_SALT_PLUME */  #endif /* ALLOW_SALT_PLUME */
705    
706  #ifdef ALLOW_RBCS  #ifdef ALLOW_RBCS
# Line 708  c     ENDIF Line 708  c     ENDIF
708          CALL RBCS_ADD_TENDENCY(          CALL RBCS_ADD_TENDENCY(
709       U                 gT_arr,       U                 gT_arr,
710       I                 k, bi, bj, 1,       I                 k, bi, bj, 1,
711       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
712        ENDIF        ENDIF
713  #endif /* ALLOW_RBCS */  #endif /* ALLOW_RBCS */
714    
# Line 717  c     ENDIF Line 717  c     ENDIF
717          CALL OBCS_SPONGE_T(          CALL OBCS_SPONGE_T(
718       U                   gT_arr,       U                   gT_arr,
719       I                   iMin,iMax,jMin,jMax, k, bi,bj,       I                   iMin,iMax,jMin,jMax, k, bi,bj,
720       I                   myTime, 0, myThid )       I                   myTime, myIter, myThid )
721        ENDIF        ENDIF
722  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
723    
# Line 725  c     ENDIF Line 725  c     ENDIF
725        IF ( useBBL ) CALL BBL_TENDENCY_APPLY_T(        IF ( useBBL ) CALL BBL_TENDENCY_APPLY_T(
726       U                       gT_arr,       U                       gT_arr,
727       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
728       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
729  #endif /* ALLOW_BBL */  #endif /* ALLOW_BBL */
730    
731  #ifdef ALLOW_MYPACKAGE  #ifdef ALLOW_MYPACKAGE
# Line 733  c     ENDIF Line 733  c     ENDIF
733          CALL MYPACKAGE_TENDENCY_APPLY_T(          CALL MYPACKAGE_TENDENCY_APPLY_T(
734       U                 gT_arr,       U                 gT_arr,
735       I                 iMin,iMax,jMin,jMax, k, bi,bj,       I                 iMin,iMax,jMin,jMax, k, bi,bj,
736       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
737        ENDIF        ENDIF
738  #endif /* ALLOW_MYPACKAGE */  #endif /* ALLOW_MYPACKAGE */
739    
# Line 834  C--   Forcing term Line 834  C--   Forcing term
834        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_S(        IF ( useAIM ) CALL AIM_TENDENCY_APPLY_S(
835       U                       gS_arr,       U                       gS_arr,
836       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
837       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
838  #endif /* ALLOW_AIM */  #endif /* ALLOW_AIM */
839    
840  #ifdef ALLOW_ATM_PHYS  #ifdef ALLOW_ATM_PHYS
841        IF ( useAtm_Phys ) CALL ATM_PHYS_TENDENCY_APPLY_S(        IF ( useAtm_Phys ) CALL ATM_PHYS_TENDENCY_APPLY_S(
842       U                       gS_arr,       U                       gS_arr,
843       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
844       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
845  #endif /* ALLOW_ATM_PHYS */  #endif /* ALLOW_ATM_PHYS */
846    
847  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
848        IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_S(        IF ( useFIZHI ) CALL FIZHI_TENDENCY_APPLY_S(
849       U                       gS_arr,       U                       gS_arr,
850       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
851       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
852  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
853    
854  #ifdef ALLOW_ADDFLUID  #ifdef ALLOW_ADDFLUID
# Line 917  C     Ocean: Add salinity surface forcin Line 917  C     Ocean: Add salinity surface forcin
917       &     CALL SHELFICE_FORCING_S(       &     CALL SHELFICE_FORCING_S(
918       U                   gS_arr,       U                   gS_arr,
919       I                   iMin,iMax,jMin,jMax, k, bi,bj,       I                   iMin,iMax,jMin,jMax, k, bi,bj,
920       I                   myTime, 0, myThid )       I                   myTime, myIter, myThid )
921  #endif /* ALLOW_SHELFICE */  #endif /* ALLOW_SHELFICE */
922    
923  #ifdef ALLOW_ICEFRONT  #ifdef ALLOW_ICEFRONT
924        IF ( useICEFRONT )        IF ( useICEFRONT )
925       &     CALL ICEFRONT_TENDENCY_APPLY_S(       &     CALL ICEFRONT_TENDENCY_APPLY_S(
926       U                   gS_arr,       U                   gS_arr,
927       I                   k, bi, bj, myTime, 0, myThid )       I                   k, bi, bj, myTime, myIter, myThid )
928  #endif /* ALLOW_ICEFRONT */  #endif /* ALLOW_ICEFRONT */
929    
930  #ifdef ALLOW_SALT_PLUME  #ifdef ALLOW_SALT_PLUME
# Line 932  C     Ocean: Add salinity surface forcin Line 932  C     Ocean: Add salinity surface forcin
932       &     CALL SALT_PLUME_TENDENCY_APPLY_S(       &     CALL SALT_PLUME_TENDENCY_APPLY_S(
933       U                     gS_arr,       U                     gS_arr,
934       I                     iMin,iMax,jMin,jMax, k, bi,bj,       I                     iMin,iMax,jMin,jMax, k, bi,bj,
935       I                     myTime, 0, myThid )       I                     myTime, myIter, myThid )
936  #endif /* ALLOW_SALT_PLUME */  #endif /* ALLOW_SALT_PLUME */
937    
938  #ifdef ALLOW_RBCS  #ifdef ALLOW_RBCS
# Line 940  C     Ocean: Add salinity surface forcin Line 940  C     Ocean: Add salinity surface forcin
940          CALL RBCS_ADD_TENDENCY(          CALL RBCS_ADD_TENDENCY(
941       U                 gS_arr,       U                 gS_arr,
942       I                 k, bi, bj, 2,       I                 k, bi, bj, 2,
943       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
944        ENDIF        ENDIF
945  #endif /* ALLOW_RBCS */  #endif /* ALLOW_RBCS */
946    
# Line 949  C     Ocean: Add salinity surface forcin Line 949  C     Ocean: Add salinity surface forcin
949          CALL OBCS_SPONGE_S(          CALL OBCS_SPONGE_S(
950       U                   gS_arr,       U                   gS_arr,
951       I                   iMin,iMax,jMin,jMax, k, bi,bj,       I                   iMin,iMax,jMin,jMax, k, bi,bj,
952       I                   myTime, 0, myThid )       I                   myTime, myIter, myThid )
953        ENDIF        ENDIF
954  #endif /* ALLOW_OBCS */  #endif /* ALLOW_OBCS */
955    
# Line 957  C     Ocean: Add salinity surface forcin Line 957  C     Ocean: Add salinity surface forcin
957        IF ( useBBL ) CALL BBL_TENDENCY_APPLY_S(        IF ( useBBL ) CALL BBL_TENDENCY_APPLY_S(
958       U                       gS_arr,       U                       gS_arr,
959       I                       iMin,iMax,jMin,jMax, k, bi,bj,       I                       iMin,iMax,jMin,jMax, k, bi,bj,
960       I                       myTime, 0, myThid )       I                       myTime, myIter, myThid )
961  #endif /* ALLOW_BBL */  #endif /* ALLOW_BBL */
962    
963  #ifdef ALLOW_MYPACKAGE  #ifdef ALLOW_MYPACKAGE
# Line 965  C     Ocean: Add salinity surface forcin Line 965  C     Ocean: Add salinity surface forcin
965          CALL MYPACKAGE_TENDENCY_APPLY_S(          CALL MYPACKAGE_TENDENCY_APPLY_S(
966       U                 gS_arr,       U                 gS_arr,
967       I                 iMin,iMax,jMin,jMax, k, bi,bj,       I                 iMin,iMax,jMin,jMax, k, bi,bj,
968       I                 myTime, 0, myThid )       I                 myTime, myIter, myThid )
969        ENDIF        ENDIF
970  #endif /* ALLOW_MYPACKAGE */  #endif /* ALLOW_MYPACKAGE */
971    
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22