- Timestamp:
- 2009-04-06T19:45:45+02:00 (15 years ago)
- Location:
- branches/dev_004_VVL/NEMO/OPA_SRC/DOM
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/dev_004_VVL/NEMO/OPA_SRC/DOM/dom_oce.F90
r1362 r1385 114 114 #if defined key_vvl 115 115 LOGICAL, PUBLIC, PARAMETER :: lk_vvl = .TRUE. !: variable grid flag 116 !! All coordinates 117 !! --------------- 118 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: gdep3w_1 !: depth of T-points (sum of e3w) (m) 119 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: gdept_1, gdepw_1 !: analytical depth at T-W points (m) 120 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: e3v_1 , e3f_1 !: analytical vertical scale factors at V--F 121 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: e3t_1 , e3u_1 !: T--U points (m) 122 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: e3vw_1 !: analytical vertical scale factors at VW-- 123 REAL(wp), PUBLIC, DIMENSION(jpi,jpj,jpk) :: e3w_1 , e3uw_1 !: W--UW points (m) 116 124 #else 117 125 LOGICAL, PUBLIC, PARAMETER :: lk_vvl = .FALSE. !: fixed grid flag -
branches/dev_004_VVL/NEMO/OPA_SRC/DOM/domvvl.F90
r1380 r1385 78 78 #endif 79 79 80 fsdept _0(:,:,:) = gdept (:,:,:)81 fsdepw _0(:,:,:) = gdepw (:,:,:)82 fsde3w _0(:,:,:) = gdep3w(:,:,:)83 fse3t _0(:,:,:) = e3t (:,:,:)84 fse3u _0(:,:,:) = e3u (:,:,:)85 fse3v _0(:,:,:) = e3v (:,:,:)86 fse3f _0(:,:,:) = e3f (:,:,:)87 fse3w _0(:,:,:) = e3w (:,:,:)88 fse3uw _0(:,:,:) = e3uw (:,:,:)89 fse3vw _0(:,:,:) = e3vw (:,:,:)80 fsdept(:,:,:) = gdept (:,:,:) 81 fsdepw(:,:,:) = gdepw (:,:,:) 82 fsde3w(:,:,:) = gdep3w(:,:,:) 83 fse3t (:,:,:) = e3t (:,:,:) 84 fse3u (:,:,:) = e3u (:,:,:) 85 fse3v (:,:,:) = e3v (:,:,:) 86 fse3f (:,:,:) = e3f (:,:,:) 87 fse3w (:,:,:) = e3w (:,:,:) 88 fse3uw(:,:,:) = e3uw (:,:,:) 89 fse3vw(:,:,:) = e3vw (:,:,:) 90 90 91 91 ! mu computation -
branches/dev_004_VVL/NEMO/OPA_SRC/DOM/domzgr_substitute.h90
r1362 r1385 8 8 !! 3.1 ! 2009-02 (G. Madec, M. Leclair) pure z* coordinate 9 9 !!---------------------------------------------------------------------- 10 #if defined key_zco11 # define fsdept(i,j,k) gdept_0(k)12 # define fsdepw(i,j,k) gdepw_0(k)13 # define fsde3w(i,j,k) gdepw_0(k)14 # define fse3t(i,j,k) e3t_0(k)15 # define fse3u(i,j,k) e3t_0(k)16 # define fse3v(i,j,k) e3t_0(k)17 # define fse3f(i,j,k) e3t_0(k)18 # define fse3w(i,j,k) e3w_0(k)19 # define fse3uw(i,j,k) e3w_0(k)20 # define fse3vw(i,j,k) e3w_0(k)21 #else22 # define fsdept(i,j,k) gdept(i,j,k)23 # define fsdepw(i,j,k) gdepw(i,j,k)24 # define fsde3w(i,j,k) gdep3w(i,j,k)25 # define fse3t(i,j,k) e3t(i,j,k)26 # define fse3u(i,j,k) e3u(i,j,k)27 # define fse3v(i,j,k) e3v(i,j,k)28 # define fse3f(i,j,k) e3f(i,j,k)29 # define fse3w(i,j,k) e3w(i,j,k)30 # define fse3uw(i,j,k) e3uw(i,j,k)31 # define fse3vw(i,j,k) e3vw(i,j,k)32 #endif33 10 #if defined key_zco 34 11 ! reference for pure z-coordinate (1D - no i,j and time dependency) … … 57 34 #endif 58 35 #if defined key_vvl 59 ! s* or z*-coordinate (3D + time dependency) 36 ! s* or z*-coordinate (3D + time dependency) + use of additional now arrays (..._1) 37 # define fsdept(i,j,k) gdept_1(i,j,k) 38 # define fsdepw(i,j,k) gdepw_1(i,j,k) 39 # define fsde3w(i,j,k) gdep3w_1(i,j,k) 40 # define fse3t(i,j,k) e3t_1(i,j,k) 41 # define fse3u(i,j,k) e3u_1(i,j,k) 42 # define fse3v(i,j,k) e3v_1(i,j,k) 43 # define fse3f(i,j,k) e3f_1(i,j,k) 44 # define fse3w(i,j,k) e3w_1(i,j,k) 45 # define fse3uw(i,j,k) e3uw_1(i,j,k) 46 # define fse3vw(i,j,k) e3vw_1(i,j,k) 47 60 48 # define fsdept_b(i,j,k) (fsdept_0(i,j,k)*(1+sshb(i,j)*mut(i,j,k))) 61 49 # define fsdepw_b(i,j,k) (fsdepw_0(i,j,k)*(1+sshb(i,j)*mut(i,j,k))) … … 91 79 # define fse3vw_a(i,j,k) (fse3vw_0(i,j,k)*(1+sshv_a(i,j)*muv(i,j,k))) 92 80 #else 93 ! time dependent coordinate (1D or 3D) 94 ! s- or zps-coordinate (3D no time dependency) 81 ! z- or s-coordinate (1D or 3D + no time dependency) use reference in all cases 82 # define fsdept(i,j,k) fsdept_0(i,j,k) 83 # define fsdepw(i,j,k) fsdepw_0(i,j,k) 84 # define fsde3w(i,j,k) fsde3w_0(i,j,k) 85 # define fse3t(i,j,k) fse3t_0(i,j,k) 86 # define fse3u(i,j,k) fse3u_0(i,j,k) 87 # define fse3v(i,j,k) fse3v_0(i,j,k) 88 # define fse3f(i,j,k) fse3f_0(i,j,k) 89 # define fse3w(i,j,k) fse3w_0(i,j,k) 90 # define fse3uw(i,j,k) fse3uw_0(i,j,k) 91 # define fse3vw(i,j,k) fse3vw_0(i,j,k) 92 95 93 # define fsdept_b(i,j,k) fsdept_0(i,j,k) 96 94 # define fsdepw_b(i,j,k) fsdepw_0(i,j,k)
Note: See TracChangeset
for help on using the changeset viewer.