Ignore:
Timestamp:
2018-06-13T17:44:54+02:00 (2 years ago)
Author:
clem
Message:

add an option in ICEDYN test case to be able to run without coriolis effect

File:
1 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/ICEDYN/MY_SRC/usrdef_hgr.F90

    r9789 r9794  
    1616   USE par_oce         ! ocean space and time domain 
    1717   USE phycst          ! physical constants 
    18    USE usrdef_nam, ONLY: rn_dx, rn_dy, rn_ppgphi0   ! horizontal resolution in meters 
    19    !                                                  and reference latitude 
     18   USE usrdef_nam, ONLY: rn_dx, rn_dy, ln_corio, rn_ppgphi0   ! horizontal resolution in meters 
     19   !                                                            coriolis and reference latitude 
    2020   USE in_out_manager  ! I/O manager 
    2121   USE lib_mpp         ! MPP library 
     
    147147      kff = 1                       !  indicate not to compute Coriolis parameter afterward 
    148148      ! 
    149       zbeta = 2._wp * omega * COS( rad * rn_ppgphi0 ) / ra 
    150       zf0   = 2._wp * omega * SIN( rad * rn_ppgphi0 ) 
    151       pff_f(:,:) = zf0 + zbeta * pphif(:,:) * 1.e+3 
    152       pff_t(:,:) = zf0 + zbeta * pphit(:,:) * 1.e+3 
     149      IF( ln_corio ) THEN 
     150         zbeta = 2._wp * omega * COS( rad * rn_ppgphi0 ) / ra 
     151         zf0   = 2._wp * omega * SIN( rad * rn_ppgphi0 ) 
     152         pff_f(:,:) = zf0 + zbeta * pphif(:,:) * 1.e+3 
     153         pff_t(:,:) = zf0 + zbeta * pphit(:,:) * 1.e+3 
     154      ELSE 
     155         pff_f(:,:) = 0. 
     156         pff_t(:,:) = 0. 
     157      ENDIF 
    153158      ! 
    154159   END SUBROUTINE usr_def_hgr 
Note: See TracChangeset for help on using the changeset viewer.