- Timestamp:
- 2015-10-31T08:40:45+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5836_NOC3_vvl_by_default/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90
r5836 r5845 52 52 53 53 !! * Substitutions 54 # include "domzgr_substitute.h90"55 54 # include "vectopt_loop_substitute.h90" 56 55 !!---------------------------------------------------------------------- … … 124 123 ENDIF 125 124 ! 126 ! - ML - sshn could be modified by istate_eel, so that initialization of fse3t_b is done here125 ! - ML - sshn could be modified by istate_eel, so that initialization of e3t_b is done here 127 126 IF( lk_vvl ) THEN 128 127 DO jk = 1, jpk 129 fse3t_b(:,:,jk) = fse3t_n(:,:,jk)128 e3t_b(:,:,jk) = e3t_n(:,:,jk) 130 129 END DO 131 130 ENDIF … … 155 154 DO jj = 1, jpj 156 155 DO ji = 1, jpi 157 un_b(ji,jj) = un_b(ji,jj) + fse3u_n(ji,jj,jk) * un(ji,jj,jk) * umask(ji,jj,jk)158 vn_b(ji,jj) = vn_b(ji,jj) + fse3v_n(ji,jj,jk) * vn(ji,jj,jk) * vmask(ji,jj,jk)156 un_b(ji,jj) = un_b(ji,jj) + e3u_n(ji,jj,jk) * un(ji,jj,jk) * umask(ji,jj,jk) 157 vn_b(ji,jj) = vn_b(ji,jj) + e3v_n(ji,jj,jk) * vn(ji,jj,jk) * vmask(ji,jj,jk) 159 158 ! 160 ub_b(ji,jj) = ub_b(ji,jj) + fse3u_b(ji,jj,jk) * ub(ji,jj,jk) * umask(ji,jj,jk)161 vb_b(ji,jj) = vb_b(ji,jj) + fse3v_b(ji,jj,jk) * vb(ji,jj,jk) * vmask(ji,jj,jk)159 ub_b(ji,jj) = ub_b(ji,jj) + e3u_b(ji,jj,jk) * ub(ji,jj,jk) * umask(ji,jj,jk) 160 vb_b(ji,jj) = vb_b(ji,jj) + e3v_b(ji,jj,jk) * vb(ji,jj,jk) * vmask(ji,jj,jk) 162 161 END DO 163 162 END DO 164 163 END DO 165 164 ! 166 un_b(:,:) = un_b(:,:) * hur(:,:)167 vn_b(:,:) = vn_b(:,:) * hvr(:,:)168 ! 169 ub_b(:,:) = ub_b(:,:) * hur_b(:,:)170 vb_b(:,:) = vb_b(:,:) * hvr_b(:,:)165 un_b(:,:) = un_b(:,:) * r1_hu_n (:,:) 166 vn_b(:,:) = vn_b(:,:) * r1_hv_n (:,:) 167 ! 168 ub_b(:,:) = ub_b(:,:) * r1_hu_b(:,:) 169 vb_b(:,:) = vb_b(:,:) * r1_hv_b(:,:) 171 170 ! 172 171 ! … … 197 196 ! 198 197 DO jk = 1, jpk 199 tsn(:,:,jk,jp_tem) = ( ( ( 7.5 - 0. * ABS( gphit(:,:) )/30. ) * ( 1.-TANH(( fsdept(:,:,jk)-80.)/30.) ) &200 & + 10. * ( 5000. - fsdept(:,:,jk) ) /5000.) ) * tmask(:,:,jk)198 tsn(:,:,jk,jp_tem) = ( ( ( 7.5 - 0. * ABS( gphit(:,:) )/30. ) * ( 1.-TANH((gdept_n(:,:,jk)-80.)/30.) ) & 199 & + 10. * ( 5000. - gdept_n(:,:,jk) ) /5000.) ) * tmask(:,:,jk) 201 200 tsb(:,:,jk,jp_tem) = tsn(:,:,jk,jp_tem) 202 201 END DO … … 251 250 ! 252 251 DO jk = 1, jpk 253 tsn(:,:,jk,jp_tem) = ( zt2 + zt1 * exp( - fsdept(:,:,jk) / 1000 ) ) * tmask(:,:,jk)252 tsn(:,:,jk,jp_tem) = ( zt2 + zt1 * exp( - gdept_n(:,:,jk) / 1000 ) ) * tmask(:,:,jk) 254 253 tsb(:,:,jk,jp_tem) = tsn(:,:,jk,jp_tem) 255 254 END DO 256 !257 IF(lwp) CALL prizre( tsn(:,:,:,jp_tem), jpi , jpj , jpk , jpj/2 , &258 & 1 , jpi , 5 , 1 , jpk , &259 & 1 , 1. , numout )260 255 ! 261 256 ! set salinity field to a constant value … … 327 322 tsn(:,:,:,jp_tem) = tsb(:,:,:,jp_tem) ! set nox temperature to tb 328 323 ! 329 IF(lwp) CALL prizre( tsn(:,:,:,jp_tem), jpi , jpj , jpk , jpj/2 , &330 & 1 , jpi , 5 , 1 , jpk , &331 & 1 , 1. , numout )332 !333 324 ! set salinity field to a constant value 334 325 ! -------------------------------------- … … 376 367 DO jj = 1, jpj 377 368 DO ji = 1, jpi 378 tsn(ji,jj,jk,jp_tem) = ( 16. - 12. * TANH( ( fsdept(ji,jj,jk) - 400) / 700 ) ) &379 & * (-TANH( (500- fsdept(ji,jj,jk)) / 150 ) + 1) / 2 &380 & + ( 15. * ( 1. - TANH( ( fsdept(ji,jj,jk)-50.) / 1500.) ) &381 & - 1.4 * TANH(( fsdept(ji,jj,jk)-100.) / 100.) &382 & + 7. * (1500. - fsdept(ji,jj,jk)) / 1500. ) &383 & * (-TANH( ( fsdept(ji,jj,jk) - 500) / 150) + 1) / 2369 tsn(ji,jj,jk,jp_tem) = ( 16. - 12. * TANH( (gdept_n(ji,jj,jk) - 400) / 700 ) ) & 370 & * (-TANH( (500-gdept_n(ji,jj,jk)) / 150 ) + 1) / 2 & 371 & + ( 15. * ( 1. - TANH( (gdept_n(ji,jj,jk)-50.) / 1500.) ) & 372 & - 1.4 * TANH((gdept_n(ji,jj,jk)-100.) / 100.) & 373 & + 7. * (1500. - gdept_n(ji,jj,jk)) / 1500. ) & 374 & * (-TANH( (gdept_n(ji,jj,jk) - 500) / 150) + 1) / 2 384 375 tsn(ji,jj,jk,jp_tem) = tsn(ji,jj,jk,jp_tem) * tmask(ji,jj,jk) 385 376 tsb(ji,jj,jk,jp_tem) = tsn(ji,jj,jk,jp_tem) 386 377 387 tsn(ji,jj,jk,jp_sal) = ( 36.25 - 1.13 * TANH( ( fsdept(ji,jj,jk) - 305) / 460 ) ) &388 & * (-TANH((500 - fsdept(ji,jj,jk)) / 150) + 1) / 2 &389 & + ( 35.55 + 1.25 * (5000. - fsdept(ji,jj,jk)) / 5000. &390 & - 1.62 * TANH( ( fsdept(ji,jj,jk) - 60. ) / 650. ) &391 & + 0.2 * TANH( ( fsdept(ji,jj,jk) - 35. ) / 100. ) &392 & + 0.2 * TANH( ( fsdept(ji,jj,jk) - 1000.) / 5000.) ) &393 & * (-TANH(( fsdept(ji,jj,jk) - 500) / 150) + 1) / 2378 tsn(ji,jj,jk,jp_sal) = ( 36.25 - 1.13 * TANH( (gdept_n(ji,jj,jk) - 305) / 460 ) ) & 379 & * (-TANH((500 - gdept_n(ji,jj,jk)) / 150) + 1) / 2 & 380 & + ( 35.55 + 1.25 * (5000. - gdept_n(ji,jj,jk)) / 5000. & 381 & - 1.62 * TANH( (gdept_n(ji,jj,jk) - 60. ) / 650. ) & 382 & + 0.2 * TANH( (gdept_n(ji,jj,jk) - 35. ) / 100. ) & 383 & + 0.2 * TANH( (gdept_n(ji,jj,jk) - 1000.) / 5000.) ) & 384 & * (-TANH((gdept_n(ji,jj,jk) - 500) / 150) + 1) / 2 394 385 tsn(ji,jj,jk,jp_sal) = tsn(ji,jj,jk,jp_sal) * tmask(ji,jj,jk) 395 386 tsb(ji,jj,jk,jp_sal) = tsn(ji,jj,jk,jp_sal) … … 466 457 zalfg = 0.5 * grav * rau0 467 458 468 zprn(:,:,1) = zalfg * fse3w(:,:,1) * ( 1 + rhd(:,:,1) ) ! Surface value459 zprn(:,:,1) = zalfg * e3w_n(:,:,1) * ( 1 + rhd(:,:,1) ) ! Surface value 469 460 470 461 DO jk = 2, jpkm1 ! Vertical integration from the surface 471 462 zprn(:,:,jk) = zprn(:,:,jk-1) & 472 & + zalfg * fse3w(:,:,jk) * ( 2. + rhd(:,:,jk) + rhd(:,:,jk-1) )463 & + zalfg * e3w_n(:,:,jk) * ( 2. + rhd(:,:,jk) + rhd(:,:,jk-1) ) 473 464 END DO 474 465
Note: See TracChangeset
for help on using the changeset viewer.