/[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.53 by gforget, Fri May 30 02:45:43 2008 UTC revision 1.57 by dimitri, Mon Apr 26 20:37:02 2010 UTC
# Line 287  C--   Forcing term Line 287  C--   Forcing term
287       &                      myTime, myThid )       &                      myTime, myThid )
288  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
289    
290    #ifdef ALLOW_ADDFLUID
291          IF ( selectAddFluid.NE.0 .AND. temp_addMass.NE.UNSET_RL ) THEN
292           IF ( ( selectAddFluid.GE.1 .AND. nonlinFreeSurf.GT.0 )
293         &      .OR. convertFW2Salt.EQ.-1. _d 0 ) THEN
294             DO j=1,sNy
295              DO i=1,sNx
296                gT(i,j,kLev,bi,bj) = gT(i,j,kLev,bi,bj)
297         &        + addMass(i,j,kLev,bi,bj)*mass2rUnit
298         &          *( temp_addMass - theta(i,j,kLev,bi,bj) )
299         &          *recip_rA(i,j,bi,bj)
300         &          *recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj)
301    C    &          *recip_deepFac2C(kLev)*recip_rhoFacC(kLev)
302              ENDDO
303             ENDDO
304           ELSE
305             DO j=1,sNy
306              DO i=1,sNx
307                gT(i,j,kLev,bi,bj) = gT(i,j,kLev,bi,bj)
308         &        + addMass(i,j,kLev,bi,bj)*mass2rUnit
309         &          *( temp_addMass - tRef(kLev) )
310         &          *recip_rA(i,j,bi,bj)
311         &          *recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj)
312    C    &          *recip_deepFac2C(kLev)*recip_rhoFacC(kLev)
313              ENDDO
314             ENDDO
315           ENDIF
316          ENDIF
317    #endif /* ALLOW_ADDFLUID */
318    
319  C     Add heat in top-layer  C     Add heat in top-layer
320        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
321         DO j=1,sNy         DO j=1,sNy
# Line 298  C     Add heat in top-layer Line 327  C     Add heat in top-layer
327         ENDDO         ENDDO
328        ENDIF        ENDIF
329    
 cph#ifndef ALLOW_AUTODIFF_TAMC  
 cph I didnt put this ifndef here.  
330        IF (linFSConserveTr) THEN        IF (linFSConserveTr) THEN
331         DO j=1,sNy         DO j=1,sNy
332          DO i=1,sNx          DO i=1,sNx
# Line 310  cph I didnt put this ifndef here. Line 337  cph I didnt put this ifndef here.
337          ENDDO          ENDDO
338         ENDDO         ENDDO
339        ENDIF        ENDIF
 cph#endif /* ndfef ALLOW_AUTODIFF_TAMC */  
340    
341  #ifdef ALLOW_SHELFICE  #ifdef ALLOW_SHELFICE
342        IF ( useShelfIce )        IF ( useShelfIce )
# Line 319  cph#endif /* ndfef ALLOW_AUTODIFF_TAMC * Line 345  cph#endif /* ndfef ALLOW_AUTODIFF_TAMC *
345       I     myTime, myThid )       I     myTime, myThid )
346  #endif /* ALLOW_SHELFICE */  #endif /* ALLOW_SHELFICE */
347    
348    #ifdef ALLOW_ICEFRONT
349          IF ( useICEFRONT )
350         &     CALL ICEFRONT_TENDENCY_APPLY_T(
351         &     bi,bj, kLev, myTime, myThid )
352    #endif /* ALLOW_ICEFRONT */
353    
354  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
355  C Penetrating SW radiation  C Penetrating SW radiation
356  c     IF ( usePenetratingSW ) THEN  c     IF ( usePenetratingSW ) THEN
# Line 439  C--   Forcing term Line 471  C--   Forcing term
471       &                      myTime, myThid )       &                      myTime, myThid )
472  #endif /* ALLOW_FIZHI */  #endif /* ALLOW_FIZHI */
473    
474    #ifdef ALLOW_ADDFLUID
475          IF ( selectAddFluid.NE.0 .AND. salt_addMass.NE.UNSET_RL ) THEN
476           IF ( ( selectAddFluid.GE.1 .AND. nonlinFreeSurf.GT.0 )
477         &      .OR. convertFW2Salt.EQ.-1. _d 0 ) THEN
478             DO j=1,sNy
479              DO i=1,sNx
480                gS(i,j,kLev,bi,bj) = gS(i,j,kLev,bi,bj)
481         &        + addMass(i,j,kLev,bi,bj)*mass2rUnit
482         &          *( salt_addMass - salt(i,j,kLev,bi,bj) )
483         &          *recip_rA(i,j,bi,bj)
484         &          *recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj)
485    C    &          *recip_deepFac2C(kLev)*recip_rhoFacC(kLev)
486              ENDDO
487             ENDDO
488           ELSE
489             DO j=1,sNy
490              DO i=1,sNx
491                gS(i,j,kLev,bi,bj) = gS(i,j,kLev,bi,bj)
492         &        + addMass(i,j,kLev,bi,bj)*mass2rUnit
493         &          *( salt_addMass - sRef(kLev) )
494         &          *recip_rA(i,j,bi,bj)
495         &          *recip_drF(kLev)*_recip_hFacC(i,j,kLev,bi,bj)
496    C    &          *recip_deepFac2C(kLev)*recip_rhoFacC(kLev)
497              ENDDO
498             ENDDO
499           ENDIF
500          ENDIF
501    #endif /* ALLOW_ADDFLUID */
502    
503  C     Add fresh-water in top-layer  C     Add fresh-water in top-layer
504        IF ( kLev .EQ. kSurface ) THEN        IF ( kLev .EQ. kSurface ) THEN
505         DO j=1,sNy         DO j=1,sNy
# Line 450  C     Add fresh-water in top-layer Line 511  C     Add fresh-water in top-layer
511         ENDDO         ENDDO
512        ENDIF        ENDIF
513    
 cph#ifndef ALLOW_AUTODIFF_TAMC  
 cph I didnt put this ifndef here.  
514        IF (linFSConserveTr) THEN        IF (linFSConserveTr) THEN
515         DO j=1,sNy         DO j=1,sNy
516          DO i=1,sNx          DO i=1,sNx
# Line 462  cph I didnt put this ifndef here. Line 521  cph I didnt put this ifndef here.
521          ENDDO          ENDDO
522         ENDDO         ENDDO
523        ENDIF        ENDIF
 cph#endif /* ndfef ALLOW_AUTODIFF_TAMC */  
524    
525  #ifdef ALLOW_SHELFICE  #ifdef ALLOW_SHELFICE
526        IF ( useShelfIce )        IF ( useShelfIce )
# Line 471  cph#endif /* ndfef ALLOW_AUTODIFF_TAMC * Line 529  cph#endif /* ndfef ALLOW_AUTODIFF_TAMC *
529       I     myTime, myThid )       I     myTime, myThid )
530  #endif /* ALLOW_SHELFICE */  #endif /* ALLOW_SHELFICE */
531    
532    #ifdef ALLOW_ICEFRONT
533          IF ( useICEFRONT )
534         &     CALL ICEFRONT_TENDENCY_APPLY_S(
535         &     bi,bj, kLev, myTime, myThid )
536    #endif /* ALLOW_ICEFRONT */
537    
538  #ifdef ALLOW_SALT_PLUME  #ifdef ALLOW_SALT_PLUME
539        IF ( useSALT_PLUME )        IF ( useSALT_PLUME )
540       &     CALL SALT_PLUME_TENDENCY_APPLY_S(       &     CALL SALT_PLUME_TENDENCY_APPLY_S(
Legend:
Removed from v.1.53  
changed lines
  Added in v.1.57
  ViewVC Help
Powered by ViewVC 1.1.22