- Timestamp:
- 2013-11-20T17:28:04+01:00 (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90
r3680 r4292 183 183 CALL iom_rstput( 0, 0, inum4, 'esigw', esigw ) 184 184 ! 185 CALL iom_rstput( 0, 0, inum4, 'e3t ', e3t )! ! scale factors186 CALL iom_rstput( 0, 0, inum4, 'e3u ', e3u)187 CALL iom_rstput( 0, 0, inum4, 'e3v ', e3v)188 CALL iom_rstput( 0, 0, inum4, 'e3w ', e3w)185 CALL iom_rstput( 0, 0, inum4, 'e3t_0', e3t_0 ) ! ! scale factors 186 CALL iom_rstput( 0, 0, inum4, 'e3u_0', e3u_0 ) 187 CALL iom_rstput( 0, 0, inum4, 'e3v_0', e3v_0 ) 188 CALL iom_rstput( 0, 0, inum4, 'e3w_0', e3w_0 ) 189 189 CALL iom_rstput( 0, 0, inum4, 'rx1', rx1 ) ! ! Max. grid stiffness ratio 190 190 ! 191 CALL iom_rstput( 0, 0, inum4, 'gdept ' , gdept )! ! stretched system192 CALL iom_rstput( 0, 0, inum4, 'gdepw ' , gdepw)191 CALL iom_rstput( 0, 0, inum4, 'gdept_1d' , gdept_1d ) ! ! stretched system 192 CALL iom_rstput( 0, 0, inum4, 'gdepw_1d' , gdepw_1d ) 193 193 ENDIF 194 194 … … 196 196 ! 197 197 IF( nmsh <= 6 ) THEN ! ! 3D vertical scale factors 198 CALL iom_rstput( 0, 0, inum4, 'e3t ', e3t)199 CALL iom_rstput( 0, 0, inum4, 'e3u ', e3u)200 CALL iom_rstput( 0, 0, inum4, 'e3v ', e3v)201 CALL iom_rstput( 0, 0, inum4, 'e3w ', e3w)198 CALL iom_rstput( 0, 0, inum4, 'e3t_0', e3t_0 ) 199 CALL iom_rstput( 0, 0, inum4, 'e3u_0', e3u_0 ) 200 CALL iom_rstput( 0, 0, inum4, 'e3v_0', e3v_0 ) 201 CALL iom_rstput( 0, 0, inum4, 'e3w_0', e3w_0 ) 202 202 ELSE ! ! 2D masked bottom ocean scale factors 203 203 DO jj = 1,jpj 204 204 DO ji = 1,jpi 205 e3tp(ji,jj) = e3t (ji,jj,mbkt(ji,jj)) * tmask(ji,jj,1)206 e3wp(ji,jj) = e3w (ji,jj,mbkt(ji,jj)) * tmask(ji,jj,1)205 e3tp(ji,jj) = e3t_0(ji,jj,mbkt(ji,jj)) * tmask(ji,jj,1) 206 e3wp(ji,jj) = e3w_0(ji,jj,mbkt(ji,jj)) * tmask(ji,jj,1) 207 207 END DO 208 208 END DO … … 212 212 ! 213 213 IF( nmsh <= 3 ) THEN ! ! 3D depth 214 CALL iom_rstput( 0, 0, inum4, 'gdept ', gdept, ktype = jp_r4 )214 CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0, ktype = jp_r4 ) 215 215 DO jk = 1,jpk 216 216 DO jj = 1, jpjm1 217 217 DO ji = 1, fs_jpim1 ! vector opt. 218 zdepu(ji,jj,jk) = MIN( gdept (ji,jj,jk) , gdept(ji+1,jj ,jk) )219 zdepv(ji,jj,jk) = MIN( gdept (ji,jj,jk) , gdept(ji ,jj+1,jk) )218 zdepu(ji,jj,jk) = MIN( gdept_0(ji,jj,jk) , gdept_0(ji+1,jj ,jk) ) 219 zdepv(ji,jj,jk) = MIN( gdept_0(ji,jj,jk) , gdept_0(ji ,jj+1,jk) ) 220 220 END DO 221 221 END DO … … 224 224 CALL iom_rstput( 0, 0, inum4, 'gdepu', zdepu, ktype = jp_r4 ) 225 225 CALL iom_rstput( 0, 0, inum4, 'gdepv', zdepv, ktype = jp_r4 ) 226 CALL iom_rstput( 0, 0, inum4, 'gdepw ', gdepw, ktype = jp_r4 )226 CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r4 ) 227 227 ELSE ! ! 2D bottom depth 228 228 DO jj = 1,jpj 229 229 DO ji = 1,jpi 230 zprt(ji,jj) = gdept (ji,jj,mbkt(ji,jj) ) * tmask(ji,jj,1)231 zprw(ji,jj) = gdepw (ji,jj,mbkt(ji,jj)+1) * tmask(ji,jj,1)230 zprt(ji,jj) = gdept_0(ji,jj,mbkt(ji,jj) ) * tmask(ji,jj,1) 231 zprw(ji,jj) = gdepw_0(ji,jj,mbkt(ji,jj)+1) * tmask(ji,jj,1) 232 232 END DO 233 233 END DO … … 236 236 ENDIF 237 237 ! 238 CALL iom_rstput( 0, 0, inum4, 'gdept_ 0', gdept_0 )! ! reference z-coord.239 CALL iom_rstput( 0, 0, inum4, 'gdepw_ 0', gdepw_0)240 CALL iom_rstput( 0, 0, inum4, 'e3t_ 0' , e3t_0)241 CALL iom_rstput( 0, 0, inum4, 'e3w_ 0' , e3w_0)238 CALL iom_rstput( 0, 0, inum4, 'gdept_1d', gdept_1d ) ! ! reference z-coord. 239 CALL iom_rstput( 0, 0, inum4, 'gdepw_1d', gdepw_1d ) 240 CALL iom_rstput( 0, 0, inum4, 'e3t_1d' , e3t_1d ) 241 CALL iom_rstput( 0, 0, inum4, 'e3w_1d' , e3w_1d ) 242 242 ENDIF 243 243 244 244 IF( ln_zco ) THEN 245 245 ! ! z-coordinate - full steps 246 CALL iom_rstput( 0, 0, inum4, 'gdept_ 0', gdept_0 )! ! depth247 CALL iom_rstput( 0, 0, inum4, 'gdepw_ 0', gdepw_0)248 CALL iom_rstput( 0, 0, inum4, 'e3t_ 0' , e3t_0 )! ! scale factors249 CALL iom_rstput( 0, 0, inum4, 'e3w_ 0' , e3w_0)246 CALL iom_rstput( 0, 0, inum4, 'gdept_1d', gdept_1d ) ! ! depth 247 CALL iom_rstput( 0, 0, inum4, 'gdepw_1d', gdepw_1d ) 248 CALL iom_rstput( 0, 0, inum4, 'e3t_1d' , e3t_1d ) ! ! scale factors 249 CALL iom_rstput( 0, 0, inum4, 'e3w_1d' , e3w_1d ) 250 250 ENDIF 251 251 ! ! ============================
Note: See TracChangeset
for help on using the changeset viewer.