New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 7069 for branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90 – NEMO

Ignore:
Timestamp:
2016-10-21T19:24:28+02:00 (7 years ago)
Author:
clem
Message:

agrif+lim3 update + trunk update

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  
    451451         e1v(:,:) =  ze1     ;      e2v(:,:) = ze1 
    452452         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(:,:) 
    453495 
    454496      CASE DEFAULT 
Note: See TracChangeset for help on using the changeset viewer.