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/DOM/domhgr.F90 – NEMO

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

rewriting of ice rheology to conserve energy

File:
1 edited

Legend:

Unmodified
Added
Removed

  • 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.