Changeset 12928 for NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/OBS/obs_readmdt.F90
- Timestamp:
- 2020-05-14T21:46:00+02:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
- Property svn:externals
-
old new 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@HEAD sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/OBS/obs_readmdt.F90
r10425 r12928 25 25 & tmask, tmask_i, e1e2t, gphit, glamt 26 26 USE obs_const, ONLY : obfillflt ! Fillvalue 27 USE oce , ONLY : ssh n! Model variables27 USE oce , ONLY : ssh ! Model variables 28 28 29 29 IMPLICIT NONE … … 37 37 REAL(wp), PUBLIC :: rn_mdtcutoff = 65.0_wp ! MDT cutoff for computed correction 38 38 39 !! * Substitutions 40 # include "do_loop_substitute.h90" 39 41 !!---------------------------------------------------------------------- 40 42 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 44 46 CONTAINS 45 47 46 SUBROUTINE obs_rea_mdt( sladata, k2dint )48 SUBROUTINE obs_rea_mdt( sladata, k2dint, Kmm ) 47 49 !!--------------------------------------------------------------------- 48 50 !! … … 59 61 TYPE(obs_surf), INTENT(inout) :: sladata ! SLA data 60 62 INTEGER , INTENT(in) :: k2dint ! ? 63 INTEGER , INTENT(in) :: Kmm ! ? 61 64 ! 62 65 CHARACTER(LEN=12), PARAMETER :: cpname = 'obs_rea_mdt' … … 106 109 ! Remove the offset between the MDT used with the sla and the model MDT 107 110 IF( nn_msshc == 1 .OR. nn_msshc == 2 ) & 108 & CALL obs_offset_mdt( jpi, jpj, z_mdt, zfill )111 & CALL obs_offset_mdt( jpi, jpj, z_mdt, zfill, Kmm ) 109 112 110 113 ! Intepolate the MDT already on the model grid at the observation point … … 169 172 170 173 171 SUBROUTINE obs_offset_mdt( kpi, kpj, mdt, zfill )174 SUBROUTINE obs_offset_mdt( kpi, kpj, mdt, zfill, Kmm ) 172 175 !!--------------------------------------------------------------------- 173 176 !! … … 183 186 !!---------------------------------------------------------------------- 184 187 INTEGER, INTENT(IN) :: kpi, kpj 188 INTEGER, INTENT(IN) :: Kmm 185 189 REAL(wp), DIMENSION(kpi,kpj), INTENT(INOUT) :: mdt ! MDT used on the model grid 186 190 REAL(wp) , INTENT(IN ) :: zfill … … 211 215 zeta2 = 0.0 212 216 213 DO jj = 1, jpj 214 DO ji = 1, jpi 215 zdxdy = e1e2t(ji,jj) * zpromsk(ji,jj) 216 zarea = zarea + zdxdy 217 zeta1 = zeta1 + mdt(ji,jj) * zdxdy 218 zeta2 = zeta2 + sshn (ji,jj) * zdxdy 219 END DO 220 END DO 217 DO_2D_11_11 218 zdxdy = e1e2t(ji,jj) * zpromsk(ji,jj) 219 zarea = zarea + zdxdy 220 zeta1 = zeta1 + mdt(ji,jj) * zdxdy 221 zeta2 = zeta2 + ssh(ji,jj,Kmm) * zdxdy 222 END_2D 221 223 222 224 CALL mpp_sum( 'obs_readmdt', zeta1 )
Note: See TracChangeset
for help on using the changeset viewer.