Changeset 15365 for branches/UKMO
- Timestamp:
- 2021-10-14T09:32:07+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/UKMO/dev_r5518_obs_oper_update_DepthBug/NEMOGCM/NEMO/OPA_SRC/OBS/obsinter_z1d.h90
r15356 r15365 63 63 z1dp = ( pobsdep(jdep) - pdep(kkco(jdep)-1) ) 64 64 65 ! Where ob is at or above top model T-point, use model top rather than extrapolate 66 IF ( pobsdep(jdep) <= pdep(kkco(jdep)-1) ) z1dp = 0.0_wp 65 ! Where both levels are masked, return a fill value 66 IF ( pobsmask(kkco(jdep)-1) == 0.0_wp ) .and. (pobsmask(kkco(jdep)) == 0.0_wp) THEN 67 pobs(jdep) = 99999. 68 ELSE 67 69 68 ! Where ob is at or below bottom model T-point, use model bottom rather than extrapolate 69 IF ( pobsdep(jdep) >= pdep(kkco(jdep)) ) z1dm = 0.0_wp 70 71 ! Where deeper level is masked, only use higher level 72 IF ( pobsmask(kkco(jdep)) == 0.0_wp ) z1dp = 0.0_wp 73 74 zsum = z1dm + z1dp 70 ! Where upper level is masked (e.g., under ice cavity), only use deeper level 71 IF ( pobsmask(kkco(jdep)-1) == 0.0_wp ) z1dm = 0.0_wp 72 73 ! Where deeper level is masked (e.g., near sea bed), only use upper level 74 IF ( pobsmask(kkco(jdep)) == 0.0_wp ) z1dp = 0.0_wp 75 76 ! Where ob is at or above upper level model T-point, use upper model level 77 ! rather than extrapolate, except where that level is masked 78 IF ( pobsdep(jdep) <= pdep(kkco(jdep)-1) ) .and. & 79 & ( pobsmask(kkco(jdep)-1) == 0.0_wp ) ) z1dp = 0.0_wp 80 81 ! Where ob is at or below deeper level model T-point, use deeper model level 82 ! rather than extrapolate, except where that level is masked 83 IF ( pobsdep(jdep) >= pdep(kkco(jdep)) ) .and. & 84 & ( pobsmask(kkco(jdep)) == 0.0_wp ) ) z1dm = 0.0_wp 85 86 87 zsum = z1dm + z1dp 75 88 76 IF ( zsum > 0.0_wp ) THEN77 89 IF ( k1dint == 0 ) THEN 78 90 … … 97 109 98 110 ENDIF 99 ELSE ! take value directly from the higher model level100 pobs(jdep) = pobsk(kkco(jdep)-1)101 ENDIF102 111 END DO 103 112
Note: See TracChangeset
for help on using the changeset viewer.