New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8143 for branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90 – NEMO

Ignore:
Timestamp:
2017-06-06T15:55:44+02:00 (7 years ago)
Author:
gm
Message:

#1880 (HPC-09) - step-7: top/bottom drag computed at T-points, zdfbfr.F90 replaced by zdfdrg.F90 + changes in namelist

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/DIA/diawri.F90

    r8093 r8143  
    6161   USE diurnal_bulk    ! diurnal warm layer 
    6262   USE cool_skin       ! Cool skin 
    63    USE wrk_nemo        ! working array 
    6463 
    6564   IMPLICIT NONE 
     
    183182         DO jj = 2, jpjm1 
    184183            DO ji = fs_2, fs_jpim1   ! vector opt. 
    185 !!gm old 
    186 !!gm BUG  missing x 0.5 
    187                zztmpx = (  bfrua(ji  ,jj) * un(ji  ,jj,mbku(ji  ,jj))  & 
    188                       &  + bfrua(ji-1,jj) * un(ji-1,jj,mbku(ji-1,jj))  )       
    189                zztmpy = (  bfrva(ji,  jj) * vn(ji,jj  ,mbkv(ji,jj  ))  & 
    190                       &  + bfrva(ji,jj-1) * vn(ji,jj-1,mbkv(ji,jj-1))  )  
    191                z2d(ji,jj) = rau0 * SQRT( zztmpx * zztmpx + zztmpy * zztmpy ) * tmask(ji,jj,1)  
    192 !!gm 
    193184               zztmp2 = (  ( rCdU_bot(ji+1,jj)+rCdU_bot(ji  ,jj) ) * un(ji  ,jj,mbku(ji  ,jj))  )**2   & 
    194185                  &   + (  ( rCdU_bot(ji  ,jj)+rCdU_bot(ji-1,jj) ) * un(ji-1,jj,mbku(ji-1,jj))  )**2   & 
     
    196187                  &   + (  ( rCdU_bot(ji,jj  )+rCdU_bot(ji,jj-1) ) * vn(ji,jj-1,mbkv(ji,jj-1))  )**2 
    197188               z2d(ji,jj) = zztmp * SQRT( zztmp2 ) * tmask(ji,jj,1)  
    198 !!gm new end 
    199189               ! 
    200             ENDDO 
    201          ENDDO 
     190            END DO 
     191         END DO 
    202192         CALL lbc_lnk( z2d, 'T', 1. ) 
    203193         CALL iom_put( "taubot", z2d )            
     
    449439      REAL(wp) ::   zsto, zout, zmax, zjulian                ! local scalars 
    450440      ! 
    451       REAL(wp), POINTER, DIMENSION(:,:)   :: zw2d       ! 2D workspace 
    452       REAL(wp), POINTER, DIMENSION(:,:,:) :: zw3d       ! 3D workspace 
     441      REAL(wp), DIMENSION(jpi,jpj)   :: zw2d       ! 2D workspace 
     442      REAL(wp), DIMENSION(jpi,jpj,jpk) :: zw3d       ! 3D workspace 
    453443      !!---------------------------------------------------------------------- 
    454444      !  
    455445      IF( nn_timing == 1 )   CALL timing_start('dia_wri') 
    456446      ! 
    457                              CALL wrk_alloc( jpi,jpj      , zw2d ) 
    458       IF( .NOT.ln_linssh )   CALL wrk_alloc( jpi,jpj,jpk  , zw3d ) 
    459       ! 
    460       ! Output the initial state and forcings 
    461       IF( ninist == 1 ) THEN                        
     447      IF( ninist == 1 ) THEN     !==  Output the initial state and forcings  ==! 
    462448         CALL dia_wri_state( 'output.init', kt ) 
    463449         ninist = 0 
     
    467453      ! ----------------- 
    468454 
    469       ! local variable for debugging 
    470       ll_print = .FALSE. 
     455      ll_print = .FALSE.                  ! local variable for debugging 
    471456      ll_print = ll_print .AND. lwp 
    472457 
     
    891876      ENDIF 
    892877      ! 
    893                              CALL wrk_dealloc( jpi , jpj        , zw2d ) 
    894       IF( .NOT.ln_linssh )   CALL wrk_dealloc( jpi , jpj , jpk  , zw3d ) 
    895       ! 
    896878      IF( nn_timing == 1 )   CALL timing_stop('dia_wri') 
    897879      ! 
Note: See TracChangeset for help on using the changeset viewer.