- Timestamp:
- 2015-04-07T10:28:07+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_v3_4_STABLE_2012/NEMOGCM/NEMO/OPA_SRC/DOM/istate.F90
r5169 r5196 113 113 ELSEIF( cp_cfg == 'gyre' ) THEN 114 114 CALL istate_gyre ! GYRE configuration : start from pre-defined T-S fields 115 ELSEIF( cp_cfg == 'seabass' ) THEN 116 CALL istate_seabass 115 117 ELSEIF( ln_tsd_init ) THEN ! Initial T-S fields read in files 116 118 CALL dta_tsd( nit000, tsb ) ! read 3D T and S data at nit000 … … 412 414 413 415 416 417 SUBROUTINE istate_seabass 418 !!---------------------------------------------------------------------- 419 !! *** ROUTINE istate_seabass *** 420 !! 421 !! ** Purpose : Initialization of the dynamics and tracers for seabass 422 !! configuration (double gyre) 423 !! 424 !! ** Method : - set temperature field following Chassignet and Gent, JPO 425 !! 21, pp1290-1299, 1991, and the law 426 !! rho/rho0=1-2.e-4(T-T0) 427 !! - set salinity field constant 428 !! 429 !!---------------------------------------------------------------------- 430 !! * Local variables 431 INTEGER :: ji, jj, jk ! dummy loop indices 432 REAL(wp) :: zsal = 35.5 433 !!---------------------------------------------------------------------- 434 435 IF(lwp) WRITE(numout,*) 436 IF(lwp) WRITE(numout,*) 'istate_seabass : initial analytical T and constant S profiles ' 437 IF(lwp) WRITE(numout,*) '~~~~~~~~~~~' 438 439 DO jk = 1, jpk 440 DO jj = 1, jpj 441 DO ji = 1, jpi 442 tsn(ji,jj,jk,jp_tem) = ( 25.+5.9e-5*800./9.81/2.e-4* & 443 & (exp(-gdept_0(jk)/800.)-1.)) * tmask(ji,jj,jk) 444 tsb(ji,jj,jk,jp_tem) = tsn(ji,jj,jk, jp_tem) 445 END DO 446 END DO 447 END DO 448 449 tsn(:,:,:,jp_sal) = zsal * tmask(:,:,:) 450 tsb(:,:,:,jp_sal) = tsn(:,:,:,jp_sal) 451 452 IF(lwp) THEN 453 WRITE(numout,*) 454 WRITE(numout,*) ' Initial temperature and salinity profiles:' 455 WRITE(numout, "(9x,' level gdept temperature salinity ')" ) 456 WRITE(numout, "(10x, i4, 3f10.2)" ) ( jk, gdept_0(jk), tsn(2,2,jk,jp_tem), tsn(2,2,jk,jp_sal), jk = 1, jpk ) 457 ENDIF 458 459 460 END SUBROUTINE istate_seabass 461 462 463 414 464 SUBROUTINE istate_uvg 415 465 !!----------------------------------------------------------------------
Note: See TracChangeset
for help on using the changeset viewer.