Changeset 1323
- Timestamp:
- 2009-02-20T10:54:40+01:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OFF_SRC/dtadyn.F90
r1291 r1323 23 23 USE zdfddm ! vertical physics: double diffusion 24 24 USE lbclnk ! ocean lateral boundary conditions (or mpp link) 25 USE zpshde 25 26 USE lib_mpp ! distributed memory computing library 26 27 … … 143 144 144 145 !! * Local declarations 145 146 147 146 NAMELIST/namdyn/ ndtadyn, ndtatot, nsptint, & 148 147 & nficdyn, lperdyn … … 198 197 199 198 !! * Local declarations 200 INTEGER :: iper, iperm1, iswap199 INTEGER :: ji,jj,iper, iperm1, iswap 201 200 202 201 REAL(wp) :: zpdtan, zpdtpe, zdemi, zt … … 297 296 CALL eos( tn, sn, rhd, rhop ) ! Time-filtered in situ density 298 297 CALL bn2( tn, sn, rn2 ) ! before Brunt-Vaisala frequency 298 IF( ln_zps ) & 299 CALL zps_hde( kt, tn , sn , rhd, & ! Partial steps: before Horizontal DErivative 300 gtu, gsu, gru, & ! of t, s, rd at the bottom ocean level 301 gtv, gsv, grv ) 299 302 CALL zdf_mxl( kt ) ! mixed layer depth 300 303 CALL ldf_slp( kt, rhd, rn2 ) 301 304 302 uslpdta (:,:,:,2)=uslp(:,:,:)303 vslpdta (:,:,:,2)=vslp(:,:,:)304 wslpidta(:,:,:,2) =wslpi(:,:,:)305 wslpjdta(:,:,:,2) =wslpj(:,:,:)305 uslpdta (:,:,:,2) = uslp(:,:,:) 306 vslpdta (:,:,:,2) = vslp(:,:,:) 307 wslpidta(:,:,:,2) = wslpi(:,:,:) 308 wslpjdta(:,:,:,2) = wslpj(:,:,:) 306 309 #endif 307 310 ! … … 369 372 CALL eos( tn, sn, rhd, rhop ) ! Time-filtered in situ density 370 373 CALL bn2( tn, sn, rn2 ) ! before Brunt-Vaisala frequency 374 IF( ln_zps ) & 375 CALL zps_hde( kt, tn , sn , rhd, & ! Partial steps: before Horizontal DErivative 376 gtu, gsu, gru, & ! of t, s, rd at the bottom ocean level 377 gtv, gsv, grv ) 371 378 CALL zdf_mxl( kt ) ! mixed layer depth 372 379 CALL ldf_slp( kt, rhd, rn2 ) … … 466 473 CALL eos( tn, sn, rhd, rhop ) ! Time-filtered in situ density 467 474 CALL bn2( tn, sn, rn2 ) ! before Brunt-Vaisala frequency 475 IF( ln_zps ) & 476 CALL zps_hde( kt, tn , sn , rhd, & ! Partial steps: before Horizontal DErivative 477 gtu, gsu, gru, & ! of t, s, rd at the bottom ocean level 478 gtv, gsv, grv ) 468 479 CALL zdf_mxl( kt ) ! mixed layer depth 469 480 CALL ldf_slp( kt, rhd, rn2 ) … … 551 562 wndm(:,:) = flx(:,:,jpwind) 552 563 fr_i(:,:) = flx(:,:,jpice) 553 emp(:,:) = flx(:,:,jpemp)564 emp(:,:) = flx(:,:,jpemp) 554 565 emps(:,:) = emp(:,:) 555 qsr(:,:) = flx(:,:,jpqsr)566 qsr(:,:) = flx(:,:,jpqsr) 556 567 557 568 END IF … … 670 681 zemp, zqsr, zmld, zice, zwspd 671 682 #if defined key_trcbbl_dif || defined key_trcbbl_adv 672 REAL(wp), DIMENSION(jpi,jpj) :: & 673 zbblx, zbbly 683 REAL(wp), DIMENSION(jpi,jpj) :: zbblx, zbbly 674 684 #endif 675 685 676 686 #if ! defined key_off_degrad 677 678 687 # if defined key_traldf_c2d 679 REAL(wp), DIMENSION(jpi,jpj) :: & 680 zahtw 688 REAL(wp), DIMENSION(jpi,jpj) :: zahtw 681 689 # if defined key_trcldf_eiv 682 REAL(wp), DIMENSION(jpi,jpj) :: & 683 zaeiw 690 REAL(wp), DIMENSION(jpi,jpj) :: zaeiw 684 691 # endif 685 692 # endif 686 693 687 694 #else 688 689 695 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & 690 696 zahtu, zahtv, zahtw ! Lateral diffusivity … … 707 713 clname_v = 'dyna_grid_V.nc', & 708 714 clname_w = 'dyna_grid_W.nc' 709 ! 715 !--------------------------------------------------------------- 710 716 ! 0. Initialization 717 711 718 ! cas d'un fichier non periodique : on utilise deux fois le premier et 712 719 ! le dernier champ temporel
Note: See TracChangeset
for help on using the changeset viewer.