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 6796 for branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC – NEMO

Ignore:
Timestamp:
2016-07-06T17:35:54+02:00 (8 years ago)
Author:
clem
Message:

rewriting of ice rheology to conserve energy

Location:
branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM
Files:
2 edited

Legend:

Unmodified
Added
Removed

  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM/dom_oce.F90

    r5123 r6796  
    168168   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::  e1e2t          !: surface at t-point (m2) 
    169169   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::  ff             !: coriolis factor (2.*omega*sin(yphi) ) (s-1) 
     170   REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) ::  ff_t           !: clem: coriolis factor at T-point (s-1) 
    170171 
    171172   !!---------------------------------------------------------------------- 
    … …  
    350351         &      glamv(jpi,jpj) , gphiv(jpi,jpj) , e1v(jpi,jpj) , e2v(jpi,jpj) , r1_e1v(jpi,jpj) , r1_e2v(jpi,jpj) ,   &   
    351352         &      glamf(jpi,jpj) , gphif(jpi,jpj) , e1f(jpi,jpj) , e2f(jpi,jpj) , r1_e1f(jpi,jpj) , r1_e2f(jpi,jpj) ,   & 
    352          &      e1e2t(jpi,jpj) , ff  (jpi,jpj) , STAT=ierr(3) )      
     353         &      e1e2t(jpi,jpj) , ff_t (jpi,jpj) , ff (jpi,jpj) , STAT=ierr(3) )      
    353354         ! 
    354355      ALLOCATE( gdep3w_0(jpi,jpj,jpk) , e3v_0(jpi,jpj,jpk) , e3f_0 (jpi,jpj,jpk) ,                         & 

  • branches/2016/dev_v3_6_STABLE_r6506_AGRIF_LIM3/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

    r6204 r6796  
    528528      CASE ( 0, 1, 4 )               ! mesh on the sphere 
    529529 
    530          ff(:,:) = 2. * omega * SIN( rad * gphif(:,:) )  
     530         ff  (:,:) = 2. * omega * SIN( rad * gphif(:,:) )  
     531         ff_t(:,:) = 2. * omega * SIN( rad * gphit(:,:) )  ! clem: coriolis at T-point  
    531532 
    532533      CASE ( 2 )                     ! f-plane at ppgphi0  
    533534 
    534          ff(:,:) = 2. * omega * SIN( rad * ppgphi0 ) 
     535         ff  (:,:) = 2. * omega * SIN( rad * ppgphi0 ) 
     536         ff_t(:,:) = 2. * omega * SIN( rad * ppgphi0 )  ! clem: coriolis at T-point 
    535537 
    536538         IF(lwp) WRITE(numout,*) '          f-plane: Coriolis parameter = constant = ', ff(1,1) 
    … …  
    551553         zf0     = 2. * omega * SIN( rad * zphi0 )                              ! compute f0 1st point south 
    552554 
    553          ff(:,:) = ( zf0  + zbeta * gphif(:,:) * 1.e+3 )                        ! f = f0 +beta* y ( y=0 at south) 
     555         ff  (:,:) = ( zf0  + zbeta * gphif(:,:) * 1.e+3 )                        ! f = f0 +beta* y ( y=0 at south) 
     556         ff_t(:,:) = ( zf0  + zbeta * gphit(:,:) * 1.e+3 )                        ! clem: coriolis at T-point 
    554557          
    555558         IF(lwp) THEN 
    … …  
    572575         zf0   = 2. * omega * SIN( rad * zphi0 )                            ! compute f0 1st point south 
    573576 
    574          ff(:,:) = ( zf0 + zbeta * ABS( gphif(:,:) - zphi0 ) * rad * ra )   ! f = f0 +beta* y ( y=0 at south) 
     577         ff  (:,:) = ( zf0 + zbeta * ABS( gphif(:,:) - zphi0 ) * rad * ra )   ! f = f0 +beta* y ( y=0 at south) 
     578         ff_t(:,:) = ( zf0 + zbeta * ABS( gphit(:,:) - zphi0 ) * rad * ra )   ! clem: coriolis at T-point 
    575579 
    576580         IF(lwp) THEN 
Note: See TracChangeset for help on using the changeset viewer.