/[MITgcm]/MITgcm_contrib/natl_12/code/external_forcing.F
ViewVC logotype

Diff of /MITgcm_contrib/natl_12/code/external_forcing.F

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

revision 1.1 by cnh, Tue Aug 5 21:22:43 2003 UTC revision 1.2 by cnh, Tue Aug 5 21:35:30 2003 UTC
# Line 69  C     Add windstress momentum impulse in Line 69  C     Add windstress momentum impulse in
69         DO j=jMin,jMax         DO j=jMin,jMax
70          DO i=iMin,iMax          DO i=iMin,iMax
71           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)           gU(i,j,kLev,bi,bj) = gU(i,j,kLev,bi,bj)
72       &   +foFacMom*surfaceTendencyU(i,j,bi,bj)*0.0100D0       &   +foFacMom*surfaceTendencyU(i,j,bi,bj)
73       &   *_maskW(i,j,kLev,bi,bj)       &   *_maskW(i,j,kLev,bi,bj)
74          ENDDO          ENDDO
75         ENDDO         ENDDO
# Line 179  C     Add windstress momentum impulse in Line 179  C     Add windstress momentum impulse in
179         DO j=jMin,jMax         DO j=jMin,jMax
180          DO i=iMin,iMax          DO i=iMin,iMax
181           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)           gV(i,j,kLev,bi,bj) = gV(i,j,kLev,bi,bj)
182       &   +foFacMom*surfaceTendencyV(i,j,bi,bj)*0.0100D0       &   +foFacMom*surfaceTendencyV(i,j,bi,bj)
183       &   *_maskS(i,j,kLev,bi,bj)       &   *_maskS(i,j,kLev,bi,bj)
184          ENDDO          ENDDO
185         ENDDO         ENDDO
# Line 299  C     Add heat in top-layer Line 299  C     Add heat in top-layer
299         ENDDO         ENDDO
300        ENDIF        ENDIF
301    
302    C--   Forcing term
303    C     Add heat in top-layer ( 90 day climatalogical average relaxation )
304          IF ( kLev .EQ. kSurface ) THEN
305           curRecipTau=90./86400.
306           DO j=jMin,jMax
307            DO i=iMin,iMax
308             gT(i,j,kLev,bi,bj)=gT(i,j,kLev,bi,bj)
309         &     +maskC(i,j,kLev,bi,bj)*(
310         &      -curRecipTau*(theta(i,j,Klev,bi,bj)-thetaRef(i,j,kLev,bi,bj))
311         &      )
312            ENDDO
313           ENDDO
314          ENDIF
315    
316  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
317  C Penetrating SW radiation  C Penetrating SW radiation
318        swfracb(1)=abs(rF(klev))        swfracb(1)=abs(rF(klev))
# Line 426  C     Add fresh-water in top-layer Line 440  C     Add fresh-water in top-layer
440          ENDDO          ENDDO
441         ENDDO         ENDDO
442        ENDIF        ENDIF
443    
444    C--   Forcing term
445    C     Add freshening/salt in top-layer ( 90 day climatalogical average relaxation )
446          IF ( kLev .EQ. kSurface ) THEN
447           curRecipTau=90./86400.
448           DO j=jMin,jMax
449            DO i=iMin,iMax
450             gS(i,j,kLev,bi,bj)=gS(i,j,kLev,bi,bj)
451         &     +maskC(i,j,kLev,bi,bj)*(
452         &      -curRecipTau*(salt(i,j,Klev,bi,bj)-saltRef(i,j,kLev,bi,bj))
453         &      )
454            ENDDO
455           ENDDO
456          ENDIF
457    
458  #if (defined (ALLOW_OBCS) && defined (ALLOW_OBCS_SPONGE))  #if (defined (ALLOW_OBCS) && defined (ALLOW_OBCS_SPONGE))
459        IF (useOBCS) THEN        IF (useOBCS) THEN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22