Changeset 1323


Ignore:
Timestamp:
2009-02-20T10:54:40+01:00 (12 years ago)
Author:
cetlod
Message:

computation horizontal derivative of density in dtadyn.F90, see ticket:347

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OFF_SRC/dtadyn.F90

    r1291 r1323  
    2323   USE zdfddm          ! vertical  physics: double diffusion 
    2424   USE lbclnk          ! ocean lateral boundary conditions (or mpp link) 
     25   USE zpshde 
    2526   USE lib_mpp         ! distributed memory computing library 
    2627 
     
    143144 
    144145      !! * Local declarations 
    145  
    146  
    147146      NAMELIST/namdyn/ ndtadyn, ndtatot, nsptint,            &  
    148147          &                nficdyn, lperdyn 
     
    198197 
    199198      !! * Local declarations 
    200       INTEGER ::   iper, iperm1, iswap    
     199      INTEGER ::   ji,jj,iper, iperm1, iswap    
    201200 
    202201      REAL(wp) :: zpdtan, zpdtpe, zdemi, zt 
     
    297296            CALL eos( tn, sn, rhd, rhop )   ! Time-filtered in situ density  
    298297            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 ) 
    299302            CALL zdf_mxl( kt )              ! mixed layer depth 
    300303            CALL ldf_slp( kt, rhd, rn2 ) 
    301304 
    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(:,:,:) 
    306309#endif 
    307310       ! 
     
    369372            CALL eos( tn, sn, rhd, rhop )   ! Time-filtered in situ density 
    370373            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 ) 
    371378            CALL zdf_mxl( kt )              ! mixed layer depth 
    372379            CALL ldf_slp( kt, rhd, rn2 ) 
     
    466473            CALL eos( tn, sn, rhd, rhop )   ! Time-filtered in situ density 
    467474            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 ) 
    468479            CALL zdf_mxl( kt )              ! mixed layer depth 
    469480            CALL ldf_slp( kt, rhd, rn2 ) 
     
    551562                    wndm(:,:) = flx(:,:,jpwind) 
    552563                    fr_i(:,:) = flx(:,:,jpice) 
    553                     emp(:,:) = flx(:,:,jpemp) 
     564                    emp(:,:)  = flx(:,:,jpemp) 
    554565                    emps(:,:) = emp(:,:) 
    555                     qsr(:,:) = flx(:,:,jpqsr) 
     566                    qsr(:,:)  = flx(:,:,jpqsr) 
    556567 
    557568          END IF  
     
    670681         zemp, zqsr, zmld, zice, zwspd 
    671682#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 
    674684#endif 
    675685 
    676686#if ! defined key_off_degrad 
    677  
    678687#  if defined key_traldf_c2d 
    679       REAL(wp), DIMENSION(jpi,jpj) ::   & 
    680          zahtw 
     688      REAL(wp), DIMENSION(jpi,jpj) :: zahtw  
    681689#   if defined key_trcldf_eiv 
    682       REAL(wp), DIMENSION(jpi,jpj) ::   & 
    683          zaeiw 
     690      REAL(wp), DIMENSION(jpi,jpj) :: zaeiw  
    684691#   endif 
    685692#  endif 
    686693 
    687694#else 
    688  
    689695   REAL(wp), DIMENSION(jpi,jpj,jpk) ::   & 
    690696      zahtu, zahtv, zahtw  !  Lateral diffusivity 
     
    707713         clname_v = 'dyna_grid_V.nc', & 
    708714         clname_w = 'dyna_grid_W.nc' 
    709       ! 
     715      !--------------------------------------------------------------- 
    710716      ! 0. Initialization 
     717       
    711718      ! cas d'un fichier non periodique : on utilise deux fois le premier et 
    712719      ! le dernier champ temporel 
Note: See TracChangeset for help on using the changeset viewer.