New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 2576 for trunk/NEMOGCM/NEMO/OPA_SRC – NEMO

Ignore:
Timestamp:
2011-02-02T18:50:35+01:00 (13 years ago)
Author:
djlea
Message:

Bug fix to velocity observation operator

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/OBS/obs_oper.F90

    r2287 r2576  
    10271027      REAL(KIND=wp) :: zdaystp 
    10281028      REAL(KIND=wp), DIMENSION(kpk) :: & 
    1029          & zobsmask, & 
    1030          & zobsk,    & 
     1029         & zobsmasku, & 
     1030         & zobsmaskv, & 
     1031         & zobsmask,  & 
     1032         & zobsk,     & 
    10311033         & zobs2k 
    10321034      REAL(KIND=wp), DIMENSION(2,2,kpk) :: & 
     
    11601162         zphi = prodatqc%rphi(jobs) 
    11611163 
     1164         ! Initialize observation masks 
     1165 
     1166         zobsmasku(:) = 0.0 
     1167         zobsmaskv(:) = 0.0 
     1168          
    11621169         ! Horizontal weights and vertical mask 
    11631170 
     
    11661173            CALL obs_int_h2d_init( kpk, kpk, k2dint, zlam, zphi,     & 
    11671174               &                   zglamu(:,:,iobs), zgphiu(:,:,iobs), & 
    1168                &                   zumask(:,:,:,iobs), zweigu, zobsmask ) 
     1175               &                   zumask(:,:,:,iobs), zweigu, zobsmasku ) 
    11691176 
    11701177         ENDIF 
     
    11751182            CALL obs_int_h2d_init( kpk, kpk, k2dint, zlam, zphi,     & 
    11761183               &                   zglamv(:,:,iobs), zgphiv(:,:,iobs), & 
    1177                &                   zvmask(:,:,:,iobs), zweigv, zobsmask ) 
     1184               &                   zvmask(:,:,:,iobs), zweigv, zobsmasku ) 
    11781185 
    11791186         ENDIF 
    11801187 
     1188         ! Ensure that the vertical mask on u and v are consistent. 
     1189 
     1190         zobsmask(:) = MIN( zobsmasku(:), zobsmaskv(:) ) 
    11811191 
    11821192         IF ( prodatqc%npvend(jobs,1) > 0 ) THEN 
Note: See TracChangeset for help on using the changeset viewer.