Changeset 4366


Ignore:
Timestamp:
2014-01-22T15:57:03+01:00 (7 years ago)
Author:
jchanut
Message:

correct vvl initialization, ticket #1227

Location:
branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM
Files:
3 edited

Legend:

Unmodified
Added
Removed
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domain.F90

    r4292 r4366  
    9090      ! 
    9191      IF( lk_c1d         )   CALL cor_c1d      ! 1D configuration: Coriolis set at T-point 
    92       ! 
    93       ! - ML - Used in dom_vvl_sf_nxt and lateral diffusion routines 
    94       !        but could be usefull in many other routines 
    95       e12t    (:,:) = e1t(:,:) * e2t(:,:) 
    96       e12u    (:,:) = e1u(:,:) * e2u(:,:) 
    97       e12v    (:,:) = e1v(:,:) * e2v(:,:) 
    98       e12f    (:,:) = e1f(:,:) * e2f(:,:) 
    99       r1_e12t (:,:) = 1._wp    / e12t(:,:) 
    100       r1_e12u (:,:) = 1._wp    / e12u(:,:) 
    101       r1_e12v (:,:) = 1._wp    / e12v(:,:) 
    102       r1_e12f (:,:) = 1._wp    / e12f(:,:) 
    103       re2u_e1u(:,:) = e2u(:,:) / e1u(:,:) 
    104       re1v_e2v(:,:) = e1v(:,:) / e2v(:,:) 
    10592      ! 
    10693      hu(:,:) = 0._wp                          ! Ocean depth at U- and V-points 
     
    460447   !!====================================================================== 
    461448END MODULE domain 
     449 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r4147 r4366  
    458458      ! -------------- 
    459459      e1e2t(:,:) = e1t(:,:) * e2t(:,:) 
    460  
     460     
     461      ! Useful shortcuts (JC: note the duplicated e2e2t array ! Need some cleaning) 
     462      ! --------------------------------------------------------------------------- 
     463      e12t    (:,:) = e1t(:,:) * e2t(:,:) 
     464      e12u    (:,:) = e1u(:,:) * e2u(:,:) 
     465      e12v    (:,:) = e1v(:,:) * e2v(:,:) 
     466      e12f    (:,:) = e1f(:,:) * e2f(:,:) 
     467      r1_e12t (:,:) = 1._wp    / e12t(:,:) 
     468      r1_e12u (:,:) = 1._wp    / e12u(:,:) 
     469      r1_e12v (:,:) = 1._wp    / e12v(:,:) 
     470      r1_e12f (:,:) = 1._wp    / e12f(:,:) 
     471      re2u_e1u(:,:) = e2u(:,:) / e1u(:,:) 
     472      re1v_e2v(:,:) = e1v(:,:) / e2v(:,:) 
    461473 
    462474      ! Control printing : Grid informations (if not restart) 
  • branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

    r4341 r4366  
    772772         IF( ln_rstart ) THEN                   !* Read the restart file 
    773773            CALL rst_read_open                  !  open the restart file if necessary 
     774            CALL iom_get( numror, jpdom_autoglo, 'sshn'   , sshn    ) 
     775            ! 
    774776            id1 = iom_varid( numror, 'fse3t_b', ldstop = .FALSE. ) 
    775777            id2 = iom_varid( numror, 'fse3t_n', ldstop = .FALSE. ) 
     
    823825            fse3t_b(:,:,:) = e3t_0(:,:,:) 
    824826            fse3t_n(:,:,:) = e3t_0(:,:,:) 
     827            sshn(:,:) = 0.0_wp 
    825828            IF( ln_vvl_ztilde .OR. ln_vvl_layer) THEN 
    826829               tilde_e3t_b(:,:,:) = 0.0_wp 
Note: See TracChangeset for help on using the changeset viewer.