Changeset 9794 for NEMO/trunk/tests/ICEDYN/MY_SRC
- Timestamp:
- 2018-06-13T17:44:54+02:00 (6 years ago)
- Location:
- NEMO/trunk/tests/ICEDYN/MY_SRC
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/ICEDYN/MY_SRC/usrdef_hgr.F90
r9789 r9794 16 16 USE par_oce ! ocean space and time domain 17 17 USE phycst ! physical constants 18 USE usrdef_nam, ONLY: rn_dx, rn_dy, rn_ppgphi0 ! horizontal resolution in meters19 ! and reference latitude18 USE usrdef_nam, ONLY: rn_dx, rn_dy, ln_corio, rn_ppgphi0 ! horizontal resolution in meters 19 ! coriolis and reference latitude 20 20 USE in_out_manager ! I/O manager 21 21 USE lib_mpp ! MPP library … … 147 147 kff = 1 ! indicate not to compute Coriolis parameter afterward 148 148 ! 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 153 158 ! 154 159 END SUBROUTINE usr_def_hgr -
NEMO/trunk/tests/ICEDYN/MY_SRC/usrdef_nam.F90
r9789 r9794 31 31 REAL(wp), PUBLIC :: rn_dy ! resolution in meters defining the horizontal domain size 32 32 REAL(wp), PUBLIC :: rn_ppgphi0 ! reference latitude for beta-plane 33 LOGICAL , PUBLIC :: ln_corio ! set coriolis at 0 (ln_corio=F) or not 33 34 34 35 !!---------------------------------------------------------------------- … … 60 61 REAL(wp):: zlx, zly ! Local scalars 61 62 !! 62 NAMELIST/namusr_def/ ln_zco, rn_dx, rn_dy, rn_ppgphi063 NAMELIST/namusr_def/ ln_zco, rn_dx, rn_dy, ln_corio, rn_ppgphi0 63 64 !!---------------------------------------------------------------------- 64 65 ! … … 111 112 WRITE(ldtxt(ii),*) ' jpjglo = ', kpj ; ii = ii + 1 112 113 WRITE(ldtxt(ii),*) ' jpkglo = ', kpk ; ii = ii + 1 114 WRITE(ldtxt(ii),*) ' Coriolis:', ln_corio ; ii = ii + 1 113 115 ! 114 116 ! ! Set the lateral boundary condition of the global domain
Note: See TracChangeset
for help on using the changeset viewer.