/[MITgcm]/MITgcm/pkg/ptracers/ptracers_forcing_surf.F
ViewVC logotype

Diff of /MITgcm/pkg/ptracers/ptracers_forcing_surf.F

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

revision 1.10 by jmc, Sun Sep 5 22:28:14 2010 UTC revision 1.11 by jmc, Thu Mar 8 17:13:31 2012 UTC
# Line 28  C !USES: =============================== Line 28  C !USES: ===============================
28  #include "FFIELDS.h"  #include "FFIELDS.h"
29  #include "PTRACERS_SIZE.h"  #include "PTRACERS_SIZE.h"
30  #include "PTRACERS_PARAMS.h"  #include "PTRACERS_PARAMS.h"
31    #include "PTRACERS_START.h"
32  #include "PTRACERS_FIELDS.h"  #include "PTRACERS_FIELDS.h"
33    
34  C !INPUT PARAMETERS: ===================================================  C !INPUT PARAMETERS: ===================================================
# Line 62  CEOP Line 63  CEOP
63    
64  C Example of how to add forcing at the surface  C Example of how to add forcing at the surface
65        DO iTrc=1,PTRACERS_numInUse        DO iTrc=1,PTRACERS_numInUse
66    c       IF ( PTRACERS_StepFwd(iTrc) ) THEN
67            DO j = jMin, jMax            DO j = jMin, jMax
68             DO i = iMin, iMax             DO i = iMin, iMax
69               surfaceForcingPTr(i,j,bi,bj,iTrc) =               surfaceForcingPTr(i,j,bi,bj,iTrc) = 0. _d 0
70       &               0. _d 0  c    &                        + surfaceForcingS(i,j,bi,bj)
 c    &               surfaceForcingS(i,j,bi,bj)  
71             ENDDO             ENDDO
72            ENDDO            ENDDO
73    c       ENDIF
74        ENDDO        ENDDO
75    
76  C--   Option to convert Salt-relaxation into additional EmP contribution  C--   Option to convert Salt-relaxation into additional EmP contribution
# Line 122  c-  NonLin_FrSurf and RealFreshWaterFlux Line 124  c-  NonLin_FrSurf and RealFreshWaterFlux
124  c   the water column height ; temp., salt, (tracer) flux associated  c   the water column height ; temp., salt, (tracer) flux associated
125  c   with this input/output of water is added here to the surface tendency.  c   with this input/output of water is added here to the surface tendency.
126  c  c
127           IF (PTRACERS_EvPrRn(iTrc).NE.UNSET_RL) THEN           IF ( PTRACERS_StepFwd(iTrc) .AND.
128         &        PTRACERS_EvPrRn(iTrc).NE.UNSET_RL ) THEN
129            DO j = jMin, jMax            DO j = jMin, jMax
130             DO i = iMin, iMax             DO i = iMin, iMax
131               surfaceForcingPTr(i,j,bi,bj,iTrc) =               surfaceForcingPTr(i,j,bi,bj,iTrc) =
# Line 150  C-    use local surface tracer field to Line 153  C-    use local surface tracer field to
153    
154          DO iTrc=1,PTRACERS_numInUse          DO iTrc=1,PTRACERS_numInUse
155    
156           IF (PTRACERS_EvPrRn(iTrc).NE.UNSET_RL) THEN           IF ( PTRACERS_StepFwd(iTrc) .AND.
157         &        PTRACERS_EvPrRn(iTrc).NE.UNSET_RL ) THEN
158  C        account for Rain/Evap tracer content (PTRACERS_EvPrRn) using  C        account for Rain/Evap tracer content (PTRACERS_EvPrRn) using
159  C        local surface tracer  C        local surface tracer
160            DO j = jMin, jMax            DO j = jMin, jMax
# Line 171  C-    use uniform tracer value to calcul Line 175  C-    use uniform tracer value to calcul
175    
176          DO iTrc=1,PTRACERS_numInUse          DO iTrc=1,PTRACERS_numInUse
177    
178           IF (PTRACERS_EvPrRn(iTrc).NE.UNSET_RL) THEN           IF ( PTRACERS_StepFwd(iTrc) .AND.
179         &        PTRACERS_EvPrRn(iTrc).NE.UNSET_RL ) THEN
180  C     account for Rain/Evap tracer content (PTRACERS_EvPrRn) assuming uniform  C     account for Rain/Evap tracer content (PTRACERS_EvPrRn) assuming uniform
181  C     surface tracer (=PTRACERS_ref)  C     surface tracer (=PTRACERS_ref)
182            DO j = jMin, jMax            DO j = jMin, jMax
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22