Changeset 9794


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

Location:
NEMO/trunk/tests/ICEDYN
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • NEMO/trunk/tests/ICEDYN/EXPREF/1_namelist_cfg

    r9789 r9794  
    2222   rn_dx       =    1000.   ! horizontal resolution in meters 
    2323   rn_dy       =    1000.   ! horizontal resolution in meters 
    24    rn_ppgphi0  =     70.    ! Reference latitude [degrees] 
     24   ln_corio    =    .false. ! set coriolis to 0 (false) or not (true) 
     25      rn_ppgphi0  =  70.    ! Reference latitude [degrees] 
    2526/ 
    2627!----------------------------------------------------------------------- 
  • NEMO/trunk/tests/ICEDYN/EXPREF/namelist_cfg

    r9789 r9794  
    2222   rn_dx       =    3000.   ! horizontal resolution in meters 
    2323   rn_dy       =    3000.   ! horizontal resolution in meters 
    24    rn_ppgphi0  =     70.    ! Reference latitude [degrees] 
     24   ln_corio    =    .false. ! set coriolis=0 (false) or not (true) 
     25      rn_ppgphi0  =  70.    ! Reference latitude [degrees] 
    2526/ 
    2627!----------------------------------------------------------------------- 
  • 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 
  • NEMO/trunk/tests/ICEDYN/MY_SRC/usrdef_nam.F90

    r9789 r9794  
    3131   REAL(wp), PUBLIC ::   rn_dy      ! resolution in meters defining the horizontal domain size 
    3232   REAL(wp), PUBLIC ::   rn_ppgphi0 ! reference latitude for beta-plane  
     33   LOGICAL , PUBLIC ::   ln_corio   ! set coriolis at 0 (ln_corio=F) or not  
    3334 
    3435   !!---------------------------------------------------------------------- 
     
    6061      REAL(wp)::   zlx, zly  ! Local scalars 
    6162      !! 
    62       NAMELIST/namusr_def/ ln_zco, rn_dx, rn_dy, rn_ppgphi0 
     63      NAMELIST/namusr_def/ ln_zco, rn_dx, rn_dy, ln_corio, rn_ppgphi0 
    6364      !!---------------------------------------------------------------------- 
    6465      ! 
     
    111112      WRITE(ldtxt(ii),*) '                                               jpjglo = ', kpj                ;   ii = ii + 1 
    112113      WRITE(ldtxt(ii),*) '                                               jpkglo = ', kpk                ;   ii = ii + 1 
     114      WRITE(ldtxt(ii),*) '         Coriolis:', ln_corio                                                 ;   ii = ii + 1 
    113115      ! 
    114116      !                             ! Set the lateral boundary condition of the global domain 
Note: See TracChangeset for help on using the changeset viewer.