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 4370 for branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90 – NEMO

Ignore:
Timestamp:
2014-01-23T18:13:16+01:00 (10 years ago)
Author:
jchanut
Message:

Time split cleaning / AMM12 restart / BDY and Agrif. See tickets #1228, #1227, #1225 and #1133

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r4366 r4370  
    182182      END DO 
    183183 
     184      ! Before depth and Inverse of the local depth of the water column at u- and v- points 
     185      ! ----------------------------------------------------------------------------------- 
     186      hu_b(:,:) = 0. 
     187      hv_b(:,:) = 0. 
     188      DO jk = 1, jpkm1 
     189         hu_b(:,:) = hu_b(:,:) + fse3u_b(:,:,jk) * umask(:,:,jk) 
     190         hv_b(:,:) = hv_b(:,:) + fse3v_b(:,:,jk) * vmask(:,:,jk) 
     191      END DO 
     192      hur_b(:,:) = umask(:,:,1) / ( hu_b(:,:) + 1. - umask(:,:,1) ) 
     193      hvr_b(:,:) = vmask(:,:,1) / ( hv_b(:,:) + 1. - vmask(:,:,1) ) 
     194 
    184195      ! Restoring frequencies for z_tilde coordinate 
    185196      ! ============================================ 
     
    514525      CALL dom_vvl_interpol( fse3t_a(:,:,:), fse3v_a(:,:,:), 'V' ) 
    515526 
     527      ! *********************************** ! 
     528      ! After depths at u- v points         ! 
     529      ! *********************************** ! 
     530 
     531      hu_a(:,:) = 0._wp                        ! Ocean depth at U-points 
     532      hv_a(:,:) = 0._wp                        ! Ocean depth at V-points 
     533      DO jk = 1, jpkm1 
     534         hu_a(:,:) = hu_a(:,:) + fse3u_a(:,:,jk) * umask(:,:,jk) 
     535         hv_a(:,:) = hv_a(:,:) + fse3v_a(:,:,jk) * vmask(:,:,jk) 
     536      END DO 
     537      !                                        ! Inverse of the local depth 
     538      hur_a(:,:) = 1._wp / ( hu_a(:,:) + 1._wp - umask(:,:,1) ) * umask(:,:,1) 
     539      hvr_a(:,:) = 1._wp / ( hv_a(:,:) + 1._wp - vmask(:,:,1) ) * vmask(:,:,1) 
     540 
    516541      CALL wrk_dealloc( jpi, jpj, zht, z_scale, zwu, zwv, zhdiv ) 
    517542      CALL wrk_dealloc( jpi, jpj, jpk, ze3t                     ) 
     
    584609      ! -------------------------------------- 
    585610      ! - ML - fse3u_b and fse3v_b are allready computed in dynnxt 
     611      ! - JC - hu_b, hv_b, hur_b, hvr_b also 
    586612      CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n (:,:,:), 'F'  ) 
    587613      ! Vertical scale factor interpolations 
     
    604630      ! Local depth and Inverse of the local depth of the water column at u- and v- points 
    605631      ! ---------------------------------------------------------------------------------- 
    606       hu(:,:) = 0. 
    607       hv(:,:) = 0. 
    608       DO jk = 1, jpk 
    609          hu(:,:) = hu(:,:) + fse3u_n(:,:,jk) * umask(:,:,jk) 
    610          hv(:,:) = hv(:,:) + fse3v_n(:,:,jk) * vmask(:,:,jk) 
     632      hu (:,:) = hu_a (:,:) 
     633      hv (:,:) = hv_a (:,:) 
     634 
     635      ! Inverse of the local depth 
     636      hur(:,:) = hur_a(:,:) 
     637      hvr(:,:) = hvr_a(:,:) 
     638 
     639      ! Local depth of the water column at t- points 
     640      ! -------------------------------------------- 
     641      ht(:,:) = 0. 
     642      DO jk = 1, jpkm1 
     643         ht(:,:) = ht(:,:) + fse3t_n(:,:,jk) * tmask(:,:,jk) 
    611644      END DO 
    612       ! Inverse of the local depth 
    613       hur(:,:) = umask(:,:,1) / ( hu(:,:) + 1. - umask(:,:,1) ) 
    614       hvr(:,:) = vmask(:,:,1) / ( hv(:,:) + 1. - vmask(:,:,1) ) 
    615645 
    616646      ! Write outputs 
     
    13221352   !!====================================================================== 
    13231353END MODULE domvvl 
     1354 
     1355 
     1356 
Note: See TracChangeset for help on using the changeset viewer.