- Timestamp:
- 2011-07-20T17:34:39+02:00 (13 years ago)
- Location:
- branches/2011/dev_r2787_NOCS_NEPTUNE/NEMOGCM
- Files:
-
- 3 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2011/dev_r2787_NOCS_NEPTUNE/NEMOGCM/ARCH/arch-ALTIX_NAUTILUS4.fcm
r2364 r2808 22 22 # Note use of -Bstatic because the library root directories are not accessible to the back-end compute nodes 23 23 %NCDF_LIB -L%HDF5_HOME/lib -L%NCDF_HOME/lib -Bstatic -lnetcdf -lhdf5_fortran -lhdf5_hl -lhdf5 -Bdynamic -lz 24 %FC mpif9024 %FC ifort 25 25 %FCFLAGS -r8 -O3 -xT -ip -vec-report0 26 26 %FFLAGS -r8 -O3 -xT -ip -vec-report0 27 %LD mpif9027 %LD ifort 28 28 %FPPFLAGS -P -C -traditional 29 %LDFLAGS 29 %LDFLAGS -lmpi 30 30 %AR ar 31 31 %ARFLAGS -r -
branches/2011/dev_r2787_NOCS_NEPTUNE/NEMOGCM/CONFIG/ORCA2_LIM/EXP00/namelist
r2795 r2808 866 866 &nam_dynnept ! Neptune effect (simplified: lateral and vertical diffusions removed) 867 867 !----------------------------------------------------------------------- 868 ! Suggested lengthscale values are those of Eby & Holloway (1994) fo a coarse model 868 869 ln_neptsimp = .false. ! yes/no use simplified neptune 869 870 870 871 ln_smoothtopo = .false. ! yes/no smooth tsu, tsv 871 rn_tslse = 1.2e4 ! value of L at the equator872 rn_tslsp = 3.0e3 ! value of L at the pole873 / 872 rn_tslse = 1.2e4 ! value of lengthscale L at the equator 873 rn_tslsp = 3.0e3 ! value of lengthscale L at the pole 874 / -
branches/2011/dev_r2787_NOCS_NEPTUNE/NEMOGCM/NEMO/OPA_SRC/DYN/dynnept.F90
r2795 r2808 2 2 !!====================================================================== 3 3 !! *** MODULE dynnept *** 4 !! Ocean dynamics: topostressas proposed by Greg Holloway,4 !! Ocean dynamics: Neptune effect as proposed by Greg Holloway, 5 5 !! recoded version of simplest case (u*, v* only) 6 6 !!====================================================================== … … 19 19 !! smooth_topo2 : 20 20 !!---------------------------------------------------------------------- 21 #if defined key_dynnept || defined key_esopa22 21 USE oce ! ocean dynamics and tracers 23 22 USE dom_oce ! ocean space and time domain … … 33 32 34 33 !! * Routine accessibility 35 PUBLIC dyn_nept_init 36 PUBLIC dyn_cor_topo 34 PUBLIC dyn_nept_init ! routine called by nemogcm.F90 35 PUBLIC dyn_cor_topo ! routine called by step.F90 37 36 !! dynnept_alloc() is called only by dyn_nept_init, within this module 38 37 !! div_cur_nept_init is called only by dyn_nept_init, within this module … … 40 39 41 40 !! * Shared module variables 42 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:) :: zunep, zvnep ! topostressu and v43 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: zhdivnep ! hor. div for Neptune vel.44 REAL(wp), ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: zmrotnep ! curl for Neptune vel.41 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:) :: zunep, zvnep ! Neptune u and v 42 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: zhdivnep ! hor. div for Neptune vel. 43 REAL(wp), PUBLIC, ALLOCATABLE, SAVE, DIMENSION(:,:,:) :: zmrotnep ! curl for Neptune vel. 45 44 46 45 47 46 !! * Namelist nam_dynnept variables 48 LOGICAL :: ln_neptsimp = .FALSE. ! yes/no use cortopo49 50 LOGICAL :: ln_smoothtopo = .FALSE. ! yes/no smooth zunep, zvnep51 REAL(wp) :: rn_tslse = 1.2e4 ! value of L at the equator52 REAL(wp) :: rn_tslsp = 3.0e3 ! value of L at the pole47 LOGICAL, PUBLIC :: ln_neptsimp = .FALSE. ! yes/no simplified neptune 48 49 LOGICAL :: ln_smoothtopo = .FALSE. ! yes/no smooth zunep, zvnep 50 REAL(wp) :: rn_tslse = 1.2e4 ! value of L at the equator 51 REAL(wp) :: rn_tslsp = 3.0e3 ! value of L at the pole 53 52 54 53 !! * Module variables … … 87 86 !! ** History : 1.0 ! 07-05 (Zeliang Wang) Original code for zunep, zvnep 88 87 !! 1.1 ! 07-06 (Michael Dunphy) namelist and initialisation 89 !! 2.0 ! 2011-0 6(Jeff Blundell, NOCS)88 !! 2.0 ! 2011-07 (Jeff Blundell, NOCS) 90 89 !! ! Simplified form for temporally invariant u*, v* 91 90 !! ! Horizontal and vertical diffusivity formulations removed … … 112 111 !!---------------------------------------------------------------------- 113 112 114 ! Define the topostressparameters115 ! ================================ 113 ! Define the (simplified) Neptune parameters 114 ! ========================================== 116 115 117 116 !! WRITE(numout,*) ' start dynnept namelist' 118 117 !! CALL FLUSH(numout) 119 REWIND( numnam ) ! Read Namelist nam_dynnept: Topostress118 REWIND( numnam ) ! Read Namelist nam_dynnept: Simplified Neptune 120 119 READ ( numnam, nam_dynnept ) 121 120 !! WRITE(numout,*) ' dynnept namelist done' … … 124 123 IF(lwp) THEN ! Control print 125 124 WRITE(numout,*) 126 WRITE(numout,*) 'dyn_nept_init : Topostressmodule enabled'125 WRITE(numout,*) 'dyn_nept_init : Simplified Neptune module enabled' 127 126 WRITE(numout,*) '~~~~~~~~~~~~~' 128 127 WRITE(numout,*) ' --> Reading namelist nam_dynnept parameters:' … … 190 189 IF( ln_smoothtopo ) THEN 191 190 CALL smooth_topo2( htn, ht, .TRUE. ) 192 !! overwrites ht with a smoothed ver ion of htn191 !! overwrites ht with a smoothed version of htn 193 192 ELSE 194 193 ht(:,:) = htn(:,:) … … 212 211 IF( ln_smoothtopo ) THEN 213 212 CALL smooth_topo2( hu, hu_n, .TRUE. ) 214 !! overwrites hu_n with a smoothed ver ion of hu213 !! overwrites hu_n with a smoothed version of hu 215 214 ELSE 216 215 hu_n(:,:) = hu(:,:) … … 229 228 IF( ln_smoothtopo ) THEN 230 229 CALL smooth_topo2( hv, hv_n, .TRUE. ) 231 !! overwrites hv_n with a smoothed ver ion of hv230 !! overwrites hv_n with a smoothed version of hv 232 231 ELSE 233 232 hv_n(:,:) = hv(:,:) … … 548 547 END SUBROUTINE smooth_topo2 549 548 550 #else551 !!------------------------------------------------------------------------552 !! Dummy module : No Topostress553 !!------------------------------------------------------------------------554 CONTAINS555 SUBROUTINE dyn_nept_init ! Dummy routine556 WRITE(*,*) 'dyn_nept_init: You should not have seen this print! error?'557 END SUBROUTINE dyn_nept_init558 559 SUBROUTINE dyn_cor_topo( kt ) ! Dummy routine560 INTEGER, INTENT(in) :: kt561 WRITE(*,*) 'dyn_cor_topo: You should not have seen this print! error?', kt562 END SUBROUTINE dyn_cor_topo563 564 SUBROUTINE div_cur_nept_init ! Dummy routine565 WRITE(*,*) 'div_cur_nept_init: You should not have seen this print! error?'566 END SUBROUTINE div_cur_nept_init567 #endif568 !!======================================================================569 549 END MODULE dynnept
Note: See TracChangeset
for help on using the changeset viewer.