/[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.4 by cnh, Thu Aug 7 13:03:31 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 86  C     tau = 1 day, 5days, 10days, 20days Line 86  C     tau = 1 day, 5days, 10days, 20days
86        recip_tauSp(5)   = 1./(86400.*60.)        recip_tauSp(5)   = 1./(86400.*60.)
87        jSouthBndy = 5        jSouthBndy = 5
88        jNorthBndy = ny-5+1        jNorthBndy = ny-5+1
89        DO j=jMin,jMax        DO j=1,sNy
90         DO i=iMin,iMax         DO i=iMin,iMax
91          jG = myYGlobalLo+(bj-1)*sNy+j-1          jG = myYGlobalLo+(bj-1)*sNy+j-1
92          jFromNBndy = jNorthBndy-jG          jFromNBndy = jNorthBndy-jG
# 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 204  C     tau = 1 day, 5days, 10days, 20days Line 204  C     tau = 1 day, 5days, 10days, 20days
204        recip_tauSp(5)   = 1./(86400.*60.)        recip_tauSp(5)   = 1./(86400.*60.)
205        jSouthBndy = 5        jSouthBndy = 5
206        jNorthBndy = ny-5+1        jNorthBndy = ny-5+1
207        DO j=jMin,jMax        DO j=1,sNy
208         DO i=iMin,iMax         DO i=iMin,iMax
209          jG = myYGlobalLo+(bj-1)*sNy+j-1          jG = myYGlobalLo+(bj-1)*sNy+j-1
210          jFromNBndy = jNorthBndy-jG          jFromNBndy = jNorthBndy-jG
# 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=1./(86400.*90.)
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 335  C     tau = 1 day, 5days, 10days, 20days Line 349  C     tau = 1 day, 5days, 10days, 20days
349        recip_tauSp(5)   = 1./(86400.*60.)        recip_tauSp(5)   = 1./(86400.*60.)
350        jSouthBndy = 5        jSouthBndy = 5
351        jNorthBndy = ny-5+1        jNorthBndy = ny-5+1
352        DO j=jMin,jMax        DO j=1,sNy
353         DO i=iMin,iMax         DO i=iMin,iMax
354          jG = myYGlobalLo+(bj-1)*sNy+j-1          jG = myYGlobalLo+(bj-1)*sNy+j-1
355          jFromNBndy = jNorthBndy-jG          jFromNBndy = jNorthBndy-jG
# Line 427  C     Add fresh-water in top-layer Line 441  C     Add fresh-water in top-layer
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=1./(86400.*90.)
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
460         CALL OBCS_SPONGE_S(         CALL OBCS_SPONGE_S(
# Line 446  C     tau = 1 day, 5days, 10days, 20days Line 474  C     tau = 1 day, 5days, 10days, 20days
474        recip_tauSp(5)   = 1./(86400.*60.)        recip_tauSp(5)   = 1./(86400.*60.)
475        jSouthBndy = 5        jSouthBndy = 5
476        jNorthBndy = ny-5+1        jNorthBndy = ny-5+1
477        DO j=jMin,jMax        DO j=1,sNy
478         DO i=iMin,iMax         DO i=iMin,iMax
479          jG = myYGlobalLo+(bj-1)*sNy+j-1          jG = myYGlobalLo+(bj-1)*sNy+j-1
480          jFromNBndy = jNorthBndy-jG          jFromNBndy = jNorthBndy-jG
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22