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 11639 for branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/OPA_SRC/OBS/obsinter_z1d.h90 – NEMO

Ignore:
Timestamp:
2019-10-02T13:13:03+02:00 (4 years ago)
Author:
rrenshaw
Message:

regional mean and transport diagnostics added for reanalysis

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/OPA_SRC/OBS/obsinter_z1d.h90

    r8058 r11639  
    2626      !!      ! 06-10 (A. Weaver) Cleanup 
    2727      !!      ! 07-01 (K. Mogensen) Use profile rather than single level 
     28      !!      ! 19-07 (R. Renshaw) Avoid division by zero 
    2829      !!--------------------------------------------------------------------- 
    2930 
     
    6263         z1dm = ( pdep(kkco(jdep)) - pobsdep(jdep)      ) 
    6364         z1dp = ( pobsdep(jdep)    - pdep(kkco(jdep)-1) ) 
     65!        Where ob is below model bottom, use model bottom rather than extrapolate 
     66         IF ( pdep(kkco(jdep)) <= pobsdep(jdep) ) z1dm = 0.0_wp 
     67!        Where lower level is missing, use higher level 
    6468         IF ( pobsmask(kkco(jdep)) == 0.0_wp ) z1dp = 0.0_wp 
    6569 
    6670         zsum = z1dm + z1dp 
     71         ! if pobsmask==0 and model level depth==observed depth, we get zsum=0 
     72         IF ( zsum > 0.0_wp ) THEN 
    6773          
    6874         IF ( k1dint == 0 ) THEN 
     
    8894 
    8995         ENDIF 
     96 
     97         ELSE  ! take value directly from the higher model level 
     98            pobs(jdep)  = pobsk(kkco(jdep)-1) 
     99         ENDIF 
     100 
    90101      END DO 
    91102 
Note: See TracChangeset for help on using the changeset viewer.