- Timestamp:
- 2017-03-17T08:46:30+01:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90
r7217 r7806 23 23 USE dom_oce ! domain: ocean 24 24 USE sbc_oce ! surface boundary condition: ocean 25 USE trc_oce ! shared ocean-passive tracers variables 25 26 USE phycst ! physical constants 26 27 USE closea ! closed seas … … 97 98 END DO 98 99 ! 99 IF( lk_vvl ) CALL dom_vvl_init ! Vertical variable mesh 100 ! 101 IF( lk_c1d ) CALL cor_c1d ! 1D configuration: Coriolis set at T-point 102 ! 103 ! 104 hu(:,:) = 0._wp ! Ocean depth at U-points 105 hv(:,:) = 0._wp ! Ocean depth at V-points 106 ht(:,:) = 0._wp ! Ocean depth at T-points 107 DO jk = 1, jpkm1 108 hu(:,:) = hu(:,:) + fse3u_n(:,:,jk) * umask(:,:,jk) 109 hv(:,:) = hv(:,:) + fse3v_n(:,:,jk) * vmask(:,:,jk) 110 ht(:,:) = ht(:,:) + fse3t_n(:,:,jk) * tmask(:,:,jk) 111 END DO 112 ! ! Inverse of the local depth 113 hur(:,:) = 1._wp / ( hu(:,:) + 1._wp - umask_i(:,:) ) * umask_i(:,:) 114 hvr(:,:) = 1._wp / ( hv(:,:) + 1._wp - vmask_i(:,:) ) * vmask_i(:,:) 100 IF( lk_c1d ) CALL cor_c1d ! 1D configuration: Coriolis set at T-point 101 ! 102 IF( .NOT.lk_offline ) THEN 103 ! 104 IF( lk_vvl ) CALL dom_vvl_init ! Vertical variable mesh 105 ! 106 hu(:,:) = 0._wp ! Ocean depth at U-points 107 hv(:,:) = 0._wp ! Ocean depth at V-points 108 ht(:,:) = 0._wp ! Ocean depth at T-points 109 DO jk = 1, jpkm1 110 hu(:,:) = hu(:,:) + fse3u_n(:,:,jk) * umask(:,:,jk) 111 hv(:,:) = hv(:,:) + fse3v_n(:,:,jk) * vmask(:,:,jk) 112 ht(:,:) = ht(:,:) + fse3t_n(:,:,jk) * tmask(:,:,jk) 113 END DO 114 ! ! Inverse of the local depth 115 hur(:,:) = 1._wp / ( hu(:,:) + 1._wp - umask_i(:,:) ) * umask_i(:,:) 116 hvr(:,:) = 1._wp / ( hv(:,:) + 1._wp - vmask_i(:,:) ) * vmask_i(:,:) 117 ! 118 ENDIF 115 119 116 120 CALL dom_stp ! time step
Note: See TracChangeset
for help on using the changeset viewer.