/[MITgcm]/MITgcm/pkg/opps/opps_calc.F
ViewVC logotype

Diff of /MITgcm/pkg/opps/opps_calc.F

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

revision 1.7 by jmc, Tue Mar 16 00:16:50 2010 UTC revision 1.8 by jmc, Tue Aug 24 14:58:48 2010 UTC
# Line 139  CML     &      Dd(k),Pd(k,1),Pd(k,2) Line 139  CML     &      Dd(k),Pd(k,1),Pd(k,2)
139  CML compute vertical velocity at cell centers from GCM velocity  CML compute vertical velocity at cell centers from GCM velocity
140         Wd(k)= - .5*(wVel(K)+wVel(K+1))         Wd(k)= - .5*(wVel(K)+wVel(K+1))
141  CML(  CML(
142  CML     avoid division by zero  CML    avoid division by zero
143  CML       IF (Wd(K) .EQ. 0.D0) Wd(K) = 2.23e-16  CML       IF (Wd(K) .EQ. 0.D0) Wd(K) = 2.23e-16
144  CML)  CML)
145  c  c
# Line 351  c Line 351  c
351           DO ktr = 1, nTracerInUse           DO ktr = 1, nTracerInUse
352            Pda(k2,ktr) = Pd(k2,ktr)            Pda(k2,ktr) = Pd(k2,ktr)
353            Paa(k2,ktr) = tracerEnv(k2+1,ktr)            Paa(k2,ktr) = tracerEnv(k2+1,ktr)
354           ENDDO             ENDDO
355  c  c
356          enddo          enddo
357          dt = min(dt,dtts)          dt = min(dt,dtts)
# Line 365  c Line 365  c
365  #endif /* ALLOW_OPPS_DEBUG */  #endif /* ALLOW_OPPS_DEBUG */
366          DO ktr=1,nTracerInUse          DO ktr=1,nTracerInUse
367           Pda(maxdepth,ktr) = Pd(maxdepth,ktr)           Pda(maxdepth,ktr) = Pd(maxdepth,ktr)
368          ENDDO            ENDDO
369  C  C
370          kmx = maxdepth-1          kmx = maxdepth-1
371          do nn=1,ntime          do nn=1,ntime
# Line 375  C Line 375  C
375           DO ktr = 1,nTracerInUse           DO ktr = 1,nTracerInUse
376            tracerEnv(k,ktr) =  tracerEnv(k,ktr)-            tracerEnv(k,ktr) =  tracerEnv(k,ktr)-
377       &        (mda(k)*(Pda(k,ktr)-Paa(k,ktr)))*dt*recip_drF(k)       &        (mda(k)*(Pda(k,ktr)-Paa(k,ktr)))*dt*recip_drF(k)
378           ENDDO             ENDDO
379  c  c
380  c now do inner points if there are any  c now do inner points if there are any
381  c  c
# Line 411  c Line 411  c
411            DO k2=1,kmx            DO k2=1,kmx
412             paa(k2,ktr) = tracerEnv(k2+1,ktr)             paa(k2,ktr) = tracerEnv(k2+1,ktr)
413            ENDDO            ENDDO
414           ENDDO             ENDDO
415  c  c
416  c end loop on number of time integration steps  c end loop on number of time integration steps
417  c  c
Legend:
Removed from v.1.7  
changed lines
  Added in v.1.8
  ViewVC Help
Powered by ViewVC 1.1.22