Changeset 6689
- Timestamp:
- 2016-06-13T16:24:43+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domwri.F90
r5836 r6689 71 71 ! 72 72 ! ! horizontal mesh (inum3) 73 CALL iom_rstput( 0, 0, inum0, 'glamt', glamt, ktype = jp_r 4) ! ! latitude74 CALL iom_rstput( 0, 0, inum0, 'glamu', glamu, ktype = jp_r 4)75 CALL iom_rstput( 0, 0, inum0, 'glamv', glamv, ktype = jp_r 4)76 CALL iom_rstput( 0, 0, inum0, 'glamf', glamf, ktype = jp_r 4)77 78 CALL iom_rstput( 0, 0, inum0, 'gphit', gphit, ktype = jp_r 4) ! ! longitude79 CALL iom_rstput( 0, 0, inum0, 'gphiu', gphiu, ktype = jp_r 4)80 CALL iom_rstput( 0, 0, inum0, 'gphiv', gphiv, ktype = jp_r 4)81 CALL iom_rstput( 0, 0, inum0, 'gphif', gphif, ktype = jp_r 4)73 CALL iom_rstput( 0, 0, inum0, 'glamt', glamt, ktype = jp_r8 ) ! ! latitude 74 CALL iom_rstput( 0, 0, inum0, 'glamu', glamu, ktype = jp_r8 ) 75 CALL iom_rstput( 0, 0, inum0, 'glamv', glamv, ktype = jp_r8 ) 76 CALL iom_rstput( 0, 0, inum0, 'glamf', glamf, ktype = jp_r8 ) 77 78 CALL iom_rstput( 0, 0, inum0, 'gphit', gphit, ktype = jp_r8 ) ! ! longitude 79 CALL iom_rstput( 0, 0, inum0, 'gphiu', gphiu, ktype = jp_r8 ) 80 CALL iom_rstput( 0, 0, inum0, 'gphiv', gphiv, ktype = jp_r8 ) 81 CALL iom_rstput( 0, 0, inum0, 'gphif', gphif, ktype = jp_r8 ) 82 82 83 83 CALL iom_rstput( 0, 0, inum0, 'e1t', e1t, ktype = jp_r8 ) ! ! e1 scale factors … … 229 229 230 230 ! ! horizontal mesh (inum3) 231 CALL iom_rstput( 0, 0, inum3, 'glamt', glamt, ktype = jp_r 4) ! ! latitude232 CALL iom_rstput( 0, 0, inum3, 'glamu', glamu, ktype = jp_r 4)233 CALL iom_rstput( 0, 0, inum3, 'glamv', glamv, ktype = jp_r 4)234 CALL iom_rstput( 0, 0, inum3, 'glamf', glamf, ktype = jp_r 4)235 236 CALL iom_rstput( 0, 0, inum3, 'gphit', gphit, ktype = jp_r 4) ! ! longitude237 CALL iom_rstput( 0, 0, inum3, 'gphiu', gphiu, ktype = jp_r 4)238 CALL iom_rstput( 0, 0, inum3, 'gphiv', gphiv, ktype = jp_r 4)239 CALL iom_rstput( 0, 0, inum3, 'gphif', gphif, ktype = jp_r 4)231 CALL iom_rstput( 0, 0, inum3, 'glamt', glamt, ktype = jp_r8 ) ! ! latitude 232 CALL iom_rstput( 0, 0, inum3, 'glamu', glamu, ktype = jp_r8 ) 233 CALL iom_rstput( 0, 0, inum3, 'glamv', glamv, ktype = jp_r8 ) 234 CALL iom_rstput( 0, 0, inum3, 'glamf', glamf, ktype = jp_r8 ) 235 236 CALL iom_rstput( 0, 0, inum3, 'gphit', gphit, ktype = jp_r8 ) ! ! longitude 237 CALL iom_rstput( 0, 0, inum3, 'gphiu', gphiu, ktype = jp_r8 ) 238 CALL iom_rstput( 0, 0, inum3, 'gphiv', gphiv, ktype = jp_r8 ) 239 CALL iom_rstput( 0, 0, inum3, 'gphif', gphif, ktype = jp_r8 ) 240 240 241 241 CALL iom_rstput( 0, 0, inum3, 'e1t', e1t, ktype = jp_r8 ) ! ! e1 scale factors … … 257 257 CALL iom_rstput( 0, 0, inum4, 'misf', zprt, ktype = jp_i2 ) ! ! nb of ocean T-points 258 258 zprt(:,:) = ssmask(:,:) * REAL( risfdep(:,:) , wp ) 259 CALL iom_rstput( 0, 0, inum4, 'isfdraft', zprt, ktype = jp_r 4) ! ! nb of ocean T-points259 CALL iom_rstput( 0, 0, inum4, 'isfdraft', zprt, ktype = jp_r8 ) ! ! nb of ocean T-points 260 260 261 261 IF( ln_sco ) THEN ! s-coordinate … … 279 279 CALL iom_rstput( 0, 0, inum4, 'gdept_1d' , gdept_1d ) ! ! stretched system 280 280 CALL iom_rstput( 0, 0, inum4, 'gdepw_1d' , gdepw_1d ) 281 CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0, ktype = jp_r 4)282 CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r 4)281 CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0, ktype = jp_r8 ) 282 CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r8 ) 283 283 ENDIF 284 284 … … 302 302 ! 303 303 IF( nmsh <= 3 ) THEN ! ! 3D depth 304 CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0, ktype = jp_r 4)304 CALL iom_rstput( 0, 0, inum4, 'gdept_0', gdept_0, ktype = jp_r8 ) 305 305 DO jk = 1,jpk 306 306 DO jj = 1, jpjm1 … … 312 312 END DO 313 313 CALL lbc_lnk( zdepu, 'U', 1. ) ; CALL lbc_lnk( zdepv, 'V', 1. ) 314 CALL iom_rstput( 0, 0, inum4, 'gdepu', zdepu, ktype = jp_r 4)315 CALL iom_rstput( 0, 0, inum4, 'gdepv', zdepv, ktype = jp_r 4)316 CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r 4)314 CALL iom_rstput( 0, 0, inum4, 'gdepu', zdepu, ktype = jp_r8 ) 315 CALL iom_rstput( 0, 0, inum4, 'gdepv', zdepv, ktype = jp_r8 ) 316 CALL iom_rstput( 0, 0, inum4, 'gdepw_0', gdepw_0, ktype = jp_r8 ) 317 317 ELSE ! ! 2D bottom depth 318 318 DO jj = 1,jpj … … 322 322 END DO 323 323 END DO 324 CALL iom_rstput( 0, 0, inum4, 'hdept', zprt, ktype = jp_r 4)325 CALL iom_rstput( 0, 0, inum4, 'hdepw', zprw, ktype = jp_r 4)324 CALL iom_rstput( 0, 0, inum4, 'hdept', zprt, ktype = jp_r8 ) 325 CALL iom_rstput( 0, 0, inum4, 'hdepw', zprw, ktype = jp_r8 ) 326 326 ENDIF 327 327 !
Note: See TracChangeset
for help on using the changeset viewer.