/[MITgcm]/MITgcm/pkg/ecco/cost_gencost_sshv4.F
ViewVC logotype

Diff of /MITgcm/pkg/ecco/cost_gencost_sshv4.F

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

revision 1.2 by gforget, Mon Oct 25 20:42:31 2010 UTC revision 1.3 by gforget, Wed May 18 03:15:38 2011 UTC
# Line 359  cgf ==================================== Line 359  cgf ====================================
359       &                       ladinit, optimcycle, mythid,       &                       ladinit, optimcycle, mythid,
360       &                       xx_psbar_mean_dummy )       &                       xx_psbar_mean_dummy )
361    
362    #ifndef ALLOW_PSBAR_MEAN
363            CALL REMOVE_MEAN_RL( 1, psbar, maskInC, maskInC, rA, drF,
364         &        'psbar', myTime, myThid )
365    #endif
366    
367  #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION  #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION
368        call cost_sla_read( topexfile, topexstartdate, topexperiod,        call cost_sla_read( topexfile, topexstartdate, topexperiod,
369       &                topexintercept, topexslope,       &                topexintercept, topexslope,
# Line 433  c--   END loop over records for the firs Line 438  c--   END loop over records for the firs
438              do j = jmin,jmax              do j = jmin,jmax
439                do i = imin,imax                do i = imin,imax
440  #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION  #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION
441                 if ( mean_psMtpobs_NUM(i,j,bi,bj) .NE. 0. ) then                 if ( ( mean_psMtpobs_NUM(i,j,bi,bj) .NE. 0. ).AND.
442         &              ( abs(YC(i,j,bi,bj)) .LE. 66. ) ) then
443                    mean_psMssh_all(i,j,bi,bj) =                    mean_psMssh_all(i,j,bi,bj) =
444       &                 mean_psMssh_all(i,j,bi,bj) +       &                 mean_psMssh_all(i,j,bi,bj) +
445       &                 mean_psMtpobs(i,j,bi,bj)       &                 mean_psMtpobs(i,j,bi,bj)
# Line 447  c--   END loop over records for the firs Line 453  c--   END loop over records for the firs
453                 endif                 endif
454  #endif  #endif
455  #ifdef ALLOW_SSH_ERSANOM_COST_CONTRIBUTION  #ifdef ALLOW_SSH_ERSANOM_COST_CONTRIBUTION
456                 if ( mean_psMersobs_NUM(i,j,bi,bj) .NE. 0. ) then                 if ( ( mean_psMersobs_NUM(i,j,bi,bj) .NE. 0. ).AND.
457         &              ( abs(YC(i,j,bi,bj)) .LE. 66. ) ) then
458                    mean_psMssh_all(i,j,bi,bj) =                    mean_psMssh_all(i,j,bi,bj) =
459       &                 mean_psMssh_all(i,j,bi,bj) +       &                 mean_psMssh_all(i,j,bi,bj) +
460       &                 mean_psMersobs(i,j,bi,bj)       &                 mean_psMersobs(i,j,bi,bj)
# Line 461  c--   END loop over records for the firs Line 468  c--   END loop over records for the firs
468                 endif                 endif
469  #endif  #endif
470  #ifdef ALLOW_SSH_GFOANOM_COST_CONTRIBUTION  #ifdef ALLOW_SSH_GFOANOM_COST_CONTRIBUTION
471                 if ( mean_psMgfoobs_NUM(i,j,bi,bj) .NE. 0. ) then                 if ( ( mean_psMgfoobs_NUM(i,j,bi,bj) .NE. 0. ).AND.
472         &              ( abs(YC(i,j,bi,bj)) .LE. 66. ) ) then
473                    mean_psMssh_all(i,j,bi,bj) =                    mean_psMssh_all(i,j,bi,bj) =
474       &                 mean_psMssh_all(i,j,bi,bj) +       &                 mean_psMssh_all(i,j,bi,bj) +
475       &                 mean_psMgfoobs(i,j,bi,bj)       &                 mean_psMgfoobs(i,j,bi,bj)
# Line 663  c -------------------------------------- Line 671  c --------------------------------------
671       &                       ladinit, optimcycle, mythid,       &                       ladinit, optimcycle, mythid,
672       &                       xx_psbar_mean_dummy )       &                       xx_psbar_mean_dummy )
673    
674    #ifndef ALLOW_PSBAR_MEAN
675            CALL REMOVE_MEAN_RL( 1, psbar, maskInC, maskInC, rA, drF,
676         &        'psbar', myTime, myThid )
677    #endif
678    
679  #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION  #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION
680        call cost_sla_read( topexfile, topexstartdate, topexperiod,        call cost_sla_read( topexfile, topexstartdate, topexperiod,
681       &                topexintercept, topexslope,       &                topexintercept, topexslope,
# Line 770  c -------------------------------------- Line 783  c --------------------------------------
783                 endif                 endif
784  #endif  #endif
785                 if ( ( anom_psMslaobs_NUM(i,j,bi,bj) .NE. 0. ).AND.                 if ( ( anom_psMslaobs_NUM(i,j,bi,bj) .NE. 0. ).AND.
786       &              ( maskc(i,j,1,bi,bj) .NE. 0. ) ) then       &              ( maskc(i,j,1,bi,bj) .NE. 0. ) .AND.
787         &              ( abs(YC(i,j,bi,bj)) .LE. 66. ) ) then
788                    anom_psMslaobs(i,j,bi,bj) =                    anom_psMslaobs(i,j,bi,bj) =
789       &                 anom_psMslaobs(i,j,bi,bj) /       &                 anom_psMslaobs(i,j,bi,bj) /
790       &                 anom_psMslaobs_NUM(i,j,bi,bj)       &                 anom_psMslaobs_NUM(i,j,bi,bj)
# Line 854  cgf ==================================== Line 868  cgf ====================================
868       &                       ladinit, optimcycle, mythid,       &                       ladinit, optimcycle, mythid,
869       &                       xx_psbar_mean_dummy )       &                       xx_psbar_mean_dummy )
870    
871    #ifndef ALLOW_PSBAR_MEAN
872            CALL REMOVE_MEAN_RL( 1, psbar, maskInC, maskInC, rA, drF,
873         &        'psbar', myTime, myThid )
874    #endif
875    
876  #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION  #ifdef ALLOW_SSH_TPANOM_COST_CONTRIBUTION
877          call cost_readtopex( irec, mythid )          call cost_readtopex( irec, mythid )
878  #endif  #endif
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22