/[MITgcm]/MITgcm/verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F
ViewVC logotype

Diff of /MITgcm/verification/tutorial_global_oce_latlon/code/ptracers_forcing_surf.F

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

revision 1.2 by dfer, Mon May 14 19:48:36 2007 UTC revision 1.3 by jmc, Mon Oct 15 15:29:31 2007 UTC
# Line 60  C Example of how to add forcing at the s Line 60  C Example of how to add forcing at the s
60             DO i = iMin, iMax             DO i = iMin, iMax
61  C     Restoring to zero at the surface on a 10-day timescale  C     Restoring to zero at the surface on a 10-day timescale
62               surfaceForcingPtr(i,j,bi,bj,iTrc) =               surfaceForcingPtr(i,j,bi,bj,iTrc) =
63       &        + 1. _d 0 / (10. _d 0 * 86400. _d 0)       &        + 1. _d 0 / (10. _d 0 * 86400. _d 0)
64       &                  * ( 0. _d 0 - pTracer(i,j,ks,bi,bj,iTrc) )       &                  * ( 0. _d 0 - pTracer(i,j,ks,bi,bj,iTrc) )
65       &                  * drF(ks) * _hFacC(i,j,ks,bi,bj)       &                  * drF(ks) * _hFacC(i,j,ks,bi,bj)
66             ENDDO             ENDDO
67            ENDDO            ENDDO
# Line 84  c Line 84  c
84       &          surfaceForcingPtr(i,j,bi,bj,iTrc)       &          surfaceForcingPtr(i,j,bi,bj,iTrc)
85       &        + PmEpR(i,j,bi,bj)       &        + PmEpR(i,j,bi,bj)
86       &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )       &          *( PTRACERS_EvPrRn(iTrc) - pTracer(i,j,ks,bi,bj,iTrc) )
87       &          *convertEmP2rUnit       &          *mass2rUnit
88             ENDDO             ENDDO
89            ENDDO            ENDDO
90           ENDIF           ENDIF
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22