Changeset 11639 for branches/UKMO/AMM15_v3_6_STABLE_package_reanalysis4/NEMOGCM/NEMO/OPA_SRC/OBS/obsinter_z1d.h90
- Timestamp:
- 2019-10-02T13:13:03+02:00 (4 years ago)
- 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 26 26 !! ! 06-10 (A. Weaver) Cleanup 27 27 !! ! 07-01 (K. Mogensen) Use profile rather than single level 28 !! ! 19-07 (R. Renshaw) Avoid division by zero 28 29 !!--------------------------------------------------------------------- 29 30 … … 62 63 z1dm = ( pdep(kkco(jdep)) - pobsdep(jdep) ) 63 64 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 64 68 IF ( pobsmask(kkco(jdep)) == 0.0_wp ) z1dp = 0.0_wp 65 69 66 70 zsum = z1dm + z1dp 71 ! if pobsmask==0 and model level depth==observed depth, we get zsum=0 72 IF ( zsum > 0.0_wp ) THEN 67 73 68 74 IF ( k1dint == 0 ) THEN … … 88 94 89 95 ENDIF 96 97 ELSE ! take value directly from the higher model level 98 pobs(jdep) = pobsk(kkco(jdep)-1) 99 ENDIF 100 90 101 END DO 91 102
Note: See TracChangeset
for help on using the changeset viewer.