- Timestamp:
- 2020-01-27T15:31:53+01:00 (4 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/branches/2019/dev_r11943_MERGE_2019/src/OCE/STO/stopar.F90
r12336 r12340 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) … … 827 829 REAL(KIND=8) :: gran ! Gaussian random number (forced KIND=8 as in kiss_gaussian) 828 830 829 DO jj = 1, jpj 830 DO ji = 1, jpi 831 CALL kiss_gaussian( gran ) 832 psto(ji,jj) = gran 833 END DO 834 END DO 831 DO_2D_11_11 832 CALL kiss_gaussian( gran ) 833 psto(ji,jj) = gran 834 END_2D 835 835 836 836 END SUBROUTINE sto_par_white … … 847 847 INTEGER :: ji, jj 848 848 849 DO jj = 2, jpjm1 850 DO ji = 2, jpim1 851 psto(ji,jj) = 0.5_wp * psto(ji,jj) + 0.125_wp * & 852 & ( psto(ji-1,jj) + psto(ji+1,jj) + & 853 & psto(ji,jj-1) + psto(ji,jj+1) ) 854 END DO 855 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 856 854 857 855 END SUBROUTINE sto_par_flt
Note: See TracChangeset
for help on using the changeset viewer.