Changeset 7069 for branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90
- Timestamp:
- 2016-10-21T19:24:28+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90
r6796 r7069 451 451 e1v(:,:) = ze1 ; e2v(:,:) = ze1 452 452 e1f(:,:) = ze1 ; e2f(:,:) = ze1 453 454 CASE ( 6 ) ! clem: f-plane with irregular grid-spacing 455 456 IF(lwp) WRITE(numout,*) 457 IF(lwp) WRITE(numout,*) ' f-plane with irregular grid-spacing (+- 10%)' 458 IF(lwp) WRITE(numout,*) ' the max is given by ppe1_m and ppe2_m' 459 460 ! Position coordinates (in kilometers) 461 ! ========== 462 glam0 = 0._wp 463 gphi0 = 0._wp 464 465 DO jj = 1, jpj 466 DO ji = 1, jpi 467 zti = FLOAT( ji - 1 + nimpp - 1 ) ; ztj = FLOAT( jj - 1 + njmpp - 1 ) 468 zui = FLOAT( ji - 1 + nimpp - 1 ) + 0.5 ; zuj = FLOAT( jj - 1 + njmpp - 1 ) 469 zvi = FLOAT( ji - 1 + nimpp - 1 ) ; zvj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5 470 zfi = FLOAT( ji - 1 + nimpp - 1 ) + 0.5 ; zfj = FLOAT( jj - 1 + njmpp - 1 ) + 0.5 471 472 glamt(ji,jj) = glam0 + ppe1_m * 1.e-5 * zti 473 glamu(ji,jj) = glam0 + ppe1_m * 1.e-5 * zui 474 glamv(ji,jj) = glam0 + ppe1_m * 1.e-5 * zvi 475 glamf(ji,jj) = glam0 + ppe1_m * 1.e-5 * zfi 476 477 gphit(ji,jj) = gphi0 + ppe2_m * 1.e-5 * ztj 478 gphiu(ji,jj) = gphi0 + ppe2_m * 1.e-5 * zuj 479 gphiv(ji,jj) = gphi0 + ppe2_m * 1.e-5 * zvj 480 gphif(ji,jj) = gphi0 + ppe2_m * 1.e-5 * zfj 481 END DO 482 END DO 483 484 ! Horizontal scale factors (in meters) 485 ! ====== 486 DO jj = 1, jpj 487 DO ji = 1, jpi 488 e1t(ji,jj) = ppe1_m * EXP(-0.8/REAL(jpiglo**2)*(mi0(ji)-REAL(jpiglo+1)*0.5)**2) 489 e2t(ji,jj) = ppe2_m * EXP(-0.8/REAL(jpjglo**2)*(mj0(jj)-REAL(jpjglo+1)*0.5)**2) 490 END DO 491 END DO 492 e1u(:,:) = e1t(:,:) ; e2u(:,:) = e2t(:,:) 493 e1v(:,:) = e1t(:,:) ; e2v(:,:) = e2t(:,:) 494 e1f(:,:) = e1t(:,:) ; e2f(:,:) = e2t(:,:) 453 495 454 496 CASE DEFAULT
Note: See TracChangeset
for help on using the changeset viewer.