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 14475 for NEMO/trunk/src – NEMO

Changeset 14475 for NEMO/trunk/src


Ignore:
Timestamp:
2021-02-17T17:19:28+01:00 (3 years ago)
Author:
techene
Message:

#2619 bug correction add u/vtau iomput call as done in step.F90

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/src/OCE/DYN/dynatf_qco.F90

    r14224 r14475  
    4343   USE isf_oce   , ONLY: ln_isf     ! ice shelf 
    4444   USE isfdynatf , ONLY: isf_dynatf ! ice shelf volume filter correction subroutine 
     45   USE zdfdrg    , ONLY: ln_drgice_imp, rCdU_top 
    4546   ! 
    4647   USE in_out_manager ! I/O manager 
     
    101102      REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zue, zve 
    102103      REAL(wp), ALLOCATABLE, DIMENSION(:,:,:) ::   zua, zva 
     104      REAL(wp), ALLOCATABLE, DIMENSION(:,:)   ::   zutau, zvtau 
    103105      !!---------------------------------------------------------------------- 
    104106      ! 
     
    239241      ENDIF 
    240242      ! 
     243      IF ( iom_use("utau") ) THEN 
     244         IF ( ln_drgice_imp.OR.ln_isfcav ) THEN 
     245            ALLOCATE(zutau(jpi,jpj)) 
     246            DO_2D( 0, 0, 0, 0 ) 
     247               jk = miku(ji,jj) 
     248               zutau(ji,jj) = utau(ji,jj) + 0.5_wp * rho0 * ( rCdU_top(ji+1,jj)+rCdU_top(ji,jj) ) * puu(ji,jj,jk,Kaa) 
     249            END_2D 
     250            CALL iom_put(  "utau", zutau(:,:) ) 
     251            DEALLOCATE(zutau) 
     252         ELSE 
     253            CALL iom_put(  "utau", utau(:,:) ) 
     254         ENDIF 
     255      ENDIF 
     256      ! 
     257      IF ( iom_use("vtau") ) THEN 
     258         IF ( ln_drgice_imp.OR.ln_isfcav ) THEN 
     259            ALLOCATE(zvtau(jpi,jpj)) 
     260            DO_2D( 0, 0, 0, 0 ) 
     261               jk = mikv(ji,jj) 
     262               zvtau(ji,jj) = vtau(ji,jj) + 0.5_wp * rho0 * ( rCdU_top(ji,jj+1)+rCdU_top(ji,jj) ) * pvv(ji,jj,jk,Kaa) 
     263            END_2D 
     264            CALL iom_put(  "vtau", zvtau(:,:) ) 
     265            DEALLOCATE(zvtau) 
     266         ELSE 
     267            CALL iom_put(  "vtau", vtau(:,:) ) 
     268         ENDIF 
     269      ENDIF 
     270      ! 
    241271      IF(sn_cfctl%l_prtctl)   CALL prt_ctl( tab3d_1=puu(:,:,:,Kaa), clinfo1=' nxt  - puu(:,:,:,Kaa): ', mask1=umask,   & 
    242272         &                                  tab3d_2=pvv(:,:,:,Kaa), clinfo2=' pvv(:,:,:,Kaa): '       , mask2=vmask ) 
Note: See TracChangeset for help on using the changeset viewer.