- Timestamp:
- 2014-01-23T18:13:16+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90
r4366 r4370 182 182 END DO 183 183 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 184 195 ! Restoring frequencies for z_tilde coordinate 185 196 ! ============================================ … … 514 525 CALL dom_vvl_interpol( fse3t_a(:,:,:), fse3v_a(:,:,:), 'V' ) 515 526 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 516 541 CALL wrk_dealloc( jpi, jpj, zht, z_scale, zwu, zwv, zhdiv ) 517 542 CALL wrk_dealloc( jpi, jpj, jpk, ze3t ) … … 584 609 ! -------------------------------------- 585 610 ! - ML - fse3u_b and fse3v_b are allready computed in dynnxt 611 ! - JC - hu_b, hv_b, hur_b, hvr_b also 586 612 CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n (:,:,:), 'F' ) 587 613 ! Vertical scale factor interpolations … … 604 630 ! Local depth and Inverse of the local depth of the water column at u- and v- points 605 631 ! ---------------------------------------------------------------------------------- 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) 611 644 END DO 612 ! Inverse of the local depth613 hur(:,:) = umask(:,:,1) / ( hu(:,:) + 1. - umask(:,:,1) )614 hvr(:,:) = vmask(:,:,1) / ( hv(:,:) + 1. - vmask(:,:,1) )615 645 616 646 ! Write outputs … … 1322 1352 !!====================================================================== 1323 1353 END MODULE domvvl 1354 1355 1356
Note: See TracChangeset
for help on using the changeset viewer.