Changeset 13189 for NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/ISOMIP+/MY_SRC/dtatsd.F90
- Timestamp:
- 2020-07-01T11:27:25+02:00 (4 years ago)
- Location:
- NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback
- Property svn:externals
-
old new 8 8 9 9 # SETTE 10 ^/utils/CI/sette@ HEADsette10 ^/utils/CI/sette@12931 sette
-
- Property svn:externals
-
NEMO/branches/2020/dev_r12472_ASINTER-05_Masson_CurrentFeedback/tests/ISOMIP+/MY_SRC/dtatsd.F90
r12077 r13189 36 36 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_tsddmp ! structure of input SST (file informations, fields read) 37 37 38 !! * Substitutions 39 # include "do_loop_substitute.h90" 38 40 !!---------------------------------------------------------------------- 39 41 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 67 69 ierr0 = 0 ; ierr1 = 0 ; ierr2 = 0 ; ierr3 = 0 68 70 ! 69 REWIND( numnam_ref ) ! Namelist namtsd in reference namelist :70 71 READ ( numnam_ref, namtsd, IOSTAT = ios, ERR = 901) 71 72 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namtsd in reference namelist' ) 72 REWIND( numnam_cfg ) ! Namelist namtsd in configuration namelist : Parameters of the run73 73 READ ( numnam_cfg, namtsd, IOSTAT = ios, ERR = 902 ) 74 74 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namtsd in configuration namelist' ) … … 191 191 ENDIF 192 192 ! 193 DO jj = 1, jpj ! vertical interpolation of T & S 194 DO ji = 1, jpi 195 DO jk = 1, jpk ! determines the intepolated T-S profiles at each (i,j) points 196 zl = gdept_0(ji,jj,jk) 197 IF( zl < gdept_1d(1 ) ) THEN ! above the first level of data 198 ztp(jk) = ptsd(ji,jj,1 ,jp_tem) 199 zsp(jk) = ptsd(ji,jj,1 ,jp_sal) 200 ELSEIF( zl > gdept_1d(jpk) ) THEN ! below the last level of data 201 ztp(jk) = ptsd(ji,jj,jpkm1,jp_tem) 202 zsp(jk) = ptsd(ji,jj,jpkm1,jp_sal) 203 ELSE ! inbetween : vertical interpolation between jkk & jkk+1 204 DO jkk = 1, jpkm1 ! when gdept(jkk) < zl < gdept(jkk+1) 205 IF( (zl-gdept_1d(jkk)) * (zl-gdept_1d(jkk+1)) <= 0._wp ) THEN 206 zi = ( zl - gdept_1d(jkk) ) / (gdept_1d(jkk+1)-gdept_1d(jkk)) 207 ztp(jk) = ptsd(ji,jj,jkk,jp_tem) + ( ptsd(ji,jj,jkk+1,jp_tem) - ptsd(ji,jj,jkk,jp_tem) ) * zi 208 zsp(jk) = ptsd(ji,jj,jkk,jp_sal) + ( ptsd(ji,jj,jkk+1,jp_sal) - ptsd(ji,jj,jkk,jp_sal) ) * zi 209 ENDIF 210 END DO 211 ENDIF 212 END DO 213 DO jk = 1, jpkm1 214 ptsd(ji,jj,jk,jp_tem) = ztp(jk) * tmask(ji,jj,jk) ! mask required for mixed zps-s-coord 215 ptsd(ji,jj,jk,jp_sal) = zsp(jk) * tmask(ji,jj,jk) 216 END DO 217 ptsd(ji,jj,jpk,jp_tem) = 0._wp 218 ptsd(ji,jj,jpk,jp_sal) = 0._wp 193 DO_2D_11_11 194 DO jk = 1, jpk ! determines the intepolated T-S profiles at each (i,j) points 195 zl = gdept_0(ji,jj,jk) 196 IF( zl < gdept_1d(1 ) ) THEN ! above the first level of data 197 ztp(jk) = ptsd(ji,jj,1 ,jp_tem) 198 zsp(jk) = ptsd(ji,jj,1 ,jp_sal) 199 ELSEIF( zl > gdept_1d(jpk) ) THEN ! below the last level of data 200 ztp(jk) = ptsd(ji,jj,jpkm1,jp_tem) 201 zsp(jk) = ptsd(ji,jj,jpkm1,jp_sal) 202 ELSE ! inbetween : vertical interpolation between jkk & jkk+1 203 DO jkk = 1, jpkm1 ! when gdept(jkk) < zl < gdept(jkk+1) 204 IF( (zl-gdept_1d(jkk)) * (zl-gdept_1d(jkk+1)) <= 0._wp ) THEN 205 zi = ( zl - gdept_1d(jkk) ) / (gdept_1d(jkk+1)-gdept_1d(jkk)) 206 ztp(jk) = ptsd(ji,jj,jkk,jp_tem) + ( ptsd(ji,jj,jkk+1,jp_tem) - ptsd(ji,jj,jkk,jp_tem) ) * zi 207 zsp(jk) = ptsd(ji,jj,jkk,jp_sal) + ( ptsd(ji,jj,jkk+1,jp_sal) - ptsd(ji,jj,jkk,jp_sal) ) * zi 208 ENDIF 209 END DO 210 ENDIF 219 211 END DO 220 END DO 212 DO jk = 1, jpkm1 213 ptsd(ji,jj,jk,jp_tem) = ztp(jk) * tmask(ji,jj,jk) ! mask required for mixed zps-s-coord 214 ptsd(ji,jj,jk,jp_sal) = zsp(jk) * tmask(ji,jj,jk) 215 END DO 216 ptsd(ji,jj,jpk,jp_tem) = 0._wp 217 ptsd(ji,jj,jpk,jp_sal) = 0._wp 218 END_2D 221 219 ! 222 220 ELSE !== z- or zps- coordinate ==! … … 226 224 ! 227 225 IF( ln_zps ) THEN ! zps-coordinate (partial steps) interpolation at the last ocean level 228 DO jj = 1, jpj 229 DO ji = 1, jpi 230 ik = mbkt(ji,jj) 231 IF( ik > 1 ) THEN 232 zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 233 ptsd(ji,jj,ik,jp_tem) = (1.-zl) * ptsd(ji,jj,ik,jp_tem) + zl * ptsd(ji,jj,ik-1,jp_tem) 234 ptsd(ji,jj,ik,jp_sal) = (1.-zl) * ptsd(ji,jj,ik,jp_sal) + zl * ptsd(ji,jj,ik-1,jp_sal) 235 ENDIF 236 ik = mikt(ji,jj) 237 IF( ik > 1 ) THEN 238 zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 239 ptsd(ji,jj,ik,jp_tem) = (1.-zl) * ptsd(ji,jj,ik,jp_tem) + zl * ptsd(ji,jj,ik+1,jp_tem) 240 ptsd(ji,jj,ik,jp_sal) = (1.-zl) * ptsd(ji,jj,ik,jp_sal) + zl * ptsd(ji,jj,ik+1,jp_sal) 241 END IF 242 END DO 243 END DO 226 DO_2D_11_11 227 ik = mbkt(ji,jj) 228 IF( ik > 1 ) THEN 229 zl = ( gdept_1d(ik) - gdept_0(ji,jj,ik) ) / ( gdept_1d(ik) - gdept_1d(ik-1) ) 230 ptsd(ji,jj,ik,jp_tem) = (1.-zl) * ptsd(ji,jj,ik,jp_tem) + zl * ptsd(ji,jj,ik-1,jp_tem) 231 ptsd(ji,jj,ik,jp_sal) = (1.-zl) * ptsd(ji,jj,ik,jp_sal) + zl * ptsd(ji,jj,ik-1,jp_sal) 232 ENDIF 233 ik = mikt(ji,jj) 234 IF( ik > 1 ) THEN 235 zl = ( gdept_0(ji,jj,ik) - gdept_1d(ik) ) / ( gdept_1d(ik+1) - gdept_1d(ik) ) 236 ptsd(ji,jj,ik,jp_tem) = (1.-zl) * ptsd(ji,jj,ik,jp_tem) + zl * ptsd(ji,jj,ik+1,jp_tem) 237 ptsd(ji,jj,ik,jp_sal) = (1.-zl) * ptsd(ji,jj,ik,jp_sal) + zl * ptsd(ji,jj,ik+1,jp_sal) 238 END IF 239 END_2D 244 240 ENDIF 245 241 !
Note: See TracChangeset
for help on using the changeset viewer.