- Timestamp:
- 2015-05-18T17:19:50+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5056_CMCC4_simplification/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90
r4998 r5282 156 156 ! Vertical scale factor interpolations 157 157 ! ------------------------------------ 158 CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n 158 CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n(:,:,:), 'W' ) 159 159 CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3uw_n(:,:,:), 'UW' ) 160 160 CALL dom_vvl_interpol( fse3v_n(:,:,:), fse3vw_n(:,:,:), 'VW' ) 161 CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3w_b 161 CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3w_b(:,:,:), 'W' ) 162 162 CALL dom_vvl_interpol( fse3u_b(:,:,:), fse3uw_b(:,:,:), 'UW' ) 163 163 CALL dom_vvl_interpol( fse3v_b(:,:,:), fse3vw_b(:,:,:), 'VW' ) … … 627 627 ! - ML - fse3u_b and fse3v_b are allready computed in dynnxt 628 628 ! - JC - hu_b, hv_b, hur_b, hvr_b also 629 CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n 629 CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3f_n(:,:,:), 'F' ) 630 630 ! Vertical scale factor interpolations 631 631 ! ------------------------------------ 632 CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n 632 CALL dom_vvl_interpol( fse3t_n(:,:,:), fse3w_n(:,:,:), 'W' ) 633 633 CALL dom_vvl_interpol( fse3u_n(:,:,:), fse3uw_n(:,:,:), 'UW' ) 634 634 CALL dom_vvl_interpol( fse3v_n(:,:,:), fse3vw_n(:,:,:), 'VW' ) 635 CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3w_b 635 CALL dom_vvl_interpol( fse3t_b(:,:,:), fse3w_b(:,:,:), 'W' ) 636 636 CALL dom_vvl_interpol( fse3u_b(:,:,:), fse3uw_b(:,:,:), 'UW' ) 637 637 CALL dom_vvl_interpol( fse3v_b(:,:,:), fse3vw_b(:,:,:), 'VW' ) … … 663 663 ! Local depth and Inverse of the local depth of the water column at u- and v- points 664 664 ! ---------------------------------------------------------------------------------- 665 hu (:,:) = hu_a(:,:)666 hv (:,:) = hv_a(:,:)665 hu(:,:) = hu_a(:,:) 666 hv(:,:) = hv_a(:,:) 667 667 668 668 ! Inverse of the local depth … … 680 680 ! ============= 681 681 z_e3t_def(:,:,:) = ( ( fse3t_n(:,:,:) - e3t_0(:,:,:) ) / e3t_0(:,:,:) * 100 * tmask(:,:,:) ) ** 2 682 CALL iom_put( "cellthc" , fse3t_n 683 CALL iom_put( "tpt_dep" , fsde3w_n 682 CALL iom_put( "cellthc" , fse3t_n(:,:,:) ) 683 CALL iom_put( "tpt_dep" , fsde3w_n(:,:,:) ) 684 684 CALL iom_put( "e3tdef" , z_e3t_def(:,:,:) ) 685 685 … … 1057 1057 ! ! ===================== 1058 1058 !! acc 1059 IF( nn_cla == 0 ) THEN1060 !1061 1059 ii0 = 139 ; ii1 = 140 ! Gibraltar Strait (e2u was modified) 1062 1060 ij0 = 102 ; ij1 = 102 … … 1105 1103 END DO 1106 1104 END DO 1107 ENDIF1108 1105 1109 1106 ii0 = 145 ; ii1 = 146 ! Danish Straits (e2u was modified)
Note: See TracChangeset
for help on using the changeset viewer.