Changeset 497 for trunk/NEMO/OFF_SRC/DOM/domrea_dimg.h90
- Timestamp:
- 2006-09-12T13:03:53+02:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OFF_SRC/DOM/domrea_dimg.h90
r343 r497 71 71 irecl8 = jpi*jpj*wp 72 72 ios1 = 0 ; ios2 = 0 73 IF (l k_zps ) ios1= 174 IF (l k_sco ) ios2= 173 IF (ln_zps ) ios1= 1 74 IF (ln_sco ) ios2= 1 75 75 76 76 SELECT CASE (nmsh ) … … 262 262 mbathy=zprt*tmask(:,:,1)+1 263 263 264 #if defined key_s_coord 265 ! 266 ! hbat 267 irec = irecv(inum4) + (narea - 1 ) 268 READ(inum4,REC=irec) hbatt(:,:) 269 irecv(inum4) = irecv(inum4) + jpnij 270 271 irec = irecv(inum4) + (narea - 1 ) 272 READ(inum4,REC=irec) hbatu(:,:) 273 irecv(inum4) = irecv(inum4) + jpnij 274 275 irec = irecv(inum4) + (narea - 1 ) 276 READ(inum4,REC=irec) hbatv(:,:) 277 irecv(inum4) = irecv(inum4) + jpnij 278 279 irec = irecv(inum4) + (narea - 1 ) 280 READ(inum4,REC=irec) hbatf(:,:) 281 irecv(inum4) = irecv(inum4) + jpnij 282 ! 283 ! gsig and esig ( as vectors of jpk per record ) 284 irec = irecv(inum4) + (narea - 1 ) 285 READ(inum4,REC=irec) gsigt(:) 286 irecv(inum4) = irecv(inum4) + jpnij 287 288 irec = irecv(inum4) + (narea - 1 ) 289 READ(inum4,REC=irec) gsigw(:) 290 irecv(inum4) = irecv(inum4) + jpnij 291 292 irec = irecv(inum4) + (narea - 1 ) 293 READ(inum4,REC=irec) gsi3w(:) 294 irecv(inum4) = irecv(inum4) + jpnij 295 296 irec = irecv(inum4) + (narea - 1 ) 297 READ(inum4,REC=irec) esigt(:) 298 irecv(inum4) = irecv(inum4) + jpnij 299 300 irec = irecv(inum4) + (narea - 1 ) 301 READ(inum4,REC=irec) esigw(:) 302 irecv(inum4) = irecv(inum4) + jpnij 303 304 # elif defined key_partial_steps 305 ! 306 ! hdep 307 irec = irecv(inum4) + (narea - 1 ) 308 READ(inum4,REC=irec) hdept(:,:) 309 irecv(inum4) = irecv(inum4) + jpnij 310 311 irec = irecv(inum4) + (narea - 1 ) 312 READ(inum4,REC=irec) hdepw(:,:) 313 irecv(inum4) = irecv(inum4) + jpnij 314 ! 315 ! e3t_ps (3D) 316 DO jk=1,jpk 317 irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 318 READ(inum4,REC=irec) e3t_ps(:,:,jk) 319 END DO 320 irecv(inum4) = irecv(inum4) + jpk * jpnij 321 322 ! e3u_ps e3v_ps e3w_ps (3D) 323 DO jk=1,jpk 324 irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 325 READ(inum4,REC=irec) e3u_ps(:,:,jk) 326 END DO 327 irecv(inum4) = irecv(inum4) + jpk * jpnij 328 329 DO jk=1,jpk 330 irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 331 READ(inum4,REC=irec) e3v_ps(:,:,jk) 332 END DO 333 irecv(inum4) = irecv(inum4) + jpk * jpnij 334 335 DO jk=1,jpk 336 irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 337 READ(inum4,REC=irec) e3w_ps(:,:,jk) 338 END DO 339 irecv(inum4) = irecv(inum4) + jpk * jpnij 340 ! 341 ! 342 ! gdep 343 irec = irecv(inum4) + (narea - 1 ) 344 READ(inum4,REC=irec) gdept(:) 345 irecv(inum4) = irecv(inum4) + jpnij 346 347 irec = irecv(inum4) + (narea - 1 ) 348 READ(inum4,REC=irec) gdepw(:) 349 irecv(inum4) = irecv(inum4) + jpnij 350 ! 351 ! e3 352 irec = irecv(inum4) + (narea - 1 ) 353 READ(inum4,REC=irec) e3t(:) 354 irecv(inum4) = irecv(inum4) + jpnij 355 356 irec = irecv(inum4) + (narea - 1 ) 357 READ(inum4,REC=irec) e3w(:) 264 #if ! defined key_zco 265 IF( ln_sco ) THEN ! s-coordinate 266 ! 267 ! hbat 268 irec = irecv(inum4) + (narea - 1 ) 269 READ(inum4,REC=irec) hbatt(:,:) 270 irecv(inum4) = irecv(inum4) + jpnij 271 272 irec = irecv(inum4) + (narea - 1 ) 273 READ(inum4,REC=irec) hbatu(:,:) 274 irecv(inum4) = irecv(inum4) + jpnij 275 276 irec = irecv(inum4) + (narea - 1 ) 277 READ(inum4,REC=irec) hbatv(:,:) 278 irecv(inum4) = irecv(inum4) + jpnij 279 280 irec = irecv(inum4) + (narea - 1 ) 281 READ(inum4,REC=irec) hbatf(:,:) 282 irecv(inum4) = irecv(inum4) + jpnij 283 ! 284 ! gsig and esig ( as vectors of jpk per record ) 285 irec = irecv(inum4) + (narea - 1 ) 286 READ(inum4,REC=irec) gsigt(:) 287 irecv(inum4) = irecv(inum4) + jpnij 288 289 irec = irecv(inum4) + (narea - 1 ) 290 READ(inum4,REC=irec) gsigw(:) 291 irecv(inum4) = irecv(inum4) + jpnij 292 293 irec = irecv(inum4) + (narea - 1 ) 294 READ(inum4,REC=irec) gsi3w(:) 295 irecv(inum4) = irecv(inum4) + jpnij 296 297 irec = irecv(inum4) + (narea - 1 ) 298 READ(inum4,REC=irec) esigt(:) 299 irecv(inum4) = irecv(inum4) + jpnij 300 301 irec = irecv(inum4) + (narea - 1 ) 302 READ(inum4,REC=irec) esigw(:) 303 irecv(inum4) = irecv(inum4) + jpnij 304 305 ! e3 (3D) 306 DO jk=1,jpk 307 irec = irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 308 READ(inum4,REC=irec) e3t(:,:,jk) 309 ENDIF 310 irecv(inum4) = irecv(inum4) + jpk * jpnij 311 312 DO jk=1,jpk 313 irec = irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 314 READ(inum4,REC=irec) e3u(:,:,jk) 315 ENDIF 316 irecv(inum4) = irecv(inum4) + jpk * jpnij 317 318 DO jk=1,jpk 319 irec = irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 320 READ(inum4,REC=irec) e3v(:,:,jk) 321 ENDIF 322 irecv(inum4) = irecv(inum4) + jpk * jpnij 323 324 DO jk=1,jpk 325 irec = irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 326 READ(inum4,REC=irec) e3w(:,:,jk) 327 ENDIF 328 irecv(inum4) = irecv(inum4) + jpk * jpnij 329 330 ! gdep 331 irec = irecv(inum4) + (narea - 1 ) 332 READ(inum4,REC=irec) gdept_0(:) 333 irecv(inum4) = irecv(inum4) + jpnij 334 335 irec = irecv(inum4) + (narea - 1 ) 336 READ(inum4,REC=irec) gdepw_0(:) 337 irecv(inum4) = irecv(inum4) + jpnij 338 ! 339 ENDIF 340 341 IF( ln_zps ) THEN 342 ! 343 ! hdep 344 irec = irecv(inum4) + (narea - 1 ) 345 READ(inum4,REC=irec) hdept(:,:) 346 irecv(inum4) = irecv(inum4) + jpnij 347 348 irec = irecv(inum4) + (narea - 1 ) 349 READ(inum4,REC=irec) hdepw(:,:) 350 irecv(inum4) = irecv(inum4) + jpnij 351 ! 352 ! e3t_ps (3D) 353 DO jk=1,jpk 354 irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 355 READ(inum4,REC=irec) e3t(:,:,jk) 356 END DO 357 irecv(inum4) = irecv(inum4) + jpk * jpnij 358 359 ! e3u_ps e3v_ps e3w_ps (3D) 360 DO jk=1,jpk 361 irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 362 READ(inum4,REC=irec) e3u(:,:,jk) 363 END DO 364 irecv(inum4) = irecv(inum4) + jpk * jpnij 365 366 DO jk=1,jpk 367 irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 368 READ(inum4,REC=irec) e3v(:,:,jk) 369 END DO 370 irecv(inum4) = irecv(inum4) + jpk * jpnij 371 372 DO jk=1,jpk 373 irec=irecv(inum4) + (jk -1 ) + jpk * (narea - 1 ) 374 READ(inum4,REC=irec) e3w(:,:,jk) 375 END DO 376 irecv(inum4) = irecv(inum4) + jpk * jpnij 377 ! 378 ! 379 ! gdep 380 irec = irecv(inum4) + (narea - 1 ) 381 READ(inum4,REC=irec) gdept_0(:) 382 irecv(inum4) = irecv(inum4) + jpnij 383 384 irec = irecv(inum4) + (narea - 1 ) 385 READ(inum4,REC=irec) gdepw_0(:) 386 irecv(inum4) = irecv(inum4) + jpnij 387 ! 388 ! e3 389 irec = irecv(inum4) + (narea - 1 ) 390 READ(inum4,REC=irec) e3t_0(:) 391 irecv(inum4) = irecv(inum4) + jpnij 392 393 irec = irecv(inum4) + (narea - 1 ) 394 READ(inum4,REC=irec) e3w_0(:) 395 ENDIF 358 396 #else 359 !360 ! gdep361 irec = irecv(inum4) + (narea - 1 )362 READ(inum4,REC=irec) gdept(:)363 irecv(inum4) = irecv(inum4) + jpnij364 365 irec = irecv(inum4) + (narea - 1 )366 READ(inum4,REC=irec) gdepw(:)367 irecv(inum4) = irecv(inum4) + jpnij368 !369 ! e3370 irec = irecv(inum4) + (narea - 1 )371 READ(inum4,REC=irec) e3t(:)372 irecv(inum4) = irecv(inum4) + jpnij373 374 irec = irecv(inum4) + (narea - 1 )375 READ(inum4,REC=irec) e3w(:)376 irecv(inum4) = irecv(inum4) + jpnij397 ! 398 ! gdep 399 irec = irecv(inum4) + (narea - 1 ) 400 READ(inum4,REC=irec) gdept_0(:) 401 irecv(inum4) = irecv(inum4) + jpnij 402 403 irec = irecv(inum4) + (narea - 1 ) 404 READ(inum4,REC=irec) gdepw_0(:) 405 irecv(inum4) = irecv(inum4) + jpnij 406 ! 407 ! e3 408 irec = irecv(inum4) + (narea - 1 ) 409 READ(inum4,REC=irec) e3t_0(:) 410 irecv(inum4) = irecv(inum4) + jpnij 411 412 irec = irecv(inum4) + (narea - 1 ) 413 READ(inum4,REC=irec) e3w_0(:) 414 irecv(inum4) = irecv(inum4) + jpnij 377 415 ! 378 416 #endif
Note: See TracChangeset
for help on using the changeset viewer.