- Timestamp:
- 2020-05-14T21:46:00+02:00 (4 years ago)
- Location:
- NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser
- Property svn:externals
-
old new 6 6 ^/vendors/FCM@HEAD ext/FCM 7 7 ^/vendors/IOIPSL@HEAD ext/IOIPSL 8 9 # SETTE 10 ^/utils/CI/sette@HEAD sette
-
- Property svn:externals
-
NEMO/branches/2019/dev_r11078_OSMOSIS_IMMERSE_Nurser/src/OCE/STO/stopar.F90
r12178 r12928 112 112 INTEGER :: nn_trc_ord = 1 ! order of autoregressive processes 113 113 114 !! * Substitutions 115 # include "do_loop_substitute.h90" 114 116 !!---------------------------------------------------------------------- 115 117 !! NEMO/OCE 4.0 , NEMO Consortium (2018) … … 261 263 262 264 ! Read namsto namelist : stochastic parameterization 263 REWIND( numnam_ref ) ! Namelist namdyn_adv in reference namelist : Momentum advection scheme264 265 READ ( numnam_ref, namsto, IOSTAT = ios, ERR = 901) 265 266 901 IF( ios /= 0 ) CALL ctl_nam ( ios , 'namsto in reference namelist' ) 266 267 267 REWIND( numnam_cfg ) ! Namelist namdyn_adv in configuration namelist : Momentum advection scheme268 268 READ ( numnam_cfg, namsto, IOSTAT = ios, ERR = 902 ) 269 269 902 IF( ios > 0 ) CALL ctl_nam ( ios , 'namsto in configuration namelist' ) … … 829 829 REAL(KIND=8) :: gran ! Gaussian random number (forced KIND=8 as in kiss_gaussian) 830 830 831 DO jj = 1, jpj 832 DO ji = 1, jpi 833 CALL kiss_gaussian( gran ) 834 psto(ji,jj) = gran 835 END DO 836 END DO 831 DO_2D_11_11 832 CALL kiss_gaussian( gran ) 833 psto(ji,jj) = gran 834 END_2D 837 835 838 836 END SUBROUTINE sto_par_white … … 849 847 INTEGER :: ji, jj 850 848 851 DO jj = 2, jpj-1 852 DO ji = 2, jpi-1 853 psto(ji,jj) = 0.5_wp * psto(ji,jj) + 0.125_wp * & 854 & ( psto(ji-1,jj) + psto(ji+1,jj) + & 855 & psto(ji,jj-1) + psto(ji,jj+1) ) 856 END DO 857 END DO 849 DO_2D_00_00 850 psto(ji,jj) = 0.5_wp * psto(ji,jj) + 0.125_wp * & 851 & ( psto(ji-1,jj) + psto(ji+1,jj) + & 852 & psto(ji,jj-1) + psto(ji,jj+1) ) 853 END_2D 858 854 859 855 END SUBROUTINE sto_par_flt
Note: See TracChangeset
for help on using the changeset viewer.