Ignore:
Timestamp:
2018-04-27T19:39:31+02:00 (3 years ago)
Author:
mathiot
Message:

fix to OBS to deal with a change in the halo at the wrap points (Dan Lea)

Location:
branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/OBS
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/OBS/mpp_map.F90

    r4245 r9523  
    1212   USE par_kind, ONLY :   wp            ! Precision variables 
    1313   USE par_oce , ONLY :   jpi, jpj      ! Ocean parameters 
    14    USE dom_oce , ONLY :   mig, mjg, nldi, nlei, nldj, nlej, narea   ! Ocean space and time domain variables 
     14   USE dom_oce , ONLY :   mig, mjg, nldi, nlei, nldj, nlej, nlci, nlcj, narea   ! Ocean space and time domain variables 
    1515#if defined key_mpp_mpi 
    1616   USE lib_mpp, ONLY :   mpi_comm_opa   ! MPP library 
     
    6464      imppmap(:,:) = 0 
    6565 
    66       ! Setup local grid points 
    67       imppmap(mig(nldi):mig(nlei),mjg(nldj):mjg(nlej)) = narea  
     66!      ! Setup local grid points 
     67      imppmap(mig(1):mig(nlci),mjg(1):mjg(nlcj)) = narea 
    6868       
    6969      ! Get global data 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/OBS/obs_grid.F90

    r6140 r9523  
    129129            IF ( cdgrid == 'T' ) THEN 
    130130               CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 
    131                   &                             nldi, nlei,nldj,  nlej,   & 
     131                  &                             1, nlci, 1, nlcj,         & 
    132132                  &                             nproc, jpnij,             & 
    133133                  &                             glamt, gphit, tmask,      & 
     
    136136            ELSEIF ( cdgrid == 'U' ) THEN 
    137137               CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 
    138                   &                             nldi, nlei,nldj,  nlej,   & 
     138                  &                             1, nlci, 1, nlcj,         & 
    139139                  &                             nproc, jpnij,             & 
    140140                  &                             glamu, gphiu, umask,      & 
     
    143143            ELSEIF ( cdgrid == 'V' ) THEN 
    144144               CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 
    145                   &                             nldi, nlei,nldj,  nlej,   & 
     145                  &                             1, nlci, 1, nlcj,         & 
    146146                  &                             nproc, jpnij,             & 
    147147                  &                             glamv, gphiv, vmask,      & 
     
    150150            ELSEIF ( cdgrid == 'F' ) THEN 
    151151               CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo, & 
    152                   &                             nldi, nlei,nldj,  nlej,   & 
     152                  &                             1, nlci, 1, nlcj,         & 
    153153                  &                             nproc, jpnij,             & 
    154154                  &                             glamf, gphif, fmask,      & 
     
    279279         zmskg(:,:) = -1.e+10 
    280280         ! Add various grids here. 
    281          DO jj = nldj, nlej 
    282             DO ji = nldi, nlei 
     281         DO jj = 1, nlcj 
     282            DO ji = 1, nlci 
    283283               zlamg(mig(ji),mjg(jj)) = glamt(ji,jj) 
    284284               zphig(mig(ji),mjg(jj)) = gphit(ji,jj) 
     
    816816             
    817817            CALL obs_grd_bruteforce( jpi, jpj, jpiglo, jpjglo,  & 
    818                &                     nldi, nlei,nldj,  nlej,    & 
     818               &                     1, nlci, 1, nlcj,          & 
    819819               &                     nproc, jpnij,              & 
    820820               &                     glamt, gphit, tmask,       & 
  • branches/2017/dev_merge_2017/NEMOGCM/NEMO/OPA_SRC/OBS/obs_inter_sup.F90

    r9125 r9523  
    365365         CALL ctl_stop( 'Error in obs_int_comm_3d_local', & 
    366366            &           'Point outside local domain' ) 
    367           
     367  
    368368      ENDIF 
    369369 
Note: See TracChangeset for help on using the changeset viewer.