- Timestamp:
- 2016-07-19T10:38:35+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/NERC/dev_r5549_BDY_ZEROGRAD/NEMOGCM/NEMO/OPA_SRC/DOM/dtatsd.F90
r5215 r6808 35 35 TYPE(FLD), ALLOCATABLE, DIMENSION(:) :: sf_tsd ! structure of input SST (file informations, fields read) 36 36 37 !! * Substitutions38 # include "domzgr_substitute.h90"39 37 !!---------------------------------------------------------------------- 40 38 !! NEMO/OPA 3.3 , NEMO Consortium (2010) … … 174 172 END DO 175 173 END DO 176 IF( nn_cla == 1 ) THEN ! Cross Land advection 177 il0 = 138 ; il1 = 138 ! set T & S profile at Gibraltar Strait 178 ij0 = 101 ; ij1 = 102 179 ii0 = 139 ; ii1 = 139 180 DO jl = mi0(il0), mi1(il1) 181 DO jj = mj0(ij0), mj1(ij1) 182 DO ji = mi0(ii0), mi1(ii1) 183 sf_tsd(jp_tem)%fnow(ji,jj,:) = sf_tsd(jp_tem)%fnow(jl,jj,:) 184 sf_tsd(jp_sal)%fnow(ji,jj,:) = sf_tsd(jp_sal)%fnow(jl,jj,:) 185 END DO 186 END DO 187 END DO 188 il0 = 164 ; il1 = 164 ! set T & S profile at Bab el Mandeb Strait 189 ij0 = 87 ; ij1 = 88 190 ii0 = 161 ; ii1 = 163 191 DO jl = mi0(il0), mi1(il1) 192 DO jj = mj0(ij0), mj1(ij1) 193 DO ji = mi0(ii0), mi1(ii1) 194 sf_tsd(jp_tem)%fnow(ji,jj,:) = sf_tsd(jp_tem)%fnow(jl,jj,:) 195 sf_tsd(jp_sal)%fnow(ji,jj,:) = sf_tsd(jp_sal)%fnow(jl,jj,:) 196 END DO 197 END DO 198 END DO 199 ELSE ! No Cross Land advection 200 ij0 = 87 ; ij1 = 96 ! Reduced temperature in Red Sea 201 ii0 = 148 ; ii1 = 160 202 sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 4:10 ) = 7.0_wp 203 sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 11:13 ) = 6.5_wp 204 sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 14:20 ) = 6.0_wp 205 ENDIF 174 ij0 = 87 ; ij1 = 96 ! Reduced temperature in Red Sea 175 ii0 = 148 ; ii1 = 160 176 sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 4:10 ) = 7.0_wp 177 sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 11:13 ) = 6.5_wp 178 sf_tsd(jp_tem)%fnow( mi0(ii0):mi1(ii1) , mj0(ij0):mj1(ij1) , 14:20 ) = 6.0_wp 206 179 ENDIF 207 180 ! … … 275 248 ENDIF 276 249 ! 277 IF( lwp .AND. kt == nit000 ) THEN278 WRITE(numout,*) ' temperature Levitus '279 WRITE(numout,*)280 WRITE(numout,*)' level = 1'281 CALL prihre( ptsd(:,:,1 ,jp_tem), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )282 WRITE(numout,*)' level = ', jpk/2283 CALL prihre( ptsd(:,:,jpk/2,jp_tem), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )284 WRITE(numout,*)' level = ', jpkm1285 CALL prihre( ptsd(:,:,jpkm1,jp_tem), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )286 WRITE(numout,*)287 WRITE(numout,*) ' salinity Levitus '288 WRITE(numout,*)289 WRITE(numout,*)' level = 1'290 CALL prihre( ptsd(:,:,1 ,jp_sal), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )291 WRITE(numout,*)' level = ', jpk/2292 CALL prihre( ptsd(:,:,jpk/2,jp_sal), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )293 WRITE(numout,*)' level = ', jpkm1294 CALL prihre( ptsd(:,:,jpkm1,jp_sal), jpi, jpj, 1, jpi, 20, 1, jpj, 20, 1., numout )295 WRITE(numout,*)296 ENDIF297 !298 250 IF( .NOT.ln_tsd_tradmp ) THEN !== deallocate T & S structure ==! 299 251 ! (data used only for initialisation)
Note: See TracChangeset
for help on using the changeset viewer.