Changeset 6285


Ignore:
Timestamp:
2016-02-04T07:23:48+01:00 (6 years ago)
Author:
kingr
Message:

Updated obs_prep.F90 to deal with cases where bottom level of obs could be masked incorrectly due to the varying depths between template profiles.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/dev_r4650_general_vert_coord_obsoper/NEMOGCM/NEMO/OPA_SRC/OBS/obs_prep.F90

    r6273 r6285  
    2626      & obs_mpp_sum_integers 
    2727   USE obs_inter_sup      ! Interpolation support 
    28    USE obs_oper           ! Observation operators 
    2928   USE lib_mpp, ONLY : & 
    3029      & ctl_warn, ctl_stop 
     30   USE dom_oce, ONLY : &   ! Geographical information 
     31      & glamt, glamu, glamv,    & 
     32      & gphit, gphiu, gphiv,    & 
     33      & gdept_1d,             & 
     34      & tmask, umask, vmask,  & 
     35      & nproc, rdt,    & 
     36      & gdepw_1d,      & 
     37      & gdepw_0,       & 
     38#if defined key_vvl  
     39      & gdepw_n,       & 
     40      & gdept_n,       & 
     41#endif  
     42      & ln_zco,        & 
     43      & ln_zps,        & 
     44      & lk_vvl 
    3145 
    3246   IMPLICIT NONE 
     
    7589      USE domstp              ! Domain: set the time-step 
    7690      USE par_oce             ! Ocean parameters 
    77       USE dom_oce, ONLY : &   ! Geographical information 
    78          & glamt,   & 
    79          & gphit,   & 
    80          & gdept_1d,& 
    81          & tmask,   & 
    82          & nproc 
    8391      !! * Arguments 
    8492      TYPE(obs_prof), INTENT(INOUT) :: profdata     ! Full set of profile data 
     
    359367      USE domstp              ! Domain: set the time-step 
    360368      USE par_oce             ! Ocean parameters 
    361       USE dom_oce, ONLY : &   ! Geographical information 
    362          & glamt,   & 
    363          & gphit,   & 
    364          & tmask,   & 
    365          & nproc 
    366369      !! * Arguments 
    367370      TYPE(obs_surf), INTENT(INOUT) :: sladata    ! Full set of SLA data 
     
    547550      USE domstp              ! Domain: set the time-step 
    548551      USE par_oce             ! Ocean parameters 
    549       USE dom_oce, ONLY : &   ! Geographical information 
    550          & glamt,   & 
    551          & gphit,   & 
    552          & tmask,   & 
    553          & nproc 
    554552      !! * Arguments 
    555553      TYPE(obs_surf), INTENT(INOUT) :: sstdata     ! Full set of SST data 
     
    732730      USE domstp              ! Domain: set the time-step 
    733731      USE par_oce             ! Ocean parameters 
    734       USE dom_oce, ONLY : &   ! Geographical information 
    735          & glamt,   & 
    736          & gphit,   & 
    737          & tmask,   & 
    738          & nproc 
    739732      !! * Arguments 
    740733      TYPE(obs_surf), INTENT(INOUT) :: seaicedata     ! Full set of Sea Ice data 
     
    916909      USE domstp              ! Domain: set the time-step 
    917910      USE par_oce             ! Ocean parameters 
    918       USE dom_oce, ONLY : &   ! Geographical information 
    919          & glamt, glamu, glamv,    & 
    920          & gphit, gphiu, gphiv,    & 
    921          & gdept_1d,             & 
    922          & tmask, umask, vmask,  & 
    923          & nproc 
    924911      !! * Arguments 
    925912      TYPE(obs_prof), INTENT(INOUT) :: profdata   ! Full set of profile data 
     
    12171204      !!---------------------------------------------------------------------- 
    12181205      !! * Modules used 
    1219       USE dom_oce, ONLY : &  ! Geographical information 
    1220          & rdt 
    12211206      USE phycst, ONLY : &   ! Physical constants 
    12221207         & rday,  &              
     
    17111696      !!---------------------------------------------------------------------- 
    17121697      !! * Modules used 
    1713       USE dom_oce, ONLY : &       ! Geographical information 
    1714          & gdepw_1d,      & 
    1715          & gdepw_0,       &                         
    1716 #if defined key_vvl  
    1717          & gdepw_n,       &  
    1718          & gdept_n,       &  
    1719 #endif  
    1720          & ln_zco,        &  
    1721          & ln_zps,        & 
    1722          & lk_vvl                         
    17231698  
    17241699      !! * Arguments 
     
    18491824               &  .OR. ( pobsphi(jobs) <  -90.         )       & 
    18501825               &  .OR. ( pobsphi(jobs) >   90.         )       & 
    1851                &  .OR. ( pobsdep(jobsp) < 0.0          )       & 
    1852                &  .OR. ( pobsdep(jobsp) > gdepw_1d(kpk)) ) THEN 
     1826               &  .OR. ( pobsdep(jobsp) < 0.0          )       ) THEN 
    18531827               kobsqc(jobsp) = kobsqc(jobsp) + 11 
    18541828               kosdobs = kosdobs + 1 
     
    18861860               IF ( ANY( zgdepw(1:2,1:2,kpk,jobs) < pobsdep(jobsp) ) & 
    18871861                  &     .OR. &  
    1888                   &  ( SUM( zgmsk(1:2,1:2,kobsk(jobsp)-1:kobsk(jobsp),jobs) ) &  
    1889                   &  == 0.0_wp) ) THEN  
     1862                  &  ( SUM( zgmsk(1,1,kobsk(jobsp)-1:kobsk(jobsp),jobs)) == 0.0_wp ) .OR. & 
     1863                  &  ( SUM( zgmsk(1,2,kobsk(jobsp)-1:kobsk(jobsp),jobs)) == 0.0_wp ) .OR. & 
     1864                  &  ( SUM( zgmsk(2,1,kobsk(jobsp)-1:kobsk(jobsp),jobs)) == 0.0_wp ) .OR. & 
     1865                  &  ( SUM( zgmsk(2,2,kobsk(jobsp)-1:kobsk(jobsp),jobs)) == 0.0_wp ) ) THEN  
    18901866                  kobsqc(jobsp) = kobsqc(jobsp) + 12  
    18911867                  klanobs = klanobs + 1  
Note: See TracChangeset for help on using the changeset viewer.